Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “properties closure”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗
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 NTRS1979

A kinetic approach to compressible turbulence

A kinetic theory of strong turbulence is developed by decomposing fluctuation of distribution function into scaled components, or ranks, to represent three transport processes: macroscopic evolution, transport property and relaxation. Closure is obtained by formulating a relaxation process for the approach of eddy diffusivity to equilibrium. The kinetic equation is solved and three transport functions of turbulence in adiabatic gas are derived: two transfer functions describe turbulent mode transfer, and a coupling function describes the build-up of kinetic energy at the expense of internal energy. The energy balance in the inertial subrange of homogeneous and isotropic turbulence in an adiabatic gas yields a k to the -3rd law, in agreement with numerical computations from one dimensional model of compressible turbulence.

Tchen, C. M.↗
NASA NTRS2000

North Atlantic Aerosol Radiative Impacts Based on Satellite Measurements and Aerosol Intensive Properties from TARFOX and ACE-2

We estimate the impact of North Atlantic aerosols on the net short-wave flux at the tropopause by combining satellite-derived aerosol optical depth (AOD) maps with model aerosol properties determined via closure analyses in TARFOX and ACE 2. We exclude African dust, primarily by restricting latitudes to 25-60 N. The analyses use in situ aerosol composition measurements and air- and ship-borne sun-photometer measurements of AOD spectra. The aerosol model yields computed flux sensitivities (dFlux/dAOD) that agree with measurements by airborne flux radiometers in TARFOX. Its midvisible single-scattering albedo is 0.9. which is in the range obtained from in situ measurements of scattering and absorption in both TARFOX and ACE 2. Combining satellite-derived AOD maps with the aerosol model yields maps of 24-hour average net radiative flux changes. For simultaneous AVHRR, radiance measurements exceeded the sunphotometer AODs by about 0.04. However. shipboard sunphotometer and AVHRR AODs agreed Within 0.02 for data acquired during satellite overflights on two other days. We discuss attempts to demonstrate column closure within the MBL by comparing shipboard sunphotometer AODs and values calculated from simultaneous shipboard in-situ aerosol size distribution measurements. These comparisons were mostly unsuccessful, but they illustrate the difficulties inherent in this type of closure analysis. Specifically, AODs derived from near-surface in-situ size distribution measurements are extremely sensitive to the assumed hygroscopic growth model that itself requires an assumption of particle composition as a function of height and size, to the radiosonde-measured relative humidity, and to the vertical profile of particle number. We investigate further the effects of hygroscopic particle growth within the MBL by using shipboard lidar aerosol backscatter profiles together with the sunphotometer AOD.

Russell, P. B.↗
DOE OSTI2016

A nonlinear structural subgrid-scale closure for compressible MHD. I. Derivation and energy dissipation properties

Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of nonlinearity involved. However, the dynamical ranges of these phenomena are much larger than what is computationally accessible. In large eddy simulations (LESs), the resulting limited resolution effects are addressed explicitly by introducing to the equations of motion additional terms associated with the unresolved, subgrid-scale dynamics. This renders the system unclosed. We derive a set of nonlinear structural closures for the ideal MHD LES equations with particular emphasis on the effects of compressibility. The closures are based on a gradient expansion of the finite-resolution operator [W. K. Yeo (CUP, 1993)] and require no assumptions about the nature of the flow or magnetic field. Thus, the scope of their applicability ranges from the sub- to the hyper-sonic and -Alfvénic regimes. The closures support spectral energy cascades both up and down-scale, as well as direct transfer between kinetic and magnetic resolved and unresolved energy budgets. They implicitly take into account the local geometry, and in particular, the anisotropy of the flow. Their properties are a priori validated in Paper II [P. Grete et al., Phys. Plasmas 23, 062317 (2016)] against alternative closures available in the literature with respect to a wide range of simulation data of homogeneous and isotropic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, C. G.↗
NASA NTRS1987

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler k-epsilon and k-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, Charles G.↗
DOE OSTI2025

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1989

The generalized drift flux approach: Identification of the void-drift closure law

The main characteristics and the potential advantages of generalized drift flux models are presented. In particular it is stressed that the issue on the propagation properties and on the mathematical nature (hyperbolic or not) of the model and the problem of closure are easier to tackle than in two fluid models. The problem of identifying the differential void-drift closure law inherent to generalized drift flux models is then addressed. Such a void-drift closure, based on wave properties, is proposed for bubbly flows. It involves a drift relaxation time which is of the order of 0.25 s. It is observed that, although wave properties provide essential closure validity tests, they do not represent an easily usable source of quantitative information on the closure laws.

Boure, J. A.↗
DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Hydrodynamic trait coordination and cost–benefit trade‐offs throughout the isohydric–anisohydric continuum in trees

Abstract Whole‐plant hydraulic performance depends on the integrated function of trait complexes, such as embolism resistance, stomatal closure mechanisms, and root properties. The diversity of such traits produces a wide range of response strategies to both short‐term variation of environmental conditions and long‐term changes to climate and hydrological cycles. This study aims to assess the role of different emergent hydraulic trait combinations in trees' vulnerability to drought using a quantitative modelling framework. This modelling framework may be helpful in studying the influence of plant hydraulic traits independently. It will also be used to assess how different groups of traits interact to form viable hydraulic strategies in response to different environmental conditions. We use the advanced plant hydrodynamic model, finite‐difference ecosystem‐scale tree crown hydrodynamics model version 2, to resolve plant functional traits of roots, stems, and leaves and simulate the transpiration. We define a multidimensional hydraulic “strategy space” by considering a continuum of these traits to test the outcomes of different hypothetical strategies in response to environmental conditions as observed in Northern Michigan, USA. We evaluate the degree to which simulated trees suffer hydraulic failure due to cavitation or carbon deficiency in response to reduced stomatal conductance. Our results demonstrate how the relationship between plant hydraulic strategy and hydraulic safety margin emerges from trait combinations along different tissue levels. Our findings suggest that hydrodynamic models present an exciting new possibility to define and study plant traits and hydrodynamics.

Mirfenderesgi, Golnazalsadat↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection, Part III: Results and parameter sensitivity

Two-fluid modelling has recently emerged as a promising approach to representing cumulus convection in weather and climate models. This study applies the two-fluid model described in Part II Thuburn et al., 2021b to a shallow cumulus convection case study over land (ARM). Large eddy simulation data is used to tune the majority of the closures that determine the properties of entrained and detrained air. The two-fluid model is generally able to reproduce the profiles of the mean and turbulent quantities over all stages of the diurnal cycle. As such, the initiation of shallow convection and the evolution of the cloud layer are well captured. The robustness of the two-fluid model is further verified using a steady-state test case (BOMEX), in which the cloud properties are also well modelled.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Structural and Electronic Transitions in Liquid FeO Under High Pressure

Abstract FeO represents an important end‐member for planetary interiors mineralogy. However, its properties in the liquid state under high pressure are poorly constrained. Here, in situ high‐pressure and high‐temperature X‐ray diffraction experiments, ab initio simulations, and thermodynamic calculations are combined to study the local structure and density evolution of liquid FeO under extreme conditions. Our results highlight a strong shortening of the Fe‐Fe distance, particularly pronounced between ambient pressure and ∼40 GPa, possibly related with the insulator to metal transition occurring in solid FeO over a similar pressure range. Liquid density is smoothly evolving between 60 and 150 GPa from values calculated for magnetic liquid to those calculated for non‐magnetic liquid, compatibly with a continuous spin crossover in liquid FeO. The present findings support the potential decorrelation between insulator/metal transition and the high‐spin to low‐spin continuous transition, and relate the changes in the microscopic structure with macroscopic properties, such as the closure of the Fe‐FeO miscibility gap. Finally, these results are used to construct a parameterized thermal equation of state for liquid FeO providing densities up to pressure and temperature conditions expected at the Earth's core‐mantle boundary.

Morard, Guillaume↗
DOE OSTI2020

2020 LANL Contributions to the BATS Test in WIPP

The DOE Office of Nuclear Energy (DOE-NE) repository research and development program seeks to provide a sound technical basis for multiple viable disposal options, increase confidence in the robustness of generic disposal concepts, and develop science and engineering tools needed to support disposal concept implementation. Sandia, Los Alamos, and Lawrence Berkeley National labs are conducting research into salt, which includes the Brine Availability Test in Salt (BATS) field test at the Waste Isolation Pilot Plant (WIPP), a DOE Office of Environmental Management facility. BATS is leveraging the existing infrastructure associated with the WIPP to advance science informing generic disposal concepts. BATS is an ongoing heated borehole experiment being conducted in the WIPP underground. The goal of this experiment is to reduce the uncertainty associated with spent nuclear fuel disposition in geologic salt formations. The BATS experiments is designed to accomplish this by increasing our understanding of brine migration in salt, assessing damaged zones from mining and drilling, simulating a post-closure environment, confirming salt properties, and providing data for model validation. Phased testing for generic salt repository R&D was proposed in 2015. Subsequently, borehole thermal testing was developed as a first step in this process. This document summarizes the work conducted on the BATS experiment during fiscal year 2020 with an emphasis on the contributions by Los Alamos National Lab.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1990

The pdf approach to turbulent flow

This paper provides a detailed discussion of the theory and application of probability density function (pdf) methods, which provide a complete statistical description of turbulent flow fields at a single point or a finite number of points. The basic laws governing the flow of Newtonian fluids are set up in the Eulerian and the Lagrangian frame, and the exact and linear equations for the characteristic functionals in those frames are discussed. Pdf equations in both frames are derived as Fourier transforms of the equations of the characteristic functions. Possible formulations for the nonclosed terms in the pdf equation are discussed, their properties are assessed, and closure modes for the molecular-transport and the fluctuating pressure-gradient terms are reviewed. The application of pdf methods to turbulent combustion flows, supersonic flows, and the interaction of turbulence with shock waves is discussed.

Kollmann, W.↗
DOE OSTI2015

Communication: Fragment-based Hamiltonian model of electronic charge-excitation gaps and gap closure

Capturing key electronic properties such as charge excitation gaps within models at or above the atomic scale presents an ongoing challenge to understanding molecular, nanoscale, and condensed phase systems. One strategy is to describe the system in terms of properties of interacting material fragments, but it is unclear how to accomplish this for charge-excitation and charge-transfer phenomena. Hamiltonian models such as the Hubbard model provide formal frameworks for analyzing gap properties but are couched purely in terms of states of electrons, rather than the states of the fragments at the scale of interest. The recently introduced Fragment Hamiltonian (FH) model uses fragments in different charge states as its building blocks, enabling a uniform, quantum-mechanical treatment that captures the charge-excitation gap. These gaps are preserved in terms of inter-fragment charge-transfer hopping integrals T and on-fragment parameters U{sup (FH)}. The FH model generalizes the standard Hubbard model (a single intra-band hopping integral t and on-site repulsion U) from quantum states for electrons to quantum states for fragments. We demonstrate that even for simple two-fragment and multi-fragment systems, gap closure is enabled once T exceeds the threshold set by U{sup (FH)}, thus providing new insight into the nature of metal-insulator transitions. This result is in contrast to the standard Hubbard model for 1d rings, for which Lieb and Wu proved that gap closure was impossible, regardless of the choices for t and U.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Integral equation theory of thermodynamics, pair structure, and growing static length scale in metastable hard sphere and Weeks-Chandler-Andersen fluids

Here, we employ the Ornstein-Zernike integral equation theory with the Percus-Yevick (PY) and modified-Verlet (MV) closures to study the equilibrium structural and thermodynamic properties of metastable monodisperse hard sphere and continuous repulsion Weeks-Chandler-Andersen (WCA) fluids under density and temperature conditions where the system is strongly overcompressed or supercooled, respectively. The theoretical results are compared to crystal-avoiding simulations of these dense monodisperse model one-component fluids. The equation of state (EOS) and dimensionless compressibility are computed using both the virial and compressibility routes. For hard spheres, the MV-based virial route EOS and dimensionless compressibility are in very good agreement with simulation for all packing fractions, much better than the PY analogs. The corresponding MV-based predictions for the static structure factor are also very good. The amplitude of density fluctuations on the local cage scale and in the long wavelength limit, and three technically different measures of the density correlation length, are studied with both closures. All five properties grow in a roughly exponential manner with density in the metastable regime up to packing fractions of 58% with no sign of saturation. The MV-based results are in good agreement with our crystal-avoiding simulations. Interestingly, the density dependences of long and short wavelength quantities are closely related. The MV-based theory is also quite accurate for the thermodynamics and structure of supercooled monodisperse WCA fluids. Overall our findings are also relevant as critical input to microscopic theories that relate the equilibrium pair correlation function or static structure factor to dynamical constraints, barriers, and activated relaxation in glass-forming liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Stochastic superparameterization in quasigeostrophic turbulence

In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent cascades of energy from the unresolved to the resolved scales resulting in complex patterns of waves, jets, and vortices. Conventional superparameterization simulates large scale dynamics on a coarse grid in a physical domain, and couples these dynamics to high-resolution simulations on periodic domains embedded in the coarse grid. Stochastic superparameterization replaces the nonlinear, deterministic eddy equations on periodic embedded domains by quasilinear stochastic approximations on formally infinite embedded domains. The result is a seamless algorithm which never uses a small scale grid and is far cheaper than conventional SP, but with significant success in difficult test problems. Various design choices in the algorithm are investigated in detail here, including decoupling the timescale of evolution on the embedded domains from the length of the time step used on the coarse grid, and sensitivity to certain assumed properties of the eddies (e.g. the shape of the assumed eddy energy spectrum). We present four closures based on stochastic superparameterization which elucidate the properties of the underlying framework: a ‘null hypothesis’ stochastic closure that uncouples the eddies from the mean, a stochastic closure with nonlinearly coupled eddies and mean, a nonlinear deterministic closure, and a stochastic closure based on energy conservation. The different algorithms are compared and contrasted on a stringent test suite for quasigeostrophic turbulence involving two-layer dynamics on a β-plane forced by an imposed background shear. The success of the algorithms developed here suggests that they may be fruitfully applied to more realistic situations. They are expected to be particularly useful in providing accurate and efficient stochastic parameterizations for use in ensemble-based state estimation and prediction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗