Cavitation(공동현상)

Cavitation(공동현상)

공동 현상은 유체 흐름의 압력이 매우 낮거나 온도 상승으로 인해 유체 내에서 증기 및 / 또는 가스 버블이 빠르게 진화하여 포화 압력을 높입니다. 기포의 갑작스런 출현 (및 후속 붕괴)은 비압축성 유체 내에서 압력의 급격한 변화를 일으켜 심각한 기계적 손상을 일으킬 수 있습니다. 캐비테이션에 의해 유발된 힘은 1983 년 글렌 캐년 댐의 유출로에서 발생한 손상으로 볼 수 있듯이 며칠 만에 수 피트의 암석을 침식 할 가능성이 있습니다 (Lee and Hoopes, 1996).

또한, 캐비테이션은 고압 다이 캐스팅에서 발생할 수 있습니다. 다이의 수축 및 곡선을 통한 용융 합금의 빠른 이동은 압력 강하를 빠르게하여 후속 캐비테이션을 유발할 수 있습니다. 생성 된 증기 기포는 최종 주조에서 다공성을 야기하거나, 더 악화되어, 다이를 손상시켜 주조를 오염시키고 다이 수명을 감소시킬 수 있습니다. 이러한 이유로, 캐비테이션이 발생할 가능성이 있는 영역을 이해하는 것이 중요합니다. 물리적 실험을 통해 캐비테이션을 시작하고 시각화하는 것은 어렵고, 잠재적으로 피해를 주기 때문에 공정을 시뮬레이션하는 것이 바람직합니다.

실증 사례

  • 물 및 환경 구조 내에서 손상된 캐비테이션 시뮬레이션
  • 고압 다이캐스팅 중 캐비테이션을 시뮬레이션하여 다이 손상 및 캐스팅 다공성을 유발할 수 있습니다.
  • MEMS 장치 내에서 열 기포 형성 시뮬레이션
  • 열전달 표면의 비등 거동 예측
  • 캐비테이션 역학으로 인한 혼합 예측

Cavitation(공동현상) FLOW-3D모델링

FLOW-3D의 현재 캐비테이션 모델은 열 버블 제트(Thermal bubble jets) 및 MEMS 장치를 시뮬레이션하는 데 성공적으로 사용되었습니다. FLOW-3D는 “Active”또는 “Passive”모델 옵션을 제공합니다. 능동형(Active) 모델은 기포 영역을 열고 수동형 모델은 흐름을 통해 캐비테이션 기포의 존재를 추적하고 전진하지만 기포 영역의 형성은 시작하지 않습니다.

능동형(Active) 모델은 더 큰 캐비테이션 영역이 예상되고 유동장에 영향을 미치는 경우에 가장 적합한 반면, 수동형(Passive) 모델은 작은 기포의 짧은 모양이 예상되는 시뮬레이션에 가장 적합합니다. 에너지 전송의 능동 모델과 계산을 통해 위상 변화(Phase change)도 옵션입니다. 기포는 계면(Surface)에서의 증발 또는 응축으로 인해 추가로 팽창 또는 수축 될 수 있습니다.

해석 사례

아래의 결과는 8m/s의 진입 속도, 18°의 수렴 기울기 및 8°의 발산 기울기를 가진 벤투리(Venturi) 내의 캐비테이션을 보여줍니다. 캐비테이션의 과도 거동이 잘 모델링되었으며, 모델은 22ms의 실험 결과와 비교하여 17.4ms의 캐비테이션 사이클 주기(Cycle period)를 예측합니다 (Stutz and Reboud 1997).

물 탱크를 통과하는 고속 발사체를 시뮬레이트하여 발사체의 웨이크에서 발생하는 저압 영역에서 공동화 깃털(Cavitation jets)을 보여줍니다. 발사체의 초기 속도는 600m/s입니다. 아래는 탱크에서의 움직임과 후미 캐비테이션 유체의 해석 결과입니다. 캐비테이팅 플룸(Cavitating plume)의 반경은 발사체가 감속함에 따라 좁아집니다.

참고 문헌

  • Lee, W., Hoopes, J.A., 1996, Prediction of Cavitation Damage for Spillways, Journal of Hydraulic Engineering, 122(9): 481-488.
  • Plesset, M.S., Prosperetti, A., 1977, Bubble Dynamics and Cavitation, Annual Revue of Fluid Mech, 9: 145-185.
  • Rouse, H., 1946. Elementary Mechanics of Fluids, New York: Dover Publications, Inc.
  • Stutz, B., Reboud, J.L., 1997, Experiments on unsteady cavitation, Experiments in Fluids, 22: 191-198.

[FLOW-3D 물리모델]Granular Flow / 입상유동

 Granular Flow / 입상유동

입상유동은 고상입자와 기체나 액체(예를 들어 모래와 공기 또는 모래와 물)인 유체와의 혼합물이다. 입상고체와 유체의 혼합물은 자유표면 경계에 의해 경계가 정해질수 있는 비압축성유체로 간주된다. 혼합 유체에서의 밀도변화는 초기에 존재할 수 있고 Drift-Flux 모델을 사용하여 계산되는 고체와 유체의 상대속도 때문에 유동중에 발생할 수도 있다. 자유표면에서의 가스의 방출은 고체가 입상간의 가스를 밀어내며 단단해질 때 발생할수있다. 액체의경우 고상이 단단해질 때 순수액체지역이 형성될수있다,

이 모델을 활성화하기 위해 General One fluid option Physics Granular Flow Granular Flow in Gas 또는 Physics Granular Flow Granular Flow in Liquid (Slurry)를 선택한다. 입상유동 창이 보이는데 여기서 입자의 직경 및 미시적 밀도와 유체의 밀도 및 점도가 정의되어야 한다. 필요하다면 고상의 최대 close packing 체적율 과 mechanical jamming 체적율이 각기 디폴트인 0.36과 0.61로부터 변경될 수 있다. 또한 자유표면의 an angle of repose(안식각?) 은 디폴트 값인 34도가 모델링하는 고상에 대해 맞지 않으면 변경될 수 있다.

입상유동을 선택하면 이는 자동으로 이 모델에서 사용되는 프로그램 내의 대 여섯 가지 모델을 활성화 시킨다. 또한 혼합물의 점도는 이 모델에서 계산되므로 정의하는 것이 불필요함에 주목한다. 사실 Fluids tab 에있는 어떤 유체물성도 정의할 필요가 없다.

입상물질이 격자 경계를 통해 계산 영역으로 들어오면 close packing 의 밀도보다 작은 고상율을 갖는 고상/기체 혼합물의 밀도를 정의하는 것이 중요하다. 그렇지 않으면 유동이 없을 것이다.

두 개의 보조 입력변수들이 있다. 하나는 Multiplier in threshold packing velocity 이고 다른 하나는 Multiplier in packing drag 이다. Multiplier in threshold packing velocity 는 이 속도 이상에서 packed solid material 의 유동을 결정하고 Multiplier in packing drag 는 입상이 충분히 높은 밀도로 packing될 때 유동을 정지시키는데 이용된다. 이 두 변수 모두 사용자가 입상체가 이 값을 변형시키는 응집력이나 다른 힘을 알지 못한다면 디폴트 값으로 남겨져야 한다.

또 다른 보조 입력 량은 마찰 각도인데 이는 보존각도보다 2~8도정도 크다. 마찰 각도는 액체인 유체의 경우에 중요하며 이 경우 마찰각도는 고상간의 충돌로 인한 전단유동 시 발생하는 분산압력에 영향을 미친다.

Granular flow application example: Core Blowing / 입상유동응용예제: 코어블로잉

코어블로잉 공정은 공기/모래 혼합물을 코어몰드에 고속 충진하는 것을 포함한다. FLOW-3D 는 코어블로잉을 각 모래 입자가 아닌 2상 연속체로 모델링 한다. 2상의 영향(공기/모래 결합)은 Drift Flux 모델을 사용하여 모델링 된다. 공기/모래 혼합물은 순수 공기와 선명한 경계면을 갖는 1유체로서 모델링 된다. 순수공기는 단열 기포로 나타내진다. 벤트는 밸브로 정의된다. 어떻게 이 모델이 실행되는지에 대한 더 많은 정보는 Flow Science Technical Note 88 at 테크니컬 노트notes/default.asp를 참조하라.

코어블로잉 모사(simulate)를 시작하는 단계는

  1. STL 파일로부터 관련 형상을 읽어 들여 생성하거나 Model Setup –> Meshing & Geometry 탭에있는 FLOW-3D 기초요소를 사용하여 형상을 생성한다.
  2. 다음 물리적 특성을 활성화하고 Model Setup Physics 탭에있는 변수들을 정의한다.

(a)   올바른 방향에서 중력을 정의하기 위해 Gravity and non-inertial reference frame 모델을 사용한다.

(b)   Viscosity and Turbulence 대화창에서 Viscosity and Laminar flow 를 활성화한다.

(c)    Activate the Granular Flow model.  Granular Flow 모델을 활성화한다.

  • Granular Flow in Gas 선택은 모래입자가 주위 매질보다 훨씬 밀도가 높다고 가정하는 Granular Flow 모델을 활성화한다.
  • Global vent 는 모래를 통과하는 공기의 전반적 배출을 조절한다. Global vent coefficient 는 모래 와 모래의 막힘에 의한 출구면적 감소에 따른 평균 손실을 나타내는 승수이다. 또한 모든 밸브의 외부압력과 모든 밸브 승수의 평균을 취한다. 추정치는 다음 식으로부터 계산될 수 있다.

여기서 Cv,g Global vent coefficient,  는 최대가능 고상율, L 은 공기 기포와 출구사이의 평균거리, 그리고 dgAverage grain diameter 이다.

  • Mechanical jamming volume fraction 은 모래의 체적율로 이 값 이상에서는 입상간의 상호작용에 의해 유동에 저항이 발생한다. 사용하기에 맞는 값은 0.61이다.
  • Close packing volume fraction 은 유동이 정지하게 되는 모래의 체적율을 기술한다. 체적율이 0.995(Close packing volume fraction) 를 넘게 되면 그 요소내의 속도는0으로 된다. 모래입자가 구형일 때 이는 일반적으로0.63이다.
  • Average grain diameter Grain density 는 정의되어야 하고 제조사로부터 알 수 있다. Gas density Gas viscosity 또한 정의되어야 한다. CGS 단위로 공기의 표준값은 각기 0.001225 g/cm3 와0.00017 poise 이다.
  • Multiplier in threshold packing velocity 와 Angle of repose 는 코어블로잉 모델링에는 필요하지 않다.
  • 입상 반발계수는 고체표면과 충돌 후에 모래입자가 유지하는 에너지의 양을 추정하는데 사용된다.

(d)    Density Evaluation 모델을 활성화한다. 일단 Granular Flow 가 활성화되면 First order approximation to density transport equation 이 자동적으로 가능하게 된다. 이는 모래의 전달을 계산하는데 필요하다. 더 나은 공간적 정확성을 위해 Second order monotonicity-preserving approximation to density transport equation 이 선택될 수 있다. 이는 모래의 농도가 급격히 변할 것으로 예측되는 모사(simulate)에 유용할 수 있다.

(e)   가스를 배출시키기 위해서는 Bubble and Phase Change 모델을 활성화시킨다. 이는 배출구와 밸브를 사용하기 위한 필요조건이다.

  1. 초기조건과 경계조건은 Meshing & Geometry 탭에서 추가될 수 있다. 공간 또는 기포영역의 초기조건은 이미 Adiabatic bubble 모델이 Bubble and phase change 모델에서 활성화될 때 정의된다. 경계조건은 Meshing Mesh Block 1 Boundaries 에서 정의된다. 모래는 공기압에 의해 코어상자를 통해 이동되므로 압력경계조건과 공기/모래 혼합물의 밀도가 경계에서 적용되어야 한다. S(대칭경계를뜻하는)를 갖는 적절한 경계상자를 택하면 경계대화상자가 나타날 것이다. Specified pressure 무선 버튼을 선택하고 입구압력, 유체율 1.0, 그리고 밀도를 정의한다.

  1. 단지 몇 개의 배출구만 있다면 밸브를 추가하거나 배출구가 너무 많아 수의 밸브로 추가할 수 없으면 Granular Flow Vent 로 정의된 형상을 사용한다. 밸브유동손실은 Bernoulli 의 차단 이론으로부터 유도된다. 밸브 생성에 관한 세부내용은 Valves 에서 찾아볼 수 있다. 배출구를 형상요소로 추가하기 위해 Meshing & Geometry 가지에서 별도 구성요소를 생성한다. 이는 배출구는 독자적 물성을 가지며, 형상요소는 그들의 물성과 운동에 따라 분류되어 있기 때문이다.

이렇게 모델링 될 때 배출구는 체적이 없다. 배출구가 같은 크기이면 이들은 하나의 STL 로써 또는 같은 구성요소의 기초요소를 사용하여 모델링 될 수 있다. 다른 크기라면 이들은 별도로 모델링 되어야 한다. 이들을 배출구로 정의하기 위해 Component Type drop down Granular Flow Vent 로부터 선택한다. 일단 형태가 정해지면 물성이 정의되어야 한다. Model Setup Meshing & Geometry Component Component Properties Granular Flow Vent Properties 에서 the Vent Flow Area, Diameter of Vent Channel 그리고 Vent External Pressure 를 정의한다.

See also: 또한 참조하라

  • Adiabatic Bubbles 단열기포
  • Flows with Density Variations 밀도 변화를 갖는 유동
  • Granular Flow. 입상유동

[FLOW-3D 물리모델]Condensation, Evaporation at Free Surfaces / 자유표면에서의 응축, 기화

Condensation/Evaporation at Free Surfaces자유표면에서의 응축/기화

1. Vaporization at Free Surfaces 자유표면에서의 기화

자유표면에서 발생하는 기화효과는 공간에서 정의된 일정 포화상태의 견지에서 모델링 될 수 있다. 이 모델을 활성화하기 위해 Physics>Bubble and phase change models>Constant pressure bubble with vaporization 를 선택한다. Fluids>Properties>Phase Change 에서의 Saturation Temperature 는 공간내의 기포의 포화상태를 정의한다. 기화 잠열은 Fluids>Phase change>Latent Heat of Vapor 에서 지정된다.

유체 에너지 방정식(열전달)은 이 모델(Physics>Heat Transfer)과 함께 해석되어야 한다. Fluids> Properties>Phase Change 에있는 Accommodation coefficient 에 양의 값을 정의한다. 자유 표면상의 액체의 온도가 포화 온도보다 높다면 액체는 다음과 같은 율로 증발할 것이다.

  • α 는 기화율을 조절하는 Accommodation coefficient이다. 이 값은 일반적으로 0.01에서0.1사이이며 1.0을 넘지 말아야 한다.
  • Hv 는 기화 잠열이다.
  • Asur 는 상변화를 위한 유효표면적이다.
  • kf 는 액체의 열전도도이다.
  • Tl 는 표면상 액체 온도이며
  • Tv1는일정한 기포 포화 온도이다
  • h 는 Prandtl 수로 정의된 표면에 있는 액체의 열전도에 대한 특정 길이이다.

여기서

  • xmin 는 (임의의 방향으로)계산 격자의 최소 셀 크기
  • Cv 는 일정 체적시의 기포 비열이며
  • µ1는 유체 #1의 점도이다.

각 표면 셀에서 기화하는 질량 유량은 후처리를 위해 저장되고 Analyze 에서 가시화될 수 있다.

기화는 자유 표면을 포함하는 셀들에서만 발생될 수 있다. 기포 포화온도는 일정 또는 변동압력을 갖는 모든 공간에 대해 일정하며 같다.

2. One Fluid with Thermal Bubbles 열기포를 갖는 하나의 유체

액체-증기 상변화에 의한 질량 전달은 열기포와 주위 액체 사이에 발생할 수 있다. 기포는 유체 #1 이 증기로 차 있다고 가정하고(즉, 기체 성분은 하나다.) 기포는 일정 압력, 온도, 그리고 밀도를 갖는다. 많은 기포 방울들이 있을 수 있고, 각 기포에서의 증기는 체적 변화와 열 및 질량 전달 때문에 고유한 시간에 따라 변하는 상을 갖는다. 유체 분율이0인 지정 압력의 격자 경계와 접하는 기포는 그 경계에서 정의된 기화 상태를 가질 것이다. 기화/응축모델은 Physics>Bubble and phase change models>Thermal bubbles with phase change 에서 활성화된다.

증기의 상태방정식은 이상 기체 방정식이며 절대 압력 P P = (γ − 1) · ρvapCvT 로부터 계산되는데 여기서

  • γ 는 1.285 ≤ γ ≤ 1.667값을 갖는 비열의 비율
  • T 는 절대온도
  • Cv 는 일정 체적에서의 증기의 비열
  • Cp 는 일정 압력에서의 증기의 비열
  • ρvap 는 기포 내의 증기 밀도

기포는 절대 단위로 이들의 초기 압력과 온도를 지정함으로써 초기화된다. 증기는 또한 Cavitation and Bubble Formation (Nucleation)에서 기술된 바와 같이 공동 또는 비등 과정을 통해 유체 내에서 생성될 수 있다. 증기 물성과 포화 곡선은 Fluids>Properties>Phase change 하위 메뉴에서 정의된다. 증기 압력은 사용자가 정의한 포화 곡선을 이용하여 그 지역의 유체 온도의 함수로써 계산된다. 디폴트 포화 곡선은 압력 P 와 온도 T 간의 Clausius-Clapeyron 관련식이다.

여기서

  • PV 1 TV 1(위의 물성치 목록에서 Saturation Pressure Saturation Temperature라고 쓰여있는) 는 포화곡선상의 한 점에서의 압력과 온도이다.
  • TEXPExponent for T-P Curve 로써 입력된다; 이의 값은 일반적으로
  • γ 는 증기의 비열 Gamma
  • Cv 는 일정 체적시의 기체 비열
  • Hv 는 기체의 잠열

형상 요소와 기포 내 증기간의 열전달은 Meshing & Geometry>Geometry>Component>Surface properties 의 component-void간의 열전달 계수에 의해 지정된다. 액체와 기포 내 증기와의 열전달도 마찬가지로 유체-void간의 열전달 계수에 의해 지정되어야 한다. 새로 생성된 증기기포는 heat transfer void type 1로 지정되는 것에 주목한다. Physics>Heat transfer>Fluid to solid heat transfer 가 증기 기포와 고체 요소간의 열전달을 가능하게 하기 위해 활성화되어야 한다.

상 변화는 계산 셀 내의 평균 유체 물성(밀도, 열에너지 그리고 액체분율)에 의존한다. 특히 액체와 증기의 온도는 한 요소에서 같으며, 표면의 얇은 유체막에서의 온도가 아니다. 이런 의미에서 상변화 모델은 현상학적이고 상변화율을 조절하기 위해 accommodation coefficient 의 조정이 필요하다. 1보다 큰 값은 사용되지 않아야 하는데, 이는 이 모델의 수렴이 힘들게 될 수도 있기 때문이다. 사실 일반적으로 사용되는 값들은 0.01과 0.1사이이다.

3. Two-fluid Model 두가지 유체 모델

이 모델은 증기 영역에서 모든 역학이 계산되는 것을 제외하고는 응축/기화 모델 (One Fluid with Thermal Bubbles)과 유사하다. 이 경우 압축 two-fluid 모델(비압축성 유체와 압축성 증기)은 경계면에서 발생하는 액체-증기 상변화가 가능하다. 순수 액체 지역에서의 핵 생성 또는 순수 증기 지역에서의 응축이 또한 가능하다. 유체 #1은 유체의 액상을 그리고 압축성 유체 #2(가스)는 증기를 기술한다. 표준 압축성 유동 모델에서와 같이 증기의 상태 방정식은 이상 기체 방정식, P = RF2 · ρ · T 이며 여기서.

  • RF2 는 증기의 기체상수
  • P 는 압력
  • ρ 는 기체 밀도
  • T 는 증기의 온도

two-fluid 상변화 모델은 Physics >Bubble and phase change models> Two-fluid phase change 에서 초기화되며, Fluids>Properties>Phase change 에서 양의 accommodation coefficient 를 필요로 한다. 상변화율은 직접적으로 accommodation coefficient 에 비례한다. 이 값은 절대적인 제한은 아니지만 일반적으로 0.01에서0.1사이이며 1.0을 넘지 말아야 한다. 증기 물성은 압축성 유체2의 물성으로 정의되며 증기 잠열과 포화곡선은 Fluids>Properties>Phase change 에서 정의된다. 포화 압력과 포화 온도로 정의되며 쌍으로 나타나는 압력-온도는 포화 곡선상의 한 점이어야 한다. T-P 곡선상의 지수는 온도-압력 포화관계의 지수이다. 디폴트 포화곡선은 압력 P 와 온도 T 간의 Clausius-Clapeyron 관련식이다.

여기서

  • PV 1 TV 1(위의 물성 목록에서 Saturation Pressure Saturation Temperature라고 쓰여있는)는 포화 곡선상의 한 점에서의 압력과 온도
  • TEXPExponent for T-P Curve 로써 입력된다; 이 값은 일반적으로 TV EXP = (γ − 1) CLHVCV 2 1
  •  Gamma 는 증기의 비열의 비율
  • CV 2 는 일정 체적시의 기체 비열
  • CLHV 1는 증기 잠열(단위질량당 에너지)

상변화는 유한 체적 내의 평균 유체 물성(밀도, 열에너지 그리고 액체분율)에 의존한다. 특히 액체와 증기의 온도는 한 요소에서 같다. 액체와 증기 경계면에서의 질량 전달율은 국부적 액체의 포화압력과 증기압사이의 차이에 의하여 계산된다.

Mass transfer rate

여기서

  • Pvap 는 증기압
  • Psat(T) 는 위에서 정의된 바와 같이 지역온도에서의 포화압력이다. 사용자가 필요에 따라 변경할 수 있는 subroutine PSAT.F에서 계산된다.
  • RSIZEAccommodation coefficient 이고 일반적으로 0.01과 0.1사이의 값이다.

액체와 증기경계에서 유체 질량의 단위면적당 상변화율이 계산되고, 후처리를 위해 Phase change mass flux 라고 불리는 공간변수로써 저장된다.
양의 값은 증발을 뜻한다:
음의 값은 응축.

액체 체적에서의 상변화는 Superheat temperature 를 지정함으로써 포화온도를 지나서까지 지연될 수 있다. 지역 포화온도보다 큰 Superheat temperature 의 값 때문에 증기 기포가 발생하기 전에 이 온도까지 유체 체적이 가열되는 것이 가능하다. 과열은 선택에따라 0이 아닌 벽의 거칠기를 사용함으로써 고체 벽 가까이에서 발생하지 않도록 할 수 있다.

4. Two Fluids with Non-condensable Gas / 비 응축가스를 갖는 Two Fluids

 

보통, 응축/기화 모델(two-fluid 모델)은 유체 #2가 완전히 액체의 증기상으로 이루어진다고 가정한다. 가스가 증기와 비응축가스(즉, 공기중의 수증기)의 혼합물로 구성되어 있는 경우에 Physics>Bubble and phase change>Two-fluid phase change>Noncondensable gas model 를 선택한다. two-fluid vapor 모델의 추가는 증기와 비응축가스의 기체상수들의 밀도 가중 평균 혼합물의 기체상수의 계산을 포함한다:

여기서

  • ρvap 는 계산된 거시적 증기밀도
  • ρnc 는 계산된 거시적 비응축 기체 밀도
  • RF2는 증기의 기체상수
  • RF 는 평균기체상수

그러므로, 압력은 P = RFρT 로 계산된다. 증기의 포화압력은 상변화(Two-fluid Model), 를 갖는 표준 Two-fluid 모델에서와 같은 방법으로 계산되지만, 질량 유량은 전체 가스압력을 사용하는 것과는 달리 증기의 부분압력을 이용하여 계산된다.

Mass transfer rate

여기서

  • Pvap 는 가스성 유체의 증기의 부분압력
  • Psat(T) 는 사용자가 필요에 따라 변경할 수 있는 subroutine PSAT.F에서 정의되는 Clausius-Clapeyron 방정식으로부터 계산되는 국부 온도에서의 포화압력이다.
  • RSIZEAccommodation coefficient 이고 일반적으로 0.01과 0.1사이의 값이다.

Accommodation coefficient 가 1.0의 값을 가진다면 모델은 한 시간단계에서 평형에 도달하기에 충분한 상변화를 예측하려고 시도할 것이다. 이 속도는 너무 급속해 실제 물리적조건과 비교될 수가 없다. 액체와 가스의 경계면의 경계층 내의 역학은 규모가 너무 작아 이 모델에 포함할 수 없으므로 FLOW-3D 가 정확히 이 계수 없이 상변화율을 예측하는 것은 불가능하다. .

이 모델을 이용하기 위해 Physics>Bubble and phase change models>Non-condensable gas model 의 체크상자를 선택한다. Gas constant Specific heat of the non-condensable gas 를 위한 값을 입력한다. 가스가 영역 경계에서 들어오는 곳에 각 mesh block 경계 조건 입력창에 있는 Non-condensable gas fraction 의 비응축가스의 체적율(0 과 1사이)을 지정한다. 비응축가스를 포함하는 초기 유체지역을 정의하기 위해 Meshing & Geometry>Initial>Global 를 지정한다. 이 양은 또한 각각의 초기유체 영역과 특정 지점에서 지정될 수 있다.

5. Vaporization Residue / 증발 잔류량

MAIN VARIABLES: SCALAR: IRESID, RMXSC
XPUT: IPHCHG

액체용제가 기화할 때 이에 포함되어 있는 용질은 더 농축된다. 마찬가지로 스칼라 농도변수로 모델링 된 용질도 유체문제의 자유표면에서 증발로 인해 자동적으로 농축될 것이다. 표면요소에 액체가 반보다 적게 있을 경우 농축변화가 표면요소의 두께의 반에 해당하는 지역으로 퍼져나가는 크기로 스칼라의 농축이 바로 주위의 표면요소에서도 또한 발생할 것이다.

 증발이 충분히 발생하고 용질의 농도가 커지면 표면에서 발생할 수도 있고 용질이 완전히 증발하면 표면상에 이의 잔류가 생성될 수 있다. 잔류형성은 Physics Bubbles and phase change 에서 활성화되는 Constant pressure bubbles with vaporization, 및 Thermal bubbles with phase change 모델과 함께 시뮬레이션 되어야 한다. 잔류모델은 IRESID = 1로 지정하고 용질 스칼라 ns, RMXSC(ns)를 최대 packing 밀도를 정의함으로써 활성화된다. 일단 용질이 최대 packing 밀도까지 농축되면 더 이상의 농축은 고정(움직이지 않는)된 잔류를 초래한다. 하나 이상의 스칼라 용질이 존재하면 잔류는 모든 용질 전체 잔류를 기록한다.

Note: 용질농도는 Physics Scalars 로부터 FLOW-3D‘s Scalars 모델을 이용하여 입력된다.