SEARCH · Search NASA

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗

DOE OSTI2022

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗

DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗

NASA NTRS2011

Analysis of ALTAIR 1998 Meteor Radar Data

We describe a new analysis of a set of 32 UHF meteor radar traces recorded with the 422 MHz ALTAIR radar facility in November 1998. Emphasis is on the velocity measurements, and on inferences that can be drawn from them regarding the meteor masses and mass densities. We find that the velocity vs altitude data can be fitted as quadratic functions of the path integrals of the atmospheric densities vs distance, and deceleration rates derived from those fits all show the expected behavior of increasing with decreasing altitude. We also describe a computer model of the coupled processes of collisional heating, radiative cooling, evaporative cooling and ablation, and deceleration - for meteors composed of defined mixtures of mineral constituents. For each of the cases in the data set we ran the model starting with the measured initial velocity and trajectory inclination, and with various trial values of the quantity mPs 2 (the initial mass times the mass density squared), and then compared the computed deceleration vs altitude curves vs the measured ones. In this way we arrived at the best-fit values of the mPs 2 for each of the measured meteor traces. Then further, assuming various trial values of the density Ps, we compared the computed mass vs altitude curves with similar curves for the same set of meteors determined previously from the measured radar cross sections and an electrostatic scattering model. In this way we arrived at estimates of the best-fit mass densities Ps for each of the cases. Keywords meteor ALTAIR radar analysis 1 Introduction This paper describes a new analysis of a set of 422 MHz meteor scatter radar data recorded with the ALTAIR High-Power-Large-Aperture radar facility at Kwajalein Atoll on 18 November 1998. The exceptional accuracy/precision of the ALTAIR tracking data allow us to determine quite accurate meteor trajectories, velocities and deceleration rates. The measurements and velocity/deceleration data analysis are described in Sections II and III. The main point of this paper is to use these deceleration rate data, together with results from a computer model, to determine values of the quantities mPs 2 (the meteor mass times its material density squared); and further, by combining these m s 2 values with meteor mass estimates for the same set of meteors determined separately from measured radar scattering

Zinn, J.↗

NASA NTRS2023

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.

Transition Model↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

NASA NTRS2022

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

CFD modeling↗

NASA NTRS2014

Turbulence Modeling Verification and Validation

Computational fluid dynamics (CFD) software that solves the Reynolds-averaged Navier-Stokes (RANS) equations has been in routine use for more than a quarter of a century. It is currently employed not only for basic research in fluid dynamics, but also for the analysis and design processes in many industries worldwide, including aerospace, automotive, power generation, chemical manufacturing, polymer processing, and petroleum exploration. A key feature of RANS CFD is the turbulence model. Because the RANS equations are unclosed, a model is necessary to describe the effects of the turbulence on the mean flow, through the Reynolds stress terms. The turbulence model is one of the largest sources of uncertainty in RANS CFD, and most models are known to be flawed in one way or another. Alternative methods such as direct numerical simulations (DNS) and large eddy simulations (LES) rely less on modeling and hence include more physics than RANS. In DNS all turbulent scales are resolved, and in LES the large scales are resolved and the effects of the smallest turbulence scales are modeled. However, both DNS and LES are too expensive for most routine industrial usage on today's computers. Hybrid RANS-LES, which blends RANS near walls with LES away from walls, helps to moderate the cost while still retaining some of the scale-resolving capability of LES, but for some applications it can still be too expensive. Even considering its associated uncertainties, RANS turbulence modeling has proved to be very useful for a wide variety of applications. For example, in the aerospace field, many RANS models are considered to be reliable for computing attached flows. However, existing turbulence models are known to be inaccurate for many flows involving separation. Research has been ongoing for decades in an attempt to improve turbulence models for separated and other nonequilibrium flows. When developing or improving turbulence models, both verification and validation are important steps in the process. Verification insures that the CFD code is solving the equations as intended (no errors in the implementation). This is typically done either through the method of manufactured solutions (MMS) or through careful step-by-step comparisons with other verified codes. After the verification step is concluded, validation is performed to document the ability of the turbulence model to represent different types of flow physics. Validation can involve a large number of test case comparisons with experiments, theory, or DNS. Organized workshops have proved to be valuable resources for the turbulence modeling community in its pursuit of turbulence modeling verification and validation. Workshop contributors using different CFD codes run the same cases, often according to strict guidelines, and compare results. Through these comparisons, it is often possible to (1) identify codes that have likely implementation errors, and (2) gain insight into the capabilities and shortcomings of different turbulence models to predict the flow physics associated with particular types of flows. These are valuable lessons because they help bring consistency to CFD codes by encouraging the correction of faulty programming and facilitating the adoption of better models. They also sometimes point to specific areas needed for improvement in the models. In this paper, several recent workshops are summarized primarily from the point of view of turbulence modeling verification and validation. Furthermore, the NASA Langley Turbulence Modeling Resource website is described. The purpose of this site is to provide a central location where RANS turbulence models are documented, and test cases, grids, and data are provided. The goal of this paper is to provide an abbreviated survey of turbulence modeling verification and validation efforts, summarize some of the outcomes, and give some ideas for future endeavors in this area.

Rumsey, Christopher L.↗

NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 1: Validation of a Computational Model for the Vaneless Diffuser Configuration with As-Manufactured Impeller Geometry

An investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. This multipart investigation focused on validating the computational model against experimental data. The validated model was then used to explore the effects of observed geometric differences between the as-manufactured and design-intent impellers. Lastly, the validated model was used for comparison against an experimental tip clearance study. Details regarding each of the parts are detailed below: A multipart investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. Part I of the investigation focused on describing the HECC vaneless computational model and its validation against the experimental data published in GT2023-103128. To generate the model, a verification of the manufactured impeller was performed and was confirmed to match the impeller CAD. Section data from the impeller CAD, referred to as As-Manufactured, were used in the model. Steady RANS simulations using a multiblock structured mesh were used to conduct a mesh convergence study. The resultant mesh was then used for simulations at both design-speed and off-design speed conditions. The simulations compared well to data in both 1D and 2D terms. Based on the observed results the model was deemed validated and was used in the latter two parts of the multipart investigation.

centrifugal compressor↗

DOE OSTI2023

Stakeholder-Driven Improvements to Nek5000 and NekRS

As a key part of the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program (NEAMS), the next generation of tools for nuclear reactor design are being developed. A primary metric for success of this program is the adoption of those tools by industry stakeholders. Based on feedback provided directly to NEAMS, two of the major weaknesses have been a lack of moderate-fidelity Reynolds-averaged Navier-Stokes (RANS) turbulence models and available documentation. Over the past several years, developments have been focused on addressing the needs of those stakeholders. The outcome of those efforts are described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗

NASA NTRS2009

Experimental Studies of the Aerothermal Characteristics of the Project Orion CEV heat Shield in High Speed Transitional and Turbulent Flows

An experimental program has been completed by CUBRC exploring laminar, transitional, and turbulent flows over a 7.0% scale model of the Project ORION CEV geometry. This program was executed primarily to answer questions concerning the increase in heat transfer on the windward, or "hot shoulder" of the CEV heat shield from laminar to turbulent flow. To answer these questions CUBRC constructed and instrumented a 14.0 inch diameter Project ORION CEV model and ran a range of Reynolds numbers based on diameter from 1.0 to over 40 million at a Mach number of 8.0. These Reynolds numbers were selected to cover laminar to turbulent heating data on the "hot shoulder". Data obtained during these runs will be used to guide design decisions as they apply to heat shield thickness and extent. Several experiments at higher enthalpies were achieved to obtain data for code validation with real gas effects and transition. CUBRC also performed computation studies of these experiments to aid in the data reduction process and study turbulence modeling.

Wadhams, T.P.↗

NASA NTRS2015

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗

DOE OSTI2025

Surface heat flux and its association with the MJO in the tropical western Pacific using ARM observations

The Madden-Julian oscillation (MJO) is a major atmospheric phenomenon in the tropics that moves eastward every 20 to 100 days. It brings heavy rain and strong winds, influencing extreme weather events far beyond the tropics – including flooding, hurricanes, tornadoes and heavy snow in the United States. Therefore, an improved understanding of the MJO is critical for enhancing weather forecasts and supporting better decision-making for communities, emergency managers, and the private sector. Past studies, based on short-term observations or long-term model data, have emphasized the dominant role of atmospheric humidity in driving the MJO. Our research, using long-term observations (2000-2014) from three U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) sites in the tropical western Pacific (Manus, Nauru, and Darwin) confirms those past findings. However, we also show that atmospheric temperature, especially in the mid and upper troposphere, also plays a key role in the MJO’s evolution from its quiet (suppressed) to convective (active) phase. These findings provide valuable insights into how the MJO evolves and may help evaluate and improve weather and climate models. Our study also examined how the exchange of heat between the Earth's surface and the atmosphere, called surface heat flux, interacts with the MJO. This flux includes components such as latent heat (related to evaporation), sensible heat (related to temperature difference between surface and atmosphere), and radiation. We found that modulation in MJO convection is well connected to surface heat flux over the tropical western Pacific. One surface heat flux component often overlooked is the sensible heat flux due to precipitation (Q P ). Since falling raindrops are typically cooler than the land or ocean surface, they can cool the surface and affect local weather. To test the impact of Q P on convection during MJO, we incorporated it into a weather model, and ran simulations for two periods: April 2009 (when the MJO was active) and June 2006 (when it was not). Including Q P reduced errors in simulating the daily peak of rainfall—by 83% during the strong MJO and 23% during the inactive phase. It also improved the timing of peak daily rainfall and reduced the overall precipitation error by about 10%. These improvements were especially noticeable during periods of heavy rainfall. Our results suggest that even small heat exchanges from rainfall can play a significant role in shaping local weather. Accounting for these effects, particularly in regions like the tropical islands, can lead to more accurate simulations of the MJO-associated precipitation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

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↗

DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS2001

The Numerical Analysis of a Turbulent Compressible Jet

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

DeBonis, James R.↗

NASA NTRS2014

Assessment of Higher-Order RANS Closures in a Decelerated Planar Wall-Bounded Turbulent Flow

A reference DNS database is presented, which includes third- and fourth-order moment budgets for unstrained and strained planar channel flow. Existing RANS closure models for third- and fourth-order terms are surveyed, and new model ideas are introduced. The various models are then compared with the DNS data term by term using a priori testing of the higher-order budgets of turbulence transport, velocity-pressure-gradient, and dissipation for both the unstrained and strained databases. Generally, the models for the velocity-pressure-gradient terms are most in need of improvement.

Jeyapaul, Elbert↗