LES of a Mach 2 Expansion Compression Flow
Presentation covers recent LES simulations of which have similar flow conditions to recent TriSonic Wind Tunnel experiments.
SEARCH · Engineering Papers
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.
Presentation covers recent LES simulations of which have similar flow conditions to recent TriSonic Wind Tunnel experiments.
Brief Description Shown is the supporting information for Horne et al. (2026). These files include all outputs from the WRF model simulations and the observational datasets used for comparison in the study. Scripts are provided so users can recreate the manuscript's figures using the provided observational and modeling data. For more information regarding the study, please contact the primary author of the associated manuscript, Jason Horne. Horne, J. P., Pan, Y., Davis, K. J., Waugh, D., Ahlswede, B. J., Prince, N. E. (2026). Simulating near-surface environments in urban neighborhoods using WRF-LES: A case study of classic atmospheric boundary layer (ABL) during a heatwave event JAMES. (to be submitted)
Brief Description Shown is the supporting information for Horne et al. (2026). These files include all outputs from the WRF model simulations and the observational datasets used for comparison in the study. Scripts are provided so users can recreate the manuscript's figures using the provided observational and modeling data. For more information regarding the study, please contact the primary author of the associated manuscript, Jason Horne. Horne, J. P., Pan, Y., Davis, K. J., Waugh, D., Ahlswede, B. J., Prince, N. E. (2026). Simulating near-surface environments in urban neighborhoods using WRF-LES: A case study of classic atmospheric boundary layer (ABL) during a heatwave event JAMES.
Several thermoelectric generators ranging in output power from 0.5 to 155 W have been completed or are undergoing testing at JPL. These generators represent a wide range of technologies, using Bi2Te3, PbTe and SiGe thermoelectric materials. Several of these generators are of a developmental type, such as HPG S/N2, and others are representative of Transit and Multi-Hundred Watt (MHW) Technology. Representative flight performance data of LES 8/9 and Voyager RTG's are presented and compared with the DEGRA computer program based on the data observed from tests of SiGe couples, modules and MHW generators.
The principal objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. A summary of work accomplished during the last six months is presented.
The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Sep. 1993 - 1 Sep. 1994, we have focused our efforts on two research problems: (1) developments of 'algebraic' moment closures for statistical descriptions of nonpremixed reacting systems, and (2) assessments of the Dirichlet frequency in presumed scalar probability density function (PDF) methods in stochastic description of turbulent reacting flows. This report provides a complete description of our efforts during this past year as supported by the NASA Langley Research Center under Grant NAG1-1122.
Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.
The extensive experimental study of dispersed two-phase turbulent flow in a vertical channel has been performed in Eaton's research group in the Mechanical Engineering Department at Stanford University. In Wang & Squires (1996), this study motivated the validation of LES approach with Lagrangian tracking of round particles governed by drag forces. While the computed velocity of the flow have been predicted relatively well, the computed particle velocity differed strongly from the measured one. Using Monte Carlo simulation of inter-particle collisions, the computation of Yamamoto et al. (2001) was specifically performed to model Eaton's experiment. The results of Yamamoto et al. (2001) improved the particle velocity distribution. At the same time, Vance & Squires (2002) mentioned that the stochastic simualtion of inter-particle collisions is too expensive, requiring significantly more CPU resources than one needs for the gas flow computation. Therefore, the need comes to account for the inter-particle collisions in a simpler and still effective way. To present such a model in the framework of LES/Lagrangian particle approach, and to compare the calculated results with Eaton's measurement and modeling of Yamamoto is the main objective of the present paper.
The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.
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.
This paper summarizes the reacting results of simulating a bluff body stabilized flame experiment of Volvo Validation Rig using a releasable edition of the National Combustion Code (NCC). The turbulence models selected to investigate the configuration are the sub-grid scaled kinetic energy coupled large eddy simulation (K-LES) and the time-filtered Navier-Stokes (TFNS) simulation. The turbulence chemistry interaction used is linear eddy mixing (LEM).
A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and Lattice-Boltzmann methods on a complex high-lift configuration across a wide range of angles of attack. Preliminary analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles of attack due to large pressure gradients present on the outboard slat and main element for𝛼 >17◦(corrected for free-air).While some scatter is reported in pitching moments at the low-angles of attack (𝛼 <11◦) - likely caused by differences in flap separation possibly due to low Reynolds number effects - excellent agreement is observed between the submissions near the𝐶𝐿,max state. Objective superiority over steady state RANS can be seen in terms of lack of excess outboard separation; a majority of the good quality WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. Differences in the onset of stall mechanism in the free-air configuration for𝛼 >20◦is reported with two distinct topologies observed. Topology A is characterized by the onset of corner-flow separation which progressively grows to produce a pitch break in free-air with an angle of attack increase from𝛼=20.55◦to𝛼=21.47◦+𝜖where𝜖is a small perturbation (varying between submissions and likely to be within±0.3◦). Topology Bis characterized by boundary layer weakness emanating from the inboard side of the wing-pylon juncture substantially larger than any weakness in the wing-body juncture (due to corner-flow);submissions in this category do not show any tendency for occurrence of a pitch break for the free-air configuration within the vicinity of𝛼=21.47◦. The in-tunnel simulations submitted by3 participants using different discretizations, grids and closure models show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of all 3 WMLES submissions over a single steady state RANS submission to the workshop.
A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and LB methods on a complex high-lift configuration across a wide range of angles-of-attack. Analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles-of-attack due to large pressure gradients present on the outboard slat and main element for 𝛼 >17◦(corrected for free-air). While some scatter is reported in pitching moment coefficient at the low-angles of attack (𝛼 <11◦) - excellent agreement is observed between submissions near the 𝐶𝐿,max state. Objective superiority of WMLES methods over RANS can be seen in terms of lack of excess outboard separation; a majority of the WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. The in-tunnel simulations show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of the WMLES submissions over RANS.
An assessment of a Hybrid RANS/LES (HRLES) approach for𝐶𝐿,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both the free air and the wind tunnel configuration of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and to two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES was shown to address some of the known shortcomings in RANS methods and prevent inboard and outboard flow separation particularly in the region of𝐶𝐿,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. HRLES predicts a weak pitch break at the highest angle-of-attack due to onset of wing-root corner flow separation whereas the free-air corrected experiment values indicate an occurrence of a much stronger pitch break. The improvements of HRLES over a URANS approach has been objectively shown by computing a set of solutions with the same grid, same numerics and time-step size and comparing the solutions. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly under performed RANS. For the in tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.
In this study, the Open National Combustion code (OpenNCC) is applied to simulate the airflow inside the high-pressure turbine (HPT) of the energy efficient engine (EEE). The main objective of this study is to validate the proposed methodology of a three-dimensional unsteady calculation (LES) of the HPT with a moving mesh capability to capture a relative motion at the stator and rotor interface. In our previous study, the adiabatic wall assumption was made at all the solid surfaces of HPT (i.e., the cooling airflows were not included.) To relax this assumption, we implement the surface source team approach, in which we impose the source term at a specific area of a cooling airflow exit at the solid surfaces by specifying an injection angle, temperature, turbulent intensity, and mass flowrate of each cooling airflow. The beneficial feature of this approach is that it does not require a mesh for each hole. For the validation of the model, the data for EEE model is obtained from the General Electric (GE) test campaign, and the full scale warm-air rig test condition representing a test point close to the integrated core/low spool design condition is considered.
Wall-modeled large-eddy simulation (LES) of a swept wing with leading-edge ice build up is performed using three mesh paradigms and associated flow solvers and compared to experimental results. The study focuses on an 8.9% scale model of the CRM65 swept wing featuring both high-fidelity and smooth leading-edge ice shapes. The assessment is conducted using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework using the three actively supported meshing paradigms in LAVA: structured curvilinear overset, unstructured, andCartesian. For the iced configurations, these paradigms use hybrid body-fitted/source-term immersed-boundary, body-fitted, and ghost-cell immersed-boundary strategies, respectively.The unstructured and Cartesian mesh paradigms are particularly attractive for complex ice shapes as they avoid the manual mesh generation effort associated with the curvilinear approach. For the high-fidelity ice shape, good agreement with the experiment can be obtained with all three strategies; however, the curvilinear method is particularly sensitive to the source-term immersed-boundary timescale and span wise mesh resolution, and the Cartesian implementation is sensitive to the choice of numerical flux. For the smooth ice, larger discrepancies are observed across all methods. However, a mesh refinement study guided by flow visualizations leads to an improved comparison with the experiment that is particularly pronounced for the unstructured mesh paradigm.
As computing technology continues to improve, simulating unsteady aerodynamic environments on launch vehicle geometries becomes more viable. For NASA’s Space Launch System, the flowfield is characterized by large regions of unsteady turbulent wall-bounded flow. Simulating these unsteady wall-bounded environments is of interest to the SLS Aerodynamics team. Wall-modeled LES methods are currently the best approach for simulating these environments in engineering applications. Less-expensive computationally than fully-resolved methods while maintaining many of the benefits in simulating unsteady environments. Before applying these methods blindly to large geometries, best practices can first be applied and understood on smaller canonical problems.
This report presents a short summary of the activities of the "Standard Model" working group for the "Physics at TeV Colliders" workshop (Les Houches, France, 12-30 June, 2023).