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 235 records · Page 13

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗

Spectra of Magnetic Turbulence in a Relativistic Plasma

Abstract We present a phenomenological and numerical study of strong Alfvénic turbulence in a magnetically dominated collisionless relativistic plasma with a strong background magnetic field. In contrast with the nonrelativistic case, the energy in such turbulence is contained in magnetic and electric fluctuations. We argue that such turbulence is analogous to turbulence in a strongly magnetized nonrelativistic plasma in the regime of broken quasi-neutrality. Our 2D particle-in-cell numerical simulations of turbulence in a relativistic pair plasma find that the spectrum of the total energy has the scaling k −3/2 , while the difference between the magnetic and electric energies, the so-called residual energy, has the scaling k −2.4 . The electric and magnetic fluctuations at scale ℓ exhibit dynamic alignment with the alignment angle scaling close to cos ϕ ℓ ∝ ℓ 1 / 4 . At scales smaller than the (relativistic) plasma inertial scale, the energy spectrum of relativistic inertial Alfvén turbulence steepens to k −3.5 .

79 ASTRONOMY AND ASTROPHYSICS↗

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas

The mechanisms that generate “seed” magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, here we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ B ≃0.4⁢u rms /L, where L/2⁢π is the driving scale and u rms is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k int ≃12⁢π/L. Our results show that turbulence—induced by, e.g., the gravitational buildup of galaxies and galaxy clusters—can magnetize collisionless plasmas with large-scale near-equipartition fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of Low-Altitude Atmospheric Turbulence Models for Aircraft Aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.

atmospheric turbulence↗

Insensitivity of a turbulent laser-plasma dynamo to initial conditions

It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetized, laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one arising from two pre-magnetized plasma jets whose creation is identical save for the addition of a strong external magnetic field imposed by a pulsed magnetic field generator. We investigate the differences between the two turbulent systems using a Thomson-scattering diagnostic, x-ray self-emission imaging, and proton radiography. The Thomson-scattering spectra and x-ray images suggest that the external magnetic field has a limited effect on the plasma dynamics in the experiment. Although the external magnetic field induces collimation of the flows in the colliding plasma jets and although the initial strengths of the magnetic fields arising from the interaction between the colliding jets are significantly larger as a result of the external field, the energies and morphologies of the stochastic magnetic fields post-amplification are indistinguishable. We conclude that, for turbulent laser-plasmas with supercritical magnetic Reynolds numbers, the dynamo-amplified magnetic fields are determined by the turbulent dynamics rather than the seed fields or modest changes in the initial flow dynamics of the plasma, a finding consistent with theoretical expectations and simulations of turbulent dynamos.

47 OTHER INSTRUMENTATION↗

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Transducer Resolution Effect on Pressure Fluctuations Beneath Hypersonic Turbulent Boundary Layers

The size of a pressure transducer is known to affect the accuracy of measurements of wall-pressure fluctuations beneath a turbulent boundary layer because of spatial averaging over the sensing area of the transducer. In this paper, the effect of finite transducer size is investigated by applying spatial averaging or wavenumber filters to a database of hypersonic wall pressure generated from a direct numerical simulation (DNS) that simulates the turbulent portion of the boundary layer over a sharp 7° half-angle cone at nominally Mach 8. Here, a good comparison between the DNS and the experiment in the Sandia Hypersonic Wind Tunnel at Mach 8 is achieved after spatial averaging is applied to the DNS data over an area similar to the sensing area of the transducer. The study shows that a finite sensor size similar to that of the PCB132 transducer can cause significant attenuation in the root-mean-square and power spectral density (PSD) of wall-pressure fluctuations, and the attenuation effect is identical between cone and flat plate configurations at the same friction Reynolds number. The Corcos theory is found to successfully compensate for the attenuated high-frequency components of the wall-pressure PSD.

42 ENGINEERING↗

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗

Potential quantum advantage for simulation of fluid dynamics

Numerical simulation of turbulent fluid dynamics needs to either parametrize turbulence—which introduces large uncertainties—or explicitly resolve the smallest scales—which is prohibitively expensive. Here, we provide evidence through analytic bounds and numerical studies that a potential quantum speedup can be achieved to simulate fluid dynamics using quantum computing. Specifically, we provide a lattice Boltzmann formulation of fluid dynamics for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems. This is achieved via a combination of reformulating the Navier-Stokes nonlinearity (u·$\triangledown$u) to lattice-Boltzmann nonlinearity (u 2 ) and accurately linearizing the dynamical equations, which effectively trades nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this, we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size compared with polynomial scaling in the best known classical algorithms. In this paper, we suggest that a quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.

42 ENGINEERING↗

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

AEflow (Autoencoder fluid flow compression network) [SWR-22-29]

As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.

King, Ryan↗

Validation of gyrokinetic simulations in NSTX and projections for high-k turbulence measurements in NSTX-U

An extensive validation effort performed for a modest-beta NSTX NBI-heated H-mode discharge predicts that electron thermal transport can be entirely explained by electron-scale turbulence fluctuations driven by the electron temperature gradient mode (ETG), both in conditions of strong and weak ETG turbulence drive. Thermal power-balance estimates computed by TRANSP as well as the shape of the high-k density fluctuation wavenumber spectrum and the fluctuation level ratio between strongly driven and weakly driven ETG-turbulence conditions can be matched by nonlinear gyrokinetic simulations and a synthetic diagnostic for high-k scattering. Linear gyrokinetic simulations suggest that the ion-scale instability in the weak ETG condition is close to the critical threshold for the kinetic ballooning mode instability, and nonlinear ion-scale gyrokinetic simulations show that turbulence might be in a state reminiscent of a Dimits' shift regime, opening speculation on the role that ion-scale turbulence might play for the weak ETG condition. A simulation that matched all experimental constraints is chosen to project high-k turbulence spectra in NSTX-U, revealing that the new high-k system [R. Barchfeld et al., Rev. Sci. Instrum. 89, 10C114 (2018)] should be sensitive to density fluctuations from radially elongated streamer structures. In this study, two schemes are designed to characterize the radial and poloidal wavenumber dependence of the density fluctuation wavenumber power spectrum around the streamer peak, suggesting future high-k fluctuation measurements could be sensitive to an asymmetry in the k r spectrum introduced due to the presence of strong background flow shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal chains and entrainment in cumulus updrafts, Part 2: Analysis of idealized simulations

Research has suggested that the structure of deep convection often consists of a series of rising thermals, or “thermal chain”, which contrasts with existing conceptual models that are used to construct cumulus parameterizations. In this work, simplified theoretical expressions for updraft properties obtained in Part 1 of this study are used to develop a hypothesis explaining why this structure occurs. In this hypothesis, cumulus updraft structure is strongly influenced by organized entrainment below the updraft’s vertical velocity maximum. In a dry environment, this enhanced entrainment can locally reduce condensation rates and increase evaporation, thus eroding buoyancy. For moderate-to-large initial cloud radius R, this breaks up the updraft into a succession of discrete pulses of rising motion (i.e., a thermal chain). For small R, this leads to the structure of a single, isolated rising thermal. In contrast, moist environments are hypothesized to favor plume-like updrafts for moderate-to-large R. In a series of axisymmetric numerical cloud simulations, R and environmental relative humidity (RH) are systematically varied to test this hypothesis. Vertical profiles of fractional entrainment rate, passive tracer concentration, buoyancy, and vertical velocity from these runs agree well with vertical profiles calculated from the theoretical expressions in Part 1. Analysis of the simulations supports the hypothesized dependency of updraft structure on R and RH, that is, whether it consists of an isolated thermal, a thermal chain, or a plume, and the role of organized entrainment in driving this dependency. Additional 3-dimensional (3-D) turbulent cloud simulations are analyzed, and the behavior of these 3-D runs is qualitatively consistent with the theoretical expressions and axisymmetric simulations.

54 ENVIRONMENTAL SCIENCES↗

High-Fidelity Simulations of a Mach 9 Shock-Turbulent Boundary Layer Interaction

High-fidelity large-eddy simulations (LES) and wall-modeled large-eddy simulations (WMLES) are performed on the interaction of a Mach 9 turbulent boundary layer with a 34-degree compression ramp. Our primary interest is the sensitivity of mean and unsteady wall quantities to grid density. Additionally within the context of WMLES, we seek to explore solution sensitivity to the wall-model exchange location.

42 ENGINEERING↗

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗