전역형 유한차분 기반 유체‑구조 결합 시뮬레이션 방법

초록

본 논문은 고정된 직교 격자 위에서 유체와 변형체의 운동 방정식을 모두 유한차분으로 풀어내는 전역형(Eulerian) 접근법을 제시한다. 다상 흐름에 널리 쓰이는 부피‑오브‑플루이드(VOF) 기법으로 복합 기하를 기술하고, 고체 변형은 라이트 카우시‑그린 변형 텐서를 Eulerian 프레임에서 시간에 따라 업데이트한다. 이를 통해 비선형 Mooney‑Rivlin 고체의 constitutive 관계를 직접 적용할 수 있다. 다양한 검증 사례와 수치 정확도 분석을 통해 제안 방법의 유효성을 입증한다.

상세 요약

이 연구는 전통적인 유체‑구조 상호작용(FSI) 해법이 갖는 격자 재구성 혹은 라그랑지안-오일러리안 혼합 접근법의 복잡성을 회피하고자, 전역형(Eulerian) 프레임에서 모든 물리량을 정의한다는 근본적인 설계를 채택하였다. 핵심은 두 가지 기술적 혁신에 있다. 첫째, 다중 물질 영역을 구분하기 위해 Hirt와 Nichols가 제안한 부피‑오브‑플루이드(VOF) 방법을 차용하였다. VOF는 각 셀에 물질 상을 나타내는 색상 함수(φ)를 정의하고, φ의 보존식 ∂φ/∂t + u·∇φ = 0 를 유한차분으로 풀어 물·고체 경계의 이동을 정확히 추적한다. 기존 VOF는 기체‑액체 인터페이스에 주로 쓰였지만, 여기서는 고체와 액체 사이의 강체 경계에도 적용함으로써 격자 재구성 없이 복합 기하를 표현한다는 점이 독창적이다.

둘째, 고체 변형을 라이트 카우시‑그린 변형 텐서(B = FFᵀ) 형태로 Eulerian 좌표계에 저장한다. 전통적인 라그랑지안 방식에서는 물질 점을 따라 변형 텐서를 직접 업데이트하지만, 이 논문은 B를 격자 기반 변수로 두고 물질 점이 이동함에 따라 B를 대류 방정식 ∂B/∂t + u·∇B = B·∇u + (∇u)ᵀ·B 로 전개한다. 이렇게 하면 고체 내부의 비선형 변형을 실시간으로 계산하면서도 격자 변형 없이 동일한 차분 스킴을 사용할 수 있다. 특히 Mooney‑Rivlin 같은 비선형 초탄성 모델의 응력 텐서 σ_s = -pI + 2∂W/∂I₁ B + 2∂W/∂I₂ B⁻¹ 형태를 B에 직접 대입해 계산함으로써 물성식 구현이 간단해진다.

수치 구현 측면에서는 2차 중앙 차분과 2차 Runge‑Kutta(또는 Adams‑Bashforth) 시간 적분을 결합해 일관된 정확도를 확보한다. 압력-속도 결합을 위해 압력 보정(Projection) 방법을 사용하고, 고체와 유체 사이의 경계 조건(속도 연속·응력 연속)은 VOF와 B 텐서가 공유하는 셀에서 강제적으로 만족시킨다.

검증 사례로는 (1) 단순 진동판 문제에서 해석 해와 수치 해의 2차 수렴을 확인, (2) 대변형 비틀림 플레이트에서 Mooney‑Rivlin 고체와 뉴턴 유체의 상호작용을 재현, (3) 복잡한 다중 고체‑유체 혼합 흐름에서 물‑고체 인터페이스가 정확히 유지되는지를 실험 데이터와 비교하였다. 특히 격자 해상도를 증가시켰을 때 오류가 L2 노름 기준으로 2차 수렴함을 보이며, 전통적인 ALE(Arbitrary Lagrangian‑Eulerian) 방식 대비 격자 재구성 비용이 현저히 낮음을 입증한다.

이 방법의 한계는 고체가 큰 변형을 겪을 경우 B 텐서가 수치적으로 불안정해질 가능성이 있다는 점이다. 이를 완화하기 위해 텐서 정규화 혹은 로그 변형 텐서(log‑B) 방식이 제안될 수 있다. 또한, VOF 기반 경계 추적은 인터페이스가 매우 얇은 경우(예: 얇은 고체 막)에서 수치 확산이 발생할 수 있어, 고해상도 재구성 혹은 고차 VOF 스키마가 필요하다.

전반적으로 이 논문은 전역형 유한차분 기반 FSI 해법을 체계화하고, 다상 흐름에 쓰이던 VOF와 라이트 카우시‑그린 텐서 업데이트를 결합함으로써 격자 재구성 없이 복합 유체‑구조 문제를 효율적으로 해결할 수 있음을 보여준다. 향후 3차원 확장, 비등방성 고체 모델, 그리고 고속 충돌·충격 문제에 적용한다면, 전산 유체역학·구조역학 분야에서 중요한 도구가 될 전망이다.