| 유체-구조 연성 Fluid-Structure Interaction | |
|---|---|
| 약칭 | FSI |
| 분야 | 전산유체역학 × 구조해석 (다물리) |
| 다루는 것 | 유체 하중 ↔ 구조 변형의 상호작용 |
| 대표 격자 기법 | ALE, 침지경계법(IBM) |
1. 개요[편집]
유체-구조 연성(Fluid-Structure Interaction, FSI)은 흐르는 유체가 고체 구조에 힘을 가해 변형시키고, 그 변형된 구조가 다시 유체의 흐름을 바꾸는 — 서로가 서로를 계속 밀고 당기는 상호작용을 해석하는 다물리(multiphysics) 시뮬레이션 분야다. 즉 전산유체역학과 구조해석을 한 판에 얹어 동시에 푸는 문제다.1
바람에 펄럭이는 깃발, 혈류에 부풀었다 오므라드는 혈관과 심장 판막, 착륙하는 여객기 날개의 휨, 다리를 무너뜨린 그 유명한 타코마 내로스 교량의 자기여기진동 — 전부 FSI다. 유체와 구조 어느 한쪽만 풀어서는 답이 안 나온다. 한쪽이 바뀌면 다른 쪽이 바뀌고, 그게 또 처음 쪽을 바꾸기 때문이다. 이 순환 고리를 어떻게 수치적으로 닫아주느냐가 FSI 해석의 전부라고 해도 과언이 아니다.
2. 단방향 vs 양방향 연성[편집]
FSI를 처음 만나면 반드시 정하고 들어가야 하는 갈림길이 있다.
2.1. 단방향 연성 (One-way)[편집]
유체 → 구조로 힘만 넘긴다. CFD로 유동장을 풀어 압력·전단 하중을 구하고, 그 하중을 구조해석에 경계조건으로 넘겨 변형을 계산한다. 단, 구조의 변형이 유동을 바꾸는 되먹임은 무시한다. 구조가 아주 뻣뻣하거나 변형이 유동장에 비해 미미할 때만 정당화된다. 계산이 싸고 간단하다는 게 장점.2
2.2. 양방향 연성 (Two-way)[편집]
유체 ↔ 구조로 힘과 변형을 매 시간스텝(또는 반복)마다 주고받는다. 구조가 변형되면 유체 도메인의 형상 자체가 바뀌므로 격자를 다시 움직여야 하고, 바뀐 유동이 다시 하중을 갱신한다. 물리적으로 정확하지만 비싸고 까다롭다. 유연한 구조, 큰 변형, 공진 근처에서는 선택의 여지가 없다.
정리하면: 구조가 유동을 안 바꾼다고 확신하면 단방향, 아니면 얌전히 양방향으로 가야 한다. “단방향으로 대충 하고 싶다”는 유혹은 강하지만, 자기여기진동(flutter) 같은 현상은 오직 양방향에서만 나타난다.
3. 강결합 vs 약결합[편집]
양방향으로 가기로 했다면, 다음 갈림길은 유체 솔버와 구조 솔버를 얼마나 단단히 묶느냐다.
- 약결합(weak / loosely coupled, 분할형 partitioned): 한 시간스텝 안에서 유체를 한 번, 구조를 한 번 순차적으로 풀고 넘어간다. 기존의 성숙한 CFD 코드와 구조 코드를 각자 두고 인터페이스로만 데이터를 주고받으면 되니 구현이 쉽고 재활용성이 좋다. 대신 인터페이스에서 힘-변형 평형이 시간스텝 내에서 완벽히 맞진 않아 오차·불안정이 생길 수 있다.
- 강결합(strong / tightly coupled): 한 시간스텝 안에서 유체-구조를 여러 번 왕복(subiteration)하며 인터페이스 평형이 수렴할 때까지 반복하거나, 아예 하나의 거대한 연립방정식(monolithic)으로 통째로 푼다. 안정적이고 정확하지만 비싸고 코드가 복잡하다.3
모놀리식(monolithic)은 유체·구조 미지수를 한 행렬에 넣어 뉴턴-랩슨법으로 동시에 푸는 극단적 강결합이다. 견고하지만 만들기 지옥이라, 산업 현장에서는 분할형 강결합(subiteration 반복)이 현실적 타협점으로 많이 쓰인다.
4. 부가질량 불안정성 — FSI의 흑막[편집]
약결합 FSI를 처음 짜본 사람은 거의 반드시 이 벽에 부딪힌다. 분명 유체도 구조도 각각 잘 도는데, 둘을 순차로 묶으면 시간스텝을 아무리 줄여도 계산이 발산해버리는 현상이다. 이게 그 악명 높은 **부가질량 불안정성(added-mass instability)**이다.4
원인은 이렇다. 구조가 유체 속에서 가속되려면 주변 유체도 함께 밀어내야 하는데, 이 “끌려오는 유체의 관성”이 마치 구조에 질량이 더 붙은 것처럼 작용한다. 이것이 부가질량이다. 문제는 가벼운 구조가 무겁고 비압축인 유체(예: 물)를 밀 때 부가질량이 구조 자체 질량을 압도하는 경우다. 이때 순차적 약결합 방식은 인터페이스에서 힘과 변형이 서로를 폭주시키며 발산한다. 놀랍게도 시간스텝을 줄여도 안정화가 안 되고 오히려 나빠지는데, 이는 CFL류의 조건이 아니라 밀도비에서 오는 구조적 불안정이기 때문이다.
혈류-혈관처럼 유체와 구조의 밀도가 비슷한 생체역학 FSI가 특히 이 지옥에 잘 빠진다. 해법은 강결합 subiteration에 완화(under-relaxation, Aitken 가속 등)를 걸거나, 아예 모놀리식으로 가는 것이다. 공기처럼 가벼운 유체를 다루는 항공 FSI에서는 상대적으로 덜 심각하다.
5. ALE — 움직이는 격자를 다루는 법[편집]
양방향 FSI에서 구조가 변형되면 유체 도메인의 경계가 움직인다. 그럼 유체 격자를 어떻게 처리할까? 여기서 등장하는 것이 ALE(Arbitrary Lagrangian-Eulerian, 임의 라그랑주-오일러) 기술이다.
유체는 보통 공간 고정 격자(오일러 관점)로 풀고, 구조는 물질을 따라가는 격자(라그랑주 관점)로 푼다. ALE는 이 둘의 절충으로, 격자를 유체 물질과도 구조와도 무관하게 임의로 움직일 수 있게 한다. 경계는 구조를 따라 움직이되, 내부 격자는 뒤틀림을 줄이도록 부드럽게 재배치(mesh motion/morphing)한다. 덕분에 경계 근처는 라그랑주처럼 물질을 정확히 추적하고, 내부는 격자 품질을 유지한다.
물론 격자가 너무 크게 뭉개지면 재격자화(remeshing)를 해야 하고, 여기서 CFD 실무자의 정신력이 갈려나간다. 대안으로 격자를 아예 안 움직이고 고체를 유동장에 “담가버리는” **침지경계법(Immersed Boundary Method, IBM)**도 널리 쓰인다. 큰 변형·복잡 형상에 강하지만 경계에서의 정확도 확보가 숙제다.
6. 어디에 쓰나[편집]
- 항공우주: 날개 플러터 예측, 낙하산 전개, 로켓 노즐 진동. 자기여기진동은 기체를 순식간에 부숴놓기에 반드시 잡아야 한다.
- 생체역학: 심장 판막, 동맥류, 인공심장. 밀도비가 1에 가까워 부가질량 지옥의 본진.
- 토목/에너지: 교량 내풍 안정성, 풍력·조류 터빈 블레이드, 송전선 갤러핑. 다상유동이 얽히면 난이도가 한 단계 더 오른다.
- 일상: 깃발, 돛, 나뭇잎 — 눈에 보이는 펄럭임은 죄다 FSI다.
7. 관련 문서[편집]
8. Footnotes[편집]
-
FSI는 다물리(multiphysics) 연성해석의 대표 격이다. 열-구조, 전자기-열 같은 다른 연성 문제와 방법론(강/약결합, 분할/모놀리식)은 공유하지만, 유체 도메인이 변형되며 격자가 움직인다는 점이 FSI를 유독 까다롭게 만든다. ↩
-
단방향은 사실 “연성”이라 부르기 민망할 정도로 반쪽짜리다. 그래도 뻣뻣한 구조물의 정적 하중 평가에는 충분하고 압도적으로 싸서, 초기 설계 단계에서 애용된다. 문제는 사람들이 양방향이 필요한 곳에도 단방향을 우겨넣는다는 것. ↩
-
강결합 subiteration이 수렴 안 하면 완화계수를 낮추는데, 너무 낮추면 수렴이 기어가고 너무 높으면 발산한다. Aitken의 동적 완화(dynamic under-relaxation)가 이 줄타기를 자동화해주는 대표적 기법이다. ↩
-
Causin, Gerbeau, Nobile(2005) 등이 이 불안정을 이론적으로 규명했다. 요지는 “유체/구조 밀도비가 임계값을 넘으면 명시적 분할 방식은 시간스텝과 무관하게 무조건 발산한다”는 것. 이걸 모르고 며칠을 태운 대학원생이 전 세계에 셀 수 없이 많다. ↩