2025-11-10T03:08:08.537072

Global stability for compressible isentropic Navier-Stokes equations in 3D bounded domains with Navier-slip boundary conditions

Liu, Wu, Zhong

We investigate the global stability of large solutions to the compressible isentropic Navier-Stokes equations in a three-dimensional (3D) bounded domain with Navier-slip boundary conditions. It is shown that the strong solutions converge to an equilibrium state exponentially in the $L^2$-norm provided the density is essentially uniform-in-time bounded from above. Moreover, we obtain that the density converges to its equilibrium state exponentially in the $L^\infty$-norm if additionally the initial density is bounded away from zero. Furthermore, we derive that the vacuum states will not vanish for any time provided vacuum appears (even at a point) initially. This is the first result concerning the global stability for large strong solutions of compressible Navier-Stokes equations with vacuum in 3D general bounded domains.

academic

3차원 유계 영역에서 Navier-슬립 경계조건을 갖는 압축성 등엔트로피 Navier-Stokes 방정식의 전역 안정성

기본 정보

논문 ID: 2504.17136
제목: Global stability for compressible isentropic Navier-Stokes equations in 3D bounded domains with Navier-slip boundary conditions
저자: Yang Liu¹, Guochun Wu², Xin Zhong³
- ¹장춘사범대학교 수학학원
- ²샤먼이공대학교 수학통계학원
- ³남서대학교 수학통계학원
분류: math.AP (수학 분석 - 편미분방정식)
제출 시간: 2025년 4월
논문 링크: https://arxiv.org/abs/2504.17136v4

초록

본 논문은 3차원 유계 영역에서 Navier-슬립 경계조건을 갖는 압축성 등엔트로피 Navier-Stokes 방정식의 대해(large solution)의 전역 안정성을 연구한다. 밀도가 본질적으로 균일하게 유계인 조건 하에서, 강해가 $L^2$ 노름 의미에서 평형상태로 지수적으로 수렴함을 증명했다. 더욱이, 초기 밀도가 양의 하한을 가질 때, 밀도는 $L^∞$ 노름에서도 평형상태로 지수적으로 수렴한다. 또한 초기 시점에 진공(비록 한 점에서만이라도)이 존재하면, 진공 상태는 임의의 시간에 소멸하지 않음을 증명했다. 이는 3차원 일반 유계 영역에서 진공을 포함하는 압축성 Navier-Stokes 방정식의 대강해의 전역 안정성에 관한 첫 번째 결과이다.

연구 배경 및 동기

문제의 중요성

압축성 Navier-Stokes 방정식은 점성 압축성 유체의 운동을 기술하며, 유체역학의 기본 방정식계이다. 그 해의 장시간 거동과 안정성을 연구하는 것은 유체 운동 법칙을 이해하는 데 중요한 의미를 갖는다.

기존 연구의 한계

경계조건 제약: 이전 연구는 주로 전체 공간 $\mathbb{R}^3$ 또는 원환면 $T^3$ 에 집중되었으며, 유계 영역에 대한 연구는 상대적으로 부족하다
경계조건 유형: 대부분의 연구는 무슬립 경계조건을 고려했으며, Navier-슬립 경계조건에 대한 연구는 상대적으로 미흡하다
해의 크기 제약: 많은 결과가 초기 데이터가 작을 것을 요구하며, 대해의 안정성 연구가 충분하지 않다
진공 상태: 진공을 포함하는 경우는 수학적 분석의 복잡성을 증가시킨다

연구 동기

Navier-슬립 경계조건은 물리적으로 거친 경계, 다공성 경계 등 특정 실제 상황에 더 부합한다
3차원 유계 영역의 기하학적 복잡성은 새로운 수학적 기법을 필요로 한다
대해의 전역 안정성 분석은 중요한 이론적 가치를 갖는다

핵심 기여

최초 결과: 3차원 일반 유계 영역에서 진공을 포함하는 압축성 Navier-Stokes 방정식의 대강해 전역 안정성에 관한 최초 결과
지수 수렴성: 적절한 조건 하에서 해의 지수 수렴성을 증명하고 구체적인 감쇠율을 제시
진공 지속성: 진공 상태의 지속성을 증명 - 일단 진공이 나타나면 소멸하지 않음
기술적 혁신: Navier-슬립 경계조건을 처리하는 새로운 기술 방법 개발

방법론 상세 설명

문제 정의

3차원 유계 영역 $\Omega \subset \mathbb{R}^3$ 에서의 압축성 등엔트로피 Navier-Stokes 방정식계를 고려한다:

$\begin{cases} \rho_t + \text{div}(\rho u) = 0 \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \mu \Delta u + (\mu + \lambda)\nabla \text{div} u \end{cases}$

여기서 $P = \rho^\gamma$ ( $\gamma > 1$ )이며, Navier-슬립 경계조건을 갖는다: $u \cdot n = 0, \quad \text{curl} u \times n = 0 \quad \text{on} \partial\Omega$

강해의 정의

$T > 0$ 과 어떤 $q \in (3,6]$ 에 대해, $(ρ,u)$ 가 강해라 함은 다음을 만족할 때이다:

$0 \leq \rho \in C([0,T]; W^{1,q})$ , $\rho_t \in C([0,T]; L^q)$
$u \in C([0,T]; H^2) \cap L^2(0,T; W^{2,q})$
$\sqrt{\rho}u_t \in L^∞(0,T; L^2)$ , $u_t \in L^2(0,T; H^1)$

주요 기술 방법

1. 에너지-소산 부등식

Lyapunov 범함수 구성: $M_1(t) = D_1 \int \left(\frac{1}{2}\rho|u|^2 + G(\rho)\right)dx - \int \rho u \cdot B[\rho - \bar{\rho}_0]dx$

여기서 $G(\rho)$ 는 상대 엔트로피이고, $B$ 는 Bogovskii 연산자이다.

2. 유효 점성 플럭스 분석

유효 점성 플럭스 정의: $F = (2\mu + \lambda)\text{div} u - (P - \bar{P})$

이것이 만족하는 타원 방정식 활용: $\begin{cases} \Delta F = \text{div}(\rho \dot{u}) & \text{in} \Omega \\ \frac{\partial F}{\partial n} = \rho \dot{u} \cdot n & \text{on} \partial\Omega \end{cases}$

3. 와도 추정

와도 $\omega = \nabla \times u$ 에 대해, 경계조건 $(\omega \times n)|_{\partial\Omega} = 0$ 과 div-curl 시스템 이론을 활용하여: $\|\nabla u\|_{L^q} \leq C(\|\text{div} u\|_{L^q} + \|\omega\|_{L^q})$

기술적 혁신점

경계 적분 처리: 항등식 $u \cdot \nabla u \cdot n = -u \cdot \nabla n \cdot u$ 를 통해 경계의 비선형 항 처리
계층적 에너지 추정: 3단계 에너지 추정 전략 채택:
- 1단계: $\|(\rho - \bar{\rho}_0, \sqrt{\rho}u)\|_{L^2}$ 의 지수 감쇠
- 2단계: $\|\nabla u\|_{L^2}$ 의 지수 감쇠
- 3단계: $\|\sqrt{\rho}\dot{u}\|_{L^2}$ 의 지수 감쇠
입자 궤적 방법: 특성선 방법을 이용한 밀도 하한 분석

주요 결과

정리 1.1 (전역 안정성)

초기 데이터가 $K = \|\rho_0 - \bar{\rho}_0\|_{L^2} + \|\sqrt{\rho_0}u_0\|_{L^2} + \|\nabla u_0\|_{H^1} < \infty$ 를 만족하고, 전역 강해가 존재하여 $\sup_{t \geq 0}\|\rho(\cdot,t)\|_{L^∞} \leq \hat{\rho}$ 를 만족한다고 가정하자. 그러면 양의 상수 $C_1, \eta_1$ 이 존재하여:

$\|(\rho - \bar{\rho}_0, \sqrt{\rho}u, \nabla u, \sqrt{\rho}\dot{u})(\cdot,t)\|_{L^2} \leq C_1 e^{-\eta_1 t}$

추가로 $\inf_{x \in \Omega}\rho_0(x) \geq \rho_* > 0$ 이면, 양의 상수 $C_2, \eta_2$ 가 존재하여: $\|\rho - \bar{\rho}_0\|_{L^∞} \leq C_2 e^{-\eta_2 t}$

정리 1.2 (진공 지속성)

정리 1.1의 조건 하에서, $\inf_{x \in \Omega}\rho_0(x) = 0$ 이면, 임의의 $t \geq 0$ 에 대해: $\inf_{x \in \Omega}\rho(x,t) = 0$

증명 전략

1단계: 기본 에너지 추정

승수법과 Bogovskii 연산자를 통해 에너지-소산 부등식 수립: $\frac{d}{dt}M_1(t) + \frac{M_1(t)}{D_1} + \frac{\|\text{div} u\|_{L^2}^2 + \|\text{curl} u\|_{L^2}^2}{D_1} \leq 0$

2단계: 기울기 추정

동량 방정식과 유효 점성 플럭스의 성질을 이용하여: $\frac{d}{dt}\Psi(t) + \frac{1}{2}\|\sqrt{\rho}\dot{u}\|_{L^2}^2 \leq C(\|\nabla u\|_{L^2}^2 + \|\nabla u\|_{L^2}^6 + \|\rho - \bar{\rho}_0\|_{L^2}^2)$

3단계: 고차 추정

개선된 에너지 추정을 통해 $\|\sqrt{\rho}\dot{u}\|_{L^2}$ 의 시간 감쇠 제어.

4단계: $L^∞$ 추정

특성선 방법과 밀도의 하한 추정을 이용하여 $L^∞$ 노름의 지수 감쇠 증명.

기술적 난제 및 혁신

주요 기술적 난제

경계항 처리: Navier-슬립 경계조건은 경계에서 속도가 0이 아니므로 경계 적분의 특수 처리 필요
비선형 추정: $\int_{\partial\Omega} F_t(u \cdot \nabla n \cdot u)dS$ 등 경계 비선형항 처리
진공 영역: 밀도가 0일 수 있어 분석의 복잡성 증가