[논문 리뷰] A homogenized limit for the 2D Euler equations in a perforated domain

Matthieu Hillairet, Christophe Lacave외 1명|arXiv (Cornell University)|2019. 07. 09.

Advanced Mathematical Modeling in Engineering참고 문헌 28인용 수 2

한 줄 요약

이 논문은 구멍 크기와 구멍 간격의 비율이 양수인 상수로 수렴하는 임계 영역에서 다공성 영역 내 2차원 오일러 방정식에 대한 최초의 균질화 극한을 확립한다. 이 경우 다공성 매질의 체적 분율은 0이 아니며 작다. 반사 방법을 사용하여, 발산-회전 문제에 대한 정밀한 추정을 도출함으로써, 비산란성과 속도를 연결하는 효과적 행렬 항을 포함하는 균질화된 오일러 유형의 시스템을 유도한다. 이는 이전의 희박한 또는 조밀한 극한에서 다루지 못한, 비영이지만 작은 다공성 비율 영역에 대해 유효하다.

ABSTRACT

We study the motion of an ideal incompressible fluid in a perforated domain. The porous medium is composed of inclusions of size $a$ separated by distances $ ilde d$ and the fluid fills the exterior. We analyse the asymptotic behavior of the fluid when $(a, ilde d) o (0,0)$. If the inclusions are distributed on the unit square, this issue is studied recently when $\frac{ ilde d}a$ tends to zero or infinity, leaving aside the critical case where the volume fraction of the porous medium is below its possible maximal value but non-zero. In this paper, we provide the first result in this regime. In contrast with former results, we obtain an Euler type equation where a homogenized term appears in the elliptic problem relating the velocity and the vorticity. Our analysis is based on the so-called method of reflections whose convergence provides novel estimates on the solutions to the div-curl problem which is involved in the 2D-Euler equations.

연구 동기 및 목표

• 구멍 크기 a와 구멍 간격 ˜d가 모두 0으로 수렴할 때 2차원 비압축성 오일러 유동의 점근적 행동을 분석하는 것. 이때 ˜d/a → ¯k > 0 이다.
• 이전 연구들이 희박한 영역(˜d/a → ∞)과 조밀한 영역(˜d/a → 0)만 고려한 이론적 간극을 메우는 것.
• 다공성 매질이 작은 비율이지만 0이 아닌 체적 분율을 가지는 임계 영역에서 오일러 방정식의 균질화 시스템을 유도하는 것.
• 반사 방법을 사용하여 다공성 영역에서의 발산-회전 문제에 대한 새로운 추정을 확립함으로써, 균질화 극한으로의 수렴을 가능하게 하는 것.

제안 방법

• 분석은 두 단계로 나뉜다: 첫째, 속도와 비산란성 사이의 관계를 스트림 함수 ψN를 통해 기술하는 다공성 영역 FN 내에서의 발산-회전 문제에 대한 정밀한 균일 추정을 도출하는 것.
• 반사 방법을 사용하여 근사 해를 구성하고, 진짜 해 ψN과 균질화된 해 ψc 사이의 오차 한계를 도출하는 것.
• 균질화된 문제의 형태는 R²에서 div[(I₂ + kMK)∇ψc] = f 이며, 여기서 MK는 구멍 형상 K에 따라 결정되는 행렬이고, k는 체적 분율의 극한이다.
• 핵심적인 신기술은 가중치 노름과 비산란성 기울기의 로그 소보레프 유형의 추정을 사용하여, 진짜 속도 uN과 균질화된 속도 uc 사이의 차이에 대한 이방성 L∞ 추정을 유도하는 것이다.
• 비산란성 수송 방정식에 에너지 추정과 그론발라 타입의 추론을 적용하여 비산란성 및 속도 장의 수렴을 확립하는 것.
• 시간에 대한 해의 존재 시간을 균일하게 연장하기 위해 부스팅 기법을 사용하여, 고정된 시간 T까지의 수렴이 유지됨을 보장하는 것.

실험 결과

연구 질문

• RQ1 다공성 매질이 작은 비율이지만 0이 아닌 체적 분율을 가지는 경우, 즉 임계 영역 ˜d/a → ¯k > 0 에서 2차원 오일러 방정식의 균질화 극한은 무엇인가?
• RQ2 비영인 다공성 매질이 오일러 설정에서 거시적 유체 역학에 어떤 영향을 미치는가? 특히 효과적 점성도나 이방성의 측면에서 설명할 수 있는가?
• RQ3 반사 방법이 임계 영역에서 진짜 해와 균질화된 해 사이의 오차를 제어하기에 충분한 추정을 제공할 수 있는가?
• RQ4 균질화 극한에서 유체 운동을 지배하는 효과적 방정식의 구조는 무엇이며, 고전적 오일러 법칙이나 다르시 법칙과 어떻게 다를까?
• RQ5 다공성 영역 내의 비산란성 및 속도 장이 어떻게 그 균질화된 대응 장으로 수렴하는가? 이 수렴에 대한 균일한 추정은 무엇인가?

주요 결과

• 임계 영역에서 2차원 오일러 방정식의 균질화 극한은 수정된 타원형 문제로 기술된다: R²에서 div[(I₂ + kMK)∇ψc] = f이며, 여기서 MK는 구멍 형상 K에 따라 결정되는 행렬이다.
• 단위 원판의 경우 MK = 2I₂ 이므로, 효과적 연산자는 div[3I₂∇ψc] = f 가 되며, 이는 고전적 라플라스 연산자에 상당한 이방성 보정이 있음을 시사한다.
• 속도 장 uN이 균질화된 장 uc로 수렴하는 것은 L∞(OT)에서 이루어지며, 수렴 속도는 F(N,k) = (a/d)^{3−η} + ||µ−k||_{W^{−1,p}}^{(1−η)(p+2)/p} + ||µ−k||_{W^{−1,p}}^{1/2} + ||k||_{L∞}^2 으로 유계이다.
• 비산란성 기울기는 로그 그론발 추정을 통해 시간에 대해 균일하게 유계이므로, 오일러 방정식의 비선형 항을 제어하는 데 필수적이다.
• 다공성 영역 내 유체 입자의 궤적 XN(t,x)은 [0,T]에서 균질화된 궤적 Xc(t,x)로 균일하게 수렴하며, 오차는 O(F(N,k)) 이다.
• 초기 비산란성이 컴actsupport를 가지며 L∞ 노름이 유계인 경우에 결과가 성립하며, 체적 분율 k가 작고 FV(ε₀)에 속할 경우, 고정된 시간 T까지의 시간에 대해 수렴은 균일하게 유지된다.