Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak turbulence”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Aeroacoustics of Turbulent High-Speed Jets

Aeroacoustic noise generation in a supersonic round jet is studied to understand in particular the effect of turbulence structure on the noise without numerically compromising the turbulence itself. This means that direct numerical simulations (DNS's) are needed. In order to use DNS at high enough Reynolds numbers to get sufficient turbulence structure we have decided to solve the temporal jet problem, using periodicity in the direction of the jet axis. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. Therefore in order to answer some questions about the turbulence we will partially compromise the overall structure of the jet. The first section of chapter 1 describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. In the second section we present preliminary work done using a TVD numerical scheme on a CM5. This work is only two-dimensional (plane) but shows very interesting results, including weak shock waves. However this is a nonviscous computation and the method resolves the shocks by adding extra numerical dissipation where the gradients are large. One wonders whether the extra dissipation would influence small turbulent structures like small intense vortices. The second chapter is an extensive discussion of preliminary numerical work using the spectral method to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which are solved in O(N) steps. A very accurate highly resolved DNS of a turbulent jet flow is expected.

Rao, Ram Mohan↗
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2026

On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf

Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.

FOS: Physical sciences↗
NASA NTRS1988

On the hazard of hydrogen explosions at space shuttle launch pads

This report was prepared in support of efforts to assess the hazard of accidental explosions of unburned hydrogen at space shuttle launch pads. It begins with presentation of fundamental detonation theory and a review of relevent experiments. A scenario for a catastrophic explosion at a KSC launch pad and a list of necessary conditions contributing to it is proposed with a view to identifying those conditions which, if blocked, would prevent a catastrophe. The balance of the report is devoted to juxtaposition of reassuring and disquieting facts, presentation of a set of recommendations that ignition of hydrogen-air mixtures by weak ignition sources in unconfined geometries may produce a detonation, provided the effective flame area in the initial fireball is rapidly increased by turbulent mixing. Another conclusion is that detonability limits can be different from and narrower than flammability limits only if one restricts the rate of work that can be done on a flammable gas by mechanical agencies acting on its boundaries.

Russell, John M.↗
NASA NTRS1990

Improvements to a nonequilibrium algebraic turbulence model

It has been noted that while the nonequilibrium turbulence model of Johnson and King (1985, 1987) performed significantly better than alternative methods, differences between predicted and observed shock locations for certain weak interactions are produced due to a defficiency in the model's inner eddy viscosity formulation. A novel formulation for the model is presented which removes this deficiency, while satisfying the law of the wall for adverse pressure-gradient conditions better than either the original formulation or mixing-length theory.

Johnson, D. A.↗
NASA NTRS2017

A Comparison of Three Random Number Generators for Aircraft Dynamic Modeling Applications

Three random number generators, which produce Gaussian white noise sequences, were compared to assess their suitability in aircraft dynamic modeling applications. The first generator considered was the MATLAB (registered) implementation of the Mersenne-Twister algorithm. The second generator was a website called Random.org, which processes atmospheric noise measured using radios to create the random numbers. The third generator was based on synthesis of the Fourier series, where the random number sequences are constructed from prescribed amplitude and phase spectra. A total of 200 sequences, each having 601 random numbers, for each generator were collected and analyzed in terms of the mean, variance, normality, autocorrelation, and power spectral density. These sequences were then applied to two problems in aircraft dynamic modeling, namely estimating stability and control derivatives from simulated onboard sensor data, and simulating flight in atmospheric turbulence. In general, each random number generator had good performance and is well-suited for aircraft dynamic modeling applications. Specific strengths and weaknesses of each generator are discussed. For Monte Carlo simulation, the Fourier synthesis method is recommended because it most accurately and consistently approximated Gaussian white noise and can be implemented with reasonable computational effort.

Grauer, Jared A.↗
NASA NTRS1975

Interaction between a shock wave and a turbulent boundary layer in transonic flow

Interaction between a shock wave and an unseparated turbulent boundary layer is considered. The method of matched asymptotic expansions is used, with solutions valid in the double limit as Reynolds number tends to infinity and Mach number tends to unity. The shock is weak enough that interaction effects can be considered as perturbations to the undisturbed flow; the case considered is that where the sonic line is near the outer edge of the boundary layer. It is shown that, with order estimates for Reynolds stress perturbations, the induced wall pressure distribution can be calculated using only the two outer interaction regions, independent of a specific closure condition and that this solution is in fact a turbulent free interaction solution. A detailed analysis of the inner regions, for which an eddy viscosity model for the Reynolds shear stress is used, provides a description of the variations in velocity, temperature and density near and at the wall.

Adamson, T. C., Jr.↗
NASA NTRS2022

Constraining Arctic Climate Projections of Wintertime Warming With Surface Turbulent Flux Observations and Representation of Surface-Atmosphere Coupling

The drivers of rapid Arctic climate change—record sea ice loss, warming SSTs, and a lengthening of the sea ice melt season—compel us to understand how this complex system operates and use this knowledge to enhance Arctic predictability. Changing energy flows sparked by sea ice decline, spotlight atmosphere-surface coupling processes as central to Arctic system function and its climate change response. Despite this, the representation of surface turbulent flux parameterizations in models has not kept pace with our understanding. The large uncertainty in Arctic climate change projections, the central role of atmosphere-surface coupling, and the large discrepancy in model representation of surface turbulent fluxes indicates that these processes may serve as useful observational constraints on projected Arctic climate change. This possibility requires an evaluation of surface turbulent fluxes and their sensitivity to controlling factors (surface-air temperature and moisture differences, sea ice, and winds) within contemporary climate models (here Coupled Model Intercomparison Project 6). The influence of individual controlling factors and their interactions is diagnosed using a multi-linear regression approach. This evaluation is done for four sea ice loss regimes, determined from observational sea ice loss trends, to control for the confounding effects of natural variability between models and observations. The comparisons between satellite- and model-derived surface turbulent fluxes illustrate that while models capture the general sensitivity of surface turbulent fluxes to declining sea ice and to surface-air gradients of temperature and moisture, substantial mean state biases exist. Specifically, the central Arctic is too weak of a heat sink to the winter atmosphere compared to observations, with implications to the simulated atmospheric circulation variability and thermodynamic profiles. Models were found to be about 50% more efficient at turning an air-sea temperature gradient anomaly into a sensible heat flux anomaly relative to observations. Further, the influence of sea ice concentration on the sensible heat flux is underestimated in models compared to observations. The opposite is found for the latent heat flux variability in models; where the latent heat flux is too sensitive to a sea ice concentration anomaly. Lastly, the results suggest that present-day trends in sea ice retreat regions may serve as suitable observational constraints of projected Arctic warming.

turbulent fluxes↗
NASA NTRS1986

Ion flow at comet Halley

The three-dimensional positive ion analyzer aboard the Giotto spacecraft has been used to study the interaction between protons and alpha-particles in the solar wind and positive ions from comet Halley. Although the first impression of the overall structure is that the plasma flow evolves smoothly as the nucleus is approached, three sharp transitions of relatively small amplitude can be identified on both the inbound and outbound legs of the trajectory. The outermost one, at about one million km from the nucleus, appears to be a multiple crossing of a weak bow shock. The innermost one, at 80,000 km, is the boundary where the flowing plasma becomes depleted. On a microscopic scale, the turbulence created by the interaction between the two ion populations extends to a distance of several million km from the nucleus. At Giotto's closest approach to the nucleus, the plasma produced around the spacecraft by dust and gas impacts was much more energetic than had been expected.

Johnstone, A.↗
DOE OSTI2023

Electrostatic turbulence in EAST plasmas with internal transport barrier

Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Are Supercells Resistant to Entrainment because of Their Rotation?

This research investigates a hypothesis posed by previous authors, which argues that the helical nature of the flow in supercell updrafts makes them more resistant to entrainment than nonsupercellular updrafts because of the suppressed turbulence in purely helical flows. It was further supposed that this entrainment resistance contributes to the steadiness and longevity of supercell updrafts. A series of idealized large-eddy simulations were run to address this idea, wherein the deep-layer shear and hodograph shape were varied, resulting in supercells in the strongly sheared runs, nonsupercells in the weakly sheared runs, and variations in the percentage of streamwise vorticity in updrafts among runs. Fourier energy spectrum analyses show well-developed inertial subranges in all simulations, which suggests that the percentages of streamwise and crosswise vorticity have little effect on turbulence in convective environments. Additional analyses find little evidence of updraft-scale centrifugally stable flow within updrafts, which has also been hypothesized to limit horizontal mass flux across supercell updrafts. Results suggest that supercells do have smaller fractional entrainment rates than nonsupercells, but these differences are consistent with theoretical dependencies of entrainment on updraft width, and with supercells being wider than nonsupercells. Thus, while supercells do experience reduced fractional entrainment rates and entrainment-driven dilution, this advantage is primarily attributable to increased supercell updraft width relative to ordinary convection, and has little to do with updraft helicity and rotation.

Peters, John M.↗
DOE OSTI2020

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Benchmarking of three DWM-based wake models at below-rated wind speeds

Wind turbine wake models are essential tools for predicting power losses and structural loads in wind farms. Among these, the dynamic wake meandering (DWM) model, included as a recommended approach in the International Electrotechnical Commission design standard, is a widely used engineering-fidelity method that balances accuracy and computational cost. This study compares the performance of three DWM-based wake model implementations (from the Technical University of Denmark, the National Renewable Energy Laboratory, and the Institute for Energy Technology) under below-rated wind speed conditions. Model predictions of wake flow, power output, and structural loads for a four-turbine row are evaluated across different ambient turbulence levels and wind-direction misalignments and compared against high-fidelity large-eddy simulation results. All three models captured the overall wake evolution and mean turbine performance with reasonable accuracy; their predicted time-averaged thrust and power were typically within 5 %–10 % of the large-eddy simulation benchmark. However, notable differences emerged in wake structure and unsteady load predictions, with discrepancies increasing for turbines further downstream. These differences highlight the importance of modelling choices such as wake summation and turbulence treatment, which strongly influence power-deficit and fatigue-load predictions. Comparison with large-eddy simulations reveals each approach's strengths and weaknesses, indicating where improvements are needed. Overall, the findings point to specific refinements for DWM models to improve their fidelity, ultimately enabling more robust wake predictions for wind farm design and operation.

17 WIND ENERGY↗
NASA NTRS1984

Mixing and the dynamics of the deep chlorophyll maximum in Lake Tahoe

Chlorophyll-temperature profiles were measured across Lake Tahoe about every 10 days from April through July 1980. Analysis of the 123 profiles and associated productivity and nutrient data identified three important processes in the formation and dynamics of the deep chlorophyll maximum (DCM): turbulent diffusion, nutrient supply rate, and light availability. Seasonal variation in these three processes resulted in three regimes: a diffusion-dominated regime with a weak DCM, a variable-mixing regime with a pronounced, nutrient supply-dominated DCM, and a stable regime with a deep, moderate light availability-dominated DCM. The transition between the first two regimes occurred in about 10 days, the transition between the last two more gradually over about 3 weeks. The degree of spatial variability of the DCM was highest in the second regime and lowest in the third. These data indicate that the DCM in Lake Tahoe is constant in neither time nor space.

Abbott, M. R.↗
NASA NTRS1992

Numerical investigation of tail buffet on F-18 aircraft

Numerical investigation of vortex induced tail buffet is conducted on the F-18 aircraft at high angles of attack. The Reynolds-averaged Navier-Stokes equations are integrated using a time-accurate, implicit procedure. A generalized overset zonal grid scheme is used to decompose the computational space around the complete aircraft with faired-over inlet. A weak coupling between the aerodynamics and structures is assumed to compute the structural oscillation of the flexible vertical tail. Time-accurate computations of the turbulent flow around the F-18 aircraft at 30 degrees angle of attack show the surface and off-surface flowfield details, including the unsteadiness created by the vortex burst and its interaction with the vertical twin tail which causes the tail buffet. The effect of installing a LEX fence on modifying the vortex structure upstream of the tail is also examined.

Rizk, Yehia M.↗
NASA NTRS1984

Investigations of flowfields found in typical combustor geometries

Studies are concerned with experimental and theoretical research on 2-D axisymmetric geometries under low speed, nonreacting, turbulent, swirling flow conditions. The flow enters the test section and proceeds into a larger chamber (the linear expansion ratio D/d = 2, 1.5 and 1) via a sudden or gradual expansion (side wall angle alpha = 90 and 45 degrees). A weak or strong nozzle (of area ratio A/a = 2 and 4) may be positioned downstream at x/D = 2 to form a contraction exit to the test section. Inlet swirl vanes are adjustable to a variety of vane angles with values of theta = 0, 38, 45, 60 and 70 degrees being emphasized. The objective is to determine the effect of these parameters on isothermal flow field patterns, time mean velocities and turbulence quantities, and to establish an improved simulation in the form of a computer prediction code equipped with a suitable turbulence model. The goal of the on going research is to perform experiments and complementary computations with the idea of doing the necessary type of research that will yield improved calculation capability. This involves performing experiments where time mean turbulence quantities are measured and taking input conditions and running an existing prediction code for a variety of test cases so as to compare predictions against experiment.

Lilley, D. G.↗
NASA NTRS1995

Saturn's central flash from the 3 July 1989 occultation of 28 Sgr

We present observations of Saturn's central flash obtained from Palomar and McDonald Observatories during the 3 July 1989 occultation of 28 Sgr. As the star passed close to the geometric center of Saturn's shadow, the focusing of the incident starlight by the planet's atmosphere formed multiple stellar images along the limb which were detected in infrared images obtained at wavelengths of 3.9 (Palomar) and 2.1 microns (McDonald). These are the first reported observations of a central flash due to Saturn, and the first of any planet in which the signal flash from each stellar image could be determined separately, permitting a comparison of both intensity and position for each image with model predictions. Four separate flashes were observed from each observatory, corresponding to points on the limb where the starlight passed through the Cassini Division and the relatively transparent C Ring, with maximum brightness reaching 1-2% of the unocculted stellar intensity. The timing of the flashes is quite sensitive to the shape of Saturn's limb, which depends in turn on the planet's zonal gravity harmonics and on the zonal wind profile in the lower stratosphere, near the 2.5-mbar pressure level. The locations of the images along the limb, as well as the timing, shapes, and amplitudes of the individual flash light curves, are well matched by a smoothed model based on a tropospheric zonal wind profile obtained from tracking cloud features in Voyager images and the Saturn ring optical depth profile obtained from the Voyager Photopolarimeter experiment. The smoothing required to give the best match to the data exceeds that attributable to the finite angular extent of the occulted star and may be due to refractive scattering by turbulence or wave structure in Saturn's atmosphere. There is no evidence of significant atmospheric absorption at the observed wavelengths, which correspond to spectral regions of weak methane absorption.

Nicholson, Philip D.↗
NASA NTRS1983

Properties of echo spectra observed by MST radars

Turbulent scatter and Fresnel reflection are the fundamental echoing mechanisms to interpret the signals observed by Mesosphere-stratosphere-troposphere (MST) radars. Turbulent scattered echoes provide information about the turbulence structure and mean flow of the atmosphere. Observational results with VHF MST radars, however, show the importance of Fresnel reflection due to the infinite gradient of reflectivity at the edges of a scattering layer. This condition is excluded for the weak fluctuation models but it is still possible to include the observed aspect sensitivity by assuming an anisotropic structure of fluctuations. Another explanation of the aspect sensitivity observed by MST radars is advanced. Spectral estimates by the widely used periodogram were related to a four-dimensional spectrum of atmospheric fluctuations with anisotropic structure. Effects of the radar system such as antenna beam width, beam direction and Fast Fourier Transformations (FFT) data length were discussed for the anisotropic turbulent atmosphere. Echo parameters were also estimated.

Wakasugi, K.↗