[논문 리뷰] Homogenization of hyperbolic equations with periodic coefficients

Mark Dorodnyi, T. A. Suslina|arXiv (Cornell University)|2016. 06. 19.

Advanced Mathematical Modeling in Engineering참고 문헌 11인용 수 3

한 줄 요약

이 논문은 $L_2(\mathbb{R}^d;\mathbb{C}^n)$에서 자기수반인 두 번째 차수의 주기적 미분 연산자에 대해, 코시안의 제곱근을 취한 연산자에 대한 연산자 오차 추정을 수립한다. 초점은 계수가 진동하는 특성을 가진 쌍곡형 방정식에 맞춰져 있다. 스펙트럼 이론과 플로케-블로흐 분해를 활용하여, $\cos(\mathcal{A}_\varepsilon^{1/2}\tau)$에 대한 $(H^s \to L_2)$-노름 근사치를 $O(\varepsilon)$ 오차로 도출하였으며, 이는 음향학 및 탄성계에 적용 가능하다.

ABSTRACT

In $L_2(\mathbb{R}^d;\mathbb{C}^n)$ we consider selfadjoint strongly elliptic second order differential operators ${\mathcal A}_\varepsilon$ with periodic coefficients depending on ${\mathbf x}/ \varepsilon$, $\varepsilon>0$. We study the behavior of the operator cosine $\cos( {\mathcal A}^{1/2}_\varepsilon τ)$, $τ\in \mathbb{R}$, for small $\varepsilon$. Approximations for this operator in the $(H^s o L_2)$-operator norm with a suitable $s$ are obtained. The results are used to study the behavior of the solution ${\mathbf v}_\varepsilon$ of the Cauchy problem for the hyperbolic equation $\partial^2_τ{\mathbf v}_\varepsilon = - \mathcal{A}_\varepsilon {\mathbf v}_\varepsilon +\mathbf{F}$. General results are applied to the acoustics equation and the system of elasticity theory.

연구 동기 및 목표

• 주기적 계수를 가진 쌍곡형 방정식에 대해 $\varepsilon \to 0$일 때 연산자余弦 $\cos(\mathcal{A}_\varepsilon^{1/2}\tau)$의 점근적 행동을 분석하는 것.
• 원래 연산자 $\cos(\mathcal{A}_\varepsilon^{1/2}\tau)$의 동질화된 대응체에 대한 $(H^s \to L_2)$-위상에서의 연산자 노름 오차 추정을 유도하는 것.
• 초기값 문제 $\partial_\tau^2 \mathbf{v}_\varepsilon = -\mathcal{A}_\varepsilon \mathbf{v}_\varepsilon + \mathbf{F}$에 대해 해의 수렴성을 보장하는 것.
• 음향 방정식과 선형 탄성계와 같은 물리적 모델에 이론을 적용하여 검증하는 것.
• 관련된 플로케-블로흐 연산자의 스펙트럼 분석을 통해 오차 추정의 정밀도를 확립하는 것.

제안 방법

• quasi-모멘타 $\boldsymbol{\theta}$에 의해 매개화된 파형 연산자 가족으로 $\mathcal{A}_\varepsilon$를 분해하기 위해 플로케-블로흐 이론에 기반한 스펙트럼 접근법을 사용한다.
• 연산자를 동질화에 적합한 형태로 표현하기 위해 $\mathcal{A}_\varepsilon = f^\varepsilon^* b(\mathbf{D})^* g^\varepsilon b(\mathbf{D}) f^\varepsilon$의 인수분해를 활용한다.
• 해석적 섭동 이론을 적용하여 해석적 연산자와 반군의 수렴을 연구하고, 함수 해석학을 통해 코시안 연산자로 확장한다.
• 원래 연산자와 효과적 연산자 간의 차이를 $L_2$ 및 소볼레프 노름에서 분석함으로써 오차 추정을 도출하며, 일阶 보정항을 포함한다.
• 상수 효과 행렬 $g^0$를 체적 평균을 통해 계산한 효과적 연산자 $\mathcal{A}^0 = b(\mathbf{D})^* g^0 b(\mathbf{D})$를 사용한다.
• 블로흐 파형 연산자 $\widehat{N}(\boldsymbol{\theta})$의 스펙트럼을 분석하고, 임계점 $\boldsymbol{\theta}^{(j)}$에서의 비영인 고유값을 검증함으로써 추정의 정밀도를 확립한다.

실험 결과

연구 질문

• RQ1 주기적 계수를 가진 쌍곡형 방정식에서 $\varepsilon \to 0$일 때, $(H^s \to L_2)$-연산자 노름에서 $\cos(\mathcal{A}_\varepsilon^{1/2}\tau)$의 수렴 속도는 어떠한가?
• RQ2 소볼레프 공간 $H^s$에서의 초기 데이터의 정(regularity)은 쌍곡형余弦 연산자에 대한 오차 추정에 어떻게 영향을 미치는가?
• RQ3 타원형 및 방산형 문제에 대해 개발된 연산자 이론 기반의 동질화 프레임워크를 주기적 계수를 가진 쌍곡형 방정식으로 확장할 수 있는가?
• RQ4 $(H^s \to L_2)$-노름에서 쌍곡형余弦 연산자에 대한 $O(\varepsilon)$ 오차 추정의 정밀도는 어떠한가?
• RQ5 이러한 결과는 주기적 매질에서의 음향 방정식과 선형 탄성계와 같은 구체적인 물리계에 어떻게 적용되는가?

주요 결과

• $\cos(\mathcal{A}_\varepsilon^{1/2}\tau)$는 $0 \leq s \leq 3/2$일 때 $(H^s \to L_2)$-연산자 노름에서 오차 $O(\varepsilon)$로 $\cos(\mathcal{A}^{0,1/2}\tau)$로 수렴한다.
• 초기값 문제 $\partial_\tau^2 \mathbf{v}_\varepsilon = -\mathcal{A}_\varepsilon \mathbf{v}_\varepsilon + \mathbf{F}$에 대해, 해 $\mathbf{v}_\varepsilon$는 오차 $\|\mathbf{v}_\varepsilon - \mathbf{v}_0\|_{L_2} \leq \varepsilon^{2s/3} \left( \widehat{\mathfrak{C}}_3(s;\tau) \|\boldsymbol{\phi}\|_{H^s} + \widehat{\mathfrak{C}}_4(s;\tau) \|\boldsymbol{\psi}\|_{H^s} \right)$로 동질화된 해 $\mathbf{v}_0$로 수렴한다.
• 오차 추정은 정밀도가 확보되어 있다: 블로흐 파형 연산자 $\widehat{N}(\boldsymbol{\theta}^{(j)})$의 비영인 고유값 $\pm \widehat{\mu} \approx \pm 0.0985$는 $O(\varepsilon)$ 수렴 속도를 향상시킬 수 없음을 확인한다.
• 힐 몸체(일관된 전단 모odulus를 가진 등방성 탄성체)의 경우, 효과 행렬 $g^0_{\wedge}$는 평균 $\underline{g_{\wedge}}$와 일치하며, $\widehat{N}(\boldsymbol{\theta}) = 0$일 때 $L_2$-노름에서 오차 추정이 $O(\varepsilon^2)$로 향상된다.
• 근사치에 포함된 보정항은 주기적 경계값 문제의 해 $\Lambda_{22}$로부터 유도되며, $\arctan$과 복소 로그 항을 포함한 명시적 표현을 가진다.
• 결과는 음향 방정식과 선형 탄성계로 확장되었으며, 여기서 인수분해 $\mathcal{A}_\varepsilon = b(\mathbf{D})^* g^\varepsilon b(\mathbf{D})$가 직접 적용 가능하고, 동질화된 연산자는 명시적으로 계산 가능하다.