이소지오메트릭 해석을 위한 모르타르 기법의 최신 동향

초록

본 논문은 복잡한 다패치 구조를 갖는 산업용 모델에서 이소지오메트릭 해석(IgA)을 위한 효율적인 인터페이스 결합 방법으로서 모르타르 기법의 최근 연구 성과를 정리한다. Lagrange multiplier 기반의 변분적 전이, 고차 스플라인 기반 다중 물리·다중 스케일 결합, 그리고 대규모 병렬 구현에 초점을 맞추어, 접촉·유동·구조 연동 등 다양한 응용 사례와 함께 수치적 안정성, 수렴성, 그리고 구현상의 어려움을 논의한다.

상세 분석

본 논문은 이소지오메트릭 해석(IgA)에서 발생하는 다패치(다중 서브도메인) 문제를 해결하기 위한 모르타르 방법의 이론적·수치적 기반을 체계적으로 정리한다. 먼저, 전통적인 도메인 분할 기법이 비정합 메쉬 간 강제적인 점대점 결합에 의존하는 반면, 모르타르 방법은 약한 적분 조건을 통해 인터페이스를 연결함으로써 Lagrange multiplier를 이용한 변분적 일관성을 확보한다는 점을 강조한다. 특히, IgA에서 사용되는 NURBS·B‑spline 기반 고차 연속성 함수는 기존 FEM 대비 높은 스무딩 효과와 고차 미분 연산이 가능하지만, 패치 간 연결 시 연속성 보장이 어려워진다. 이를 해결하기 위해 저자들은 다음과 같은 핵심 기술을 제시한다.

1. 이중형 Lagrange multiplier 공간 설계 – 고차 스플라인 기반의 주 도메인 함수와 일치하도록 multiplier 공간을 선택함으로써 inf‑sup 조건을 만족하고, 수치적 안정성을 확보한다. 특히, bi‑orthogonal spline을 이용한 효율적인 정합(Condensation) 기법을 도입해 전역 시스템의 차원을 크게 감소시킨다.

2. 다중 차원·다중 물리 결합 – 1D‑3D, 2D‑3D 등 차원이 다른 서브도메인 간 전이 연산자를 변분적으로 정의하고, 이를 이용해 유동‑구조 상호작용(FSI), 다중 스케일(MD‑FEM) 및 접촉 문제를 하나의 프레임워크로 통합한다. 특히, 파티션‑오브‑유니티 기반의 가중치 함수를 원자 위치에 부착해 MD와 FEM 사이의 고주파 성분을 필터링함으로써 물리적 일관성을 유지한다.

3. 병렬 구현 및 인터섹션 탐색 – 대규모 HPC 환경에서 비정합 메쉬 간 인터섹션을 효율적으로 찾기 위해 k‑d 트리 기반의 계층적 검색과 space‑filling curve를 활용한 로드 밸런싱 전략을 제시한다. 이와 함께, MoonoLith 라이브러리를 통해 표면·부피 전이 연산을 완전 병렬화하고, 전역 검색의 O(N²) 복잡도를 회피한다.

4. 수렴성·패치 테스트 검증 – 기존 노드 기반 접촉 모델이 패치 테스트를 통과하지 못하고 수렴성이 저하되는 문제를 지적하고, 변분적 모르타르 접촉 조건이 이를 해결함을 수치 실험(예: 바이오프롭테식 대동맥 판막)으로 입증한다.

5. 응용 사례 – 심장 시뮬레이션에서 유동‑구조‑접촉을 동시에 다루는 복합 모델, 지구과학에서 파쇄 네트워크 흐름을 위한 비정합 메쉬 결합, 그리고 다중 스케일 재료 모델링(MD‑FEM) 등을 통해 제안된 방법론의 실용성을 보여준다.

전반적으로 저자들은 이론적 엄밀성(변분적 일관성, inf‑sup 안정성)과 실용적 효율성(고차 스플라인, 병렬 구현, 인터섹션 자동 탐색)을 동시에 만족시키는 모르타르 기반 IgA 결합 프레임워크를 제시한다. 이는 향후 복잡한 산업용 CAD 모델을 직접 해석에 활용하고, 다중 물리·다중 스케일 문제를 대규모 병렬 환경에서 해결하는 데 중요한 기반이 될 것으로 기대된다.