경계조건 페널티를 이용한 피FEM 방식의 포아송 방정식 해법

초록

본 논문은 레벨셋 기반의 비맞춤형 유한요소법인 φFEM을 변형하여, 경계조건을 페널티 항으로 강제하는 새로운 스키마를 제안한다. 고정된 격자 위에서 경계 근처 셀만 레벨셋 값을 필요로 하며, 고스트 페널티를 통해 수치적 안정성을 확보한다. H¹ 반노름에서 최적 수렴률을, L² 노름에서는 준최적 수렴률을 이론적으로 증명하고, 조건수와 계산 비용을 실험적으로 검증한다.

상세 분석

이 연구는 기존 φFEM이 전 영역에서 u=φw 형태의 해를 가정하고 경계조건을 자동으로 만족시키는 방식과 달리, 경계 근처 셀에만 레벨셋 φ를 사용하고, 추가 변수 p를 도입해 u=φp 관계를 페널티 항으로 강제한다는 점에서 차별화된다. 페널티 파라미터 γ와 고스트 페널티 계수 σ_D를 충분히 크게 잡으면, bilinear form a_h가 정의된 복합 노름 ‖(u,p)‖_h에 대해 coercivity를 확보할 수 있다. 이를 위해 저자들은 두 개의 보조 Lemma를 활용한다. 첫 번째는 φ_h p_h 를 ‖∇u_h‖와 ‖u_h−φ_h p_h‖ 로 제어하는 불등식이며, 두 번째는 경계 셀에 대한 추적 및 라플라시안 항을 포함한 불등식으로, α∈(0,1)와 β>0를 적절히 선택해 ‖∇v_h‖_ΩΓ_h 를 전체 ‖∇v_h‖_Ω_h 와 경계 항들의 조합으로 제한한다. 이러한 결과를 바탕으로 Proposition 1에서 a_h가 ‖·‖h 에 대해 coercive함을 증명하고, Theorem 1에서는 H¹ 반노름에서 |u−u_h| ≤ C h^k ‖f‖{k−1} 라는 최적 수렴률을 도출한다. L² 노름에 대해서는 O(h^{k+½})의 준최적 결과를 얻지만, 수치 실험에서는 실제 O(h^{k+1}) 수렴을 확인한다. 또한, 행렬 조건수는 O(h^{−2}) 로 기존 비맞춤형 방법과 동일하게 스케일링됨을 보인다. 구현 측면에서는 DOLFINx 기반 파이썬 코드를 사용했으며, 2차원 복잡한 간판형 도메인과 3차원 구형 도메인 두 가지 테스트에서 표준 FEM, 직접 φFEM, 제안된 페널티 φFEM(‘Dual φFEM’)을 비교한다. 결과는 세 방법 모두 H¹ 및 L² 노름에서 최적 수렴을 보이지만, 페널티 φFEM은 레벨셋을 경계 셀에만 필요로 하므로 전처리 비용이 감소하고, 시스템 크기가 약간 커짐에도 불구하고 전체 계산 시간이 직접 φFEM보다 크게 늘어나지 않는다. 조건수 역시 안정적인 범위에 머무른다. 따라서 이 논문은 경계조건을 강제하는 새로운 페널티 기반 φFEM이 이론적 안정성과 실용적 효율성을 동시에 만족함을 입증한다.