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 451 records · Page 25

A simulation-based parametric study of CLOU chemical looping reactor performance

Chemical looping with oxygen uncoupling (CLOU) is a variant on chemical looping combustion in which the oxygen carrier releases gaseous O 2 in the fuel reactor, making it well-suited for solid fuels, since the released gaseous O 2 readily reacts with solid char. This work presents several computational fluid dynamic (CFD) simulations of copper-based CLOU in a dual fluidized bed system, each with different operating conditions. The modeling predicted that coal particle sizes as large as 1000 μm did not significantly affect performance. Increased oxygen carrier copper loading resulted in an excess of gaseous oxygen in the product gas stream. Decreasing the oxygen carrier bed mass as well as reducing the air reactor fluidizing velocity did not supply enough oxygen to the fuel reactor to complete combustion of the coal. This generated a failure state in which the temperature continued to decrease in the fuel reactor from the lack of combustion, which in turn reduced the O 2 equilibrium partial pressure, further lowering the amount of combustion possible. Sufficient O 2 can be maintained in the fuel reactor by ensuring a high enough air reactor velocity and a large enough supply of oxygen carrier inventory to handle the chosen coal feed rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗

Modeling and simulation of multiphase flows

This presentation provides an overview of the National Energy Technology Laboratory’s (NETL) multiphase computational fluid dynamics codes. The highly successful Multiphase Flows with Interphase eXchanges (MFIX) suite has been used to model a wide range of applications including post-combustion carbon capture, bioreactor optimization, and bio-FCC regeneration. MFIX-Exa, a state-of-the-art CFD code, developed under DOE’s Exascale Computing Project, is built on the AMReX software framework (https://amrex-codes.github.io/) and is designed to leverage modern accelerator-based compute architectures. This presentation further reviews the underlying physical models of both MFIX and MFIX-Exa and contrasts their similarities and differences. Examples of past and present CFD simulations will illustrate how scientific computing at NETL is being used not only for scientific exploration but also for design, optimization and scale-up of multiphase flow devices.

Musser, Jordan [NETL]↗

Performance and Emissions of an SI Engine Fueled With DME-Propane Blends

Dimethyl Ether (DME) is an alternative fuel that can be produced renewably and has the potential for lower CO and NOx emissions than conventional petroleum-based fuels. Blending DME with another gaseous fuel such as propane, which has a lower knock tendency than gasoline, can allow this fuel to be leveraged on SI engines. In this study, the use of DME-propane blends on a spark ignition (SI) engine was studied via computer simulations in order to understand the impact on engine performance and emissions and to identify the knock limitations of using such fuel blends. A 2L Hyundai SI engine was modeled in GT Power and the model was validated by comparing it with computational fluid dynamics (CFD) simulation results. Starting from pure propane, DME was added incrementally until knock was observed in the engine. Results indicate that it is feasible to run propane with DME percentages up to 35% before severe knock impacts were observed. The BTE was higher, but the BSFC also increased for DME-propane blends as compared to gasoline. Here, an increase in NOx emissions was detected along with a significant decrease in CO emissions. CO 2 emissions declined for propane as compared to gasoline but increased with the addition of DME.

alternative fuels↗

Gas Phase Effects on Slosh Dynamics

Gas phase effects on slosh dynamics were quantified using computational fluid dynamics (CFD) simulation for a range of propellant and gas combinations. Historical slosh modeling using potential flow solutions typically neglects gas phase effects. Regardless, the results have been shown to compare well with slosh ground tests typically performed with water and air at standard temperature and pressure. CFD analysis reveals that as the liquid-to-gas density ratio decreases, slosh dynamics change due to the relative increase in gas inertia and thus influence on liquid motion. The result is a profound impact on slosh dynamics over certain parameter spaces particularly for liquid hydrogen. Gas phase effects on slosh dynamics should be considered in slosh models especially for liquid hydrogen propellant tanks.

Jacob M Brodnick↗

Gas Phase Effects on Slosh Dynamics

Gas phase effects on slosh dynamics were quantified using computational fluid dynamics (CFD) simulation for a range of propellant and gas combinations. Historical slosh modeling using potential flow solutions typically neglects gas phase effects. Regardless, the results have been shown to compare well with slosh ground tests typically performed with water and air at standard temperature and pressure. CFD analysis reveals that as the liquid-to-gas density ratio decreases, slosh dynamics change due to the relative increase in gas inertia and thus influence on liquid motion. The result is a profound impact on slosh dynamics over certain parameter spaces particularly for liquid hydrogen. Gas phase effects on slosh dynamics should be considered in slosh models especially for liquid hydrogen propellant tanks.

Jacob M Brodnick↗

Comparing Methods to Extract Sectional Angle of Attack From Blade-Resolved Rotor Simulations

The angle of attack is widely used to describe the aerodynamics of rotor blades. Models based on lifting lines, actuator lines, and actuator disks use it as a primary variable, assuming the aerodynamics can be represented by the flow around a collection of equivalent 2D airfoils. This paper evaluates the accuracy of existing methods that extract the angle of attack from blade-resolved CFD simulations, towards enabling a direct comparison of angle of attack between blade-resolved and blade-modeled simulations. Results are compared for a 2D airfoil for which the angle of attack is known by definition. The hovering rotor case is also considered, where we look at two different rotor blades. We discuss the factors that impact the accuracy of the 2D predictions. For the 3D rotors, we show that three-dimensional effects in the root and tip regions lead to large discrepancies between methods. This suggests that the very definition of angle of attack is unable to accurately represent the flow in these regions. Finally, based on guidelines that we establish, we propose an alternate method that uses local velocities sampled on a contour located 0.25 chord, or less, away from the airfoil surface and generally yields improved accuracy compared to most existing techniques.

ARMD↗

Dynamic Stall Suppression Using Combustion-Powered Actuation (COMPACT)

Retreating blade stall is a well-known phenomenon that limits rotorcraft speed, maneuverability, and efficiency. Airfoil dynamic stall is a simpler problem, which demonstrates many of the same flow phenomena. Combustion Powered Actuation (COMPACT) is an active flow control technology, which at the outset of this work, had been shown to mitigate static and dynamic stall at low Mach numbers. The attributes of this technology suggested strong potential for success at higher Mach numbers, but such experiments had never been conducted. The work detailed in this report documents a 3-year effort focused on assessing the effectiveness of COMPACT for dynamic stall suppression at freestream conditions up to Mach 0.5. The work done has focused on implementing COMPACT on a high-lift rotorcraft airfoil: the VR-12. This selection was made in order to ensure that any measured benefits are over and above the capabilities of state-of-the-art high-lift rotorcraft airfoils. The detailed Computational Fluid Dynamics (CFD) simulations, wind-tunnel experiments, and system-level modeling conducted have shown the following: (1) COMPACT, in its current state of development, is capable of reducing the adverse effects of deep dynamic stall at Mach numbers up to 0.4; (2) The two-dimensional (2D) CFD results trend well compared to the experiments; and (3) Implementation of the CFD results into a system-level model suggest that significant rotor-level benefits are possible.

Matalanis, Claude G.↗

Gas Phase Effects on Slosh Dynamics

Gas phase effects on slosh dynamics were quantified using computational fluid dynamics (CFD) simulation for a range of propellant and ullage gas combinations. Historical slosh modeling using potential flow solutions typically neglects gas phase effects. Regardless, the results have been shown to compare well with slosh ground tests typically performed with water and air at standard temperature and pressure. CFD analysis reveals that as the liquid-to-gas density ratio decreases, slosh dynamics change due to the relative increase in gas inertia and thus influence on liquid motion. The result is a profound impact on slosh dynamics over certain parameter spaces particularly for liquid hydrogen. Gas phase effects on slosh dynamics should be considered in slosh models especially for liquid hydrogen propellant tanks.

Jacob M Brodnick↗

Spray process of multi-component gasoline surrogate fuel under ECN Spray G conditions

As modern gasoline direct injection (GDI) engines utilize sophisticated injection strategies, a detailed understanding of the air-fuel mixing process is crucial to further improvements in engine emission and fuel economy. In this study, a comprehensive evaluation of the spray process of single-component iso-octane (IC8) and multi-component gasoline surrogate E00 (36 % n-pentane, 46 % iso-octane, and 18 % n-undecane, by volume) fuels was conducted using an Engine Combustion Network (ECN) Spray G injector. High-speed extinction, schlieren, and microscopy imaging campaigns were carried out under engine-like ambient conditions in a spray vessel. Experimental results including liquid/vapor penetration, local liquid volume fraction, droplet size, and projected liquid film on the nozzle tip were compared under ECN G1 (573 K, 3.5 kg/m 3 ), G2 (333 K, 0.5 kg/m 3 ), and G3 (333 K, 1.01 kg/m 3 ) conditions. In addition to the experiments, preferential evaporation process of the E00 fuel was elucidated by Large–Eddy Simulations (LES). The three-dimensional liquid volume fraction measurement enabled by the computed tomographic reconstruction showed substantial plume collapse for E00 under the G2 and G3 conditions having wider plume growth and plume-to-plume interaction due to the fuel high vapor pressure. Here, the CFD simulation of E00 showed an inhomogeneity in the way fuel components vaporized, with more volatile components carried downstream in the spray after the end of injection. The high vapor pressure of E00 also results in ~4 μm smaller average droplet diameter than IC8, reflecting a higher rate of initial vaporization even though the final boiling point temperature is higher. Consistent with high vapor pressure, E00 had a wider plume cone angle and enhanced interaction with the wall to cover the entire surface of the nozzle tip in a film. However, the liquid fuel underwent faster evaporation, so the final projected tip wetting area was smaller than the IC8 under the flash-boiling condition.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical Investigation of Film Coefficient Engineering Methodology for Dissociated, Chemically Reacting Boundary Layers

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flow field for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for their application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. The sensitivity of the engineering model is assessed based on the impact of the CFD assumptions as an effort to make it available to design engineers across industry.

numerical boundary layer↗

Aerodynamic loss generation in a low-pressure turbine stage with high lift blading

The present work aims to investigate the loss generation mechanisms in a highly loaded low-pressure turbine (LPT) stage. The turbine employs high-lift airfoils to reduce the number of blades and therefore the overall engine weight. Large Eddy Simulation (LES) is applied to study the detailed flow fields at two representative Reynolds numbers (150,000 and 40,000). The calculated flow fields are compared with available measured data. LES calculates the measured aerodynamic performance of the current turbine stage well. Both the measurements and the CFD simulations indicate that most of the aerodynamic loss is generated by the development of secondary flows near the rotor hub. An area of low momentum and high aerodynamic loss forms at roughly 25 % span above the hub downstream of the rotor. Formation and migration of a vortex system near the hub is the primary source of the loss generation. The hub wall vortex system in this highly loaded turbine stage is quite different from the horse-shoe vortex system in a typical stationary turbine cascade. Calculated instantaneous and time-averaged flow fields at two Reynolds numbers are examined in detail to study Reynolds number effects on the turbine stage performance. Steady and unsteady flow separation on the rotor suction side is calculated at both high and low Reynolds numbers. At the lower Reynolds number of 40,000, thicker blade boundary layers increase the area of flow separation on the rotor suction side, resulting in a 1.9 % increase in aerodynamic loss compared to the higher Reynolds number case.

low pressure turbine↗

Anti-Phase Vortex Reduction Control for Rotor Noise Suppression

An investigation is conducted in 2019 under the NASA Ames Center Innovation Funds (CIF) project entitled “Anti-Phase Vortex Reduction Control for Rotor Noise Suppression“ to experimentally validate an anti-phase rotor noise suppression concept. The first objective of the investigation is to conduct computational fluid dynamics (CFD) simulations to investigate the noise characteristics of several anti-phase rotor designs. The second objective is to conduct a series of acoustic tests of anti-phase rotors in an anechoic chamber at Pennsylvania State University (PSU) to evaluate the merit of the anti-phase rotor concept. The CFD investigation seeks to optimize the anti-phase alternating trailing edge patterns for rotor noise suppression. The design objective is to maximize the noise reduction while maintaining the aerodynamic thrust. The investigation is performed using a three-dimensional (3D) Unsteady Reynolds-Averaged Navier-Stokes (URANS) commercial solver STAR-CCM+ together with the Ffowcs-Williams and Hawkings (FW-H) formula to obtain the aerodynamic thrust and far-field noise level. An acoustic study is conducted for 13 anti-phase design candidates based on a proprietary rotor design. The CAD geometry of the rotor is furnished by PSU. These design candidates include different alternating trailing edge (TE) waveforms, TE segment lengths, TE deflection amplitudes, and transition characteristics. The best design candidate among those explored is an anti-phase rotor that has a four-period TE waveform which results in a reduction in far-field noise level of 2.1 dB in the hover condition and a reduction of 1.1 dB in the forward flight condition at 9.7 m/s. A further acoustic study is conducted for a different rotor manufactured by APC. Five APC rotor design candidates are simulated. The best anti-phase design candidate for the APC rotor results in a reduction in far-field noise level of 4.0 dB in the hover condition. An in-phase design candidate is also studied. This in-phase design provides a noise reduction of 2.5 dB. A series of acoustic experiments in the PSU anechoic chamber have been conducted in July 2019 and October 2019. Both anti-phase and in-phase rotors fabricated for the left-hand and right-hand rotations are tested. In the hover condition, all the rotors do not provide sufficient evidence of improved acoustic performance. However, the experimental data in the hover condition are deemed to be inconclusive due to the flow recirculation in the anechoic chamber caused by the rotor downwash. In the forward flight condition at 9.7 m/s, the anti-phase right-hand rotor produces a noise reduction by as much as 6.5 dB in the frequency range of 2000-4000 Hz, while the in-phase right-hand rotor produces a noise reduction by as much as 5 dB in the same frequency range. Both the anti-phase and in-phase left hand rotors offer no evidence of noise reduction. The difference in the acoustic performance of the left-hand and right-hand rotors could be explained by the location of the microphone array which is placed to the left side of the rotors. This microphone array location could create a bias in the sound pressure level in favor of the right-hand rotor. Using the average noise reduction values, the anti-phase 4H rotor could offer a noise reduction by as much as 3.25 dB while the in-phase 4I rotor could produce up to a 2.5 dB noise reduction. Both the computational and experimental results have provided sufficient evidence to support the noise suppression capability of the proposed anti-phase rotor concepts.

Rotor noise↗

Advanced Turbulence Modeling Concepts

The ZCET program developed at NASA Glenn Research Center is to study hydrogen/air injection concepts for aircraft gas turbine engines that meet conventional gas turbine performance levels and provide low levels of harmful NOx emissions. A CFD study for ZCET program has been successfully carried out. It uses the most recently enhanced National combustion code (NCC) to perform CFD simulations for two configurations of hydrogen fuel injectors (GRC- and Sandia-injector). The results can be used to assist experimental studies to provide quick mixing, low emission and high performance fuel injector designs. The work started with the configuration of the single-hole injector. The computational models were taken from the experimental designs. For example, the GRC single-hole injector consists of one air tube (0.78 inches long and 0.265 inches in diameter) and two hydrogen tubes (0.3 inches long and 0.0226 inches in diameter opposed at 180 degree). The hydrogen tubes are located 0.3 inches upstream from the exit of the air element (the inlet location for the combustor). To do the simulation, the single-hole injector is connected to a combustor model (8.16 inches long and 0.5 inches in diameter). The inlet conditions for air and hydrogen elements are defined according to actual experimental designs. Two crossing jets of hydrogen/air are simulated in detail in the injector. The cold flow, reacting flow, flame temperature, combustor pressure and possible flashback phenomena are studied. Two grid resolutions of the numerical model have been adopted. The first computational grid contains 0.52 million elements, the second one contains over 1.3 million elements. The CFD results have shown only about 5% difference between the two grid resolutions. Therefore, the CFD result obtained from the model of 1.3-million grid resolution can be considered as a grid independent numerical solution. Turbulence models built in NCC are consolidated and well tested. They can handle both coarse and fine grids near the wall. They can model the effect of anisotropy of turbulent stresses and the effect of swirling. The chemical reactions of Magnusson model and ILDM method were both used in this study.

Shih, Tsan-Hsing↗

Computational and Experimental Validation of Fractal-Fin, Dimpled Solar Heat Collector

A 3 foot x 3 foot x 3 foot aluminum solar collector was manufactured using computer numerical control. The interior of the device included six triangular dimpled fins for enhanced heat transfer. The interior vertical wall on the south side was also dimpled. The solar collector working fluid was based on water, and the collector consisted solely of passive heat transfer mechanisms (no moving parts), making it ideal for off-the-grid and rural applications. Two types of heat transfer experiments were conducted. One experiment had external flat heaters attached on the top and the front side, while the other four sides were insulated. Except for the bottom surface, the second experiment had all its exterior surfaces sprayed with black solar paint to collect as much solar heat as possible. Temperature data as a function of time was collected using 14 thermocouples spread strategically throughout the solar collector. In addition, computational fluid dynamics (CFD) simulations were conducted using the dynamic Smagorinsky large eddy simulation turbulence model. The first simulation considered that both the top and front surfaces were exposed to a fixed temperature of 313.7 K (105 °F), while the remaining four surfaces were insulated. For the second simulation, all conditions were the same, except that the temperature for both heated surfaces was raised to 350 K (170.3 °F). The two temperatures are expected to bound the solar collector operational temperature during the late- Spring, Summer, and early-Fall months. The solar collector design, experimental data, CFD output, and a discussion of five manufacturing approaches and costs are documented in this report.

14 SOLAR ENERGY↗

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

Europa Clipper Preliminary Design Review Propellant Slosh Analysis

This presentation outlines the preliminary design review (PDR) propellant slosh analysis done for the Europa Clipper Mission. It provides sample results for both high acceleration pendulum-damper models and low acceleration pendulum-spring-damper models. The high acceleration pendulum-damper models were derived from STAR-CCM+ computational fluid dynamic (CFD) simulations and the low acceleration pendulum-spring-damper models were derived from Surface Evolver models.

slosh↗

Measurements and Computational Analysis of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow Experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This report presents and discusses the results of a recent experiment with the Juncture Flow Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of the test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model angle of incidence. Laser Doppler velocimetry (LDV) measurements were made at three model angles of incidence (0° : fully attached, 1° : incipient separation, and 5° : separated flow) and for each angle, mean-flow and Reynolds-stress data were obtained on the fuselage and at several locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and include model and tunnel wall static pressures, tunnel wall boundary-layer profiles, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. A companion set of Reynolds-averaged Navier-Stokes (RANS) CFD simulations were also performed and those results were used to support the discussion of the experimental test results.

Juncture Flow↗