[논문 리뷰] Potential Singularity Formation of 3D Axisymmetric Navier-Stokes Equations with Degenerate Variable Diffusion Coefficients.

Thomas Y. Hou, De Huang|arXiv (Cornell University)|2021. 02. 12.

Navier-Stokes equation solutions참고 문헌 31인용 수 5

한 줄 요약

이 논문은 축대칭 3D 라우즈-스토크스 방정식에서 비정상적인, 변화하는 확산 계수를 가질 경우, 원점에서 유한 시간 내에 국소적으로 자가유사적인 특이점이 발생할 수 있다는 수치적 증거를 제시한다. 이 특이점은 반대칭 복소소용돌이가 두 스케일의 이동파를 생성함으로써 유도되며, 비선형적 복소소용돌이 강도 증가와 관련된 것으로 나타났다. 특이점은 비선형적 정렬 현상에 의해 발생하며, 유한 에너지를 지닌 매끄러운 초기 조건에서 양성 피드백 루프에 의해 안정화된다.

ABSTRACT

In this paper, we present strong numerical evidences that the $3$D incompressible axisymmetric Navier-Stokes equations with degenerate diffusion coefficients and smooth initial data of finite energy develop a potential finite time locally self-similar singularity at the origin. The spatial part of the degenerate diffusion coefficient is a smooth function of $r$ and $z$ independent of the solution and vanishes like $O(r^2)+O(z^2)$ near the origin. This potential singularity is induced by a potential singularity of the $3$D Euler equations. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the property that the center of the traveling wave approaches the origin at a slower rate than the rate of the collapse of the singularity. The driving mechanism for this potential singularity is due to two antisymmetric vortex dipoles that generate a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop an additional small scale characterizing the thickness of the sharp front. On the other hand, the Navier-Stokes equations with a constant diffusion coefficient regularize the two-scale solution structure and do not develop a finite time singularity for the same initial data. The initial condition is designed in such a way that it generates a positive feedback loop that enforces a strong nonlinear alignment of vortex stretching, leading to a stable locally self-similar blowup at the origin. We perform careful resolution study and asymptotic scaling analysis to provide further support of the potential finite time locally self-similar blowup.

연구 동기 및 목표

• 비정상적인 확산 계수를 가진 3D 축대칭 라우즈-스토크스 방정식이 매끄러운 초기 조건 하에서도 유한 시간 내 특이점이 발생할 수 있는지 조사하기.
• 특히 원점 부근에서 O(r²) + O(z²)의 형태로 사라지는 변수 확산 계수의 역할이 특이점 형성에 어떻게 기여하는지 이해하기.
• 반대칭 복소소용돌이가 비선형 정렬을 유도하고 잠재적 폭발을 이끄는 메커니즘을 탐구하기.
• 일정한 확산 계수를 가진 경우 정규화가 발생하여 유한 시간 폭발이 방지되는 것과의 대비를 통해 특이점 발생 메커니즘을 분석하기.
• 해상도 분석과 점점 가까워지는 스케일 분석을 통해 안정적이고 국소적으로 자가유사적인 폭발의 존재를 검증하기.

제안 방법

• 원점 부근에서 O(r²) + O(z²)의 형태로 사라지는 공간에 따라 변화하는 확산 계수를 가진 3D 비압축성 축대칭 라우즈-스토크스 방정식의 수치적 시뮬레이션.
• 복소소용돌이 강도 증가 역학에서 양성 피드백 루프를 유도할 수 있도록 설계된 유한 에너지를 지닌 매끄러운 초기 조건 사용.
• 원점에 가까워지는 두 스케일의 이동파 구조를 해상도를 높여 정확히 해석하기 위한 고해상도 수치 기법 구현.
• 잠재적 특이점의 국소적으로 자가유사성과 폭발 행동을 확인하기 위해 점점 가까워지는 스케일 분석 수행.
• 특이점 발생 메커니즘에서 점성 정규화의 역할을 분리하기 위해 3D 오일러 방정식과 결과 비교.
• 반경 방향 및 축 방향 속도장의 진화를 추적하여 반대칭 복소소용돌이에 의해 생성된 강력한 비틀림 층의 형성을 확인하기.

실험 결과

연구 질문

• RQ1 비정상적인 확산 계수를 가진 3D 축대칭 라우즈-스토크스 방정식이 매끄러운 초기 조건 하에서도 유한 시간 내에 국소적으로 자가유사적인 특이점이 발생할 수 있는가?
• RQ2 두 스케일의 이동파 구조가 이러한 특이점 형성에 어떤 역할을 하는가? 그리고 이 파동의 중심이 원점에 접근하는 속도와 특이점의 붕괴 속도는 어떻게 비교되는가?
• RQ3 반대칭 복소소용돌이는 어떻게 비선형 정렬을 유도하고 폭발 메커니즘을 유지하는 강력한 비틀림 층을 생성하는가?
• RQ4 동일한 초기 조건에서 일정한 확산 계수를 사용할 경우 왜 특이점이 발생하지 않으며, 이 경우 점성 정규화의 역할은 무엇인가?
• RQ5 점점 가까워지는 스케일 분석과 해상도 연구가 안정적이고 자가유사적인 폭발 프로파일의 존재를 어느 정도 뒷받침하는가?

주요 결과

• 비정상적인 확산 계수를 가진 3D 축대칭 라우즈-스토크스 방정식은 고해상도 시뮬레이션을 통해 원점에서 잠재적인 유한 시간 특이점이 발생함을 확인하였다.
• 특이점은 국소적으로 자가유사적이며, 중심이 원점에 더 천천히 접근하는 두 스케일의 이동파에 의해 유도된다.
• 초기 조건에 의해 강제되는 복소소용돌이 강도 증가의 양성 피드백 루프는 폭발 과정을 안정화시키고 비선형 정렬을 유지한다.
• 특이점 메커니즘은 반대칭 복소소용돌이에 의해 반경 및 축 방향 속도장에서 강력한 비틀림 층이 생성됨에 기반한다.
• 일정한 확산 계수를 가진 라우즈-스토크스 방정식은 해를 정규화하여 폭발을 방지하는 데 반해, 비정상적인 경우 특이점 형성이 가능하다.
• 점점 가까워지는 스케일 분석과 해상도 연구는 잠재적 특이점의 자가유사성과 안정성을 강력히 뒷받침한다.