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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
DOE OSTI2025

Coupled Climate Simulations With E3SM‐MMF

Simulations of the recent historical period from 1950 to 2014 are conducted with E3SM‐MMF, which uses an embedded 2D cloud resolving model that runs efficiently on GPUs in place of traditional parameterizations for cloud and turbulence. Analysis of the climate and variability reveal several aspects where E3SM‐MMF produces smaller biases compared to E3SMv2, including better agreement with the observed evolution of global mean surface temperature, although the representation of ENSO is too weak and fast. Three idealized abrupt CO 2 experiments were also conducted to assess climate sensitivity and feedbacks. These yield three estimates of effective climate sensitivity (4.38, 5.21, and 6.06 K), with a corresponding spread in the shortwave cloud feedbacks. These estimates are on the higher end of sensitivity estimates from CMIP ensembles, and the spread indicates substantial state‐dependent feedbacks. These results demonstrate how multiscale modeling framework (MMF) models can be used for climate relevant experiments and projections by leveraging modern GPU enabled computational platforms. The unique qualities of E3SM‐MMF shown in previous literature are largely still present, but various instances of reduced biases suggest that MMF models have utility in improving future projections.

E3SM↗
NASA NTRS2006

Progress in Development of an Airborne Turbulence Detection System

Aircraft encounters with turbulence are the leading cause of in-flight injuries (Tyrvanas 2003) and have occasionally resulted in passenger and crew fatalities. Most of these injuries are caused by sudden and unexpected encounters with severe turbulence in and around convective activity (Kaplan et al 2005). To alleviate this problem, the Turbulence Prediction and Warning Systems (TPAWS) element of NASA s Aviation Safety program has investigated technologies to detect and warn of hazardous in-flight turbulence. This effort has required the numerical modeling of atmospheric convection: 1) for characterizing convectively induced turbulence (CIT) environments, 2) for defining turbulence hazard metrics, and 3) as a means of providing realistic three-dimensional data sets that can be used to test and evaluate turbulence detection sensors. The data sets are being made available to industry and the FAA for certification of future airborne turbulence-detection systems (ATDS) with warning capability. Early in the TPAWS project, a radar-based ATDS was installed and flight tested on NASA s research aircraft, a B-757. This ATDS utilized new algorithms and hazard metrics that were developed for use with existing airborne predictive windshear radars, thus avoiding the installation of new hardware. This system was designed to detect and warn of hazardous CIT even in regions with weak radar reflectivity (i.e. 5-15 dBz). Results from an initial flight test of the ATDS were discussed in Hamilton and Proctor (2002a; 2002b). In companion papers (Proctor et al 2002a; 2002b), a numerical simulation of the most significant encounter from that flight test was presented. Since the presentation of these papers a second flight test has been conducted providing additional cases for examination. In this paper, we will present results from NASA s flight test and a numerical model simulation of a turbulence environment encountered on 30 April 2002. Progress leading towards FAA certification of industry built ATDS will also be discussed.

Hamilton, David W.↗
NASA NTRS1981

The preferred-mode coherent structure in the near field of an axisymmetric jet with and without excitation

Visualization and phase-locked hot-wire measurements have been used to investigate the preferred-mode coherent structure in the near field of an axisymmetric jet for different Reynolds numbers and for the two limiting asymptotic states of the exit boundary layer: laminar and fully developed turbulent. Distributions of phase-averaged longitudinal and lateral velocities, coherent vorticity, coherent and incoherent turbulence Reynolds stresses, and coherent strain rate and production have been obtained over the spatial extent of the coherent structure for a Strouhal number of 0.3. It is found that these properties are independent of the initial conditions, but weakly dependent on the Reynolds number up to about 100,000. The naturally occurring coherent structure in the unexcited jet, educed via an optimized conditional sampling scheme, shows no clear dependence on the jet Reynolds number up to 800,000.

Hussain, A. K. M. F.↗
NASA NTRS1991

Mars boundary layer simulations - Comparison with Viking lander and entry observations

Diurnal variations of wind and temperature in the lower Martian atmosphere are simulated with a boundary layer model that includes radiative heating in a dusty CO2 atmosphere, turbulence generated by convection and/or shear stresses, a surface heat budget, and time varying pressure forces due to sloping terrain. Model results for early northern summer are compared with Viking lander observations to determine the model's strengths and weaknesses, and suitability as an engineering model.

Haberle, R. M.↗
NASA NTRS1995

Radio interferometer observations of solar wind turbulence from the orbit of Helios to the solar corona

We report observations of Very Long Baseline Interferometer (VLBI) phase scintillations due to turbulence in the solar wind. The observations were made at 5.00 and 8.42 GHz with the Very Long Baseline Array (VLBA) on three dates in 1991 July and August. We observed the sources 0851 + 202 and 0735 + 178 at solar elongations ranging from 2.66 deg to 13.29 deg; the closest approach of the line of sight to the Sun ranged from 10 to 49.8 solar radii. We have also included previously unpublished 5 GHz VLBI phase scintillation measurements from 1989. These measurements probe solar wind density fluctuations on spatial scales from about 200 km to 2000 km. Our measurements are in quite good agreement with the Coles & Harmon model for the radio phase structure function, which was largely determined from observations on both shorter and larger spatial scales. Departures from the Coles & Harmon functions are attributable to day-to-day variations in the solar wind conditions. Phase scintillations at the greatest solar elongations are in very good agreement with extrapolated estimates from direct measurements made with the Helios spacecraft at slightly larger heliocentric distances. Thus there is a consistency between the in-situ spacecraft and radio sensing measurements of density turbulence. All of the VLBI data are consistent with a Kolmogorov spectrum for the density fluctuations, although at the closest elongations there may be excess power at small spatial scales. An advantage of interferometric techniques over other radio propagation measurements is that they provide a measure of the anisotropy of the irregularities. Our observations at closest approach (10 solar radii) show weak evidence for anisotropic, field-aligned density irregularities with an axial ratio of order 2. This degree of anisotropy would appear to be less than that measured at similar solar elongations but on smaller spatial scales by Armstrong and colleagues. Finally, a combination of the radio propagation data and in situ Helios measurements is used to determine the heliocentric distance dependence of the normalization coefficient of the density power spectrum, C(sup 2)(sub N). Modeling C(sup2)(sub N) varies as (R/solar radii)(sup -Kappa), we find Kappa = 3.72 +/- 0.30. Subject to a number of assumptions, this result is consistent with the conclusion of D. A. Roberts that turbulence within 0.28 astronomical units adheres to a WKB formula for turbulence amplitude as a function of heliocentric distance.

Spangler, Steven R.↗
DOE OSTI2020

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2010

The Impact of the Evolving Satellite Data Record on Reanalysis Water and Energy Fluxes During the Past 30 Years

Retrospective analyses (reanalyses) use a fixed assimilation model to take diverse observations and synthesize consistent, time-dependent fields of state variables and fluxes (e.g. temperature, moisture, momentum, turbulent and radiative fluxes). Because they offer data sets of these quantities at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Of course, one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics. Another issue, perhaps less appreciated, is the fact that the discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (~ the last 30 years). In particular we highlight the effect on reanalyses of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We examine these data sets from two perspectives. The first is the ability to capture modes of variability that have coherent spatial structure (e.g. ENSO events and near-decadal coupling to SST changes) and how these modes are contained within trends in near global averages of key quantities. Secondly, we consider diagnostics that measure the consistency in energetic scaling in the hydrologic cycle, particularly the fractional changes in column-integrated water vapor versus precipitation as they are coupled to radiative flux constraints. These results will be discussed in the context of implications for science objectives and priorities of the NASA Energy and Water Cycle Study, NEWS.

Robertson, Franklin R.↗
NASA NTRS1991

An experimental study of fluctuating pressure loads beneath swept shock/boundary-layer interactions

A database is established on the fluctuating pressure loads produced on aerodynamic surfaces beneath 3-D shock wave/boundary layer interactions. Such loads constitute a fundamental problem of critical concern to future supersonic and hypersonic flight vehicles. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins. Fin angles from 5 to 25 deg at freestream Mach numbers between 2.5 and 4 produce a variety of interaction strengths from weak to very strong. Miniature Kulite pressure transducers mounted in the flat plate were used to measure interaction-induced wall pressure fluctuations. These data will be correlated with proposed new optical data on the fluctuations of the interaction structure, especially that of the lambda-shock system and its associated high-speed jet impingement.

Settles, Gary S.↗
NASA NTRS1994

Evaluation of turbulence models in the PARC code for transonic diffuser flows

Flows through a transonic diffuser were investigated with the PARC code using five turbulence models to determine the effects of turbulence model selection on flow prediction. Three of the turbulence models were algebraic models: Thomas (the standard algebraic turbulence model in PARC), Baldwin-Lomax, and Modified Mixing Length-Thomas (MMLT). The other two models were the low Reynolds number k-epsilon models of Chien and Speziale. Three diffuser flows, referred to as the no-shock, weak-shock, and strong-shock cases, were calculated with each model to conduct the evaluation. Pressure distributions, velocity profiles, locations of shocks, and maximum Mach numbers in the duct were the flow quantities compared. Overall, the Chien k-epsilon model was the most accurate of the five models when considering results obtained for all three cases. However, the MMLT model provided solutions as accurate as the Chien model for the no-shock and the weak-shock cases, at a substantially lower computational cost (measured in CPU time required to obtain converged solutions). The strong shock flow, which included a region of shock-induced flow separation, was only predicted well by the two k-epsilon models.

Georgiadis, N. J.↗
NASA NTRS2000

Flow Simulation of Solid Rocket Motors: Sub-Scale Air Flow Simulation of Port Flows - 2

The injection-flow issuing from a porous medium in the cold-flow simulation of internal port flows in solid rocket motors is characterized by a spatial instability termed pseudoturbulence that produces a rather non-uniform (lumpy) injection-velocity profile. The objective of this study is to investigate the interaction between the injection- and the developing axial-flows. The findings show that this interaction generally weakens the lumpy injection profile and affects the subsequent development of the axial flow. The injection profile is found to depend on the material characteristics, and the ensuing pseudoturbulence is a function of the injection velocity, the axial position and the distance from the porous wall. The flow transition (from laminar to turbulent) of the axial-flow is accelerated in flows emerging from smaller pores primarily due to the higher pseudoturbulence produced by the smaller pores in comparison to that associated with larger pores. In flows with rather uniform injection-flow profiles (weak or no pseudoturbulence), the axial and transverse velocity components in the porous duct are found to satisfy the sine/cosine analytical solutions derived from inviscid assumptions. The transition results from the present study are compared with previous results from surveyed literature, and detailed flow development measurements are presented in terms of the blowing fraction, and characterizing Reynolds numbers.

Yeh, Y. P.↗
NASA NTRS2001

Quantitative Studies on the Propagation and Extinction of Near-Limit Premixed Flames under Normal- and Micro-gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion.. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Studies at high strain rates are important in quantifying and understanding the response of vigorously burning flames and determine extinction states. Studies of weakly strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds, S(sup o)(sub u), can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames are stabilized by weak strain rates. Previous studies have shown that near-limit flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In our past studies the laminar flame speeds and extinction strain rates were directly measured at ultra-low strain rates. The laminar flame speeds were measured by having a positively strained planar flame undergoing a transition to a negatively strained Bunsen flame and by measuring the propagation speed during that transition. The extinction strain rates of near-limit flames were measured in mu-g. Results obtained for CH4/air and C3H8/air mixtures are in agreement with those obtained by Maruta et al.

Egolfopoulos, F. N.↗
NASA NTRS2022

The Puzzling Structure of Solar Convection: Window into the Dynamo

The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun’s turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo’s organization and timing. Ultimately, the convective and large-scale flow structure derive from the Coriolis force. When the Coriolis effect is weak (rapid convective flow), convection exhibits little organization in its spatial structure, and the meridional flow is expected to assume a single-celled profile within each hemisphere. Convection subject to strong Coriolis forces instead organizes into compact, spiraling columnar structures and sustains meridional flows that possess multiple cells in latitude and radius. The implications for the dynamo in these two regimes are substantial, impacting the flux-transport properties of any assumed meridional flow and the convectively-driven EMF. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun’s outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.

SMD↗
NASA NTRS1984

Analytical estimates of turbulent MHD transport coefficients

Turbulent transfer rates from small-scale MHD excitations to large-scale Fourier modes are calculated algebraically, using the method of Biskamp and Welter. Three cases are considered: two-dimensional Navier-Stokes flows, two-dimensional incompressible MHD, and the weakly three-dimensional Strauss equations. In all cases, an initially large spectral gap between the small-scale and large-scale excitations is assumed, and attention focusses on the initial values of the back-transfer rates. The sign of the transfer is determined by the sign of an analytically calculable eddy viscosity and/or anomalous resistivity. We are able to confirm the results of Biskamp and Welter for the case of two-dimensional MHD, but find some differences for the case of the Strauss equations. It is argued that the Strauss equations may not exhibit an inverse cascade phenomenon for the spatially periodic case unless their initial spectra are such that the behavior is essentially that of two-dimensional MHD.

Montgomery, D.↗
NASA NTRSDate not supplied

MAPL3: A Hierarchical Framework Based on ESMF

MAPL is a hierarchical, ESMF-based, coupling framework developed to support NASA's GEOS data assimilation system, but also now supports some external applications from partner institutions. The framework continues to evolve in response novel coupling requirements and is currently undergoing a substantial re-engineering effort to improve extensibility and maintainability. Some of the major new features that will be provided in the upcoming release of MAPL3 include the following: (1) "weak" NUOPC compliance that expresses a MAPL hierarchy as single NUOPC model, (2) component extensions that represent an improvement over conventional component couplers in many contexts, (3) an expansion of the concept of import and export to include component "services" such as advection, turbulence, etc., (4) automated OpenMP threading of column-based component subtrees, (5) and (5) improved I/O scalability.

GEOS↗
NASA NTRS2007

Slow Mode Waves in the Heliospheric Plasma Sheet

We report the results of a search for waves/turbulence in the Heliospheric Plasma Sheet (HPS) surrounding the Heliospheric Current Sheet (HCS). The HPS is treated as a distinctive heliospheric structure distinguished by relatively high Beta, slow speed plasma. The data used in the investigation are from a previously published study of the thicknesses of the HPS and HCS that were obtained in January to May 2004 when Ulysses was near aphelion at 5 AU. The advantage of using these data is that the HPS is thicker at large radial distances and the spacecraft spends longer intervals inside the plasma sheet. From the study of the magnetic field and solar wind velocity components, we conclude that, if Alfven waves are present, they are weak and are dominated by variations in the field magnitude, B, and solar wind density, NP, that are anti-correlated.

heliospheric magnetic field↗
DOE OSTI2021

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories

High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear – correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view – the size of the camera focal plane for the Vera C. Rubin Observatory in Chile – using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗