Results for “solvers” — Page 73

NASA NTRSDate not supplied

Dynamic Stability Analysis Using Free-Flight CFD

Blunt body probes are the primary design for atmospheric entry. The blunted capsule shape, while desirable for the decreased forebody heating, typically experiences dynamic instability at low supersonic and transonic flight. At moderate to low supersonic conditions, the unsteady wake begins to adversely interact with the vehicle in free-flight resulting in attitude oscillations which can grow in time and result in tumbling. Understanding vehicle dynamics is thus important when considering vehicle entry design and guidance navigation and control (GN&C) considerations. Historically, dynamic stability has been carried out exclusively using experimental techniques. The wide range of experimental approaches provide reasonable coverage in Mach and dynamic pressure space for a given trajectory, however achieving flight similitude is often difficult or impossible. Advancements in Computational Fluid Dynamics (CFD) techniques allow for dynamic motion of the vehicle surface during the simulation with either forced/prescribed motion or free motion in response to fluid forces. Several examples of CFD software with the ability to simulate fluid dynamics and vehicle motion are available in literature. The Free-Flight CFD (FF-CFD) capability developed by NASA Ames within the US3D flow solver has shown extensive validation and verification (V&V) efforts from moderate to low supersonic flow across a range of vehicle architectures. The solver is able to simulate full rotational and translation motion resulting in full 6 degrees-of-freedom (DoF) motion. The validation and verification efforts the solver's ability to predict the motion of a ballistic range model as compared to experimental results. Similar to other computational capabilities in literature, FF-CFD relied on heritage data reduction methodologies and processes for reducing dynamic data into aerodynamic coefficients intended for use in design and trajectory codes. These heritage approaches, initially developed with the assumptions (and restrictions) of ground test facilities, resulted in the significant culling of the full 6-DoF dynamic data set to discrete values which are then fit using an aerodynamic model. Recent development of simulation and data reduction methodologies has resulted in an improved approach to FF-CFD simulations of vehicle dynamics which leverage the control of dynamic constraints within a simulation set-up. The presentation of this work will cover the constrained dynamic simulation approach for deriving dynamic coefficients for a Genesis capsule. The methodology will be outlined and a discussion of implications in motion restriction will be discussed. Finally, reconstructed vehicle trajectories using FF-CFD derived aerodynamic coefficients will be presented.

Joseph M Brock↗

NASA NTRSDate not supplied

Free-Flight CFD Simulations and Dynamic StabilityAnalysis of the Orion Crew Module

Dynamic stability analysis of the Multi-Purpose Crew Vehicle (MPCV) is performed us-ing the US3D flow solver and a plugin library to solve rigid-body dynamics, Free-Flight CFD(FF-CFD). Previous efforts investigated the free-flight behavior of higher speed, open-backaeroshell shapes. Three primary experimental sources are used to evaluate the predictivecapability of the FF-CFD solver in the low supersonic range (Mach≈1). First, the ballisticrange results obtained at the HFFAF facility at NASA Ames is used to verify the six degree-of-freedom (6-DoF) dynamic capability of the FF-CFD solver. Next, FF-CFD simulationsare preformed using restricted motion and the resultant trajectories are post-processed toobtain pitch damping coefficient as a function of angle-of-attack. The comparison of thepitch damping results obtained with one degree-of-freedom (1-DoF) FF-CFD compare wellwith experimental fits derived from ballistic range data. Finally, the atmospheric flight ca-pability of FF-CFD is compared to data from the Ascent Abort-2 (AA-2) flight experiment.Two simulations were performed using forced and forced-free flight to investigate surfacepressure predictive capability and free-flight aerodynamic performance through a varyingatmosphere at low speed (Mach 0.6-0.2). Surface pressure predicted with FF-CFD agreeswell with experimental trends, with slight over-prediction near the end of the trajectory.The total angle-of-attack for the free-flight portion agrees well with experimental data.

FFCFD↗

NASA NTRSDate not supplied

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

NASA NTRSDate not supplied

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

Chun Tang↗

NASA NTRSDate not supplied

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

NASA NTRSDate not supplied

Least Squares Reverse Time Migration (LSRTM) for Damage Imaging in Composite Laminates

A method for adapting least squares reverse time migration (LSRTM) for ultrasonic guided wave imaging of composite laminates is proposed in this paper. As composites become more widely used in fields such as the aerospace industry, the need for high-resolution imaging in structural health monitoring (SHM) and nondestructive evaluation (NDE) is also growing. For instance, delamination is a common problem in composite laminates, which has led to a certain degree of apprehension in the use of composite materials for load-bearing structures. Although the solver-based imaging techniques using conventional reverse time migration (RTM) methods illuminate damage with a wide range of damage-scattering effects, the resulting images do not fully define the damage regions due to the limited data acquisition aperture, sensor density, frequencies/wavelengths, and incompleteness of adjoint reconstruction. Previously, the LSRTM theory was derived by incorporating a least squares migration (LSM) algorithm into RTM, and benchmarked its high-resolution damage imaging performance for isotropic plates. To improve damage imaging in composite laminates, a process to create an ultrasonic guided wave-based LSRTM method for anisotropic materials is proposed in this paper. The derivation of the forward modeling operator and the adjoint operator is presented. Numerical case studies were conducted to show the improvement of LSRTM over RTM in mapping damage in composite plates. Multiple damage sites or damage with a complex shape were created in the numerical studies based on 1) Born approximation-based modeling or 2) reduction of the plate stiffness-based modeling. LSRTM begins by using a forward modeling process. The forward modeling accuracy for the flexural wave solver used in this study for LSRTM was also compared with the results from a three-dimensional rotated staggered grid (RSG) elastodynamic finite difference solver, previously open-sourced by NASA Langley Research Center. These benchmarking studies have shown that the LSRTM method is more effective in reducing artifacts, improving resolution, and enhancing damage model reconstruction compared to RTM in composite laminates.

Jiaze He↗

NASA NTRSDate not supplied

Comparison of Aerodynamic Analysis Tools Applied to a Propeller-Blown Wing

This paper describes initial results from an ongoing activity comparing two aerodynamic analysis tools: OVERFLOW, a Navier-Stokes flow solver, and FlightStream, a viscous surface-vorticity, unstructured panel method flow solver. Both tools are applied to analyses of an unblown wing and a propeller-blown wing, and then comparisons are performed between the predicted aerodynamic loads and moments from each tool. Computational performance and efficiencies of each tool and their usability in conceptual and preliminary air vehicle design are also discussed. A description of the workflow, underlying theoretical foundations, and best practices are provided. In comparing FlightStream and OVERFLOW, reasonable agreement between predicted performance is observed for flow regimes of interest at high angle of attack prior to stall, including blown wing aero-propulsive effects on lift augmentation and stall delay. In the post-stall regions, FlightStream appears to underestimate the effects of turbulent flow separation that contribute heavily to drag. However, as expected for a high-fidelity, time-accurate solver, the setup and computation time are substantially greater for OVERFLOW than for FlightStream.

Vivek Ahuja↗

NASA NTRSDate not supplied

Comparison of Aerodynamic Analysis Tools Applied to a Propeller-Blown Wing

This paper describes initial results from an ongoing activity comparing two aerodynamic analysis tools: OVERFLOW, a Navier-Stokes flow solver, and FlightStream, a viscous surface-vorticity, unstructured panel method flow solver. Both tools are applied to analyses of an unblown wing and a propeller-blown wing, and then comparisons are performed between the predicted aerodynamic loads and moments from each tool. Computational performance and efficiencies of each tool and their usability in conceptual and preliminary air vehicle design are also discussed. A description of the workflow, underlying theoretical foundations, and best practices are provided. In comparing FlightStream and OVERFLOW, reasonable agreement between predicted performance is observed for flow regimes of interest at high angle of attack prior to stall, including blown wing aero-propulsive effects on lift augmentation and stall delay. In the post-stall regions, FlightStream appears to underestimate the effects of turbulent flow separation that contribute heavily to drag. However, as expected for a high-fidelity, time-accurate solver, the setup and computation time are substantially greater for OVERFLOW than for FlightStream.

Vivek Ahuja↗

NASA NTRSDate not supplied

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗

NASA NTRS2022

Hypersonic Simulations with US3D using Unstructured Grids from Fidelity Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. This article outlines a workflow that applies best practices for structured and unstructured grids. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured flow solver. Results from the two codes are compared and they show excellent agreement. The unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

NASA NTRSDate not supplied

Numerical Simulations of A Conceptual MSR-EES Shoulder Recession

The present study demonstrates our in-house material response solver, Icarus's capability to simulate the ablative conditions, including pyrolysis effects due to the interaction between hypersonic boundary layers and the thermal protection system (TPS). A conceptual aeroshell shoulder design, which undergoes mission-relevant flow and material conditions, is selected for the demonstration purpose. LAURA, a structured flow solver, is used to solve flow around the shoulder at several trajectory points of a flight path. The aerothermal dataset obtained from LAURA is used to enforce boundary conditions on the aeroshell wall to simulate ablation processes. The two-layered material system is stacked with HEEET (TPS) and Aluminum (actual material). A set of base parameters, such as the angle of material orientation with respect to the flow direction, convective heat-transfer co-efficient, and aerothermal boundary conditions on the rear end of the HEEET material layer, is selected to conduct a base-case simulation. A small amount of recession was observed, indicating that the design should be fine to go through the selected flow and material conditions. Furthermore, the parameters mentioned above are individually altered and are observed to affect the shoulder design's recession compared to the base case. Verification of our Icarus setup was also carried out using a one-dimensional grid and FIAT, a one-dimensional material-response solver, to build credibility for our results. As the current work is an uncoupled fluid-material-response simulation, a local sharp mesh deformation in the recessed surface near the shoulder corner is observed.

Prakash Shrestha↗

NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

NASA NTRSDate not supplied

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

NASA NTRSDate not supplied

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

NASA NTRS2023

Simulations of Zero Boiloff, Densification, and Solidification in a Large-Scale Liquid Hydrogen Tank

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Colin Philip Mahony↗

NASA NTRSDate not supplied

Coupled Material Response Simulations of Dragonfly Capsule and DrEAM Reconstruction

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center and has been recently used to analyze the material response of the Dragonfly capsule for a variety of problems. Since the Dragonfly capsule will be instrumented in a similar manner to the Mars 2020 and MSL capsules, it is important to assess how our current simulation tools can aid in understanding Dragonfly Entry Aerosciences Measurements (DrEAM). In this presentation, the Ares multi-physics tool that couples Icarus to the flow and radiation physics of the aerothermal environment will be used to better understand how modeling sensitivities might affect environment reconstructions and our understanding of measured data. Ares couples Icarus to the US3D flow solver and NERO, a reduced-order, finite-volume radiation transport solver and uses a customized time-scale management to enable coupled simulations for a large portion of the Dragonfly trajectory. As a result of the coupling, a more accurate and consistent assessment of modelling sensitivities to aerothermal reconstruction can be assessed. For example, the radiometer is sensitive to the quantity of CN in the flow, which is a byproduct of the free-stream methane but also the pyrolysis gas products injected into the boundary layer due to ablation of the heatshield. Ares will be used to conduct a simulation of the full Dragonfly capsule for 50 seconds of the trajectory. The different aerothermal and material response models will be discussed and the key results will be presented in terms of the simulated thermocouple and radiometer measurements on the backshell.

Ablation↗

NASA NTRSDate not supplied

Validation of Gas-Liquid Sharp Interface Model in Loci-Stream for Propellant Tank Self-pressurization under Normal Gravity

Deep space missions require advances in cryogenic fluid management (CFM) for long term storage of propellants. One of the important phenomena is self-pressurization due to heat leakages into the tank. Managing self-pressurization is one of the key technologies for deep space exploration and long-term space missions. The complex interactions involving natural convection, thermal gradients, turbulence, and phase change near the gas-liquid interface cannot be modeled using reduced order or nodal analysis models, and 3-D CFD analyses are necessary to fully characterize the dynamics. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. In this paper, we discuss modeling approach for self-pressurization in situations where the liquid interface is static and validate NASA MSFC’s Loci-Stream CFD solver for predicting self-pressurization in a test tank in normal gravity. This tank was designed to capture quantifiable and accurate data for understanding various two-phase fluid phenomena and to validate modeling tools. Specifically, we demonstrate the capability of the Loci-Stream solver with a two-phase sharp-interface treatment to predict self-pressurization of a tank with an unperturbed gas-liquid interface. This validation highlights the reliability of our modeling approach and our solver to serve as a design and analysis tool for NASA’s CFM application needs.

cryogenic fluid management↗