비국소 2차원 전송 방정식 스킴의 수렴성 연구

초록

본 논문은 변위 밀도 역학을 기술하는 비국소 2차원 전송 방정식 시스템에 대해, Fejér 커널을 이용한 속도 정규화와 IMEX(반명시적) 유한 차분 스킴을 제안한다. 제안된 스킴은 이산적인 기울기 엔트로피 추정량을 보존하며, 이를 기반으로 격자 해가 연속 해로 수렴함을 증명한다. 최초의 수렴 결과이며, 수치 실험을 통해 스킴의 효율성을 확인한다.

상세 분석

이 연구는 변위 밀도(ρ⁺, ρ⁻)를 기술하는 비국소 시스템을 다루며, 시스템은 Riesz 변환을 통한 비국소 항 K∗(ρ⁺−ρ⁻)와 외부 전단 변형 a(t)·b·∇ρ 형태의 비선형 전송 방정식으로 구성된다. 주요 난점은 속도장이 Riesz 변환의 합성으로 정의돼 급격히 특이하고, 초기 데이터와 속도 모두 낮은 정규성을 가진다는 점이다. 기존 이론에서는 엔트로피 추정식
∫∂₁ρ ln ∂₁ρ dx + ∫₀ᵗ∫|R₁R₂(ρ⁺−ρ⁻)|² dx ds ≤ 초기값
을 이용해 존재성을 확보했지만, 수치 스킴이 이 추정량을 보존하지 못하면 수렴을 보장할 수 없었다.

저자들은 먼저 K를 Fejér 커널 F₂M과의 컨볼루션 σ_KM=F₂M∗K 로 정규화한다. M을 크게 하면 σ_KM→K가 되며, 정규화된 시스템은 부드러운 속도장을 제공해 고전적인 차분 연산이 가능하도록 만든다. 정규화된 초기 데이터도 동일하게 Fejér 커널을 적용해 부드럽게 만든다.

수치 스킴은 시간에 대해 IMEX(반명시적) 오일러, 공간에 대해 upwind 차분을 사용한다. 비보존형이지만, 스킴은 이산 기울기 엔트로피
Eⁿ=∑_i Δx (∂₁ρ_iⁿ+L) ln(∂₁ρ_iⁿ+L)
를 감소시키는 구조를 갖는다. 핵심은 σ_KM∗(ρ⁺−ρ⁻)의 Fourier 계수를 이용해 파르세발 등식에서 비음성 항을 얻는 것이다. 이를 통해 스킴이 L¹ 안정성, 총 변동성 제한, 그리고 시간 연속성(시간 차분에 대한 BV 추정)을 만족함을 증명한다.

다음 단계에서는 격자 크기 Δx, Δt에 대한 균일 추정량을 도출한다. 특히, M 고정 하에 Δx→0, Δt→0을 취하면 재구성된 해가 강한 L¹ 수렴을 보이며, 한계 과정에서 비국소 항은 σ_KM∗(·) 형태로 남아 연속 문제의 약해 해와 일치한다. 마지막으로 M→∞ 한계에서 σ_KM이 원래 커널 K로 수렴함을 이용해 정규화된 해가 원 시스템의 해로 수렴함을 보인다.

이러한 일련의 분석은 비국소, 비보존, 비엄격 초과파동계에 대해 최초로 전산적 수렴성을 확보한 사례이며, 엔트로피 보존을 핵심 도구로 삼은 점이 혁신적이다. 또한, Fejér 커널을 통한 정규화가 물리적 의미(고주파 차단)와 수학적 편리성(양의 Fourier 계수) 모두를 제공한다는 점에서 유용하다.