다중 스케일 유동‑구조 상호작용을 위한 단일체 고차 ALE 유한요소법

초록

본 논문은 고차 등변형 ALE 프레임워크와 IMEX‑PRK 시간 적분을 결합한 단일체(FSI) 해법인 MLH‑ALE를 제안한다. P₂ 등변형 요소와 국부적 메쉬 업데이트 전략을 통해 미세 입자와 마이크로채널 같은 다중 스케일 문제에서 높은 정확도와 효율성을 확보한다. 수치 실험에서 최적 차수 수렴을 확인했으며, 나선형 마이크로채널 내 입자 집속 시뮬레이션이 실험 결과와 일치함을 보였다.

상세 분석

MLH‑ALE는 기존 ALE 기반 FSI 해법의 두 가지 핵심 한계를 동시에 극복한다. 첫째, 다중 스케일 상황에서 전체 메쉬를 전역적으로 재생성하면 계산 비용이 급증하고 메쉬 품질이 저하된다. 저자들은 “국부적 업데이트” 알고리즘을 도입해 구조물 주변 영역만 동적으로 재메쉬하고, 나머지 영역은 기존 메쉬를 그대로 유지한다. 이 방식은 장거리 입자 이동을 정확히 추적하면서도 전체 자유도 감소 효과를 가져온다. 둘째, 고차 기하학적 정확성을 확보하기 위해 등변형 P₂ 요소와 등변형 매핑을 사용한다. 등변형 요소는 기하학적 오차를 최소화하고, Cauchy‑Green 텐서 B의 P₁ 근사를 통해 변형률을 안정적으로 계산한다. 시간 적분 측면에서는 IMEX‑PRK 스킴을 채택해 강성 있는 구조‑유체 결합을 암시적으로 처리하면서, 비강성 부분은 명시적으로 계산한다. 이는 큰 시간 스텝에서도 안정성을 유지하고, 고차 정확도를 보장한다. 또한, ALE 매핑을 고정된 기준 영역으로 전이시켜 시간‑의존적 도메인 문제를 정적 형태로 변환함으로써, 고차 시간 적분에 필요한 Jacobian 및 변형 텐서의 정확한 평가가 가능해진다. 변형 텐서의 진화 방정식을 직접 포함시켜 물성 비선형성을 자연스럽게 다루며, 압력-속도 연동을 단일 시스템으로 풀어 결합 강성을 완전히 제거한다. 수치 실험에서는 2D/3D 벤치마크에서 이론적 수렴 차수(p+1)를 만족하고, 나선형 마이크로채널에서 입자 집속 현상이 실험 데이터와 정량적으로 일치함을 보여준다. 전체적으로 MLH‑ALE는 고차 정확도, 메쉬 효율성, 그리고 다중 스케일 물리 현상의 정밀한 포착이라는 세 축을 동시에 만족시키는 혁신적인 FSI 해법이라 할 수 있다.