공간 시간 매끄러운 인공 점성법으로 비선형 보존법칙 해결

초록

본 논문은 압축성 유체의 1차원 Euler 방정식에 국한하여, 충격파·희소파·접촉불연속을 정확히 포착하기 위한 새로운 인공 점성 기법인 C‑method를 제안한다. 스칼라 반응‑확산 방정식으로 정의된 가중치 함수 C(x,t)를 이용해 점성 계수를 공간‑시간에 걸쳐 부드럽게 조절함으로써, 충격 전후에서만 국부적으로 큰 점성을 부여하고, 평탄한 영역에서는 점성을 거의 없앤다. FEM‑C, WENO‑C, WENO‑LF‑C 등 세 가지 수치 구현을 통해 전형적인 테스트 케이스(Sod, Osher‑Shu, Woodward‑Colella, Leblanc)에서 높은 차수의 정확도와 최소한의 오버슈트·진동을 확인하였다.

상세 분석

C‑method의 핵심은 기존 인공 점성 기법이 갖는 ‘전역적 점성’ 문제를 피하기 위해, 점성 계수를 동적으로 생성하는 별도의 PDE를 도입한 점이다. 구체적으로, 선형 스칼라 반응‑확산 방정식
∂ₜC – εΔC = (1/τ)·S(u) – (1/τ)·C
을 Euler 시스템에 결합한다. 여기서 S(u)는 압축성 흐름에서 발생하는 급격한 변화(예: 압력·밀도 기울기)를 정량화하는 감지 함수이며, ε와 τ는 각각 점성의 공간적 확산 정도와 시간적 반응 속도를 조절한다. 이 방정식은 C가 충격 전후에서 급격히 상승하고, 이후 확산에 의해 부드럽게 감소하도록 보장한다. 따라서 인공 점성 항 ν·C·∂ₓu는 충격 근처에서만 유의미하게 작용하고, 평탄 구간에서는 거의 0에 가까워진다.

수학적으로는 C‑method가 기존 Richtmyer‑Von Neumann 인공 점성의 정규화 버전임을 보이며, C가 L²‑bounded하고, ε→0, τ→0 한계에서 원래 보존법칙의 약해해(weak solution)와 일치함을 증명한다. 이는 점성 계수가 급격히 변하는 전통적 방법에 비해 수치적 안정성을 크게 향상시키면서도, 물리적 충격 구조를 왜곡하지 않는다는 강점을 제공한다.

세 가지 구현 방식은 각각 다른 차수와 스키마를 사용한다. FEM‑C는 연속 Galerkin 형식으로 2차 시간·공간 스키마를 적용해 고전적인 유한요소 해법과 결합한다. WENO‑C는 고차 비선형 재구성을 이용해 C‑방정식과 동시에 해결함으로써, 기존 WENO‑Riemann 해법에서 발생하는 과도한 점성 없이도 고해상도 충격 포착이 가능하다. 마지막으로 WENO‑LF‑C는 Lax‑Friedrichs 플럭스를 사용해 계산 비용을 낮추면서도 C‑가 제공하는 국부 점성 제어를 유지한다.

실험 결과는 특히 Leblanc 충격관에서 두드러진다. 기존 고차 WENO‑Riemann 해법은 압력·에너지 스케일 차이 때문에 큰 오버슈트와 비물리적 진동을 보였지만, C‑method 기반 구현은 이러한 문제를 현저히 감소시켰으며, 전반적인 L¹·L² 오차도 기존 방법보다 우수했다. 또한, 충격 전후의 접촉면에서도 비대칭적인 스무딩이 최소화되어, 접촉불연속의 정확한 위치와 강도를 유지하였다.

요약하면, C‑method는 (1) 점성 계수를 충격에만 국부화, (2) 공간‑시간에 부드럽게 변하도록 설계, (3) 다양한 수치 스키마에 쉽게 통합 가능, (4) 이론적 수렴 보장 및 실험적 정확도 향상이라는 네 가지 장점을 제공한다. 이는 고차 비선형 보존법칙 해석에 있어 인공 점성 설계의 새로운 패러다임을 제시한다.