3차원 유체와 2차원 판 상호작용 시스템의 재구성을 위한 유한요소 근사법

초록

본 논문은 3차원 비압축성 점성 유체와 2차원 Kirchhoff‑Plate가 결합된 유동‑구조 상호작용 문제를 다룬다. 판의 4차 미분 방정식을 보조 변수 도입으로 2차 연립식으로 전환하고, 인터페이스 압력의 평균값을 라그랑주 승수로 강제한다. 시간·공간 이산화에 대한 존재·안정성 이론을 구축하고, a‑priori 오차 추정식을 도출한다. 또한 고정점 기반 분할 도메인 분해 알고리즘을 제시해 구현 효율성을 검증한다. 수치 실험을 통해 이론적 수렴률과 실제 물리 문제 적용 가능성을 확인한다.

상세 분석

이 연구는 3차원 Stokes 흐름과 2차원 Kirchhoff‑Plate 사이의 강한 결합을 수치적으로 다루는 데 있어, 기존에 요구되던 H²‑연속 혹은 비정형 P₂‑Morley 요소를 회피하는 새로운 전술을 제시한다. 핵심 아이디어는 판 변위 w에 대해 보조 변수 z = −Δw를 정의함으로써 원래의 4차 방정식(Δ²w) 를 두 개의 2차 방정식(Δz = −∂ₜₜw + p, z + Δw = 0) 으로 분해하는 것이다. 이 변환은 판 변수에 대한 Sobolev 정규성을 H¹ 수준으로 낮추어, 표준 Lagrange 요소(예: P₁) 로도 충분히 근사할 수 있게 만든다.

인터페이스 결합 조건은 “유체의 수직 속도 성분 u₃와 판의 시간 미분 ∂ₜw가 일치한다”는 물리적 제약을 라그랑주 승수 g ∈ G 로 강제한다. 여기서 g는 평균값이 0인 압력 p₀의 트레이스이며, 압력 자체는 평균값 s(t)와 p₀(x,t) 로 분해한다. 이러한 처리 덕분에 압력 공간을 L²₀(Ω_f) 로 제한하면서도 압력 평균을 별도 변수 s 로 관리해 질량 보존을 정확히 만족한다.

시간 이산화는 후방 오일러(1차) 스킴을 사용하고, 공간 이산화는 다음과 같은 혼합 유한 요소 조합을 채택한다. 유체 속도 u와 압력 p₀는 스테이븐스‑퍼렐(𝑃₁–𝑃₀) 혹은 MINI 요소와 같은 inf‑sup 안정 쌍을 이용하고, 판 변수 w와 보조 변수 z는 동일 차수의 H¹‑연속 Lagrange 요소를 사용한다. 라그랑주 승수 g는 H^{-1/2}(Γ) 에 대한 적절한 트레이스 공간 G 로 선택한다.

수학적 분석에서는 먼저 연속‑시간·불연속‑공간 반정형 문제에 대해 연속성 및 coercivity 를 보이고, b_Δt 라는 결합 연산자가 만족하는 inf‑sup 조건을 레머 3.1 로 증명한다. 이를 통해 반정형 시스템(26)의 유일해 존재성을 정리 3.2 로 확보한다. 안정성 정리 3.3 은 에너지 추정식을 전개해, 시간 단계 Δt 와 물성 파라미터(ρ_f, ν_f, ρ) 에 독립적인 상수 C 로 제어되는 ‖u‖, ‖∂ₜw‖, ‖z‖, ‖g‖ 등의 노름을 얻는다.

오차 분석은 연속 문제의 정규성 가정 하에, 시간 차분 오차 O(Δt) 와 공간 근사 오차 O(h^k) (k는 선택한 Lagrange 차수) 를 결합한 a‑priori 추정식을 도출한다. 특히, 라그랑주 승수 g 에 대한 추정이 압력 평균 s 와 결합해 전체 시스템의 수렴성을 보장한다는 점이 주목할 만하다.

계산 측면에서는 고정점 기반 분할 도메인 분해 알고리즘을 설계한다. 각 반복에서 (i) 유체 서브문제를 라그랑주 승수 gⁿ⁺¹ 를 고정하고 풀고, (ii) 판 서브문제를 업데이트된 유체 압력 트레이스를 이용해 해결한다. 수렴 조건은 인터페이스 강성 및 시간 단계 크기에 의해 제한되지만, 실험에서는 비교적 큰 Δt 에도 안정적으로 수렴한다.

수치 실험은 두 가지 케이스로 구성된다. 첫 번째는 인공적으로 만든 제조 해(solution) 를 이용해 공간 차수 1,2,3 에 대한 수렴률을 확인했으며, 이론적 O(h^k) 와 O(Δt) 를 정확히 재현한다. 두 번째는 실제 물리적 파라미터(예: 혈류와 혈관벽, 마이크로플루이딕 채널) 를 적용한 시뮬레이션으로, 판 변형과 유체 속도/압력 필드가 물리적으로 일관된 결과를 보여준다. 특히, H²‑요소 없이도 판 변위의 스무스한 프로파일을 얻을 수 있음을 강조한다.

전체적으로, 본 논문은 4차 판 방정식을 2차 연립식으로 변환함으로써 FEM 구현의 복잡성을 크게 낮추고, 라그랑주 승수 기반 인터페이스 처리와 분할 알고리즘을 결합해 효율적이고 이론적으로 견고한 FSI 해법을 제공한다. 이는 고차원 유체와 얇은 구조물 간 상호작용을 다루는 다양한 공학·생물학적 응용에 바로 적용 가능할 것으로 기대된다.