Search NASA⌕ Search

SEARCH · Search NASA

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 19 records

An efficient algorithm for generating diverse microstructure sets and delineating properties closures

Properties closures delineate the theoretical objective space for materials design problems, allowing designers to make informed trade-offs between competing constraints and target properties. In this paper, we present a new algorithm called hierarchical simplex sampling (HSS) that approximates properties closures more efficiently and faithfully than traditional optimization based approaches. By construction, HSS generates samples of microstructure statistics that span the corresponding microstructure hull. As a result, we also find that HSS can be coupled with synthetic polycrystal generation software to generate diverse sets of microstructures for subsequent mesoscale simulations. Finally, by more broadly sampling the space of possible microstructures, it is anticipated that such diverse microstructure sets will expand our understanding of the influence of microstructure on macroscale effective properties and inform the construction of higher-fidelity mesoscale structure-property models.

36 MATERIALS SCIENCE↗

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗

Accelerate Nuclear Research and Development by Reducing Time and Cost Spend in the Pre-conceptual Design Phase of Advanced Reactor Experiments

The design process of every new concept, such as advanced nuclear reactors or associated experiments, starts with the pre-conceptual design phase. In this phase, the viability of a wide range of design options needs to be assessed quickly, to understand the operating envelope and its feasibility. A variety of physics models (thermal-hydraulics, neutronics, mechanical design, etc.) has to be considered at this very first design stage and optimum component sizes and materials (e.g. heat exchangers, piping, turbomachinery, coolant type, etc.) have to be chosen for a given set of boundary conditions (e.g. heat source, heat sink, flow rate, etc.). Detailed solutions such as provided by high fidelity methods like computational fluid dynamics (CFD), Monte Carlo methods, etc. and even lower fidelity tools such as system or subchannel codes, etc. are usually not used during the pre-conceptual design due to the relatively long time needed to create input models, the computational time to obtain a solution and the lack of flexibility to quickly investigate different combinations of components, individual component sizes and material properties. High fidelity tools are usually only employed in the conceptual design and later phases once a base concept has been identified during the pre-conceptual design stage. The current practice during the pre-conceptual design stage is that analysts collect the needed equations, material properties, closure laws, etc. and create ad-hoc solutions form scratch for every new problem. There clearly is a lack of a flexible scoping tool that can be used during pre-conceptional design before higher fidelity tools (as described above) come into play. To reduce user errors in ad-hoc solutions and increase fidelity and efficiency, this project aims to investigate and develop a user-friendly scoping tool to address the thermal-hydraulic designing needs during preconceptual experiment design, i.e. Thermal-hydraulic Research Universal Scoping Tool (TRUST). The success of TRUST will provide the nuclear engineers with an easy-to-use and affordable calculator for early reactor system design and optimization.

42 ENGINEERING↗

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comprehensive characterization of the irradiation effects of glassy carbon

Carbon materials have become increasingly diverse, finding applications in high-temperature and high-radiation environments. Glassy carbon, an allotrope known for its exceptional chemical inertness and desirable mechanical properties. However, understanding neutron irradiation effects in glassy carbon has proven challenging, primarily because of its unique nanopore structure. Here, this study presents a highly detailed microstructural characterization investigation of neutron-induced changes in glassy carbon, revealing how changes in nanopore structure and crystallinity impact the irradiation-induced shrinkage. Aberration-corrected scanning transmission electron microscopy (STEM) reveals pore closure that leads to material densification in the irradiated samples. Dimensional analysis combined with comparison to historical data suggest significant length shrinkage to occur. Neutron and in situ electron irradiation experiments suggest that glassy carbon transforms into so-called carbon onions, supporting the concept of shrinkage saturation. Investigating irradiation temperature effects using STEM, electron energy loss spectroscopy, x-ray diffraction, and Raman spectroscopy revealed partial amorphization at 210 °C–230 °C and preserved order in glassy carbon at 860 °C, coinciding with pore closure. Thermal property measurements were also conducted to assess the effects of densification and other changes in the atomic structure of glassy carbon. The results of this study have broad implications in the deployment of glassy carbon to nuclear environments, based around the observed changes in the thermal properties, and demonstrates the operational window for the onset of densification.

36 MATERIALS SCIENCE↗

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗

A radiation closure study of Arctic stratus cloud microphysical properties using the collocated satellite-surface data and Fu-Liou radiative transfer model

Retrievals of cloud microphysical properties based on passive satellite imagery are especially difficult over snow-covered surfaces because of the bright and cold surface. To help quantify their uncertainties, single-layered overcast liquid-phase Arctic stratus cloud microphysical properties retrieved by using the Clouds and the Earth's Radiant Energy System Edition 2 and Edition 4 (CERES Ed2 and Ed4) algorithms are compared with ground-based retrievals at the Atmospheric Radiation Measurement North Slope of Alaska (ARM NSA) site at Barrow, AK, during the period from March 2000 to December 2006. A total of 206 and 140 snow-free cases (Rsfc ≤ 0.3), and 108 and 106 snow cases (R sfc > 0.3), respectively, were selected from Terra and Aqua satellite passes over the ARM NSA site. The CERES Ed4 and Ed2 optical depth (τ) and liquid water path (LWP) retrievals from both Terra and Aqua are almost identical and have excellent agreement with ARM retrievals under snow-free and snow conditions. In order to reach a radiation closure study for both the surface and top of atmosphere (TOA) radiation budgets, the ARM precision spectral pyranometer-measured surface albedos were adjusted (63.6% and 80% of the ARM surface albedos for snow-free and snow cases, respectively) to account for the water and land components of the domain of 30 km × 30 km. Most of the radiative transfer model calculated SW↓ sfc and SW↑ TOA fluxes by using ARM and CERES cloud retrievals and the domain mean albedos as input agree with the ARM and CERES flux observations within 10 W m -2 for both snow-free and snow conditions. Sensitivity studies show that the ARM LWP and r e retrievals are less dependent on solar zenith angle (SZA), but all retrieved optical depths increase with SZA.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Aqueous stream characterization from biomass fast pyrolysis and catalytic fast pyrolysis

Here, biomass pyrolysis offers a promising means to rapidly depolymerize lignocellulosic biomass for subsequent catalytic upgrading to renewable fuels. Substantial efforts are currently ongoing to optimize pyrolysis processes including various fast pyrolysis and catalytic fast pyrolysis schemes. In all cases, complex aqueous streams are generated containing solubilized organic compounds that are not converted to target fuels or chemicals and are often slated for wastewater treatment, in turn creating an economic burden on the biorefinery. Valorization of the species in these aqueous streams, however, offers significant potential for substantially improving the economics and sustainability of thermochemical biorefineries. To that end, here we provide a thorough characterization of the aqueous streams from four pilot-scale pyrolysis processes: namely, from fast pyrolysis, fast pyrolysis with downstream fractionation, in situ catalytic fast pyrolysis, and ex situ catalytic fast pyrolysis. These configurations and processes represent characteristic pyrolysis processes undergoing intense development currently. Using a comprehensive suite of aqueous-compatible analytical techniques, we quantitatively characterize between 12 g kg -1 of organic carbon of a highly aqueous catalytic fast pyrolysis stream and up to 315 g kg -1 of organic carbon present in the fast pyrolysis aqueous streams. In all cases, the analysis ranges between 75 and 100% of mass closure. The composition and stream properties closely match the nature of pyrolysis processes, with high contents of carbohydrate-derived compounds in the fast pyrolysis aqueous phase, high acid content in nearly all streams, and mostly recalcitrant phenolics in the heavily deoxygenated ex situ catalytic fast pyrolysis stream. Overall, this work provides a detailed compositional analysis of aqueous streams from leading thermochemical processes -- analyses that are critical for subsequent development of selective valorization strategies for these waste streams.

09 BIOMASS FUELS↗

Implications of Aerosol Physicochemical Properties Including Ice Nucleation at ARM Mega Sites for Improved Understanding of Microphysical Atmospheric Cloud Processes (Final Technical Report)

Continuous, long-term measurements of atmospheric ice-nucleating particles (INPs) that influence clouds and precipitation were conducted to investigate the abundance and variability of ground-level INPs across the world. Three field campaigns were organized by the DOE Atmospheric Radiation Measurement (ARM) program, including Examining INP from Southern Great Plains (ExINP-SGP, 2019), Examining INP from Eastern North Atlantic (ExINP-ENA, 2020 – 2021), and Examining INP from North Slope of Alaska (ExINP-NSA, 2021 – 2024). Additional small-scale supporting field experiments were performed in 2019 and 2021 to collect airborne particulate matter at SGP for complementary laboratory characterization of the particles’ physical and chemical properties [Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY), 2019; ExINP-SGP II, 2021]. In these studies, the PI’s team measured INP concentration with both real-time and laboratory measurements in a wide range of freezing temperatures ($T$ from 0 °C to about –30 °C). This project elucidated spatial variability and seasonality in the abundance of immersion mode active INPs across three ARM sites using a single instrument, a Portable Ice Nucleation Experiment (PINE) chamber version 03 (PINE-03 hereafter). Collocated aerosol and meteorological data were analyzed to assess the correlation between ambient INP abundance, air mass origin region, and meteorological variability. Our findings suggest very high freezing efficiency of INPs at the NSA site across the measured temperatures (ice nucleation active surface site density, $n_s(T)$, $\approx 2 \times 10^{8} - 10^{10}$ m -2 for from –16 to –31 °C), which is a factor of 10 – 1000 times greater efficiency as compared to that found in the previous mid-latitude INP measurements in autumn using the same instrument; surprisingly high INP abundance ($\ge 1 \text{ L}^{-1}$ at –25 °C) for the examined temperatures throughout the year that PINE-03 did not measure at other sites; and high INP concentration in spring, possibly related to arctic haze episodes.

58 GEOSCIENCES↗

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗