Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

INTEGRATE – Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements

The INTEGRATE (Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements) project developed a new inverse-design capability for the aerodynamic design of wind turbine rotors using invertible neural networks. Training data was obtained from improved turbulence and transition models for RANS and hybrid RANS/LES solvers with machine-learned physics-based data-augmented corrections and then using the resulting neural-network(s) augmented RANS model to run thousands of 2-D and 3-D CFD simulations.

17 WIND ENERGY↗
DOE OSTI2023

Rapid CFD Using Machine Learning Algorithms (CRADA Final Report)

This is a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (“LLNL”) and Guardian Glass, LLC ("Guardian Glass") to develop a fast-running emulator of the reactive Computational Fluid Dynamics (“CFD”) simulations needed to understand the complex reactions and flows in the glass melting, fining, and forming subprocesses. This CRADA project is sponsored under the High-Performance Computing for Manufacturing (“HPC4Mfg”) Program of the Department of Energy’s Advanced Manufacturing Office (“AMO”) within the Energy Efficiency and Renewable Energy (“EERE”) Office.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Plasma-Assisted Pre-Chamber Ignition System for Highly Dilute Stoichiometric Heavy-Duty Natural Gas Engines (Final Technical Report)

This project explored advanced ignition technologies to significantly enhance efficiency and reduce operating costs for heavy-duty natural gas engines operating at stoichiometric conditions, while meeting ultra-low NOx emission standards. The main goal was to develop and validate a plasma-assisted pre-chamber ignition system that could deliver at least a 2% increase in brake thermal efficiency (BTE) and a 4% decrease in total cost of ownership (TCO) compared to a typical multi-cylinder engine with three-way catalyst aftertreatment, ensuring compatibility with the expected 2027 EPA/CARB regulations. In the first half of the project, the research team concentrated on developing and testing plasma-assisted pre-chamber ignition using nanosecond pulsed discharges. Extensive experiments were conducted in an optically accessible rapid-compression and expansion machine, a constant-volume chamber, and an optical single-cylinder engine. Experiments were coupled with CFD simulations. The work produced unique insights into pre-chamber flame formation, jet ignition, dilution effects, and flame quenching at pressures, temperatures, and dilution levels relevant to engines. Although plasma-assisted ignition showed promise in controlled lab settings, the research also identified fundamental and practical challenges when applying this technology to real engine conditions. Midway through the project, a crucial pivot was made, guided by three key findings. First, the power electronics required for nanosecond plasma discharges were found to be too costly for commercial use, undermining the project’s cost-of-ownership goals. Second, nanosecond plasma ignition was highly sensitive to turbulent flow in the pre-chamber, resulting in lower ignition reliability than traditional spark under engine-like conditions. Third, achieving a truly diffuse low-temperature plasma at high pressures near top dead center was not possible, reducing the anticipated chemical enhancement benefits. These results collectively suggested that continuing with plasma-assisted ignition was unlikely to meet both efficiency and cost objectives. In response, the project shifted focus to a more realistic approach: enhancing traditional spark-based pre-chamber ignition with significantly less spark energy. Using insights gained earlier in the project, the team redesigned the pre-chamber to maintain high dilution tolerance and quick combustion, even with lower ignition energy. Testing confirmed that with optimized pre-chamber design and combustion timing, a lower-energy spark could reliably ignite highly diluted stoichiometric mixtures, reduce burn time, and boost thermal efficiency. Final engine testing and techno-economic analysis verified that this revised approach successfully achieved the project goals. The optimized pre-chamber ignition system provided over a 2% increase in calculated brake thermal efficiency compared to the baseline engine. Notably, the lower ignition energy and simplified hardware reduced component stress, extended maintenance intervals, and lowered the total cost of ownership. When used with stoichiometric operation and traditional three-way aftertreatment, the system remained compatible with near-zero NOx emissions targets without increasing cost or complexity in the emissions control system. In summary, although the project deviated from its initial plasma-assisted ignition idea, the work produced a more practical and commercially viable solution. The results show that precisely optimized, low-energy pre-chamber spark ignition can significantly improve efficiency and reduce overall ownership costs for heavy-duty natural gas engines. This directly aligns with DOE goals for cleaner, more efficient, and cost-effective transportation technologies.

03 NATURAL GAS↗
DOE OSTI2021

Parametric Evaluation of RDE Inlet Performance using a Reduced-Order Model

This paper presents the results of a parametric study on the effects of inlet restriction, along with a number of other controlling parameters, using a previously developed reduced-order RDE inlet model. The results of this study will be used to inform subsequent high-fidelity CFD simulations and/or selection of experimental operating conditions which may be conducive to maximizing performance.

Bedick, Clinton↗
DOE OSTI2022

Design-Space Exploration for Inverse-Design of Wind Turbine Blades Using Data-Driven Methods

The state-of-the-practice aerodynamic design methods for wind turbine blades is typically based on Blade-Element Momentum (BEM) theory using a pre-designed frozen family of airfoils. The airfoils are themselves typically designed using panel methods. The design of the next-generation of large flexible rotors will need to capture the non-linear aerodynamics and three-dimensional flow to reduce the levelized cost of wind energy. Data-driven methods for aerodynamic design using data generated by computational fluid dynamics offer an attractive alternative to BEM-based methods that captures the non-linear aerodynamics of the component airfoils as well as the root and tip sections. In this work, we develop and demonstrate a framework to "smartly" explore the relevant design space in combination with an appropriate automated CFD pipeline to evaluate the aerodynamics of each design. The design-space exploration framework uses appropriate perturbations to the airfoil shape and induction profile from a baseline shape in combination with the inverse-design using BEM. The resulting designs are evaluated using an automated CFD pipeline using the in-house CFD solver framework "Mercury". We perform verification and validation to establish the capability of the Mercury framework to predict the aerodynamic performance of wind turbines. The CFD simulation of the perturbed blade shapes are optimized to restart from the converged baseline simulation to reduce the computational time. Finally, we demonstrate the design-space exploration technique for the design of the outboard section and the full rotor using perturbations to the shape and operating conditions of the NREL 5-MW turbine.

aerodynamic design methods↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2023

Measurements and Prediction of Ash Deposition in a Cyclone-Fired Boiler Operating under Variable Load Conditions

Measurements of ash deposition rates were made between the secondary superheater and reheater sections of a 450 MW cyclone-fired lignite boiler as the operational load varied from 33 to 100%. Significant reductions in deposition rates with a decrease in operational load were observed. To uncover the causative mechanisms behind these observations, operational data from the power plant were used to carry out computational fluid dynamic (CFD) simulations of the boiler. After ascertaining that the gas temperatures and velocities at various sections within the boiler were being represented adequately, decoupled simulations of the ash deposition process on the deposit probe were carried out using a finely resolved boundary layer mesh. Fly ash particle size distribution (PSD) and its concentration for the decoupled calculations were determined from stand-alone cyclone barrel simulations. The ash partitioning (mass %) between the fly ash and slag was found to be ~50:50, which was in line with previous field observations, and it did not vary significantly across different cyclone loads. The predicted PSD of the deposit ash was concentrated in the size range 10–30 microns, which was in agreement with cross-sectional images of the deposit obtained from the measurements. At lower loads, sharp variations in the deposition rates were predicted in the gas temperature range 950–1150 K. The particle kinetic energy—particle viscosity-based capture methodology utilized in this study in conjunction with appropriate ash compositions, ash viscosity models and gas temperature estimates can help estimate slagging propensities at different loads reasonably well in these systems.

Krishnamoorthy, Gautham (ORCID:0000000255205092)↗
DOE OSTI2023

Ducted Fuel Injection Provides Consistently Lower Soot Emissions in Sweep to Full-Load Conditions

Earlier studies have proven how ducted fuel injection (DFI) substantially reduces soot for low- and mid-load conditions in heavy-duty engines, without significant adverse effects on other emissions. Nevertheless, no comprehensive DFI study exists showing soot reductions at high- and full-load conditions. This study investigated DFI in a single-cylinder, 1.7-L, optical engine from low- to full-load conditions with a low-net-carbon fuel consisting of 80% renewable diesel and 20% biodiesel. Over the tested load range, DFI reduced engine-out soot by 38.1–63.1% compared to conventional diesel combustion (CDC). This soot reduction occurred without significant detrimental effects on other emission types. Thus, DFI reduced the severity of the soot–NO x tradeoff at all tested conditions. While DFI delivered considerable soot reductions in thepresent study, previous DFI studies at low- and mid-load conditions delivered larger soot reductions (>90%) compared to CDC operation at the same conditions. Therefore, the DFI configuration used here has been deemed nonoptimal (in terms of parameters such as the injector-spray and piston geometries), and several improvements are recommended for future studies with high-load DFI. Further, these improvements include employing better spray-duct alignment, a deeper piston bowl with a smaller injector umbrella angle, and a fuel injector that opens and closes faster. The study also suggests future research to make DFI ready for commercialization, such as metal-engine tests to ensure desirable DFI performance over an engine’s complete speed/load map. Overall, this study supports the continued development and commercialization of DFI to meet upcoming emissions regulations for heavy-duty vehicles. Specifically, multicylinder engine experiments and CFD simulations should be utilized to optimize the performance and clarify the full potential of DFI.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗
DOE OSTI2026

Investigation of Flux Spreading in a Light-Trapping, Planar-Cavity Receiver for Enclosed Solar Particle Heating

Concentrating solar thermal power (CSP) technology development has recently focused on increasing the operating temperatures to accommodate high efficiency power cycles and thermochemical processes. Inert solid particles as heat transfer media enable solar receivers to operate above 700 degrees Celsius resulting in increased system thermal efficiency compared to the conventional molten salt based CSP system. An open-cavity falling-particle solar receiver that can efficiently heat particles by direct heating from concentrated solar radiation faces challenges with large particle losses from wind and unable to support thermochemical reactions. A light-trapping, planar cavity reiver (LTPCR) where particles are indirectly heated can significantly minimize the particle losses during the operation, support thermochemical reactions, and offer scalability potential. The LTPCR features an array of vertical planar receiver/absorber panels arranged within a cavity configuration. Concentrated solar radiation from heliostats is focused onto the receiver walls, where heat is indirectly transferred to solid particles flowing inside the receiver channels. Heat transfer occurs through direct contact between the receiver panel walls and particles, and can be enhanced by fluidizing particles with air. This fluidization increases particle-wall contact and extends particle residence time, maximizing heat transfer efficiency. The unique vertical planar receiver structure originated from a near-blackbody tubular light absorber, effectively distributing the incoming solar beam spread across the panel walls and trapping light. This flux spreading effect, driven by cosine projection, converts high incident solar flux into a lower, more uniform heat flux on the panel walls. This redistribution enhances heat transfer efficiency between particle-wall or reaction gases-wall, while preventing localized overheating of the receiver panel. Indirect planar cavity solar receivers completely separate solid particles from the ambient environment that can greatly reduce the thermal losses in heated particles resulting in high efficiency at high temperatures above 700 degrees Celsius. This design ensures no particle losses to the environment during the operation while open-cavity designs can experience significant particle losses from wind. An experimental investigation was conducted to observe flux spreading on the receiver panel wall. A lab-scale prototype planar receiver, fabricated using Haynes 230 alloy, was tested under direct concentrated solar radiation using the high-flux solar furnace (HFSF) facility at NREL. The experiment was performed under normal peak radiative heat fluxes ranging from 800 to 1900 kW/m2. A temperature distribution on the panel wall was measured using a thermal imaging camera (FLIR A 6600). To prevent overheating at the receiver front tip, prism-shaped heat shields (Zircar UNIFROM C1) were placed in front of the receiver, and their influence on flux spreading was also studied. Absorbed flux distribution on the panel wall was modeled using SolTrace. The total solar power and flux distributions delivered from HFSF were determined based on the heliostat mirror optical properties, direct normal irradiance (DNI) on the on-sun testing days, peak flux measurement during the on-sun testing, and shutter/attenuator settings Due to the large incident angles of the solar beam on the panel wall, the angular optical properties of Haynes 230 alloy and Zircar heat shields were incorporated into the model. This flux distribution model was then integrated into a computational fluid dynamics (CFD) simulation to predict the receiver panel wall temperature, which was compared with the experimental measurements. Both prediction and measurements identified a temperature hotspot at the backside of the panel, indicating that the incident solar beam can fully reach to the rear of the receiver. The heat shields positioned at the front of the receiver effectively reduced the excessive temperature rise at the receiver front tip. Overall, the temperature was well distributed over the panel wall, with a minor hotspot at the back of the receiver. The model slightly overpredicted the temperature, possibly due to discrepancies in optical properties of the panel and an underprediction of thermal loss in the receiver. The advancement of the particle LTPCR offers a viable alternative to open-cavity receivers by addressing particle loss issues. Additionally, it presents a pathway for enabling solar thermochemical processes, extending CSP technology beyond power generation to fuel and chemical production.

14 SOLAR ENERGY↗
NASA NTRS1990

Supersonic flameholding by attached oblique shock waves

The object of this investigation was to explore various regimes of shock-induced combustion that may occur as a result of supersonic flow of a reactive gas mixture past a simple compression ramp or wedge. Attention was confined to combinations of wedge angles, fuel-air mixture ratios and freestream velocities that permit the resulting oblique shock or oblique detonation wave to remain attached. Computational fluid dynamic (CFD) simulations, both inviscid and viscous, were performed and analyzed in order to investigate two modes of compression ramp flameholding; namely, complete oblique detonation wave (ODW) and shock-induced combustion. In addition, the possibility of flame persistence or hysteresis was explored. It was preliminarily concluded that inclusion of viscous effects may be essential to accurately predict the occurrence of the various flameholding modes.

Fort, James A.↗
NASA NTRS1989

A turbulent boundary layer at low Reynolds number with adverse pressure gradient

The evolution of a low Re(theta) turbulent boundary layer in an adverse pressure gradient (APG) is studied for comparison with CFD simulations by Spalart (1988). A short region of favorable pressure gradient is applied first to establish a self-preserving layer which is a suitable initial condition for the simulations. The APG is then applied rapidly such that Beta(x) of about 2 for Re(theta) of about 1500.

Watmuff, J. H.↗
NASA NTRS1991

Aerothermodynamic environment predictions in support of the Aeroassist Flight Experiment

A computational fluid dynamic (CFD) simulation has been made of the flow past an Aeroassist Flight Experiment (AFE) aerobrake model with a payload and sting. The simulation involves solving the complete Navier-Stokes equations in three dimensions. Comparisons with data taken in a Mach 10 wind tunnel test are made and agreement is shown to be very good. Points of comparison include the aerobrake surface pressure and heat transfer rate, the sting surface streamlines, and the sting heat transfer rate. The dependence of the solution on the grid is also explored. Finally, predictions are made for conditions which correspond to a hypersonic shock tunnel force and moment test to be carried out later this year.

Fay, John F.↗
NASA NTRS1992

Modeling of solid/porous wall boundary conditions for the validation of computational fluid dynamics codes

A computational study has been undertaken to investigate method of modeling solid and porous wall boundary conditions in computational fluid dynamics (CFD) codes. The procedure utilizes experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution is then imposed as a boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the boundary conditions at the wind tunnel walls using only sparse experimental measurements has been investigated. Position and refinement of experimental measurement locations required to describe porous wall boundary conditions has also been considered.

Beutner, Thomas J.↗
NASA NTRS1991

Aerobrake aerothermodynamic environment predictions in support of the Aeroassist Flight Experiment

A computational fluid dynamic (CFD) simulation of the flow past an Aeroassist Flight Experiment (AFE) aerobrake under wind tunnel test conditions has been made. The Navier-Stokes equations in three dimensions were solved numerically using a finite-volume, implicit approach. Comparisons with experimental data include surface pressures and heat transfer rates and aerodynamic coefficients. Agreement with experiment is shown to be excellent. The dependence of the solution on the computational grid is explored. The present work is preliminary to simulation of the flow past the AFE under flight conditions, which at present cannot be duplicated with ground-based experimental facilities.

Fay, John F.↗
NASA NTRS1992

A summary of the computational fluid dynamics analysis of the hypersonic flights of Pegasus

A computational fluid dynamics (CFD) analysis of the hypersonic flights of the NASA's Pegasus launch vehicle was carried out using the PARC code. The flight data included heat transfer rates and sparse pressure coefficients primarily in the fillet region of the vehicle. Since the CFD simulations did not account for ablation at the vehicle surface, the test heat transfer results were obtained for the low Mach number simulations; as the Mach number was increased, the ablative effects at the vehicle surface became more apparent and contributed to the error.

Fricker, Darren↗
NASA NTRS1992

Analysis of film cooling in rocket nozzles

Progress during the reporting period is summarized. Analysis of film cooling in rocket nozzles by computational fluid dynamics (CFD) computer codes is desirable for two reasons. First, it allows prediction of resulting flow fields within the rocket nozzle, in particular the interaction of the coolant boundary layer with the main flow. This facilitates evaluation of potential cooling configurations with regard to total thrust, etc., before construction and testing of any prototype. Secondly, CFD simulation of film cooling allows for assessment of the effectiveness of the proposed cooling in limiting nozzle wall temperature rises. This latter objective is the focus of the current work. The desired objective is to use the Finite Difference Navier Stokes (FDNS) code to predict wall heat fluxes or wall temperatures in rocket nozzles. As prior work has revealed that the FDNS code is deficient in the thermal modeling of boundary conditions, the first step is to correct these deficiencies in the FDNS code. Next, these changes must be tested against available data. Finally, the code will be used to model film cooling of a particular rocket nozzle. The third task of this research, using the modified code to compute the flow of hot gases through a nozzle, is described.

Woodbury, Keith A.↗
NASA NTRS1993

Determination of solid/porous wall boundary conditions from wind tunnel data for computational fluid dynamics codes

A computational and experimental study has been undertaken to investigate methods of modeling solid and porous wall boundary conditions in CFD codes. The procedure utilizes experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution is then imposed as a boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the boundary conditions at the wind tunnel walls using only sparse experimental measurements has been investigated. Verification of the approach using computational studies has been carried out using an incompressible flow solver. The current work demonstrates this technique for low speed flows and compares the result with experimental data obtained from a heavily instrumented variable porosity test section. Position and refinement of experimental measurements required to describe porous wall boundary conditions has also been considered for application to other porous wall wind tunnels.

Beutner, Thomas J.↗