열변동을 고려한 관성 유체‑구조 상호작용의 체계적 확률적 축소

초록

본 논문은 열변동이 존재하는 미세 탄성 구조물과 유체 사이의 관성 상호작용을 다루며, (i) 유체‑구조 간 동력 전달이 강한 경우, (ii) 구조물 질량이 유체가 대체하는 질량에 비해 작을 때, (iii) 유체 응력이 빠르게 준정상 상태에 도달하는 경우를 대상으로 확률적 축소 모델을 유도한다. Stochastic Eulerian Lagrangian Method(SELM)를 기반으로 Backward Kolmogorov 방정식에 대한 특이 섭동 분석을 수행해 유효 동역학 방정식을 도출한다.

상세 요약

본 연구는 열변동이 중요한 마이크로스케일 시스템에서 유체와 탄성 구조물 간의 상호작용을 정량적으로 기술하기 위해, 기존의 관성 유체‑구조 상호작용 모델에 확률적 요소를 체계적으로 삽입한다. 저자들은 먼저 Stochastic Eulerian Lagrangian Method(SELM)를 이용해 연속체 수준에서 유체와 구조물의 운동 방정식을 결합한다. SELM은 라그라지안 입자(구조물)와 오일러적 유체 필드 사이의 상호작용을 커플링 연산자 Γ와 전이 연산자 Λ로 표현하며, 열변동은 가우시안 백색 잡음 ξ(t)와 연관된 플럭스‑디스펜스 관계에 의해 보존된다.

핵심은 세 가지 물리적 제한 조건을 동시에 고려한 확률적 축소 과정이다. 첫 번째 조건은 커플링 강도 κ가 무한대에 가까워지는 강결합(limit)으로, 이 경우 구조물과 유체의 속도 차이가 거의 사라져 상대 속도가 제로에 수렴한다. 두 번째는 구조물 질량 m이 유체가 대체하는 질량 ρV에 비해 충분히 작아, 관성 항이 무시될 수 있는 저질량(limit)이다. 세 번째는 유체의 점성 확산 시간 τν가 구조물의 동역학 시간 스케일 τs보다 훨씬 짧아, 유체 응력이 즉시 준정상 상태에 도달한다는 가정이다.

이러한 가정을 수학적으로 구현하기 위해 저자들은 확률 과정의 Backward Kolmogorov 방정식을 도입하고, 작은 매개변수 ε(=m/ρV, 1/κ, τν/τs 등)의 함수로 전개한다. 특이 섭동 해법을 적용하면, 빠른 변수(유체 속도와 압력)는 ε→0에서 고정점 분포를 형성하고, 느린 변수(구조물 위치와 변형)는 평균화된 drift와 유효 확산항을 갖는 저차원 확률 미분 방정식으로 기술된다. 특히, 유효 드리프트는 원래의 유체‑구조 커플링 항이 ΓΛ 형태로 축소된 결과이며, 유효 확산은 fluctuation‑dissipation theorem에 의해 결정된 온도와 점성에 의존한다.

결과적으로, 원래 3‑D 유체와 N‑차원 구조물의 복합 시스템은, 강결합·저질량·빠른 유체 응답이라는 세 가지 제한 하에, 구조물 위치와 변형만을 포함하는 저차원 확률 미분 방정식으로 완전히 대체될 수 있다. 이 방정식은 수치적으로 구현이 용이하고, 열변동에 의한 랜덤 포스가 정확히 보존되는 형태를 유지한다는 점에서, 마이크로플루이딕스, 바이오물리컬 시뮬레이션, 나노스케일 재료 설계 등에 직접적인 활용 가능성을 제시한다.