⚡ electrical engineering

CFD Model Predictions of Zero Boil‑Off Unvented and Vented Tank Filling in Microgravity

본 논문은 충전율, 가이저 유동 영역, 그리고 압력 진화 사이의 관계를 규명하기 위해 미세중력 탱크 충전 실험을 통해 검증된 2D 및 3D CFD 모델을 제시하며, 궁극적으로 안정적인 가이저가 더 높은 피크 압력을 생성한다는 점과 흔들리는 가이저 거동을 포착하는 데 3D 모델링이 필수적임을 입증한다.

원저자: Rebecca L. Winter, Mohammad Kassemi

게시일 2026-07-03

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Rebecca L. Winter, Mohammad Kassemi

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

우주 공간에 떠 있는 거대하고 빈 물풍선을 채우려고 한다고 상상해 보세요. 지구에서는 중력이 물을 아래로 끌어당겨서 물이 바닥에 고입니다. 하지만 우주에는 "아래"라는 개념이 없습니다. 만약 당신이 이 떠 있는 풍선 안으로 차가운 액체 줄기를 쏜다면, 그것은 그냥 가라앉는 것이 아니라 통제 불능의 예측 불가능한 분수처럼 움직일 것입니다.

이 논문은 바로 그 문제, 즉 우주에서 폭발하거나 증발하지 않도록 초저온 연료 탱크를 채우는 방법에 관한 컴퓨터 시뮬레이션 연구입니다.

연구진이 발견한 내용을 쉬운 비유를 들어 정리해 드립니다.

핵심 문제: "간헐천(Geyser)" 효과

차가운 액체를 따뜻하고 빈 탱크 안으로 쏘아 넣으면, 액체는 그냥 가만히 있지 않습니다. 마치 간헐천처럼 위로 솟구칩니다.

문제점: 차가운 액체가 탱크 내부의 따뜻한 공기와 부딪힐 때, 공기는 응축(액체로 변함)되려 하고, 동시에 액체는 끓으려고 합니다. 이 과정에서 일종의 줄다리기가 발생하며 탱크 내부의 압력을 변화시킵니다.
위험성: 압력이 너무 빠르게 높아지면 탱크가 파손될 수 있습니다. 탱크를 보호하기 위해 가스를 배출(벤팅)하려고 하면 귀중한 연료를 잃게 됩니다. 목표는 가스를 전혀 배출하지 않으면서(제로 보일 오프, Zero Boil-Off) 탱크가 터지지 않게 탱크를 채우는 것입니다.

액체의 네 가지 "춤 동작"

연구진은 슈퍼컴퓨터를 사용하여 액체가 어떻게 움직이는지 관찰했습니다. 그들은 액체를 쏘는 속도에 따라 액체가 네 가지 뚜렷한 "춤 동작"을 한다는 것을 발견했습니다.

낮은 안정적 간헐천: 액체가 약간 솟구친 뒤 그대로 유지됩니다. 마치 잔잔한 분수처럼 안정적이고 예측 가능합니다.
흔들리는 간헐천: 액체가 솟구치긴 하지만, 해파리나 취객처럼 좌우로 기울거나 흔들립니다. 여전히 안정적이긴 하지만 중심에서 벗어나 움직입니다.
물방울 분리: 액체가 솟구치다가 높이가 높아지면, 마치 정원용 호스를 너무 세게 틀었을 때처럼 개별 물방울로 부서집니다. 이 물방울들은 탱크 반대편으로 날아가 부딪힙니다.
전체 제트 기둥: 액체가 끊기지 않는 하나의 단단한 줄기가 되어 레이저 빔처럼 탱크 반대쪽 벽까지 가로질러 쏩니다.

놀라운 발견: 느린 것이 항상 안전한 것은 아니다

"탱크를 천천히 채우면 압력이 낮게 유지될 것"이라고 생각할 수도 있습니다. 하지만 컴퓨터 모델은 이것이 항상 사실은 아니라고 보여주었습니다.

비유: 정원용 호스로 양동이를 채운다고 상상해 보세요. 물을 아주 미세하게(느린 속도로) 틀면, 물이 곧게 뻗어 나가 바닥에 부딪히며 큰 열을 발생하는 물보라를 일으킬 수 있습니다. 이는 실제로 시작 단계에서 엄청난 압력 급상승을 유발할 수 있습니다.
실제 상황: 물을 조금 더 빠르게 틀면, 액체 줄기가 물방울로 부서질 수 있습니다. 이 물방울들은 공기와 더 잘 섞여서 공기를 더 효율적으로 냉각시킵니다. 이는 결과적으로 더 낮은 압력 급상승을 만들어낼 수 있습니다.
교훈: 채우는 속도를 늦추는 것이 반드시 안전을 보장하는 것은 아닙니다. 때로는 액체의 행동 방식이 달라지기 때문에 오히려 압력 급상승을 악화시킬 수도 있습니다.

"마법의 숫자" (웨버 수, Weber Number)

연구진은 **웨버 수(Weber Number)**라는 수학적 개념을 사용하여 액체가 어떤 "춤 동작"을 할지 예측하는 방법을 찾아냈습니다.

이것은 액체의 속도 제한 표지판과 같습니다.
숫자가 1.5 미만이면, 액체는 안정적입니다 (안정적 간헐천).
숫자가 1.7 초과이면, 액체는 격렬하게 움직입니다 (물방울 또는 전체 제트).
회색 지대: 1.5와 1.7 사이에서는 결과가 불확실합니다. 액체는 안정적일 수도 있고 격렬할 수도 있으며, 컴퓨터조차 속도만 보고는 이를 완벽하게 예측할 수 없었습니다. 이는 다른 숨겨진 요인들이 작용하고 있음을 시사합니다.

2D vs 3D: 평면인가 실제인가

연구진은 두 가지 방식으로 모델링을 시도했습니다.

2D (평면): 탱크의 그림자를 보는 것과 같습니다. 계산은 빠르지만 "흔들리는" 춤 동작을 놓칩니다. 2D 모델은 액체가 완벽하게 직선을 유지한다고 가정하는데, 이는 우주에서는 사실이 아닙니다.
3D (실제): 실제 탱크를 보는 것과 같습니다. 흔들리는 간헐천과 날아다니는 물방울을 관찰하기 위해서는 3D 모델이 필수적입니다. 2D 모델은 액체의 흔들림을 볼 수 없기 때문에 압력에 대해 잘못된 답을 내놓기도 합니다.

벤트(Vents)에 대하여

그들은 또한 밸브를 열어 가스를 배출하는 경우(Vented)를 살펴보았습니다.

결과: 가스를 배출하면 압력이 낮게 유지되고 온도가 일정하게 유지됩니다. 액체는 밀폐된 탱크와 거의 비슷하게 행동하지만, 위험한 압력 급상승은 발생하지 않습니다.
주의점: 가스를 배출하지 않으면 압력이 상승하고 온도가 변하며 액체의 행동 방식도 달라집니다. 따라서 "배출형(vented)" 시뮬레이션으로 "밀폐형(sealed)" 탱크에서 일어날 일을 예측할 수는 없습니다.

진공 상태의 경우

그들은 거의 진공 상태인(공기가 없는) 탱크를 채우는 것도 시뮬레이션했습니다.

시작: 액체가 진공에 처음 닿으면 즉시 격렬하게 끓어 넘치며(flash boiling) 거대한 구름을 형성합니다.
안정화: 약 30초가 지나면 압력이 충분히 쌓여서, 액체는 끓기를 멈추고 앞서 설명한 일반적인 "간헐천" 사례처럼 행동하기 시작합니다.

결론

우주에서 연료를 잃지 않고 탱크가 터지는 일 없이 안전하게 연료 탱크를 채우려면, 단순히 속도를 짐작해서는 안 됩니다. 액체가 어떻게 "춤을 출지"(간헐천, 흔들림, 또는 물방울 분리)를 정확히 이해해야 합니다. 때로는 속도를 늦추는 것이 압력을 더 높게 만들 수도 있으며, 전체적인 모습을 보기 위해서는 3D 컴퓨터 모델이 필요합니다. "마법의 숫자"(웨버 수)는 행동을 예측하는 데 도움을 주지만, 여로 예측하기 까다로운 구간이 여전히 존재합니다.

기술 요약: 미세 중력 환경에서의 제로 보일-오프(Zero Boil-Off) 무벤트 및 벤트 탱크 충전의 CFD 모델 예측

문제 정의
미세 중력 환경에서 극저온 액체 연료를 안전하고 효율적으로 저장하고 이송하는 것은 NASA의 향후 장기 우주 임무를 위한 핵심 과제입니다. 궤도 내 재보급은 달 및 화성 탐사 임무에 필수적이지만, 중력의 부재는 유체 관리를 복잡하게 만듭니다. 주요 우려 사항은 충전 중 탱크 내부의 압력 상승이며, 이는 열 누출이나 충전 과정 자체로 인해 발생할 수 있으며, 결과적으로 파이프 파손이나 연료 자가 발화로 이어질 수 있습니다. 벤팅(venting)은 압력 상승을 방지하지만 전달 효율을 감소시킵니다. 따라서 "무벤트(no-vent)"(무벤트) 충전은 선호되는 방식이지만, 미세 중력에서의 2상 유동(two-phase flow), 열 전달 및 상변화 역학에 대한 깊은 이해를 요구하는 복잡한 운용 모드입니다. 본 연구는 무벤트 및 벤트 탱크 충전 시 미세 중력에서의 과도적 압력 진화와 액체-기체 계면 거동(특히 "가이저(geyser)" 형성)을 이해하기 위한 예측 모델의 필요성을 다룹니다.

방법론
저자들은 열역학적 분석과 전산유체역학(CFD)을 결합한 다각적인 모델링 접근 방식을 채택하였습니다:

열역학 모델: 최종 정상 상태 압력을 예측하기 위해 0차원(0D) 균질 열역학 모델을 개발했습니다. 이 모델은 유체가 포화 온도 및 압력에서 균질하다고 가정하며, 시스템의 이론적 한계를 설정하기 위해 질량 및 에너지 보존 방정식을 해결합니다.
CFD 모델: ANSYS Fluent를 사용하여 VOF(Volume of Fluid) 기법, 명시적 정식화(explicit formulation), 그리고 날카로운 계면 추적(sharp interface tracking)을 적용한 2상 CFD 시뮬레이션을 수행했습니다.
- 차원: 2D 축대칭 모델과 3D 파라메트릭 모델을 모두 사용했습니다.
- 유체: PNP(Perfluoro-n-pentane)를 시뮬런트 유체로 사용하였으며, 물성은 253.18 K에서의 DIPPR 데이터베이스로부터 도출되었습니다.
- 조건: 미세 중력 조건( $9.81 \times 10^{-6}$ m/s²) 하의 원통형 탱크(높이 0.094 m, 반경 0.03 m)에서 1.0–1.65 mL/s의 충전율을 대상으로 시뮬레이션을 진행했습니다.
- 시나리오: 무벤트(no-vent) 및 벤트(측면 및 상단) 구성을 조사하였고, 다양한 입구 속도 프로파일(평평한 플러그 유동 및 포아죄에 유동)과 다양한 초기 탱크 조건(빈 탱크, 부분 사전 충전, 근진공 상태)을 검토했습니다.
- 검증: 솔루션의 정확성을 보장하기 위해 격자 독립성 및 타임 스텝 연구를 수행했습니다. CFD 결과는 드롭 타워 연구(예: Symons et al., Govindan and Dreyer)의 기존 실험 데이터와 비교하여 검증되었습니다.

주요 기여 및 결과

가이저 거동 분류: 3D CFD 모델은 입구 속도와 유동 프로파일에 따라 네 가지 뚜렷한 가이저 거동을 식별했습니다:
1. 저안정 가이저(Low Stable Geyser): 꾸준히 성장하는 작은 액체 둔덕 형태.
2. 흔들리는 가이저(Wobbly Geyser): 안정적이지만 불규칙하게 기울거나 휘어지는 액체 기둥.
3. 액적 탈락(Droplet Shedding): 액체 기둥이 끊어져 액적으로 분리되는 불안정한 영역.
4. 전체 제트 기둥(Full Jet Column): 액체가 끊어지지 않고 반대편 탱크 벽에 도달하는 불안정한 영역.
  본 연구는 2D 축대칭 모델이 비대칭적 특성 때문에 "흔들리는 가이저" 거동을 포착하지 못하고, 종종 비물리적인 목 조임(necking) 현상을 예측한다는 점을 밝혀냈습니다.
압력 진화 및 과도적 피크: 중요한 발견은 가이저 안정성과 과도적 압력 사이의 관계입니다. 낮은 충전율이 항상 낮은 압력을 생성할 것이라는 가정과 달리, 모델은 안정적인 가이저 케이스(저안정 및 흔들리는 가이저)가 불안정한 케이스보다 더 높은 과도적 압력 피크를 생성함을 보여주었습니다. 이는 안정적인 가이저의 더 작은 계면 면적이 증기 응축 효율을 감소시키기 때문입니다. 불안정한 가이저(액적 탈락/제트)는 더 넓은 표면적을 제공하여 빠른 응축을 촉진하고 압력 급증을 완화합니다.
웨버 수(Weber Number) 임계값: 본 연구는 웨버 수(입구 속도 프로파일에 따라 조정됨)가 안정적 영역과 불안정 영역 사이의 전이를 효과적으로 예측함을 확인했습니다. 1.5와 1.7 사이의 임계 임계값이 확인되었으며, 이는 이전의 실험적 결과와 일치합니다. 그러나 저자들은 이 범위 내의 불일치성을 언급하며, 웨버 수만으로는 저안정 가이저와 흔들리는 가이저를 완전히 구분할 수 없음을 지적했고, 이는 다른 변수들이 안정성에 영향을 미칠 수 있음을 시사합니다.
열역학 vs. CFD 일치성: 0D 열역학 모델은 최종 정상 상태 압력을 성공적으로 예측(열역학적 한계로 수렴)했지만, CFD 모델에서 관찰된 초기 과도적 압력 오버슈트(overshoot)는 포착하지 못했습니다. 또한 열역학 모델은 낮은 충전율이 항상 더 낮은 초기 압력을 초래할 것이라고 잘못 예측하여, 가이저에 의한 압력 급증이라는 미묘한 차이를 놓쳤습니다.
벤트 vs. 무벤트 및 초기 조건:
- 벤트 충전: 벤트된 경우의 가이저 형상은 무벤트인 경우와 거의 동일했으나, 벤트는 압력 상승과 온도 상승을 방지했습니다. 측면 벤트는 등온 조건을 형성한 반 반면, 상단 벤트는 액체에 의한 부분적 폐쇄로 인해 일부 온도 상승를 허용했습니다.
- 초기 조건: 부분적으로 사전 충전된 탱크는 빈 탱크에 비해 온도 상승이 적었습니다. 근진공 충전은 초기 플래시 증발 단계를 보였으나, 약 30초 후에는 가이저 거동이 표준 충전 사례와 유사하게 수렴했습니다.
- 중력 효과: 미세 중력과 무중력(0g) 간의 비교에서 가이저 거동에 눈에 띄는 차이가 없었으며, 이는 드롭 타워(0g) 데이터를 저지구 궤도 및 심우주 임무 성능 예측에 사용하는 것이 타당함을 입증합니다.

의의
본 논문은 특정 가이저 거동을 이해하는 것이 단순화된 열역학 모델에서는 간과되기 쉬운 요소인, 무벤트 탱크 충전 시의 과도적 압력 진화를 정확하게 예측하는 데 매우 중요하다는 점을 주장합니다. 연구는 충전율을 낮추는 것이 항상 과압 방지를 위한 효과적인 전략은 아니며, 오히려 시스템을 더 높은 압력 피크를 갖는 안정적인 가이저 영역으로 전환시킬 수 있다는 점을 강조합니다. CFD 모델을 실험 데이터와 대조하여 검증하고, 비대칭 유동에 대한 2D 축대칭 근사의 한계를 확립함으로써, 본 연구는 미래 우주 임무를 위한 극저온 유체 관리 시스템을 설계하기 위한 더욱 견고한 프레임워크를 제공합니다. 이 작업은 향후 궤도 내 재보급 작업의 안전성과 효율성을 보장하기 위해 상세한 2상 유동 물리학을 포착해야 할 필요성을 강조합니다.