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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

On recent advances and future research directions for computational fluid dynamics

This paper highlights some recent accomplishments regarding CFD numerical algorithm constructions for generation of discrete approximate solutions to classes of Reynolds-averaged Navier-Stokes equations. Following an overview of turbulent closure modeling, and development of appropriate conservation law systems, a Taylor weak-statement semi-discrete approximate solution algorithm is developed. Various forms for completion to the final linear algebra statement are cited, as are a range of candidate numerical linear algebra solution procedures. This development sequence emphasizes the key building blocks of a CFD RNS algorithm, including solution trial and test spaces, integration procedure and added numerical stability mechanisms. A range of numerical results are discussed focusing on key topics guiding future research directions.

Baker, A. J.↗
NASA NTRS1992

The electrodynamic and hydrodynamic phenomena in magnetically-levitated molten droplets. I - Steady state behavior

A mathematical formulation is given and computed results are presented describing the behavior of electromagnetically-levitated metal droplets under the conditions of microgravity. In the formulation the electromagnetic force field is calculated using a modification of the volume integral method and these results are then combined with the FIDAP code to calculate the steady state melt velocities. The specific computational results are presented for the conditions corresponding to the planned IML-2 Space Shuttle experiment, using the TEMPUS device, which has separate 'heating' and 'positioning' coils. While the computed results are necessarily specific to the input conditions, some general conclusions may be drawn from this work. These include the fact that for the planned TEMPUS experiments to positioning coils will produce only a weak melt circulation, while the heating coils are like to produce a mildly turbulent recirculating flow pattern within the samples. The computed results also allow us to assess the effect of sample size, material properties and the applied current on these phenomena.

Zong, Jin-Ho↗
NASA NTRS2010

Interstellar Magnetic Fields and Polarimetry of Dust Emission

Magnetic fields are an important ingredient in the stormy cosmos. Magnetic fields: (1) are intimately involved with winds from Active Galactic Nuclei (AGN) and stars (2) create at least some of the structures observed in the ISM (3) modulate the formation of clouds, cores, and stars within a turbulent medium (4) may be dynamically important in protostellar accretion disks (5) smooth weak shocks (C-shocks).

stormy cosmos↗
NASA NTRS1991

Direct numerical simulation and analysis of shock turbulence interaction

Two kinds of linear analysis, rapid distortion theory (RDT) and linear interaction analysis (LIA), were used to investigate the effects of a shock wave on turbulence. Direct numerical simulations of two-dimensional isotropic turbulence interaction with a normal shock were also performed. The results from RDT and LIA are in good agreement for weak shock waves, where the effects of shock front curvature and shock front unsteadiness are not significant in producing vorticity. The linear analyses predict wavenumber-dependent amplification of the upstream one-dimensional energy spectrum, leading to turbulence scale length scale decrease through the interaction. Instantaneous vorticity fields show that vortical structures are enhanced while they are compressed in the shock normal direction. Entrophy amplfication through the shock wave compares favorably with the results of linear analyses.

Lee, Sangsan↗
NASA NTRS2022

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗
NASA NTRS2015

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved.

Carlson, Jan-Renee↗
DOE OSTI2020

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Here, we use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number S limit (above ~104), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate ~0.01–0.13 is weakly dependent on S. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a superdiffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗
DOE OSTI2024

A two-dimensional numerical study of ion-acoustic turbulence

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Potential vorticity transport in weakly and strongly magnetized plasmas

Tangled magnetic fields, often coexisting with an ordered mean field, have a major impact on turbulence and momentum transport in many plasmas, including those found in the solar tachocline and magnetic confinement devices. We present a novel mean field theory of potential vorticity mixing in β-plane magnetohydrodynamic (MHD) and drift wave turbulence. Our results show that mean square stochastic fields strongly reduce Reynolds stress coherence. This decoherence of potential vorticity flux due to stochastic field scattering leads to suppression of momentum transport and zonal flow formation. A simple calculation suggests that the breaking of the shear-eddy tilting feedback loop by stochastic fields is the key underlying physics mechanism. Furthermore, a dimensionless parameter that quantifies the increment in power threshold is identified and used to assess the impact of stochastic field on the L-H transition. We discuss a model of stochastic fields as a resisto-elastic network.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

An interplanetary magnetic field ensemble at 1 AU

A method for calculation ensemble averages from magnetic field data is described. A data set comprising approximately 16 months of nearly continuous ISEE-3 magnetic field data is used in this study. Individual subintervals of this data, ranging from 15 hours to 15.6 days comprise the ensemble. The sole condition for including each subinterval in the averages is the degree to which it represents a weakly time-stationary process. Averages obtained by this method are appropriate for a turbulence description of the interplanetary medium. The ensemble average correlation length obtained from all subintervals is found to be 4.9 x 10 to the 11th cm. The average value of the variances of the magnetic field components are in the approximate ratio 8:9:10, where the third component is the local mean field direction. The correlation lengths and variances are found to have a systematic variation with subinterval duration, reflecting the important role of low-frequency fluctuations in the interplanetary medium.

Matthaeus, W. H.↗
NASA NTRS1985

Cascade properties of shear Alfven wave turbulence

Nonlinear three-wave interactions of linear normal modes are investigated for two-dimensional incompressible magnetohydrodynamics and the weakly three-dimensional Strauss equations in the case where a strong uniform background field B0 is present. In both systems the only resonant interaction affecting Alfven waves is caused by the shear of the background field plus the zero frequency components of the perturbation. It is shown that the Alfven waves are cascaded in wavenumber space by a mechanism equivalent to the resonant absorption at the Alfven resonance. For large wavenumbers perpendicular to B0, the cascade is described by Hamilton's ray equations, dk/dt = -(first-order) partial derivative of omega with respect to vector r, where omega includes the effects of the zero frequency perturbations.

Bondeson, A.↗
NASA NTRS1986

An interplanetary magnetic field ensemble at 1 AU

A method for calculation ensemble averages from magnetic field data is described. A data set comprising approximately 16 months of nearly continuous ISEE-3 magnetic field data is used in this study. Individual subintervals of this data, ranging from 15 hours to 15.6 days comprise the ensemble. The sole condition for including each subinterval in the averages is the degree to shich it represents a weakly time-stationary process. Averages obtained by this method are appropriate for a turbulence description of the interplanetary medium. The ensemble average correlation length obtained from all subintervals is found to be 4.9 x 10 to the 11th cm. The average value of the variances of the magnetic field components are in the approximate ratio 8:9:10, where the third component is the local mean field direction. The correlation lengths and variances are found to have a systematic variation with subinterval duration, reflecting the important role of low-frequency fluctuations in the interplanetary medium.

Matthaeus, W. H.↗
NASA NTRS1986

On the use of colour reflectivity plots to monitor the structure of the troposphere and stratosphere

The radar reflectivity, defined as the range squared corrected power of VHF radar echoes, can be used to monitor and study the temporal development of inversion layer, frontal boundaries and convective turbulence. From typical featurs of upward or downward motion of reflectivity structures, the advection/convection of cold and warm air can be predicted. High resolution color plots appear to be useful to trace and to study the life history of these structures, particularly their persistency, descent and ascent. These displays allow an immediate determination of the tropopause height as well as the determination of the tropopause structure. The life history of warm fronts, cold fronts, and occlusions can be traced, and these reflectivity plots allow detection of even very weak events which cannot be seen in the traditional meteorological data sets. The life history of convective turbulence, particular evolving from the planetary boundary layer, can be tracked quite easily. Its development into strong convection reaching the middle troposphere can be followed and predicted.

Rottger, J.↗
DOE OSTI2022

GPU-accelerated DNS of compressible turbulent flows

Here, this paper explores strategies to transform an existing CPU-based high-performance computational fluid dynamics solver, HyPar, for compressible flow simulations on emerging exascale heterogeneous (CPU+GPU) computing platforms. The scientific motivation for developing a GPU-enhanced version of HyPar is to simulate canonical turbulent flows at the highest resolution possible on such platforms. We show that optimizing memory operations and thread blocks results in 200x speedup of computationally intensive kernels compared with a CPU core. Using multiple GPUs and CUDA-aware MPI communication, we demonstrate both strong and weak scaling of our GPU-based HyPar implementation on the NVIDIA Volta V100 GPUs. We simulate the decay of homogeneous isotropic turbulence in a triply periodic box on grids with up to 1024 3 points (5.3 billion degrees of freedom) and on up to 1,024 GPUs. We compare the wall times for CPU-only and CPU+GPU simulations. The results presented in the paper are obtained on the Summit and Lassen supercomputers at Oak Ridge and Lawrence Livermore National Laboratories, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1975

Copernicus observations of Betelgeuse and Antares

Copernicus observations of the M-supergiants, alpha Ori and alpha Sco, are presented. The MgII h and k resonance lines are strongly in emission in both stars. The k line is highly asymmetric in both stars but the h line is symmetric. Upper limits for several other resonance lines are given for alpha Ori. The possibility is explored that the k line asymmetry is caused by overlying resonance lines of MnI and FeI formed in the cool circumstellar gas shells around these stars. Observations of the MnI 4030-4033 A lines are used to show that circumstellar shell absorption is too weak to explain the asymmetry. It is suggested that the absorption occurs in a cool turbulent region between the base of the circumstellar shell and the top of the chromosphere.

Bernat, A. P.↗
NASA NTRS1985

A multiple-scales model of the shock-cell structure of imperfectly expanded supersonic jets

The present investigation is concerned with the development of an analytical model of the quasi-periodic shock-cell structure of an imperfectly expanded supersonic jet. The investigation represents a part of a program to develop a mathematical theory of broadband shock-associated noise of supersonic jets. Tam and Tanna (1982) have suggested that this type of noise is generated by the weak interaction between the quasi-periodic shock cells and the downstream-propagating large turbulence structures in the mixing layer of the jet. In the model developed in this paper, the effect of turbulence in the mixing layer of the jet is simulated by the addition of turbulent eddy-viscosity terms to the momentum equation. Attention is given to the mean-flow profile and the numerical solution, and a comparison of the numerical results with experimental data.

Tam, C. K. W.↗