[논문 리뷰] Implementation for blow up of tornado-type solutions for complex version of 3D Navier-Stokes system

Максим Дмитриевич Арнольд, А. В. Хохлов|ArXiv.org|2008. 06. 26.

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

한 줄 요약

이 논문은 복소수 3차원 나비에-스토크스 시스템에서 토네이도 유형의 해에 대한 자가유사 폭발 시나리오를 푸리에 변환된 방정식과 순수 허수의 초기 자료를 사용하여 수치적으로 실현한다. (0,0,R)에 중심을 둔 압축된 지원을 가진 스케일링된 에르미트 함수의 선형 조합으로 초기 조건을 구성함으로써, 저자들은 임계 시간 $T_{\text{cr}}$ 근처에서 다항식 폭발 행동을 보이는 에너지 증가를 시뮬레이션한다. 관측된 증가율은 약 15–20로, 시나리오의 이론적 $1/(T_{\text{cr}}-t)^5$ 점 渐진값에 가까워지며, 시나리오의 예측과 일치한다.

ABSTRACT

We consider Cauchy problem for Fourier transformation of 3-dimensional Navier-Stokes system with zero external force. Using initial data purposed by Dong Li and Ya.G.Sinai we implement self-similar regime producing fast growing behavior of the energy of solution while time tends to critical value.

연구 동기 및 목표

• 시나리오의 이론적 폭발 시나리오를 복소수 3차원 나비에-스토크스 시스템에서 실현하기 위해 시나리오의 제안자인 시나이와 리가 제안한 바를 수치적으로 실현한다.
• 스케일링된 에르미트 함수를 사용하여 (0,0,R)를 중심으로 한 3차원 푸리에 공간에서 컴팩트한 지원을 가진 초기 자료를 구성함으로써, 자가유사 폭발 역학의 수치 시뮬레이션을 가능하게 한다.
• 유한한 임계 시간 $T_{\text{cr}}$에 가까워질수록 해의 에너지 증가율을 조사하여, 시나리오의 이론적 $1/(T_{\text{cr}}-t)^5$ 점 渐진값을 검증한다.
• 해의 폭발 행동이 메esh 크기 및 초기 조건 형태의 미세한 변화에 대해 얼마나 견고한지 검증한다.

제안 방법

• 연구는 외부 힘을 제로로 설정하고 점성 계수를 1로 설정한 푸리에 변환된 3차원 나비에-스토크스 시스템을 사용하며, 식 (1)을 진폭 매개변수 $A$에 대한 멱급수 전개를 통해 적분 방정식을 해결한다.
• 초기 자료는 $v_A(k,0) = A c_0(k)$로 구성되며, $c_0(k)$는 (0,0,R)를 중심으로 반지름 $r$인 구 안에서 컴팩트한 지원을 가지며, 최대 차수 $D=3$까지의 스케일링된 에르미트 함수의 선형 조합을 사용한다.
• 해는 $nX=nY=160$, $nZ=1440$, $dt=0.001$의 3차원 정렬 격자에서 수치적으로 시뮬레이션되며, 3차원 컨볼루션 계산에 인텔 MKL을 사용한다.
• 에너지 $\mathcal{E}(t) = \int |v(k,t)|^2 dk$는 $L^2$ 노름의 통합된 제곱근으로 계산되며, 연속된 시간 단계에서 감시된다.
• 자기유사 행동 분석을 위해 $g_p(y,s) = h_p((0,0,R)+\sqrt{R}y,s)$를 통한 재규합화 방법을 사용하며, $p \to \infty$의 극한에서 재귀 관계의 불안정한 모드에 집중한다.
• 폭발의 안정성은 시간 간격, 공간 해상도, 초기 조건 형태를 최대 25%까지 변화시켜 테스트하며, 자연스러운 재스케일링 하에 일관된 폭발 행동을 관찰한다.

실험 결과

연구 질문

• RQ1 푸리에 변환된 방정식과 컴팩트한 지원을 가진 초기 자료를 사용하여 3차원 나비에-스토크스 시스템에서 토네이도 유형의 해에 대한 자가유사 폭발 시나리오의 수치적 실현이 가능할 수 있는가?
• RQ2 이러한 시뮬레이션에서 임계 시간 $T_{\text{cr}}$에 가까워질수록 관측된 에너지 증가율 $\mathcal{E}(t)$는 얼마인가?
• RQ3 메쉬 크기 및 초기 조건 형태의 미세한 변화에 대해 폭발 행동은 얼마나 견고한가?
• RQ4 수치 시뮬레이션은 시나리오의 이론적 $1/(T_{\text{cr}}-t)^5$ 점 渐진값과 일치하는 에너지 증가를 보이는가?
• RQ5 재귀 관계 $h_p(k,s)$의 불안정한 고유 모드는 폭발 역학의 발생에 어떤 역할을 하는가?

주요 결과

• 에너지 $\mathcal{M}(t) = \sqrt{\mathcal{E}(t)}$는 $t=0.000$에 1.00000에서 시작하여 $t=0.042$에 0.3664로 감소한 후, $t=0.045$에 19.972로 급격히 증가하고, $t=0.046$에 101.20, $t=0.047$에 698.57로 증가하여 $T_{\text{cr}}$ 근처에서 급격한 증가를 보임.
• 다른 메쉬 설정(공간 간격 0.20328, $dt=0.0008$)에서도 에너지는 $t=0.0424$에 0.3047에서 시작하여 $t=0.0480$에 1027.7로 증가하여 매개변수 변화에 대해 일관된 폭발 행동을 보임.
• 관측된 에너지 증가율은 15에서 20 사이로 추정되며, 이는 이론적 $1/(T_{\text{cr}}-t)^5$ 비율보다는 높지만, 계산 제약 조건을 감안할 때 예상 범위 내에 있음.
• 시간 간격, 공간 해상도, 초기 조건 형태에 대해 최대 25%의 변화가 가해져도 폭발 행동은 안정되었으며, 자연스러운 재스케일링 하에 해의 구조가 변형됨.
• 수치적 해는 푸리에 공간에서 뚜렷한 '구름'을 보였으며, 첫 번째 구름은 $z \sim 20$에 중심을 두었고, 이후 구름들은 $z \sim 40$ 및 $z \sim 60$에 위치하여, $h_p$의 지지 집합의 자가유사 스케일링 $O(Rp)$와 일치함.
• 결과는 유한 시간 내에 폭발하는 데에 기여하는 10개의 매개변수를 가진 초기 자료의 가족 존재를 지지하며, $t \to T_{\text{cr}}$에 따라 에너지가 발산함을 보여주지만, 해상도 제약으로 인해 정확한 $1/(T_{\text{cr}}-t)^5$ 비율은 달성되지 않음.