새로운 혼합 유한요소법으로 본 캔 힐러 방정식

초록

본 논문은 캔-힐러 방정식에 대한 새로운 혼합 유한요소법을 제안한다. 혼합 변분형을 도입해 2·3 차원에서 임의 차수의 다항식으로 구현 가능하도록 하였으며, 연속성 요구를 완화한 H(divDiv) 공간을 이용한다. 이론적으로 해의 존재·유일성을 증명하고, 선형화된 전후진 시간스킴에 대한 최적 차수의 오차 추정식을 얻었다. 수치 실험을 통해 이론적 수렴률과 효율성을 확인하였다.

상세 분석

본 연구는 캔-힐러 방정식의 4차 미분 연산자를 직접 다루는 대신, 보조 변수 σ = ∇²u − ε⁻²f(u)I 를 도입해 2차 미분 연산자 두 개로 분해하는 혼합 변분형을 제시한다. σ는 대칭 텐서값 함수이며, divDiv 연산자를 적용할 수 있는 H(divDiv,Ω;S) 공간에 속한다. 저자는 Σ라는 부분공간을 정의하고, Green 항등식을 이용해 Σ에 대한 경계조건을 명시적으로 도출한다. 특히, Π_F(σ n_F)=0 및 n_FᵀDiv σ=0 와 같은 조건이 Σ에 포함됨을 보이며, 이는 기존 비정형 방법보다 구현이 간단하고 차원에 독립적인 통일된 구성을 가능하게 한다.

이후, 원래의 4차 방정식과 혼합 변분형(3.4)이 완전 동치임을 정리 3.4를 통해 증명한다. 여기서는 Nečas 부등식과 Green 항등식을 활용해 u∈H²(Ω)와 σ∈Σ를 각각 얻는다. 특히, σ∈Σ임을 보이기 위해 경계항이 소거되는 조건을 상세히 검토한다.

시간 이산화는 후진 오일러와 반명시적(선형화) 처리로 구성한다. 비선형 항 f(u) = u³ − u는 이전 시간 단계의 값으로 대체해 선형 시스템을 만든다. 공간 이산화는 차수 k≥3인 다항식 기반의 H(div,Ω;S) 및 H(divDiv,Ω;S) 적합 요소를 사용한다. Σ_h는 Π_F(τ_h n_F)=0 및