Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic closures”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Variable-moment fluid closures with Hamiltonian structure

Abstract Based on ideas due to Scovel–Weinstein, I present a general framework for constructing fluid moment closures of the Vlasov–Poisson system that exactly preserve that system’s Hamiltonian structure. Notably, the technique applies in any space dimension and produces closures involving arbitrarily-large finite collections of moments. After selecting a desired collection of moments, the Poisson bracket for the closure is uniquely determined. Therefore data-driven fluid closures can be constructed by adjusting the closure Hamiltonian for compatibility with kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

The development and preliminary application of an invariant coupled diffusion and chemistry model

In many real-world pollution chemical reaction problems, the rate of reaction problems, the rate of reaction may be greatly affected by unmixedness. An approximate closure scheme for a chemical kinetic submodel which conforms to the principles of invariant modeling and which accounts for the effects of inhomogeneous mixing over a wide range of conditions has been developed. This submodel has been coupled successfully with invariant turbulence and diffusion models, permitting calculation of two-dimensional diffusion of two reacting (isothermally) chemical species. The initial calculations indicate the ozone reactions in the wake of stratospheric aircraft will be substantially affected by the rate of diffusion of ozone into the wake, and in the early wake, by unmixedness.

Hilst, G. R.↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗
NASA NTRS1996

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗
NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗
NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of thier predictions.

Speziale, Charles G.↗
NASA NTRS1990

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗
NASA NTRS1978

Kinetic model of turbulence in an incompressible fluid

A statistical description of turbulence in an incompressible fluid obeying the Navier-Stokes equations is proposed, where pressure is regarded as a potential for the interaction between fluid elements. A scaling procedure divides a fluctuation into three ranks representing the three transport processes of macroscopic evolution, transport property, and relaxation. Closure is obtained by relaxation, and a kinetic equation is obtained for the fluctuation of the macroscopic rank of the distribution function. The solution gives the transfer function and eddy viscosity. When applied to the inertia subrange of the energy spectrum the analysis recovers the Kolmogorov law and its numerical coefficient.

Tchen, C. M.↗
NASA NTRSDate not supplied

Development of a Unified SHOC-EDMF PBL Parameterization in the NASA GEOS Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies.

SHOC-EDMF↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

A review of near-wall Reynolds-stress

The advances made in second-order near-wall turbulence closures are summarized. All closures examined are based on some form of high Reynolds number models for the Reynolds stress and the turbulent kinetic energy dissipation rate equations. Consequently, most near-wall closures proposed to data attempt to modify the high Reynolds number models for the dissipation rate equation so that the resultant models are applicable all the way to the wall. The near-wall closures are examined for their asymptotic behavior so that they can be compared with the proper near-wall behavior of the exact equations. A comparison of the closure's performance in the calculation of a low Reynolds number plane channel flow is carried out. In addition, the closures are evaluated for their ability to predict the turbulence statistics and the limiting behavior of the structure parameters compared to direct simulation data.

So, R. M. C.↗
NASA NTRS1990

Modeling of near-wall turbulence

An improved k-epsilon model and a second order closure model is presented for low Reynolds number turbulence near a wall. For the k-epsilon model, a modified form of the eddy viscosity having correct asymptotic near wall behavior is suggested, and a model for the pressure diffusion term in the turbulent kinetic energy equation is proposed. For the second order closure model, the existing models are modified for the Reynolds stress equations to have proper near wall behavior. A dissipation rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The calculations are compared with direct numerical simulations. It is shown that the present models, both the k-epsilon model and the second order closure model, perform well in predicting the behavior of the near wall turbulence. Significant improvements over previous models are obtained.

Shih, T. H.↗
NASA NTRS1997

Numerical Solutions of the Complete Navier-Stokes Equations

This report details the development of a new two-equation turbulence closure model based on the exact turbulent kinetic energy k and the variance of vorticity, zeta. The model, which is applicable to three dimensional flowfields, employs one set of model constants and does not use damping or wall functions, or geometric factors.

Robinson, David F.↗
NASA NTRS2007

Small Scale Response and Modeling of Periodically Forced Turbulence

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied using two-point closures. The frequency response of the turbulent kinetic energy and dissipation rate, and the phase shifts between production, energy and dissipation are determined as functions of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit nontrivial frequency and Reynolds number dependence that reveals a filtering effect of the energy cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on distant interactions between widely separated scales of motion. Finally, the extent to which finite dimensional models (standard two-equation models and various generalizations) can reproduce the observed behavior is discussed.

Bos, Wouter↗