[논문 리뷰] Observable Divergence Theorem: Evolution Equations for Inviscid Regularization of Shocks and Turbulence

Kamran Mohseni|arXiv (Cornell University)|2010. 10. 13.

Cosmology and Gravitation Theories인용 수 4

한 줄 요약

이 논문은 사용자가 정의한 척도 $\alpha$에서 공간 평균된 유량을 사용하여 발산 정리를 재정의하는 새로운 프레임워크인 관측 가능한 발산 정리(Observable Divergence Theorem)를 소개한다. 이는 충격파와 난류에 대한 정규화된, 점성 없는 진화 방정식 유도를 가능하게 하며, 미분 방정식을 유도하기 이전에 이 정리를 적용함으로써 점성 없는 레일리 모델을 도출하고, 점성 항이나 수치적 소산 없이도 안정적이고 진동 없는 충격파와 난류 시뮬레이션을 가능하게 한다.

ABSTRACT

The divergence theorem of Gauss plays a central role in the derivation of the governing differential equations in fluid dynamics, electrodynamics, gravitational fields, and optics. One is often interested in an evolution equation for the large scale quantities without resolving the details of the small scales. As a result, there has been a significant effort in developing time-averaged and spatially-filtered equations for large scale dynamics from the fully resolved governing differential equations. One should realize that by starting from these fully-resolved equations (e.g. the Euler or Navier-Stokes equations) to derive an averaged evolution equation one has already taken the limit of the wave-numbers approaching infinity with no regards to our observational abilities at such a limit. As a result, obtaining the evolution equations for large scale quantities (low wave-numbers) by an averaging or filtering process is done after the fact. This could explain many of the theoretical and computational difficulties with the Euler or Navier-Stokes equations. Here, a rather different approach is proposed. The averaging process in implemented before the derivation of the differential form of the transport equations. A new observable divergence concept is defined based on fluxes calculated from observable quantities at a desired averaging scale, $α$. An observable divergence theorem is then proved and applied in the derivation of the observable and regularized transport equations. We further show that the application of the observable divergence theorem to incompressible flows results in a formal derivation of the inviscid Leray turbulence model first proposed in 1934. It is argued that such a methodology in deriving fluid evolution equations removes many of the theoretical and computational difficulties in multi-scale problems such as turbulence and shocks.

연구 동기 및 목표

• 고파르모드 불안정성의 이론적 및 계산적 과제를 고전적 오일러 및 레이놀즈 평균 나비에-스토크스 방정식을 사용하여 충격파와 난류에서 해결하기 위해.
• 후처리 평균화나 필터링의 한계를 극복하기 위해, 보존 법칙의 기본 유도 과정에 정규화를 통합하기 위해.
• 새로운 관측 가능한 발산 개념을 사용하여, 제1원리에서 점성 없는 레일리 난류 모델을 공식적으로 도출하기 위해.
• 인위적 점성 또는 수치적 소산에 의존하지 않고도 다중 스케일 유체역학에서 수치적 안정성과 물리적 일관성을 보장하는 프레임워크를 개발하기 위해.

제안 방법

• 특성 길이 척도 $\alpha$를 갖는 정규화되며, 음이 아닌, 반사 대칭이며 감소하는 척도 의존성 평균 커널 $g^\alpha$를 정의한다.
• 관측 가능한 발산 정리를 도입하여, 고전적 발산 정리에서 점별 유량을 제어체적 크기 $\alpha$ 내에서 공간 평균된 장 양에 기반한 유량으로 대체한다.
• 보존 법칙의 적분형을 적용하여 관측 가능한 형태의 보존 방정식(질량, 운동량, 에너지)을 도출한다.
• 관측 가능한 발산 정리를 사용하여 관측 가능한 비점성 오일러 방정식을 공식적으로 도출하며, 이는 한계 $\alpha \to 0$에서 점성 없는 레일리 모델으로 축소된다.
• 압축성 유동에 대해, 대류 항에 평균된 속도 및 밀도 장을 사용하여 관측 가능한 오일러 방정식을 유도한다.
• 푸리에 공간에서의 가짜 스펙트럼 방법과 룬게-쿠타-페르블베르(Runge-Kutta-Fehlberg, RK45) 시간 적분을 사용하여 수치 시뮬레이션을 구현하고, 기준 충격파 및 난류 문제에서 모델을 검증한다.

실험 결과

연구 질문

• RQ1 미분 방정식 유도 이전에 발산 정리에 평균을 통합함으로써, 정규화된 유체 방정식의 공식적 유도가 가능할 수 있는가?
• RQ2 관측 가능한 발산 정리는 점성 정규화 없이도 충격파와 난류에 대해 일관되고 안정적인 진화 방정식을 도출하는가?
• RQ3 유도된 관측 가능한 오일러 방정식이 엔트로피 해와 이동 충격파와 같은 알려진 해를 재현할 수 있는가?
• RQ4 고전적 오일러 시뮬레이션과 비교했을 때, 관측 가능한 프레임워크는 충격파를 가로질러 발생하는 비물리적 진동을 어느 정도 제거하는가?
• RQ5 인위적 점성 없이 평균 척도 $\alpha$는 충격파 두께와 수치적 안정성 제어에 어떻게 기여하는가?

주요 결과

• 관측 가능한 발산 정리는 유도 과정 이전에 유량을 유한 척도 $\alpha$에서 평균함으로써, 충격파와 난류에 대한 정규화된 점성 없는 진화 방정식의 유도를 성공적으로 가능하게 한다.
• 관측 가능한 비점성 오일러 방정식은 1934년 레일리가 처음 제안한 점성 없는 레일리 모델로 공식적으로 축소되며, 이 오랜 기간 동안의 모델에 대한 제1원리적 유도를 제공한다.
• 2차원 충격-비틀림 상호작용의 수치 시뮬레이션은 점성 항이나 수치적 소산 없이도 충격파를 가로질러 비물리적 진동이 발생하지 않음을 보여준다. 이는 고해상도에서도 동일하게 유지된다.
• 충격파 두께는 평균 척도 $\alpha$에 의해 제어되며, 더 작은 $\alpha$는 더 높은 해상도가 필요하지만 물리적 구조는 유지된다.
• 관측 가능한 오일러 방정식은 고전적 오일러 방정식과 동일한 이동 충격파 및 엔트로피 해를 지님으로써, $\alpha \to 0$일 때 엔트로피 해로 수렴함을 보여준다.
• 이 방법은 인위적 점성 또는 업윈드 스킴에 의존하지 않고도 복잡한 유동, 예를 들어 충격-난류 상호작용 및 등방성 난류의 안정적이고 정확한 시뮬레이션을 가능하게 한다.