Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulent 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 469 records · Page 26

CMG-Augmented Control of a Hovering VTOL Platform

This paper describes how Control Moment Gyroscopes (CMGs) can be used for stability augmentation to a thrust vectoring system for a generic Vertical Take-Off and Landing platform. The response characteristics of the platform which uses only thrust vectoring and a second configuration which includes a single-gimbal CMG array are simulated and compared for hovering flight while subject to severe air turbulence. Simulation results demonstrate the effectiveness of a CMG array in its ability to significantly reduce the agility requirement on the thrust vectoring system. Albeit simplifying physical assumptions on a generic CMG configuration, the numerical results also suggest that reasonably sized CMGs will likely be sufficient for a small hovering vehicle.

Lim, K. B.↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed in the K-site flightweight hydrogen storage tank at NASA Glenn Research Center. Laminar and turbulent simulations together with conjugated heat transfer analysis are performed. Effects of turbulence, as well as tank wall conduction are presented and discussed. Predicted tank pressures and fluid temperatures are compared with the experimental data at two different heat loads for validating the CFD model.

phase change↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for modeling two-phase flow, as well as interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed using the K-site flightweight hydrogen storage tank at NASA Glenn Research Center1. Laminar and turbulent simulations are performed together with conjugated heat transfer analysis. Effects of turbulence, the value of accommodation coefficient used for predicting phase change rates, as well as tank wall geometry are presented and discussed. Predicted tank pressures, fluid and wall temperatures are compared with the experimental data at 49% fill level and two different heat loads for validating this CFD model.

Computational Fluid Dynamics↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for modeling two-phase flow, as well as interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model is validated against self-pressurization experiment performed using the K-site flightweight hydrogen storage tank at NASA Glenn Research Center 1 . Laminar and turbulent simulations are performed together with conjugate heat transfer analysis. Effects of turbulence, the value of accommodation coefficient used for predicting phase change rates, as well as tank wall geometry are presented and discussed. Predicted tank pressures, fluid and wall temperatures are compared with the experimental data at 49% fill level and two different heat loads for validating this CFD model.

Computational Fluid Dynamics↗

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗

Direct simulation of homogeneous turbulence at low Reynolds numbers

The results of direct simulations of anisotropic turbulent flows, wherein the mean strain causes growth in the length and velocity scales of the turbulence, are summarized. Attention was given to plane strain, axisymmetric strain, uniform shear, and uniform rotation turbulent flows. The strains were limited to four and total shears to about 10. Quantitative accuracy was obtained whereever experimental data were available. The simulation results provide a data base for investigating terms in the time-averaged Navier-Stokes equations.

Rogallo, R. S.↗

SPECTRA OF STRONG MAGNETOHYDRODYNAMIC TURBULENCE FROM HIGH-RESOLUTION SIMULATIONS

Magnetohydrodynamic (MHD) turbulence is present in a variety of solar and astrophysical environments. Solar wind fluctuations with frequencies lower than 0.1 Hz are believed to be mostly governed by Alfvénic turbulence with particle transport depending on the power spectrum and the anisotropy of such turbulence. Recently, conflicting spectral slopes for the inertial range of MHD turbulence have been reported by different groups. Spectral shapes from earlier simulations showed that MHD turbulence is less scale-local compared with hydrodynamic turbulence. This is why higher-resolution simulations, and careful and rigorous numerical analysis is especially needed for the MHD case. In this Letter, we present two groups of simulations with resolution up to 4096{sup 3}, which are numerically well-resolved and have been analyzed with an exact and well-tested method of scaling study. Our results from both simulation groups indicate that the asymptotic power spectral slope for all energy-related quantities, such as total energy and residual energy, is around –1.7, close to Kolmogorov's –5/3. This suggests that residual energy is a constant fraction of the total energy and that in the asymptotic regime of Alfvénic turbulence magnetic and kinetic spectra have the same scaling. The –1.5 slope for energy and the –2 slope for residual energy, which have been suggested earlier, are incompatible with our numerics.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulent transport in the solar nebula

This paper describes the current state of an ongoing project to simulate turbulent flow in a solar nebula, which is the flattened disk of dust and gas out of which a solar system forms. The goal of this project is to determine a model for the transport of mass and angular momentum in the nebula. The nebula flow exhibits compressibility, thermal conduction, viscosity, internal heating through viscous dissipation, a stable shear due to Keplerian rotation, and a gravitational acceleration in the vertical direction which is linear with altitude. These properties combine to give flow patterns not seen in terrestrial applications. Primordial solar systems are known to exist and are presumably undergoing an evolution similar to the early stages of our own solar system; for example, the IRAS infrared telescope has discovered such a protoplanetary system around the star Vega. Solar nebula evolution is the subject of much research in the astrophysical community. In the long run, researchers hope to gain a better understanding of planetary formation and the processes which dissipate the solar nebula with time.

Thompson, K. W.↗

The Energy Exascale Earth System Model Simulations With High Vertical Resolution in the Lower Troposphere

Abstract General circulation models (GCMs) are typically run with coarse vertical resolution. For example, the Energy Exascale Earth System Model (E3SM) has a vertical resolution of about 200 m in the boundary layer, which is far too coarse to resolve sharp gradients often found in the thermodynamic fields capping subtropical marine stratocumulus. In this article, we present a series of multiyear atmosphere only simulations of E3SM version 1 where we progressively increase the vertical resolution in the lower troposphere to scales approaching those often used in large eddy simulation (LES). We report marginal impacts in regards to the simulation of boundary layer clouds when vertical resolution is moderately increased, yet find significant positive impacts when the vertical resolution approaches that typically used in LES (∼10 m). In these experiments, there is a marked change in the simulated turbulence and thermodynamics which leads to more abundant marine stratocumulus. However, these simulations are burdened with excessive computational cost. They are also subject to degradations in overall climate metrics due to time step sensitivities and because some processes and parameterizations are sensitive to changes in the vertical resolution.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulation of a Convective Turbulence Encounter

A numerical simulation of a convective turbulence event is investigated and compared with observational data. The numerical results show severe turbulence of similar scale and intensity to that encountered during the test flight. This turbulence is associated with buoyant plumes that penetrate the upper-level thunderstorm outflow. The simulated radar reflectivity compares well with that obtained from the aircraft's onboard radar. Resolved scales of motion as small as 50 m are needed in order to accurately diagnose aircraft normal load accelerations. Given this requirement, realistic turbulence fields may be created by merging subgrid-scales of turbulence to a convective-cloud simulation. A hazard algorithm for use with model data sets is demonstrated. The algorithm diagnoses the RMS normal loads from second moments of the vertical velocity field and is independent of aircraft motion.

Proctor, Fred H.↗

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

Implementation and Validation of the BHR Turbulence Model in the FLAG Hydrocode

The BHR-2 turbulence model, developed at Los Alamos National Laboratory for variable density and compressible flows, is implemented in an Arbitrary Lagrangian-Eulerian hydrocode, FLAG. The BHR-2 formulation is discussed, with emphasis on its connection to multi-component flow formulations that underlie FLAG's treatment of multi-species flow. One-dimensional and two-dimensional validation tests are performed and compared to experiment and Eulerian simulations. Turbulence is an often studied and ubiquitous phenomenon in nature, and modeling its effects is essential in many practical applications. Specifically the behavior of turbulence in the presence of strong density gradients and compressibility is of fundamental importance in applications ranging from Inertial Confinement Fusion (ICF) [1], supernovae [2], and atmospheric flows. The BHR closure approach [3] seeks to model the physical processes at work in variable density turbulence including Kelvin-Helmholtz (KH) [4], Rayleigh-Taylor (RT) [5], and Richtmyer-Meshkov (RM) [6], driven turbulence. The effectiveness of the BHR-2 implementation has been demonstrated for variable density mixing in the KH, RT, and RM cases in an Eulerian framework [7]. The primary motivation of the present work is to implement the BHR-2 turbulence model in the Arbitrary Lagrangian-Eulerian (ALE) hydrodynamics code FLAG. The goal is not only to demonstrate results in agreement with previous Eulerian calculations, but also document behavior that arises from the underlying differences in code philosophy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of gyrokinetic simulations with measurements of electron temperature fluctuations and density-temperature phase angles on ASDEX Upgrade

Measurements of turbulent electron temperature fluctuation amplitudes, $δT_e⊥/T_e$, frequency spectra, and radial correlation lengths, $L_r(T_{e⊥}$), have been performed at ASDEX Upgrade using a newly upgraded Correlation ECE diagnostic in the range of scales $k_⊥<1.4 cm^{-1}, k_r < 3.5 cm^{-1} (k_⊥ρ_s < 0.28$ and $k_rρ_s < 0.7$). The phase angle between turbulent temperature and density fluctuations, α nT , has also been measured by using an ECE radiometer coupled to a reflectometer along the same line of sight. These quantities are used simultaneously to constrain a set of ion-scale non-linear gyrokinetic turbulence simulations of the outer core ( ρ tor = 0.75) of a low density, electron heated L-mode plasma, performed using the gyrokinetic simulation code, GENE. The ion and electron temperature gradients were scanned within uncertainties. It is found that gyrokinetic simulations are able to match simultaneously the electron and ion heat flux at this radius within the experimental uncertainties. The simulations were performed based on a reference discharge for which $δT_{e⊥}/T_e$ measurements were available, and $L_r(T_{e⊥}$) and α nT were then predicted using synthetic diagnostics prior to measurements in a repeat discharge. While temperature fluctuation amplitudes are overestimated by >50% for all simulations within the sensitivity scans performed, good quantitative agreement is found for $L_r(T_{e⊥}$) and α nT . Finally, a validation metric is used to quantify the level of agreement of individual simulations with experimental measurements, and the best agreement is found close to the experimental gradient values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resolving the mystery of transport within internal transport barriers

The Trapped Gyro-Landau Fluid (TGLF) quasi-linear model [G. M. Staebler, et al., Phys. Plasmas 12, 102508 (2005)], which is calibrated to nonlinear gyrokinetic turbulence simulations, is now able to predict the electron density, electron and ion temperatures, and ion toroidal rotation simultaneously for internal transport barrier (ITB) discharges. This is a strong validation of gyrokinetic theory of ITBs, requiring multiple instabilities responsible for transport in different channels at different scales. The mystery of transport inside the ITB is that momentum and particle transport is far above the predicted neoclassical levels in apparent contradiction with the expectation from the theory of suppression of turbulence by E×B velocity shear. The success of TGLF in predicting ITB transport is due to the inclusion of ion gyro-radius scale modes that become dominant at high E×B velocity shear and to improvements to TGLF that allow momentum transport from gyrokinetic turbulence to be faithfully modeled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SIMULATED FARADAY ROTATION MEASURES TOWARD HIGH GALACTIC LATITUDES

We study the Faraday rotation measure (RM) due to the Galactic magnetic field (GMF) toward high Galactic latitudes. The RM arises from the global, regular component as well as from the turbulent, random component of the GMF. We model the former based on observations and the latter using the data of magnetohydrodynamic turbulence simulations. For a large number of different GMF models, we produce mock RM maps around the Galactic poles and calculate various statistical quantities with the RM maps. We find that the observed medians of RMs toward the north and south Galactic poles, {approx}0.0 {+-} 0.5 rad m{sup -2} and {approx} + 6.3 {+-} 0.5 rad m{sup -2}, are difficult to explain with any of our many alternate GMF models. The standard deviation of observed RMs, {approx}9 rad m{sup -2}, is clearly larger than that of simulated RMs. The second-order structure function of observed RMs is substantially larger than that of simulated RMs, especially at small angular scales. We discuss other possible contributions to RM toward high Galactic latitudes. Besides observational errors and the intrinsic RM of background radio sources against which RM is observed, we suggest that the RM due to the intergalactic magnetic field may account for a substantial fraction of the observed RM. Finally, we note that reproducing the observed medians may require additional components or/and structures of the GMF that are not present in our models.

79 ASTRONOMY AND ASTROPHYSICS↗