거친 잡음이 포함된 확률적 Cahn‑Hilliard 방정식의 강인 사후오차 분석 및 적응 FEM 구현

초록

본 논문은 공간‑시간 백색 잡음을 적절히 정규화한 확률적 Cahn‑Hilliard 방정식에 대해, 인터페이스 폭 파라미터 ε와 잡음 정규화 파라미터에 강인한 사후오차 추정식을 제시한다. 선형·비선형 부분을 분리하고, 변형된 선형 변환과 새로운 보간 부등식을 이용해 경로별 오차를 고확률 집합에서 제어한다. 이를 바탕으로 잔차 기반 적응 메쉬 재조정 알고리즘을 설계하고, 완전 이산화된 1차 유한요소·시간 스킴에 대한 수치 실험으로 이론을 검증한다.

상세 분석

이 연구는 기존에 H⁴∩W¹,∞‑trace class 잡음만을 다루던 사후오차 분석을 넘어, 실제 물리 모델에서 등장하는 가장 거친 형태인 공간‑시간 백색 잡음(정규화 전)까지 포괄한다. 백색 잡음은 L²‑정규화가 불가능해 해석적 해가 H² 수준의 공간 정규성을 갖지 못한다는 근본적인 난점을 가지고 있다. 저자들은 (4)식에 의해 백색 잡음을 Ṽ_h에 대한 선형 보간으로 근사함으로써, 잡음이 H¹‑정규성을 갖는 수준으로 낮춘다. 이 정규화는 평균값을 보존하도록 설계돼, 질량 보존 특성을 유지한다.

주요 수학적 전략은 해를 u = ũ + û 로 분해하고, ũ가 선형 SPDE(6)를, û가 비선형 RPDE(7)를 만족하도록 하는 것이다. 선형 부분에 대해서는 기존의 변환 y = ũ – ∫₀ᵗ dW̃ 를 그대로 적용할 수 없으며, 잡음이 H¹‑정규성만을 가짐을 고려해 변형된 변환(9)을 도입한다. 이 변환은 y가 H⁻¹‑시간 미분을 갖도록 보장하고, 결과적으로 (11) 형태의 결정적 PDE와 동등하게 만든다.

비선형 RPDE에 대해서는 기존 연구와 달리 고확률 집합 Ω_{δ,ε}∩Ω_{γ,ε}∩Ω_{ε} 를 정의한다. Ω_{ε}는 선형 SPDE의 근사 오차가 일정 이하인 사건, Ω_{δ,ε}는 L^∞(0,T;H⁻¹)·L⁴(0,T;L⁴) 노름이 제한된 사건, Ω_{γ,ε}는 L^∞(0,T;L⁴) 노름이 제한된 사건이다. Lemma A.3·A.4와 새로운 보간 부등식(Lemma A.1)을 이용해 이 세 집합이 각각 확률 1‑ε 수준으로 큰 확률을 가진다는 것을 증명한다. 따라서 경로별 오차 추정식(Theorem 6.1)은 거의 모든 시뮬레이션 경로에 대해 유효하다.

수치적 측면에서는 비균일 시공간 메쉬와 가변 시간 스텝을 허용하는 완전 이산화 스킴(15)·(17)·(18)을 제시한다. 여기서 사용된 잔차 지표 η_{SPACE,i}, η_{TIME,i}는 전통적인 잔차 기반 적응 FEM 이론을 SPDE에 맞게 확장한 것으로, 각 요소·면에 대한 국소 오차를 정확히 포착한다. 특히, η_{SPACE,1}은 시간 차분에 의한 오차와 면에 대한 화살표(∇e·n) 불연속을 동시에 포함해, 4차 미분 연산자와 잡음의 비정규성을 동시에 제어한다.

주요 정리로는 Theorem 7.1이 있다. 이는 전체 확률 공간 Ω 위에서, ε와 잡음 정규화 파라미터에 독립적인 상수 C(ε,δ)·h^{α}+τ^{β} 형태의 전역 사후오차 경계를 제공한다. 이는 기존 연구가 Ω_{∞}와 같은 제한된 사건에만 적용되던 점을 크게 개선한다. 부록 B에서는 선형 4차 SPDE(6)의 완전 이산화 해에 대해 새로운 수렴률(O(h)+O(τ^{1/2}))을 제시함으로써, 선형 부분 자체도 독립적인 이론적 기여를 한다.

전체적으로, 이 논문은 (i) 거친 잡음에 대한 정규화와 변환 기법, (ii) 고확률 집합을 활용한 경로별 오차 제어, (iii) 잔차 기반 적응 메쉬 전략을 결합해, 실용적인 시뮬레이션 환경에서도 신뢰할 수 있는 오차 추정과 효율적인 계산을 가능하게 만든다.