이질 다공성 매질에서 변분다중스케일 요소와 라비아트 트루스톤 요소의 성능 비교

초록

본 연구는 이질적인 투과성(비율 1 × 10⁵)과 개방형 압력 구동 경계조건을 갖는 2차원 대규모 지하수 흐름·수송 문제에 대해, 저차 라비아트‑트루스톤(RT0) 요소와 변분다중스케일(VMS) 안정화 요소의 수치적 성능을 비교한다. 흐름장을 속도 입력으로 사용했을 때 질량 보존 위반 정도를 평가한 결과, VMS 요소는 투과성 계층에 대해 흐름이 층과 거의 평행할 때 큰 질량 생성·소실을 보였으며, 흐름이 층에 수직에 가깝거나 무작위 배향일 경우에는 국부 질량 균형 위반이 미미하였다. 이러한 결과는 VMS 요소를 사용할 때 문제의 기하·물성 특성을 사전에 검토해야 함을 시사한다.

상세 분석

본 논문은 두 가지 혼합 유한요소, 즉 저차 라비아트‑트루스톤(RT0) 요소와 변분다중스케일(VMS) 안정화 요소를 이질 다공성 매질의 흐름·수송 시뮬레이션에 적용해 성능을 정량적으로 비교하였다. 연구자는 투과성 비율이 10⁵에 달하는 고이질성 매질, 압력 구동형 개방 경계조건, 그리고 2차원 대규모 영역을 대상으로 네 가지 대표적인 테스트 케이스를 설계하였다. 첫 번째 케이스는 층상 구조가 흐름과 직각으로 배열된 경우이며, 두 번째는 흐름이 층에 거의 평행한 경우, 세 번째는 무작위로 배치된 고투과성/저투과성 블록, 네 번째는 복합적인 경계조건을 가진 실제 지형 모델이다. 각 케이스마다 RT0와 VMS 요소로 흐름장을 계산한 뒤, 이를 대류 속도로 사용해 대류‑확산 수송 방정식을 풀었다. 질량 보존 위반은 (1) 전역 질량 손실/생성, (2) 셀별 질량 불균형, (3) 수송 해석에서의 오차 증폭 세 가지 지표로 평가되었다.

결과는 다음과 같이 요약된다. RT0 요소는 연속적인 유속 보존과 셀 경계에서의 정확한 질량 흐름을 제공해 전역·국부 모두에서 거의 0에 가까운 질량 불균형을 보였다. 반면 VMS 요소는 안정화 항에 의해 인위적인 확산(수치적 인공점성)이 도입되어, 흐름이 고투과성·저투과성 경계에 평행하게 흐를 때(특히 케이스 2) 유속이 과도하게 스무딩되어 질량이 인위적으로 생성·소멸하는 현상이 관찰되었다. 이러한 현상은 VMS가 기본적인 보존 법칙을 완전히 만족시키지 못한다는 근본적인 한계와 연결된다. 그러나 흐름이 경계에 수직에 가깝거나, 매질 이질성이 무작위적으로 분포된 경우(케이스 1·3·4)에는 VMS의 인공 확산이 실제 물리적 확산과 유사하게 작용해 질량 위반이 미미했다.

또한, 퍼포먼스 측면에서 VMS 요소는 고차 자유도를 사용하지 않음에도 불구하고 행렬 조건수가 RT0보다 우수해 선형 해석기에서 수렴 속도가 빠른 장점이 있었다. 하지만 질량 보존이 중요한 장기 수송 예측(오염 물질 이동, 탄소 격리 등)에서는 이러한 속도 이점이 질량 오류와 트레이드오프 관계에 있음을 명확히 보여준다. 논문은 최종적으로 VMS 요소를 적용할 때는 (i) 흐름 방향과 층 구조의 정렬 정도를 사전 분석하고, (ii) 필요 시 보존 강화 기법(예: Lagrange multiplier, 재조정된 안정화 파라미터)과 결합해 사용해야 함을 제언한다.