조합격자 메쉬를 이용한 수송 현상의 변분 형식

초록

본 논문은 셀 복합체 위에서 정의된 조합격자 메쉬 계산(CMC) 체계를 통해 질량, 열, 전하 및 다공성 유체 흐름과 같은 수송 현상의 원시 및 혼합 변분 형식을 제시한다. 연속 미분기하와 이산 형태를 일대일 대응시키며, 물질 계수를 포함한 블록 대각 질량 행렬을 얻어 지역 소거가 가능하도록 설계하였다. 기존의 이산 외부 미분법(DEC)이나 유한요소 외부 미분법(FEEC)과는 달리 매쉬의 매끄러운 임베딩이나 원형 이중 격자를 요구하지 않는다. 규칙적인 및 불규칙적인 2·3차원 메쉬에 대한 수치 실험을 통해 해석해와의 일치를 확인하였다.

상세 분석

이 연구는 복합 재료와 다공성 매체와 같이 서로 다른 위상 차원을 가진 내부 구조가 물성에 직접적인 영향을 미치는 문제를 다루기 위해, 셀 복합체(cell complex)라는 위상학적 기반 위에 물리량을 직접 정의한다는 점에서 혁신적이다. 기존 연속체 모델은 매크로 스케일에서 평균화된 물성만을 반영해 미세 구조의 경계 효과를 손실시키는 반면, CMC는 0‑형식(전위), (D‑1)‑형식(플럭스), D‑형식(보조 전위) 등을 각각 해당 차원의 셀에 할당함으로써 물리량을 자연스럽게 구분한다. 특히, Forman의 조합 미분형식에 Hodge 별과 코미터 연산자를 메트릭 의존적으로 확장한 것이 핵심이다.

연속 미분기하 섹션에서는 외부 미분(d), 코미터(δ), Hodge 별(*)를 이용해 전위‑플럭스 관계와 보존 법칙을 깔끔히 정리한다. 이를 바탕으로 원시 변분(전위만 사용)과 혼합 변분(전위와 플럭스를 동시에 사용) 두 가지 형식을 도출한다. 혼합 형식은 연속 문제에서 잘 알려진 Dirichlet‑Neumann 혼합 형태와 동일 구조를 가지며, 물성 계수를 포함한 내적에 의해 정의된 질량 행렬이 블록 대각 형태가 된다. 이는 각 셀 수준에서 물성 차이를 반영하면서도 전역 시스템을 효율적으로 축소할 수 있게 한다.

이산화 단계에서는 셀 복합체를 ‘조합격자’(combinatorial mesh)라 명명하고, Forman subdivision을 통해 원시 셀과 그 서브디비전 사이의 동치 관계를 설정한다. 여기서 중요한 점은 셀의 볼록성(convexity)을 강제하지 않아 곡면 셀이나 비정형 메쉬도 허용한다는 것이다. 그러나 비정형성은 Hodge 별 정의에 있어 가중치 조정이 필요하고, 수치 실험에서 수렴 속도가 다소 저하되는 원인으로 작용한다.

혼합 변분을 이산화하면 연속 연산자와 정확히 대응되는 차분 연산자 dₕ, δₕ, *_h 가 얻어지며, 물성 가중치가 포함된 내적 ⟨·,·⟩ₕ에 의해 질량 행렬 Mₕ가 형성된다. Mₕ는 각 차원별 셀에 국한된 블록 대각 구조를 가지므로, 플럭스 변수에 대한 로컬 소거(local elimination)를 통해 전위만 남은 스칼라 시스템으로 변환할 수 있다. 이는 기존 DEC가 요구하는 원형 이중 격자와 달리, 임의의 셀 연결성에서도 동일한 소거 절차가 적용 가능함을 의미한다.

수치 실험에서는 2D 정규 사각형 메쉬, 3D 정규 육면체 메쉬, 그리고 비정형 Delaunay 기반 메쉬에 대해 정밀 해와 비교하였다. 질량 확산, 열 전도, 전하 전송, 다공성 유체 흐름 네 가지 물리 모델 모두에서 L² 오차가 메쉬 크기 h에 대해 1차 수렴을 보였으며, 특히 정규 메쉬에서는 이론적 수렴률에 근접하였다. 비정형 메쉬에서는 경계 셀의 형태 왜곡이 오차 상수를 증가시켰지만, 전반적인 수렴 특성은 유지되었다.

결론적으로, CMC는 위상 구조와 물성 이질성을 동시에 고려하는 수송 현상 모델링에 강력한 도구이며, 블록 대각 질량 행렬을 통한 효율적인 선형 시스템 해석, 매쉬 형태에 대한 낮은 제약조건, 그리고 연속·이산 해석 사이의 명확한 수학적 대응 관계를 제공한다. 향후 비정형 메쉬에 대한 사전 가중치 최적화, 비선형 및 시간 의존 문제 확장, 그리고 다중 물리 결합 시뮬레이션에 적용하는 연구가 기대된다.