[논문 리뷰] Homogenization of the Neumann problem for higher-order elliptic equations with periodic coefficients

T. A. Suslina|arXiv (Cornell University)|2017. 05. 21.

Advanced Mathematical Modeling in Engineering참고 문헌 16인용 수 4

한 줄 요약

이 논문은 주기적 계수를 가진 고계 타원형 뉴먼 문제의 균질화에 대해 날카운 연산자 오차 추정을 수립한다. $ L_2 $ 및 $ H^p $ 노름에서 $ (A_{N,\rho}-\zeta I)^{-1} $의 근사값을 도출하며, 오차 한계가 $ \varepsilon^{1/2}\rho_\flat(\zeta)^{1/2} + \varepsilon|\zeta+1|\rho_\flat(\zeta) $ 순서임을 보이며, 최적의 정규성 및 질량 조건 하에서 보정자 항 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $ 를 포함한다.

ABSTRACT

Let $\mathcal{O}\subset\mathbb{R}^d$ be a bounded domain of class $C^{2p}$. In $L_2(\mathcal{O};\mathbb{C}^n)$, we study a selfadjoint strongly elliptic operator $A_{N,\varepsilon}$ of order $2p$ given by the expression $b({\mathbf D})^* g({\mathbf x}/\varepsilon) b({\mathbf D})$, $\varepsilon >0$, with the Neumann boundary conditions. Here $g({\mathbf x})$ is a bounded and positive definite $(m imes m)$-matrix-valued function in ${\mathbb R}^d$, periodic with respect to some lattice; $b({\mathbf D})=\sum_{|α|=p} b_α{\mathbf D}^α$ is a differential operator of order $p$ with constant coefficients; $b_α$ are constant $(m imes n)$-matrices. It is assumed that $m\geqslant n$ and that the symbol $b({\boldsymbol ξ})$ has maximal rank for any $0 e {\boldsymbol ξ}\in {\mathbb C}^d$. We find approximations for the resolvent $\left(A_{N,\varepsilon}-ζI ight)^{-1}$ in the $L_2(\mathcal{O};\mathbb{C}^n)$-operator norm and in the norm of operators acting from $L_2(\mathcal{O};\mathbb{C}^n)$ to the Sobolev space $H^p(\mathcal{O};\mathbb{C}^n)$, with error estimates depending on $\varepsilon$ and $ζ$.

연구 동기 및 목표

• 유계 $ C^{2p} $ 영역에서 주기적 계수를 가진 자기수반, 강한 타원형 고계 미분 연산자($ 2p $차수)의 균질화를 분석하는 것.
• $ \varepsilon $가 $ 0 $으로 수렴할 때 $ L_2 $ 및 $ H^p $ 공간에서 $ (A_{N,\rho} - \zeta I)^{-1} $의 연산자 노름 오차 추정을 도출하는 것, 여기서 $ \varepsilon $는 주기성의 척도이다.
• $ H^p $ 노름에서 리솔베ント 근사의 개선을 위해 보정자 항 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $ 를 구성하는 것.
• 스펙트럼 매개변수 $ \zeta \in \mathbb{C} \setminus [c_\flat, \infty) $ 에 따라 오차가 어떻게 의존하는지 명시적으로 기술하는 것. 이는 저주파수 및 고주파수 행동을 모두 포괄한다.

제안 방법

• 스케일링 변환, 플루에트-블로흐 이론, 해석적 섭동 이론을 기반으로 한 연산자 이론적 접근을 사용한다.
• 미분 연산자는 $ A_{N,\varepsilon} = b(\mathbf{D})^* g(\mathbf{x}/\varepsilon) b(\mathbf{D}) $ 로 인수분해되며, 여기서 $ b(\mathbf{D}) $ 는 최대 질량 기호를 가진 $ p $차수의 일정 계수 연산자이다.
• 효과적 연산자 $ A_N^0 $ 는 $ g(\mathbf{x}) $ 의 주기적 구조로부터 유도된 균질화된 행렬 $ g^0 $ 을 통해 구성된다.
• 스펙트럼 분석을 통해 푸리에 도메인에서 오차 추정을 유도하며, 리솔베ント 항등식과 스펙트럼 프로젝터 및 보정자 연산자의 추정을 사용한다.
• 보정자 $ K_N^0(\zeta;\varepsilon) $ 는 명시적으로 정의되며, $ \varepsilon $ 척도에서 해의 진동 행동을 포착함을 보이며, $ \|K_N^0(\zeta;\varepsilon)\|_{L_2 \to H^p} = O(\varepsilon^{-p}) $ 임을 보인다.
• 추정은 두 가지 연산자 노름에서 증명된다: $ L_2 \to L_2 $ 및 $ L_2 \to H^p $, 추정치는 $ \varepsilon $, $ \zeta $, 그리고 스펙트럼 함수 $ \rho_\flat(\zeta) $ 에 따라 달라진다.

실험 결과

연구 질문

• RQ1 주기성 척도 $ \varepsilon \to 0 $ 일 때, $ (A_{N,\varepsilon} - \zeta I)^{-1} $ 이 $ (A_N^0 - \zeta I)^{-1} $ 으로 수렴하는 속도는 $ L_2 $-연산자 노름에서 어떻게 되는가?
• RQ2 보정자 항을 $ H^p $-노름에서 사용하여 주로 근사값을 어떻게 개선할 수 있는가?
• RQ3 오차는 복소수 평면에서 $ \zeta \in \mathbb{C} \setminus [c_\flat, \infty) $ 에 따라 어떻게 의존하는가?
• RQ4 보정자 항이 사라지거나 단순화되는 조건은 무엇이며, 이는 오차 추정에 어떻게 영향을 미치는가?
• RQ5 이전의 2차 연산자 결과를 확장하여, 주기적 계수를 가진 고계 타원형 뉴먼 문제에 대해 날카운 연산자 오차 추정을 유도할 수 있는가?

주요 결과

• 리솔베ント $ (A_{N,\varepsilon} - \zeta I)^{-1} $ 는 $ L_2(\mathcal{O};\mathbb{C}^n) \to L_2(\mathcal{O};\mathbb{C}^n) $ 노름에서 오차 $ \leq \mathfrak{C}_1 \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) $ 로 수렴하며, $ 0 < \varepsilon \leq \varepsilon_1 $ 이다.
• $ L_2(\mathcal{O};\mathbb{C}^n) \to H^p(\mathcal{O};\mathbb{C}^n) $ 노름에서 오차는 $ \widetilde{\mathfrak{C}}_{10} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) $ 로 유계이며, 보정자 항 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $ 가 포함되어 있다.
• $ g^0 = \overline{g} $ 일 때, 보정자가 사라지고 오차 추정은 $ \| \mathbf{u}_\varepsilon - \mathbf{u}_0 \|_{H^p} \leq \mathfrak{C}_{10} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $ 로 단순화된다.
• 유속 $ \mathbf{p}_\varepsilon = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon $ 에 대해 $ L_2 $-노름에서 오차는 $ \widetilde{\mathfrak{C}}_{11} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $ 로 유계이다.
• $ g^0 = \underline{g} $ 일 때, 오차 추정은 $ \| \mathbf{u}_\varepsilon - \mathbf{v}_\varepsilon^0 \|_{H^p} \leq \widehat{\mathfrak{C}}_{10} \varepsilon \left( \rho_\flat(\zeta)^{1/2} + |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $ 로 변형되며, $ \varepsilon $-의존성이 향상됨을 보여준다.
• 추정치는 날카롭고 스펙트럼 매개변수 $ \zeta $ 에 따라 명시적으로 의존하며, $ \rho_\flat(\zeta) $ 는 리솔베ント의 스펙트럼 갭 및 감쇠 성질을 캐릭터라이즈한다.