Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulence theory and modelling”

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 73 records · Page 4

A new and alternative look at nonlinear Alfvénic states

The formulation for studying nonlinear Alfvenic states, sustainable in Hall Magnetohydrodynamics (HMHD), becomes considerably simpler and more tractable when circularly polarized Beltrami vectors (the eigenstates of linear HMHD) are used as the basis functions. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange/depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfvenic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a “universal” picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: (1) the energy tends to flow from lower to higher k and (2) the higher k z (in the direction of the ambient magnetic field) components of a mode with a given k are depleted/oscillate faster—in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gain degradation and amplitude scintillation due to tropospheric turbulence

It is shown that a simple physical model is adequate for the prediction of the long term statistics of both the reduced signal levels and increased peak-to-peak fluctuations. The model is based on conventional atmospheric turbulence theory and incorporates both amplitude and angle of arrival fluctuations. This model predicts the average variance of signals observed under clear air conditions at low elevation angles on earth-space paths at 2, 7.3, 20 and 30 GHz. Design curves based on this model for gain degradation, realizable gain, amplitude fluctuation as a function of antenna aperture size, frequency, and either terrestrial path length or earth-space path elevation angle are presented.

Theobold, D. M.↗

Validation of Heat Transfer and Film Cooling Capabilities of the 3-D RANS Code TURBO

The capabilities of the 3-D unsteady RANS code TURBO have been extended to include heat transfer and film cooling applications. The results of simulations performed with the modified code are compared to experiment and to theory, where applicable. Wilcox s k-turbulence model has been implemented to close the RANS equations. Two simulations are conducted: (1) flow over a flat plate and (2) flow over an adiabatic flat plate cooled by one hole inclined at 35 to the free stream. For (1) agreement with theory is found to be excellent for heat transfer, represented by local Nusselt number, and quite good for momentum, as represented by the local skin friction coefficient. This report compares the local skin friction coefficients and Nusselt numbers on a flat plate obtained using Wilcox's k-model with the theory of Blasius. The study looks at laminar and turbulent flows over an adiabatic flat plate and over an isothermal flat plate for two different wall temperatures. It is shown that TURBO is able to accurately predict heat transfer on a flat plate. For (2) TURBO shows good qualitative agreement with film cooling experiments performed on a flat plate with one cooling hole. Quantitatively, film effectiveness is under predicted downstream of the hole.

Shyam, Vikram↗

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

Comparison of private flux region instability in conventional and super-X divertor configurations

Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kinetic modeling of neutral transport for a continuum gyrokinetic code

In this work we present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper, the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. Furthermore, it is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On Galilean invariance of mean kinetic helicity

While kinetic helicity is not Galilean invariant locally, it is known [Moffatt, J. Fluid Mech. 35, 117 (1969)] that its spatial integral quantifies the degree of knottedness of vorticity field lines. Being a topological property of the flow, mean kinetic helicity is Galilean invariant. Here, we provide a direct mathematical proof that kinetic helicity is Galilean invariant when spatially integrated over regions enclosed by vorticity surfaces, i.e., surfaces of zero vorticity flux. Here, we also discuss so-called relative kinetic helicity, which is Galilean invariant when integrated over any region in the flow.

42 ENGINEERING↗

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beam Propagation Through Atmospheric Turbulence Using an Altitude-Dependent Structure Profile with Non-Uniformly Distributed Phase Screens

For free-space optical communication systems, numerical wave optics simulations provide a useful technique for modeling turbulence-induced beam degradation when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such applications we use a split-step method modeling the turbulence as a series of random phase screens using the Hufnagel-Valley turbulence profile. We employ a space-to-ground case study to examine the irradiance and phase statistics for uniformly and non-uniformly located screens and find better agreement with theory using a non-uniform discretization minimizing the contribution of each screen to the total scintillation. In this poster, we summarize the method and the results of the case study including a comparison to layered models used in astronomical imaging applications.

Chahine, Yousef K.↗

On the Yakhot-Orszag renormalization group method for deriving turbulence statistics and models

An independent, comprehensive, critical review of the 'renormalization group' (RNG) theory of turbulence developed by Yakhot and Orszag (1986) is provided. Their basic theory for the Navier-Stokes equations is confirmed, and approximations in the scale removal procedure are discussed. The YO derivations of the velocity-derivative skewness and the transport equation for the energy dissipation rate are examined. An algebraic error in the derivation of the skewness is corrected. The corrected RNG skewness value of -0.59 is in agreement with experiments at moderate Reynolds numbers. Several problems are identified in the derivation of the energy dissipation rate equations which suggest that the derivation should be reformulated.

Smith, L. M.↗

Role of secondary instability theory and parabolized stability equations in transition modeling

In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.

El-Hady, Nabil M.↗

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modelling

ABSTRACT Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here, we introduce an error propagation formula that relaxes both limitations, allowing for non-linear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum – the latter being a non-linear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory.

79 ASTRONOMY AND ASTROPHYSICS↗