Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence simulation”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Transition to turbulence in plane channel flows

Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗
NASA NTRS1985

Parameterization of turbulence and the planetary boundary layer in the GLA Fourth Order GCM

A new scheme has been developed to model the planetary boundary layer in the GLAS Fourth Order GCM through explicit resolution of its vertical structure into two or more vertical layers. This involves packing the lowest layers of the GCM close to the ground and developing new parameterization schemes that can express the turbulent vertical fluxes of heat, momentum and moisture at the earth's surface and between the layers that are contained with the PBL region. Offline experiments indicate that the combination of the modified level 2.5 second-order turbulent closure scheme and the 'extended surface layer' similarity scheme should work well to simulate the behavior of the turbulent PBL even at the coarsest vertical resolution with which such schemes will conceivably be used in the GLA Fourth Order GCM.

Helfand, H. M.↗
NASA NTRS2011

Global Variation of Meteor Trail Plasma Turbulence

We present the first global simulations on the occurrence of meteor trail plasma irregularities. These results seek to answer the following questions: when a meteoroid disintegrates in the atmosphere will the resulting trail become plasma turbulent, what are the factors influencing the development of turbulence, and how do they vary on a global scale. Understanding meteor trail plasma turbulence is important because turbulent meteor trails are visible as non-specular trails to coherent radars, and turbulence influences the evolution of specular radar meteor trails, particularly regarding the inference of mesospheric temperatures from trail diffusion rates, and their usage for meteor burst communication. We provide evidence of the significant effect that neutral atmospheric winds and density, and ionospheric plasma density have on the variability of meteor trail evolution and the observation of nonspecular meteor trails, and demonstrate that trails are far less likely to become and remain turbulent in daylight, explaining several observational trends using non-specular and specular meteor trails.

Dyrud, L. P.↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
NASA NTRS2012

Creating a Realistic Weather Environment for Motion-Based Piloted Flight Simulation

A flight simulation environment is being enhanced to facilitate experiments that evaluate research prototypes of advanced onboard weather radar, hazard/integrity monitoring (HIM), and integrated alerting and notification (IAN) concepts in adverse weather conditions. The simulation environment uses weather data based on real weather events to support operational scenarios in a terminal area. A simulated atmospheric environment was realized by using numerical weather data sets. These were produced from the High-Resolution Rapid Refresh (HRRR) model hosted and run by the National Oceanic and Atmospheric Administration (NOAA). To align with the planned flight simulation experiment requirements, several HRRR data sets were acquired courtesy of NOAA. These data sets coincided with severe weather events at the Memphis International Airport (MEM) in Memphis, TN. In addition, representative flight tracks for approaches and departures at MEM were generated and used to develop and test simulations of (1) what onboard sensors such as the weather radar would observe; (2) what datalinks of weather information would provide; and (3) what atmospheric conditions the aircraft would experience (e.g. turbulence, winds, and icing). The simulation includes a weather radar display that provides weather and turbulence modes, derived from the modeled weather along the flight track. The radar capabilities and the pilots controls simulate current-generation commercial weather radar systems. Appropriate data-linked weather advisories (e.g., SIGMET) were derived from the HRRR weather models and provided to the pilot consistent with NextGen concepts of use for Aeronautical Information Service (AIS) and Meteorological (MET) data link products. The net result of this simulation development was the creation of an environment that supports investigations of new flight deck information systems, methods for incorporation of better weather information, and pilot interface and operational improvements for better aviation safety. This research is part of a larger effort at NASA to study the impact of the growing complexity of operations, information, and systems on crew decision-making and response effectiveness; and then to recommend methods for improving future designs.

Daniels, Taumi S.↗
NASA NTRS2005

Coherent Optical Receiver for PPM Signals under Atmospheric Turbulence

Adaptive combining of experimentally obtained heterodyned pulse position modulated (PPM) signals with pulse-to-pulse coherence in the presence of simulated spatial distortions resembling atmospheric turbulence is demonstrated. The adaptively combined PPM signals are phased up via an LMS algorithm suitably optimized to operate with PPM in the presence of additive shot-noise. A convergence analysis of the algorithm is presented, and results with both, computer simulated and experimentally obtained PPM signals are analyzed.

coherent detection↗
NASA NTRS2004

Adaptive Combining of Coherently Detected PPM Signals in the Presence of Atmospheric Turbulence

Adaptive combining of experimentally obtained heterodyned pulse position modulated (PPM) signals with pulse-to-pulse coherence, in the presence of simulated spatial distortions resembling atmospheric turbulence, is demonstrated. The adaptively combined PPM signals are phased up via an LMS algorithm suitably optimized to operate with PPM in the presence of additive shot-noise. A convergence study of the algorithm is presented, and results with both computer simulated and experimentally obtained PPM signals are presented and analyzed.

pulse position modulation (PPM)↗
NASA NTRSDate not supplied

Large-Eddy Simulation of the Variable Speed Power Turbine Cascade with Inflow Turbulence

The main goal of this is to implement a digital filtering method to generate turbulence upstream and a sub-grid model (Localized dynamic k-equation model (LDKM)) in the framework of LES which together enable reproduction of the effect of inflow turbulence on the transition seen in the VSPT experimental data at the cruise condition (incidence angle of 40o and Tu=0.5%, 5%, 10% and 15%.).

CFD↗
DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗
NASA NTRS1979

Full-scale aircraft simulation with cryogenic tunnels and status of the National Transonic Facility

The effect of thermal and caloric imperfections in cryogenic nitrogen on boundary layers was determined to indicate that in order to simulate nonadiabatic laminar or turbulent boundary layers in a cryogenic nitrogen wind tunnel, the flight enthalpy ratio, rather than the temperature ratio, should be reproduced. The absence of significant real gas effects on both viscous and inviscid flows makes it unlikely that there will be large real gas effects on the cryogenic tunnel simulation of shock boundary layer interactions or other complex flow conditions encountered in flight. Condensation effects were studied to determine the minimum usable temperature and indicated that under most circumstances free stream Mach number rather than maximum local Mach number determines the onset of condensation effects.

Kilgore, R. A.↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Studying marine stratus with large eddy simulation

Data sets from field experiments over the stratocumulus regime may include complications from larger scale variations, decoupled cloud layers, diurnal cycle, or entrainment instability, etc. On top of the already complicated turbulence-radiation-condensation processes within the cloud-topped boundary layer (CTBL), these complexities may sometimes make interpretation of the data sets difficult. To study these processes, a better understanding is needed of the basic processes involved in the prototype CTBL. For example, is cloud top radiative cooling the primary source of the turbulent kinetic energy (TKE) within the CTBL. Historically, laboratory measurements have played an important role in addressing the turbulence problems. The CTBL is a turbulent field which is probably impossible to generate in laboratories. Large eddy simulation (LES) is an alternative way of 'measuring' the turbulent structure under controlled environments, which allows the systematic examination of the basic physical processes involved. However, there are problems with the LES approach for the CTBL. The LES data need to be consistent with the observed data. The LES approach is discussed, and results are given which provide some insights into the simulated turbulent flow field. Problems with this approach for the CTBL and information from the FIRE experiment needed to justify the LES results are discussed.

Moeng, Chin-Hoh↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
NASA NTRS2005

On the Origin of the Saturnian Satellite System: Did Iapetus Form In-Situ?

Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.

Mosqueira, I.↗
NASA NTRS2015

TPSAS-NF1676L-22800-DND

This lecture will provide an overview of an advanced turbulence closure scheme used in cloud-resolving models (CRMs) and climate modeling effort with the superparameterization approach. The CRM is embedded into each atmospheric column of a climate model as a super parameterization. It replaces the traditional subgrid-scale physical parameterizations in a conventional climate model. Two sensitivity simulations with doubling of CO2 and an increase of sea surface temperature, respectively, are performed and compared to similar simulations without this advanced turbulence closure in CRM to better understand the role of low-level clouds in climate sensitivity.

Kuan-Man Xu↗
NASA NTRS1990

A subgrid-scale model based on the second-order velocity structure function

A series of tests were performed to help extend the use of a subgrid-scale model to compressible and wall-bounded flows. A priori tests were done in the case of the incompressible turbulent channel flow. They showed that a 1-D formulation of the structure-function model is more appropriate, leading to a satisfactory behavior of the model at the walls without requiring any damping function. This model is consistent with the original formulation of Metais & Lesieur (1990). In large-eddy simulations of compressible isotropic turbulence, both models performed well up to an initial rms Mach number of 0.6.

Comte, P.↗
NASA NTRS1982

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗
DOE OSTI2022

Structure of the thermal boundary layer in turbulent channel flows at transcritical conditions

Enhanced fluctuations, steep gradients, and intensified heat transfer are characteristics of wall-bounded turbulence at transcritical conditions. Although such conditions are prevalent in numerous technical applications, the structure of the thermal boundary layer under realistic density gradients and heating conditions remains poorly understood. Specifically, statistical descriptions of the temperature field in such flows are provided inconsistently using existing models. To address this issue, direct numerical simulations are performed by considering fully developed transcritical turbulent channel flow at pressure and temperature conditions that cause density changes of a factor of up to $O(20)$ between the hot and cold walls. As a consequence of the proximity of the Widom line to the hot wall, significant asymmetries are observed when comparing regions near the cold wall and near the hot wall. Previous transformations that attempt to collapse the near-wall mean temperature profiles among different cases to a single curve are examined. By addressing model deficiencies of these transformations, a formulation for an improved mean temperature transformation is proposed, with appropriate considerations for real fluid effects that involve strong variations in thermodynamic quantities. Our proposed transformation is shown to perform well in collapsing the slope of the logarithmic region to a single universal value with reduced uncertainty. Coupled with a predictive framework to estimate the non-universal shift parameter of the logarithmic region using a priori information, our transformation provides an analytic profile to model the near-wall mean temperature. These results thus provide a framework to guide the development of models for wall-bounded transcritical turbulence.

42 ENGINEERING↗