Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Sensible-Enthalpy-Based Conditional Moment Closure Model for Homogeneous Charge Compression Ignition with Temperature Inhomogeneity

Abstract This paper presents a sensible-enthalpy-based conditional moment closure (CMC) model for homogeneous charge compression ignition (HCCI) with temperature inhomogeneity. The focus of study is placed on the effects of a conditioning variable on the prediction of heat release rates and developing the mixing models for sensible enthalpy, a reacting scalar. The sensible-enthalpy-based CMC model is validated a priori with 2-D direct numerical simulation (DNS) data for HCCI, and compared with the total-enthalpy-based formulation. Sensible enthalpy is shown to be strongly correlated with temperature, which makes sensible enthalpy well suited for a conditioning variable in CMC for HCCI. Results show that sensible-enthalpy-based CMC performs well for the prediction of heat release rates for all temperature stratification levels investigated in this study, while total-enthalpy-based CMC leads to the overprediction of heat release rates when the temperature stratification level is high. While conditional fluctuations of species mass fractions are substantial, sensible-enthalpy-based CMC is found to accurately estimate heat release rates during ignition, as well as reaction rates of key species, due to substantially suppressed conditional temperature fluctuations. The mapping closure models for the probability density function (PDF) and conditional dissipation rates of sensible enthalpy are proposed. The proposed models are found to capture the characteristics of reactive scalar mixing.

42 ENGINEERING↗
DOE OSTI2018

Groups in Which the Normal Closures of Cyclic Subgroups Have Bounded Finite Hirsch–Zaitsev Rank

In this paper, we study generalized soluble groups with restriction on normal closures of cyclic subgroups. A group G is said to have finite Hirsch–Zaitsev rank if G has an ascending series whose factors are either infinite cyclic or periodic and if the number of infinite cyclic factors is finite. It is not hard to see that the number of infinite cyclic factors in each of such series is an invariant of a group G. This invariant is called the Hirsch–Zaitsev rank of G and will be denoted by r{sub hz}(G). We study the groups in which the normal closure of every cyclic subgroup has the Hirsch–Zaitsev rank at most b (b is some positive integer). For some natural restrictions we find a function k{sub 1}(b) such that r{sub hz}([G/Tor(G),G/Tor(G)]) ≤ k{sub 1}(b).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

von Kármán–Howarth and Corrsin equations closure based on Lagrangian description of the fluid motion

A new approach to obtain the closure formulas for the von Kármán–Howarth and Corrsin equations is presented, which is based on the Lagrangian representation of the fluid motion, and on the Liouville theorem associated to the kinematics of a pair of fluid particles. This kinematics is characterized by the finite scale separation vector which is assumed to be statistically independent from the velocity field. Such assumption is justified by the hypothesis of fully developed turbulence and by the property that this vector varies much more rapidly than the velocity field. This formulation leads to the closure formulas of von Kármán–Howarth and Corrsin equations in terms of longitudinal velocity and temperature correlations following a demonstration completely different with respect to the previous works. Some of the properties and the limitations of the closed equations are discussed. In particular, we show that the times of evolution of the developed kinetic energy and temperature spectra are finite quantities which depend on the initial conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2020

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2020

Selective bond formation triggered by short optical pulses: quantum dynamics of a four-center ring closure

Here we report bond formation induced by an ultrashort UV pulse. The photochemical process is described by quantum dynamics as coherent electronic and nuclear motions during the ultrashort pulse induced ring closure of norbornadiene to quadricyclane. Norbornadiene consists of two ethylene moieties connected by a rigid (CH 2 ) 3 bridge. Upon photoexcitation, two new sigma bonds are formed, resulting in the closure of a four-atom ring. As a medium-sized polyatomic molecule, norbornadiene exhibits a high density of strongly coupled electronic states from about 6 eV above the ground state. We report on inducing the formation of the new bonds using a short femtosecond UV pulse to pump a non-equilibrium electronic density in the open form that evolves towards the closed ring form. As the coherent electronic-nuclear coupled dynamics unfold, the excited states change character through non-adiabatic interactions and become valence states for the two new C–C bonds of quadricyclane. Our three-dimensional fully quantum dynamical grid simulations during the first 200 fs show that short UV pulses of different polarization initiate markedly different initial non-equilibrium electronic densities that follow different dynamical paths to the S 0 /S 1 conical intersection. They lead to different initial relative yields of quadricyclane, thereby opening the way to controlling bond-making with attopulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Closure and transport theory for high-collisionality electron-ion plasmas

Systems of algebraic equations for a high-collisionality electron-ion plasma are constructed from the general moment equations with linearized collision operators [J.-Y. Ji and E. D. Held, Phys. Plasmas 13, 102103 (2006) and J.-Y. Ji and E. D. Held, Phys. Plasmas 15, 102101 (2008)]. A systematic geometric method is invented and applied to solve the system of equations to find closure and transport relations. It is known that some closure coefficients of Braginskii [S. I. Braginskii, Reviews of Plasma Physics (Consultants Bureau, New York, 1965), Vol. 1] are in error up to 65% for some finite values of x (cyclotron frequency Multiplication-Sign electron-ion collision time) and have significant error in the large-x limit [E. M. Epperlein and M. G. Haines, Phys. Fluids 29, 1029 (1986)]. In this work, fitting formulas for electron coefficients are obtained from the 160 moment (Laguerre polynomial) solution, which converges with increasing moments for x{<=}100 and from the asymptotic solution for large x-values. The new fitting formulas are practically exact (less than 1% error) for arbitrary x and Z (the ion charge number, checked up to Z = 100). The ion coefficients for equal electron and ion temperatures are moderately modified by including the ion-electron collision operator. When the ion temperature is higher than the electron temperature, the ion-electron collision and the temperature change terms in the moment equations must be kept. The ion coefficient formulas from 3 moment (Laguerre polynomial) calculations, precise to less than 0.4% error from the convergent values, are explicitly written.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Transitional nuclei near shell closures

High spin states in Bismuth and Thallium nuclei near the Z = 82 shell closure and Cesium nuclei near the N = 82 shell closure in A = 190 and A = 130 regions, respectively, have been experimentally investigated using heavy-ion fusion evaporation reaction and by detecting the gamma rays using the Indian National Gamma Array (INGA). Interesting shape properties in these transitional nuclei have been observed. The results were compared with the neighboring nuclei in these two regions. The total Routhian surface (TRS) calculations have been performed for a better understanding of the observed properties. In mass region A = 190, a change in shape from spherical to deformed has been observd around neutron number N = 112 for the Bi (Z = 83) isotopes with proton number above the magic gap Z = 82, whereas, the shape of Tl (Z = 81) isotopes with proton number below the magic gap Z = 82 remains stable as a function of neutron number. An important transition from aplanar to planar configuration of angular momentum vectors leading to the occurance of nuclar chirality and magnetic rotation, respectively, has been proposed for the unique parity πh{sub 11/2}⊗νh{sub 11/2} configuration in Cs isotopes in the mass region A ∼ 130 around neutron number N = 79. These results are in commensurate with the TRS calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Closure of the single fluid magnetohydrodynamic equations in presence of electron cyclotron current drive

In the presence of electron cyclotron current drive (ECCD), the Ohm's law of single fluid magnetohydrodynamics is modified as E + v × B = η(J – J{sub EC}). This paper presents a new closure relation for the EC driven current density appearing in this modified Ohm's law. The new relation faithfully represents the nonlocal character of the EC driven current and its main origin in the Fisch-Boozer effect. The closure relation is validated on both an analytical solution of an approximated Fokker-Planck equation as well as on full bounce-averaged, quasi-linear Fokker-Planck code simulations of ECCD inside rotating magnetic islands. The new model contains the model put forward by Giruzzi et al. [Nucl. Fusion 39, 107 (1999)] in one of its limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A pressure consistent bridge correction of Kovalenko-Hirata closure in Ornstein-Zernike theory for Lennard-Jones fluids by apparently adjusting sigma parameter

Ornstein-Zernike (OZ) integral equation theory is known to overestimate the excess internal energy, U{sup ex}, pressure through the virial route, P{sub v}, and excess chemical potential, μ{sup ex}, for one-component Lennard-Jones (LJ) fluids under hypernetted chain (HNC) and Kovalenko-Hirata (KH) approximatons. As one of the bridge correction methods to improve the precision of these thermodynamic quantities, it was shown in our previous paper that the method to apparently adjust σ parameter in the LJ potential is effective [T. Miyata and Y. Ebato, J. Molec. Liquids. 217, 75 (2016)]. In our previous paper, we evaluated the actual variation in the σ parameter by using a fitting procedure to molecular dynamics (MD) results. In this article, we propose an alternative method to determine the actual variation in the σ parameter. The proposed method utilizes a condition that the virial and compressibility pressures coincide with each other. This method can correct OZ theory without a fitting procedure to MD results, and possesses characteristics of keeping a form of HNC and/or KH closure. We calculate the radial distribution function, pressure, excess internal energy, and excess chemical potential for one-component LJ fluids to check the performance of our proposed bridge function. We discuss the precision of these thermodynamic quantities by comparing with MD results. In addition, we also calculate a corrected gas-liquid coexistence curve based on a corrected KH-type closure and compare it with MD results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Separating divertor closure effects on divertor detachment and pedestal shape in DIII-D

Comparison between an open divertor and a more-closed divertor in DIII-D demonstrates detachment at an up to 40% lower pedestal density (n e,ped ) in the closed divertor, due to a combination of decreased fueling of the pedestal and increased dissipation in the scrape off layer (SOL) in the closed divertor, both resulting from increased neutral trapping in the divertor. Predicting whether the relationship between divertor closure and detachment will hold for an opaque SOL, in which the contribution of ionizing neutrals to fueling the pedestal is lessened, requires separating out the different mechanisms contributing to the density difference at detachment. A series of experiments on DIII-D characterizes matched discharges using various divertor configurations to isolate the effects of divertor closure. These experiments show detachment at up to 25% lower n e,sep in the closed divertor than the open, supported by simulations showing increased neutral trapping, and hence increased dissipation, in the closed divertor. A difference in n e,ped /n e,sep is also seen: for matched n e,sep , the closed divertor has up to 20% lower ne,ped, consistent with modeling showing a smaller ionization fraction inside the separatrix in this case. Understanding how these pieces fit together will help in develop of predictive models of pedestal density and detached divertors compatible with a high-performance core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On closures for reduced order models - A spectrum of first-principle to machine-learned avenues

For over a century, reduced order models (ROMs) have been a fundamental discipline of theoretical fluid mechanics. Early examples include Galerkin models inspired by the Orr–Sommerfeld stability equation and numerous vortex models, of which the von Karman vortex street is one of the most prominent. Subsequent ROMs typically relied on first principles, like mathematical Galerkin models, weakly nonlinear stability theory, and two- and three-dimensional vortex models pioneered the data-driven proper orthogonal decomposition (POD) modeling. In early POD modeling, available data were used to build an optimal basis, which was then utilized in a classical Galerkin procedure to construct the ROM, but data have made a profound impact on ROMs beyond the Galerkin expansion. In this paper, we take a modest step and illustrate the impact of data-driven modeling on one significant ROM area. Specifically, we focus on ROM closures, which are correction terms that are added to the classical ROMs in order to model the effect of the discarded ROM modes in under-resolved simulations. Through simple examples, we illustrate the main modeling principles used to construct the classical ROMs, motivate and introduce modern ROM closures, and show how data-driven modeling, artificial intelligence, and machine learning have changed the standard ROM methodology over the last two decades. Finally, we outline our vision on how the state-of-the-art data-driven modeling can continue to reshape the field of reduced order modeling.

42 ENGINEERING↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2016

Mitochondrial Hsp90 is a ligand-activated molecular chaperone coupling ATP binding to dimer closure through a coiled-coil intermediate

Heat-shock protein of 90 kDa (Hsp90) is an essential molecular chaperone that adopts different 3D structures associated with distinct nucleotide states: a wide-open, V-shaped dimer in the apo state and a twisted, N-terminally closed dimer with ATP. Although the N domain is known to mediate ATP binding, how Hsp90 senses the bound nucleotide and facilitates dimer closure remains unclear. In this paper we present atomic structures of human mitochondrial Hsp90 N (TRAP1 N ) and a composite model of intact TRAP1 revealing a previously unobserved coiled-coil dimer conformation that may precede dimer closure and is conserved in intact TRAP1 in solution. Our structure suggests that TRAP1 normally exists in an autoinhibited state with the ATP lid bound to the nucleotide-binding pocket. ATP binding displaces the ATP lid that signals the cis-bound ATP status to the neighboring subunit in a highly cooperative manner compatible with the coiled-coil intermediate state. We propose that TRAP1 is a ligand-activated molecular chaperone, which couples ATP binding to dramatic changes in local structure required for protein folding.

59 BASIC BIOLOGICAL SCIENCES↗