Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD Modeling”

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 361 records · Page 20

A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP

The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.

42 ENGINEERING↗

Thermomechanical Stress and Creep-Fatigue Analysis of a High-Temperature Prototype Receiver for Heating Particles

This work presents a three-dimensional (3D) thermomechanical model of a prototype-scale enclosed light trapping solar receiver for heating particles. Results of the thermoelastic model are used to estimate receiver lifetime under maximum flux conditions. A computational fluid dynamics (CFD) model is first developed to predict the temperature fields in a multi-panel assembly under steady operating conditions. Solar flux distributions on the receiver are obtained from the software package SolTrace and applied to the 3D thermal model. The subsequent particle heating is captured through a simplified 1D energy balance. Panel reradiation is considered through a surface-to-surface radiation model and natural convection loss to the surrounding air is captured in a representative fluid domain surrounding the receiver. The resulting temperature fields from the CFD analysis are used as inputs for a thermoelastic mechanical model with representative boundary conditions. With the resultant temperature and stress fields, a creep-fatigue damage and lifetime analysis is performed using the linear damage accumulation (LDA) theory. The Manson-Coffin formula and Larson Miller correlation are used to calculate the fatigue and creep, respectively. A maximum damage (corresponding to a 30-year service life) is defined for design assessment. The model was first developed and verified in detail by comparing with published results in the literature (temperature and stress profiles and distributions, and creep/fatigue damage fractions) for tubular solar receivers with supercritical carbon dioxide as the working fluid. It was then implemented to model a planar-cavity receiver with various design parameters. Specifically, three different design geometries are considered, and the results show that a maximum temperature of approximately 1200 K could be reached for each design with the given incident solar flux, with the main difference being the distribution of these temperatures. Preliminary resulting stresses for the small-scale prototype without design optimization vary from 20 MPa to 250 MPa for each design, with the maximum stresses occurring on the front face and concave geometry on the side of the panel. In future work, the developed methodology shown here will be applied to analyze a full-scale (50-150 MWth) receiver.

concentrated solar power↗

Aerothermal Ground Testing of Flexible Thermal Protection Systems for Hypersonic Inflatable Aerodynamic Decelerators

Flexible TPS development involves ground testing and analysis necessary to characterize performance of the FTPS candidates prior to flight testing. This paper provides an overview of the analysis and ground testing efforts performed over the last year at the NASA Langley Research Center and in the Boeing Large-Core Arc Tunnel (LCAT). In the LCAT test series, material layups were subjected to aerothermal loads commensurate with peak re-entry conditions enveloping a range of HIAD mission trajectories. The FTPS layups were tested over a heat flux range from 20 to 50 W/cm with associated surface pressures of 3 to 8 kPa. To support the testing effort a significant redesign of the existing shear (wedge) model holder from previous testing efforts was undertaken to develop a new test technique for supporting and evaluating the FTPS in the high-temperature, arc jet flow. Since the FTPS test samples typically experience a geometry change during testing, computational fluid dynamic (CFD) models of the arc jet flow field and test model were developed to support the testing effort. The CFD results were used to help determine the test conditions experienced by the test samples as the surface geometry changes. This paper includes an overview of the Boeing LCAT facility, the general approach for testing FTPS, CFD analysis methodology and results, model holder design and test methodology, and selected thermal results of several FTPS layups.

Bruce, Walter E., III↗

Parametric Study of Panel PCM–Air Heat Exchanger Designs

Heat exchangers, devices for the transfer of heat between two or more working fluids, are extensively used in cooling applications and heating applications. Heat exchangers in buildings are typically components of space-conditioning systems, as well as of water-heating applications. Heat exchangers are also sometimes used in applications that require storage and release of energy at specific times. Phase change materials (PCMs) enhance these heat-exchange processes, given their ability to melt and solidify at a fixed range of temperatures, absorbing or releasing significant amounts of latent heat. Five different configurations of PCM–air heat exchangers for thermal control in buildings are analyzed in this work. The heat exchangers were fitted with PCM encapsulated in plastic and composite pouches of various shapes, and packaged in stackable panel layers. Three-dimensional computational fluid dynamics (CFD) modeling of coupled incompressible fluid and conjugate heat transfer were performed on the designs. The phase change process was numerically modelled using the apparent heat capacity method. Steady-state CFD simulations provided quantification of pressure drop as a function of air flow velocity. Transient simulation results describe the thermal evolution of PCM in the pouches, helping to determine the best performing configuration with respect to total thermal charging time.

10 SYNTHETIC FUELS↗

Quasi 1-D Study of Pulse Detonation Rocket Engine Blowdown Gasdynamics and Performance

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional designs, but represent a challenging modeling task. A quasi 1-D, finite-rate chemistry CFD model for a PDRE is described and implemented. A parametric study of the effect of blowdown pressure ratio on the performance of several different PDRE nozzle configurations is reported.

Morris, Christopher I.↗

Pulse Detonation Rocket Engine Research at NASA Marshall

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional designs, but represent a challenging modeling task. A quasi 1-D, finite-rate chemistry CFD model for a PDRE is described and implemented. A parametric study of the effect of blowdown pressure ratio on the performance of an optimized, fixed PDRE nozzle configuration is reported. The results are compared to a steady-state rocket system using similar modeling assumptions.

Morris, Christopher I.↗

Quasi-One-Dimensional Modeling of Pulse Detonation Rocket Engines

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional designs, but represent a challenging modeling task. A quasi 1-D, finite-rate chemistry CFD model for a PDRE & described and implemented. A parametric study of the effect of blowdown pressure ratio on the performance of an optimized, fixed PDRE nozzle configuration is reported. The results are compared to a steady-state rocket system using similar modeling assumptions.

Morris, Christopher I.↗

Predicting Orion Launch Abort Acoustics

Computational fluid dynamics (CFD) scale-resolving simulations were performed to help assess the Orion launch abort system’s vibro-acoustic environment in collaboration with the Orion Loads and Dynamics team at the Johnson Space Center (JSC). The primary objectives are to establish best practices and to validate CFD predictions for surface fluctuating pressure spectra with ground, wind tunnel, and flight test data. Excellent agreement is obtained between CFD predictions by the Cartesian adaptive mesh refinement (AMR) module of the Launch, Ascent, and Vehicle Aerodynamics (LAVA) solver and measurements from a ground test representative of the start of a pad abort scenario, and with measurements from a wind tunnel test case corresponding to a transonic, high total angle of attack ascent abort scenario. A simplified (static and at rest) simulation of a pad abort flight test shows good agreement with the flight data, whereas a more costly simulation where the vehicle is accelerating relative to the ground from rest to Mach 0.4 shows consistent over-prediction of the spectra. Simulation results for a supersonic ascent abort flight test are in reasonable agreement with the measurements, but show larger excursions in regions where non-linear interactions between the plume turbulence, shocks, and turbulent boundary layer are strongest – likely due to CFD modeling assumptions, the effect of which are discussed. Overall, the CFD results demonstrate that scale-resolving simulation methodologies like the one presented can help characterize full-scale launch abort vehicle near-field noise and can complement wind tunnel, ground and flight tests by providing trends and reducing uncertainty for scenarios that cannot or are not planned to be tested.

Orion↗

RSRM Chamber Pressure Oscillations: Transit Time Models and Unsteady CFD

Space Shuttle solid rocket motor low frequency internal pressure oscillations have been observed since early testing. The same type of oscillations also are present in the redesigned solid rocket motor (RSRM). The oscillations, which occur during RSRM burn, are predominantly at the first three motor cavity longitudinal acoustic mode frequencies. Broadband flow and combustion noise provide the energy to excite these modes at low levels throughout motor burn, however, at certain times during burn the fluctuating pressure amplitude increases significantly. The increased fluctuations at these times suggests an additional excitation mechanism. The RSRM has inhibitors on the propellant forward facing surface of each motor segment. The inhibitors are in a slot at the segment field joints to prevent burning at that surface. The aft facing segment surface at a field joint slot burns and forms a cavity of time varying size. Initially the inhibitor is recessed in the field joint cavity. As propellant burns away the inhibitor begins to protrude into the bore flow. Two mechanisms (transit time models) that are considered potential pressure oscillation excitations are cavity-edge tones, and inhibitor hole-tones. Estimates of frequency variation with time of longitudinal acoustic modes, cavity edge-tones, and hole-tones compare favorably with frequencies measured during motor hot firing. It is believed that the highest oscillation amplitudes occur when vortex shedding frequencies coincide with motor longitudinal acoustic modes. A time accurate computational fluid dynamic (CFD) analysis was made to replicate the observations from motor firings and to observe the transit time mechanisms in detail. FDNS is the flow solver used to detail the time varying aspects of the flow. The fluid is approximated as a single-phase ideal gas. The CFD model was an axisymmetric representation of the RSRM at 80 seconds into burn.Deformation of the inhibitors by the internal flow was determined through an iterative structural and CFD analysis. The analysis domain ended just upstream of the nozzle throat. This is an acoustic boundary condition that caused the motor to behave as a closed-open organ pipe. This differs from the RSRM which behaves like a closed-closed organ pipe. The unsteady CFD solution shows RSRM chamber pressure oscillations predominately at the longitudinal acoustic mode frequencies of a closed-open organ pipe. Vortex shedding in the joint cavities and at the inhibitors contribute disturbances to the flow at the second longitudinal acoustic mode frequency. Further studies are planned using an analysis domain that extends downstream of the nozzle throat.

Nesman, Tom↗

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA Iced Aerodynamics and Controls Current Research

This slide presentation reviews the state of current research in the area of aerodynamics and aircraft control with ice conditions by the Aviation Safety Program, part of the Integrated Resilient Aircraft Controls Project (IRAC). Included in the presentation is a overview of the modeling efforts. The objective of the modeling is to develop experimental and computational methods to model and predict aircraft response during adverse flight conditions, including icing. The Aircraft icing modeling efforts includes the Ice-Contaminated Aerodynamics Modeling, which examines the effects of ice contamination on aircraft aerodynamics, and CFD modeling of ice-contaminated aircraft aerodynamics, and Advanced Ice Accretion Process Modeling which examines the physics of ice accretion, and works on computational modeling of ice accretions. The IRAC testbed, a Generic Transport Model (GTM) and its use in the investigation of the effects of icing on its aerodynamics is also reviewed. This has led to a more thorough understanding and models, both theoretical and empirical of icing physics and ice accretion for airframes, advanced 3D ice accretion prediction codes, CFD methods for iced aerodynamics and better understanding of aircraft iced aerodynamics and its effects on control surface effectiveness.

Addy, Gene↗

Numerical Prediction of Heat Transfer Coefficients during Xenon Tank Fill in Microgravity

NASA’s Gateway will serve as an orbital outpost to enable sustained human lunar exploration and facilitate access to destinations beyond the Moon. To maintain sustainable deep space missions, on-orbit propellant transfer and refueling must be realized. One of Gateway’s critical components, the Power and Propulsion Element (PPE), will use solar electric propulsion to perform attitude control maneuvers and orbit transfers. Xenon will be used as the propellant for the PPE and stored in a composite overwrapped pressure vessel (COPV). During the propellant refueling process, xenon will compress causing the fluid and tank walls to warm. There is concern that the bond between the tank liner and composite may degrade at elevated temperatures, along with the possibility of the tank exceeding design pressure. While on-orbit testing is being proposed, numerical models are being used to predict the thermal response and reduce risk during a refueling operation in microgravity where natural convection is diminished. Accurate numerical models may be used to inform decisions for follow-on testing, or design and operation of refueling architectures. Challenges of modeling xenon tank fill in microgravity include non-ideal gas behavior, operating near the fluid’s critical point, and lack of experimental flight data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict averaged heat transfer coefficients during tank filling. At the anticipated injection flow rates, the CFD results showed that forced convection dominates in microgravity. Initial results from a multi-node Thermal Desktop model that assumed a purely conducting fluid (Nu=1) over-predicted the CFD fluid temperature and pressure, while a similar model allowing for natural convection under-predicted. The heat transfer coefficients predicted by CFD were implemented into the multi-node model, and the COPV thermal response from the two computational approaches were compared with good agreement.

Heat Transfer↗

Solar Thermal Upper Stage Liquid Hydrogen Pressure Control Testing and Analytical Modeling

The demonstration of a unique liquid hydrogen (LH2) storage and feed system concept for solar thermal upper stage was cooperatively accomplished by a Boeing/NASA Marshall Space Flight Center team. The strategy was to balance thermodynamic venting with the engine thrusting timeline during a representative 30-day mission, thereby, assuring no vent losses. Using a 2 cubic m (71 cubic ft) LH2 tank, proof-of-concept testing consisted of an engineering checkout followed by a 30-day mission simulation. The data were used to anchor a combination of standard analyses and computational fluid dynamics (CFD) modeling. Dependence on orbital testing has been incrementally reduced as CFD codes, combined with standard modeling, continue to be challenged with test data such as this.

Olsen, A. D.↗

Low embodied energy and carbon, high lifetime silicon boules via a combined chemical vapor deposition/float zone process

This work evaluates a new process route to making float zone (Fz)-quality silicon wafers using a combination of computational fluid dynamics (CFD) modeling and technoeconomic analysis. Our analysis finds that the new process competes with Czochralski (Cz)-grown wafers on a levelized cost of energy system level. The new process also decreases embodied energy and carbon of silicon photovoltaics (PV) by ~6x circumventing the energy-costly Siemens process used in polycrystalline silicon (poly-Si) production plants to generate feedstock for Fz and Cz boules. Instead of using poly-Si from the Siemens process to feed crystallization, the new process uses the high-purity, trichlorosilane (TCS) precursor gas to grow a poly-Si feed rod in-situ during a modified Fz1,2 boule growth process. The gas-to-boule float zone process enables opportunity to produce high-purity (low metals and oxygen content), uniformly doped single crystal silicon boules and wafers with high bulk lifetimes (τ bulk > 15 ms) to enable higher efficiency cells (>27 %) with fewer known degradation mechanisms than Czochralski (Cz)-grown wafers. These benefits reduce the levelized cost of electricity (LCOE) of PV-produced electricity. Here we show the results of our CFD and chemical modeling of the process to prove feasibility and economic viability.

14 SOLAR ENERGY↗

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

Combustion Model of Supersonic Rocket Exhausts in an Entrained Flow Enclosure

This paper describes the Computation Fluid Dynamics (CFD) model developed to simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used to study the plume-induced environment due to static firing test of the Taurus II launch vehicle. The finite rate chemistry is used to model the combustion process involving rocket propellant (RP 1) and liquid oxidizer (LOX). A similar chemical reacting model is also used to simulate the mixing of rocket plume and ambient air. The model provides detailed information on the gas concentration and other flow parameters within the enclosed region thus allowing different operating scenarios to be examined in an efficient manner. It is shown that the real gas influence is significant and yields better agreement with the theory.

Vu, Bruce↗

Combustion Model of Supersonic Rocket Exhausts in an Entrained Flow Enclosure

This paper describes the Computational Fluid Dynamics (CFD) model developed to simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used to study the plume-induced environment due to static firing tests of the Taurus-II launch vehicle. The finite-rate chemistry is used to model the combustion process involving rocket propellant (RP-1) and liquid oxidizer (LOX). A similar chemical reacting model is also used to simulate the mixing of rocket plume and ambient air. The model provides detailed information on the gas concentration and other flow parameters within the enclosed region, thus allowing different operating scenarios to be examined in an efficient manner. It is shown that the real gas influence is significant and yields better agreement with the theory.

Vu, Bruce T.↗

Measurement of a radial flow profile with eddy current flow meters and deep neural networks

Eddy current flow meters (ECFMs) measure flows of conductive fluids. Recent interest in ECFMs has increased due to applications in advanced nuclear reactors. ECFMs are well suited for such applications, as they can provide non-invasive measurements of flow in fluids that are often difficult to measure. Traditionally, ECFMs are operated using an alternating current at a single frequency, limiting ECFMs to measure average fluid velocities, blockages, or voids. Here, we expand the capabilities of ECFMs by measuring the fluid radial velocity profile of liquid mercury. To accomplish this, we made several ECFM sensitivity measurements at a range of frequencies. Different frequencies vary the electromagnetic skin depth of the device. By adjusting frequencies, we probed the fluid velocity at various radial locations and constructed a flow-velocity profile. The relationship between the ECFM measurements and velocity profile is nonlinear and requires solving an inverse problem. Using electromagnetic finite-element simulations to train a deep neural network (DNN), we created a model that provides a stable general relationship between the sensitivity measurements of an ECFM and the fluid velocity profile. Using ECFM measurements of liquid mercury, our DNN model calculates a flow profile that agrees well with computational fluid dynamics (CFD) simulations. This technique has potential to improve flow monitoring for optimization, safe operation of conductive fluid loops, and/or validating complex CFD models.

47 OTHER INSTRUMENTATION↗