Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

45 records · Page 3

Validation & Verification of CFD Models for Cryogenic Fluid Management of Propellant Tanks in Space

This article describes the need and the strategy for CFD model development, validation, and verification for Cryogenic Fluid Management (CFM) of Propellant Tanks in Space. It describes the type of CFD models that must be developed to address the future needs of Space CFM. It also discusses the two classes of experiments currently used to validate the fidelity of the CFD models. These two experiment classes are: (a) the small-scale simulant fluid science experiments that are equipped with scientific diagnostics to elucidate the underlying two-phase fluid physics of the CFM processes; and (b) the large-scale cryogenic experiments that assess the engineering performance of the propellant tank for storage and transfer. The current status of the CFD model development and validation is briefly assessed by presenting examples of segregated two-phase flow problems that have been successfully modeled. The future model development directions for CFM situations involving more complex interpenetrating phases are also defined.

Cryogenic Fluid Management↗

Columbus IFHX Ammonia Leak Analysis

After the Columbus Moderate Temperature Loop (MTL) InterFace Heat eXchanger (IFHX) low temperature event of GMT 345-2013, NASA investigated relevant transient scenarios involving IFHX rupture after water freezing and subsequent thawing. NASA recommended development of a Fault Detection Isolation and Recovery (FDIR) plan that would, in the event of a heat exchanger freeze event, close the Water On/Off Valves (WOOVs) to isolate the heat exchanger and prevent ammonia from the external flow loops from spreading into the cabin. NASA performed a preliminary simplified analysis for the reference case of IFHX rupture, but for a deeper understanding TAS developed detailed SINDA-FLUINT models of the Columbus ITCS that were built and run through the SINAPS GUI. This allowed simulation of the ammonia leakage physics including the variation of environmental parameters, thus providing more accurate and specific input to the FDIR under development. The result was finalization of the IFHX WOOVs closure sequence and wait times to contain the ammonia propagation to Columbus and allow identification of the leaking IFHX. In addition, the analysis results provided reference pressure profiles to be used on console and by the Engineering as support for the telemetry data assessment in case of failure.This paper gives an overview on the issue and focuses on the analytical aspects of the multiphase fluid dynamics involved.

Tilloca, A.↗

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.

Computational Fluid Dynamics↗

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗