비음수 유지 유한요소법: 곱셈 잡음이 있는 파라볼릭 SPDE를 위한 새로운 접근

초록

본 논문은 초기값이 비음수인 경우 해가 비음수를 유지하는 파라볼릭 SPDE에 대해, 질량-라밍(lumped mass) 유한요소법을 도입하여 공간 이산화 단계에서 비음수 보존을 무조건적으로 달성한다. 제시된 반분산(semidiscrete) 스킴은 1차 정확도를 가지며, 선형 경우에는 완전 이산화(time‑stepping) 스키마까지 확장하여 비음수 보존과 수렴성을 동시에 만족한다. 수치 실험을 통해 기존 유한차분·유한요소 방법보다 음수 발생이 현저히 적은 것을 확인하였다.

상세 분석

논문은 먼저 (1) (du-\Delta u,dt = f(u),dW) 형태의 파라볼릭 SPDE를 정의하고, 초기값 (u_0\ge0) 와 (f(0)=0) 조건 하에 해가 비음수를 유지한다는 비교 원리를 정리한다. 이를 바탕으로 비음수 보존을 위한 유한요소 설계 원칙을 제시한다. 핵심은 (i) ‘아크루스(weakly acute)’ 메쉬 조건을 만족하는 삼각형(또는 사면체) 분할을 사용하고, (ii) 질량 행렬을 라밍(lumped) 처리하여 대각화된 형태로 만든다. 라밍 질량 행렬은 각 자유도에 대한 질량을 해당 노드에만 집중시켜, 변분식에서 ((\pi_h