Search NASASearch

SEARCH · Search NASA

Results for “CFD 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.

At least 73 records · Page 4

Method for CFD Simulation of Propellant Slosh in a Spherical Tank

Propellant sloshing can impart unwanted disturbances to spacecraft, especially if the spacecraft controller is driving the system at the slosh frequency. This paper describes the work performed by the authors in simulating propellant slosh in a spherical tank using computational fluid dynamics (CFD). ANSYS-CFX is the CFD package used to perform the analysis. A 42 in spherical tank is studied with various fill fractions. Results are provided for the forces on the walls and the frequency of the slosh. Snapshots of slosh animation give a qualitative understanding of the propellant slosh. The results show that maximum slosh forces occur at a tank fill fraction of 0.4 and 0.6 due to the amount of mass participating in the slosh and the room available for sloshing to occur. The slosh frequency increases as the tank fill fraction increases.

Benson, David J.

Implementation of tritium transport in a gas-liquid contactor CFD simulation of tritium extraction from lead-lithium in ANSYS fluent

Modifications to the Computational Fluid Dynamic (CFD) software ANSYS Fluent were done to quantify and characterize tritium transport in Gas-Liquid Contactors (GLCs). A double-slit, Ergun-like equation was employed for the porous media model, with Ergun coefficients validated with Sulzer’s Sulcol software. Tritium transport from PbLi within the GLC was verified against analytical models. The geometry of the CFD model was based on the MELODIE GLC experiment. The hydrodynamic CFD pressure drop results align well with SulCol estimations and fall between the predictions of the analytical Delft-Olujic and Billet & Schultes models. In terms of mass transfer efficiency, traditional mass transfer models showed a significant deviation from experimental results when using varying values of H solubility in PbLi. A saturation phenomenon occurred when utilizing high solubility values for hydrogen in PbLi. In conclusion, a modified film theory mass transfer coefficient, incorporating either the Delft-Olujic or Billet & Schultes wettability model, yielded CFD-predicted extraction efficiencies that closely matched experimental measurements.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY

Cause and Cure - Deterioration in Accuracy of CFD Simulations with Use of High-Aspect-Ratio Triangular/Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD researchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where simplex elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identifies the reason behind the difficulties in use of such high-aspect ratio simplex elements is formulated using two different approaches and presented here. Drawing insights from the analysis, a potential solution to avoid that pitfall is also provided as part of this work. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, how the gradient evaluation procedures of the CESE framework can be effectively used to produce accurate and stable results on such high-aspect ratio simplex meshes is also showcased.

CESE Method

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian

CFD Simulations of the Supersonic Inflatable Aerodynamic Decelerator (SIAD) Ballistic Range Tests

A series of ballistic range tests were performed on a scaled model of the Supersonic Flight Demonstration Test (SFDT) intended to test the Supersonic Inflatable Aerodynamic Decelerator (SIAD) geometry. The purpose of these experiments were to provide aerodynamic coefficients of the vehicle to aid in mission and vehicle design. The experimental data spans the moderate Mach number range, $3.8-2.0$, with a total angle of attack ($alpha_T$) range, $10o-20o$. These conditions are intended to span the Mach-$alpha$ space for the majority of the SFDT experiment. In an effort to validate the predictive capabilities of Computational Fluid Dynamics (CFD) for free-flight aerodynamic behavior, numerical simulations of the ballistic range experiment are performed using the unstructured finite volume Navier-Stokes solver, US3D. Comparisons to raw vehicle attitude, and post-processed aerodynamic coefficients are made between simulated results and experimental data. The resulting comparisons for both raw model attitude and derived aerodynamic coefficients show good agreement with experimental results. Additionally, near body pressure field values for each trajectory simulated are investigated. Extracted surface and wake pressure data gives further insights into dynamic flow coupling leading to a potential mechanism for dynamic instability.

CFD modeling

Computational Fluid Dynamics (CFD) Simulations of Jet Mixing in Tanks of Different Scales

For long-duration in-space storage of cryogenic propellants, an axial jet mixer is one concept for controlling tank pressure and reducing thermal stratification. Extensive ground-test data from the 1960s to the present exist for tank diameters of 10 ft or less. The design of axial jet mixers for tanks on the order of 30 ft diameter, such as those planned for the Ares V Earth Departure Stage (EDS) LH2 tank, will require scaling of available experimental data from much smaller tanks, as well designing for microgravity effects. This study will assess the ability for Computational Fluid Dynamics (CFD) to handle a change of scale of this magnitude by performing simulations of existing ground-based axial jet mixing experiments at two tank sizes differing by a factor of ten. Simulations of several axial jet configurations for an Ares V scale EDS LH2 tank during low Earth orbit (LEO) coast are evaluated and selected results are also presented. Data from jet mixing experiments performed in the 1960s by General Dynamics with water at two tank sizes (1 and 10 ft diameter) are used to evaluate CFD accuracy. Jet nozzle diameters ranged from 0.032 to 0.25 in. for the 1 ft diameter tank experiments and from 0.625 to 0.875 in. for the 10 ft diameter tank experiments. Thermally stratified layers were created in both tanks prior to turning on the jet mixer. Jet mixer efficiency was determined by monitoring the temperatures on thermocouple rakes in the tanks to time when the stratified layer was mixed out. Dye was frequently injected into the stratified tank and its penetration recorded. There were no velocities or turbulence quantities available in the experimental data. A commercially available, time accurate, multi-dimensional CFD code with free surface tracking (FLOW-3D from Flow Science, Inc.) is used for the simulations presented. Comparisons are made between computed temperatures at various axial locations in the tank at different times and those observed experimentally. The affect of various modeling parameters on the agreement obtained are assessed.

Breisacher, Kevin

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank

CFD Simulations of Tiltrotor Configurations in Hover

Navier-Stokes computational fluid dynamics calculations are presented for isolated, half-span, and full-span V-22 tiltrotor hover configurations. These computational results extend the validity of CFD hover methodology beyond conventional rotorcraft applications to tiltrotor configurations. Computed steady-state, isolated rotor performance agrees well with experimental measurements, showing little sensitivity to grid resolution. However, blade-vortex interaction flowfield details are sensitive to numerical dissipation and are more difficult to model accurately. Time-dependent, dynamic, half- and full-span installed configurations show sensitivities in performance to the tiltrotor fountain flow. As such, the full-span configuration exhibits higher rotor performance and lower airframe download than the half-span configuration. Half-span rotor installation trends match available half-span data, and airframe downloads are reasonably well predicted. Overall, the CFD solutions provide a wealth of flowfield details that can be used to analyze and improve tiltrotor aerodynamic performance.

Potsdam, Mark a.

Correlating CFD Simulation with Wind Tunnel Test for the Full-Scale UH-60A Airloads Rotor

Data from the recent UH-60A Airloads Test in the National Full-Scale Aerodynamics Complex 40- by 80- Foot Wind Tunnel at NASA Ames Research Center are presented and compared to predictions computed by a loosely coupled Computational Fluid Dynamics (CFD)/Comprehensive analysis. Primary calculations model the rotor in free-air, but initial calculations are presented including a model of the tunnel test section. The conditions studied include a speed sweep at constant lift up to an advance ratio of 0.4 and a thrust sweep at constant speed into deep stall. Predictions show reasonable agreement with measurement for integrated performance indicators such as power and propulsive but occasionally deviate significantly. Detailed analysis of sectional airloads reveals good correlation in overall trends for normal force and pitching moment but pitching moment mean often differs. Chord force is frequently plagued by mean shifts and an overprediction of drag on the advancing side. Locations of significant aerodynamic phenomena are predicted accurately although the magnitude of individual events is often missed.

Romandr, Ethan

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad/Separation Bolt Wedge Tests

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Each panel test article included a metallic separation bolt imbedded in Orion compression-pad and heatshield materials, resulting in a circular protuberance over a flat plate. The protuberances produce complex model flowfields, containing shock-shock and shock-boundary layer interactions, and multiple augmented heating regions on the test plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the nonequilibrium flowfield in the facility nozzle, test box, and flowfield over test articles, and comparisons with the measured calibration data.

CFD