Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass conservation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mass and density of individual frozen hydrometeors

A new precipitation sensor, the Differential Emissivity Imaging Disdrometer (DEID), is used to provide the first continuous measurements of the mass, diameter, and density of individual hydrometeors. The DEID consists of an infrared camera pointed at a heated aluminum plate. It exploits the contrasting thermal emissivity of water and metal to determine individual particle mass by assuming that energy is conserved during the transfer of heat from the plate to the particle during evaporation. Particle density is determined from a combination of particle mass and morphology. A Multi-Angle Snowflake Camera (MASC) was deployed alongside the DEID to provide refined imagery of particle size and shape. Broad consistency is found between derived mass–diameter and density–diameter relationships and those obtained in prior studies. However, DEID measurements show a generally weaker dependence with size for hydrometeor density and a stronger dependence for aggregate snowflake mass.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗
DOE OSTI2021

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Performance and Accuracy Implications of Parallel Split Physics-Dynamics Coupling in the Energy Exascale Earth System Atmosphere Model

Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Benchmark of the Chlorine Worth Study Experiments in Support of Chlorine Nuclear Data Validation for Nuclear Criticality Safety

The Chlorine Worth Study (CWS) was a critical experiment to address an urgent need for thermal chlorine nuclear data validation in plutonium systems. This urgent need is tied directly to plutonium recycle and recovery operations in the plutonium facility at Los Alamos National Laboratory, where exceptionally conservative criticality safety limits are used because no credit is taken for the neutron capture by chlorine. The experiment used weapons-grade plutonium metal plates clad in stainless steel, known as the PANN (plutonium aluminum no nickel) ZPPR (zero power physics reactor) plates. The plutonium was reflected and moderated by high-density polyethylene and included combinations of polyvinyl chloride (PVC) and chlorinated polyvinyl chloride (CPVC) as absorbers. The experiment and benchmark included three configurations mimicking 30 g 239 Pu/L plutonium, 300 g 239 Pu/L plutonium, and 600 g 239 Pu/L plutonium in an aqueous chloride solution. Uncertainties in the benchmark included five broad categories: (1) criticality measurement, (2) mass and density, (3) dimensions, (4) material compositions, and (5) positioning. The largest contribution to the overall uncertainties for all three cases came from the material compositions, in particular the PVC and CPVC absorber compositions. A detailed model was created to be a near match (that is within expectations of transport code users) and a simplified model was created to minimize offset dimensions and expedite modeling for code validation. Sample calculations were completed in MCNP6.3 with ENDF/B-VIII.0 and ENDF/B-VII.1 nuclear data. For the detailed and simplified models, the average difference between the computed and experimental k eff was 951 pcm. CWS will serve as the key validation experiment for nuclear criticality safety in support of aqueous chloride operations. The sensitivity to the chlorine capture cross section is orders of magnitude greater than other existing benchmarks. The current limits, as defined by nuclear criticality safety, are 520 g Pu per batch, i.e. the minimum critical mass of the Pu solution infinitely reflected by water [Criticality Handbook: Volume II, (1969)]. This extremely conservative critical mass limit does not credit any neutron capture by chlorine (in particular neutron capture by 35 Cl) and greatly impedes the throughput required for current and future operations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Insight into the autoproteolysis mechanism of the RsgI9 anti‐σ factor from Clostridium thermocellum

Abstract Clostridium thermocellum is a potential microbial platform to convert abundant plant biomass to biofuels and other renewable chemicals. It efficiently degrades lignocellulosic biomass using a surface displayed cellulosome, a megadalton sized multienzyme containing complex. The enzymatic composition and architecture of the cellulosome is controlled by several transmembrane biomass‐sensing RsgI‐type anti‐σ factors. Recent studies suggest that these factors transduce signals from the cell surface via a conserved RsgI extracellular (CRE) domain (also called a periplasmic domain) that undergoes autoproteolysis through an incompletely understood mechanism. Here we report the structure of the autoproteolyzed CRE domain from the C. thermocellum RsgI9 anti‐σ factor, revealing that the cleaved fragments forming this domain associate to form a stable α/β/α sandwich fold. Based on AlphaFold2 modeling, molecular dynamics simulations, and tandem mass spectrometry, we propose that a conserved Asn‐Pro bond in RsgI9 autoproteolyzes via a succinimide intermediate whose formation is promoted by a conserved hydrogen bond network holding the scissile peptide bond in a strained conformation. As other RsgI anti‐σ factors share sequence homology to RsgI9, they likely autoproteolyze through a similar mechanism.

Takayesu, Allen↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Contributions to the BeEST sterile neutrino search

Sterile neutrinos are hypothetical right-handed leptons that have been postulated as extensions to the Standard Model to help explain the origin of neutrino masses and the matter-antimatter asymmetry in the Universe. Sterile neutrinos with keV-scale masses are also a promising candidate for dark matter. LLNL and the Colorado School of Mines are developing an experiment to search for sterile neutrinos by measuring the \Beryllium-7 Electron capture in Superconducting Tunnel junctions (STJs)" (the \BeEST" Experiment). When 7 Be implanted into STJ radiation detectors decays by electron capture, the neutrino escapes and the energy of the recoiling 7 Li nucleus can be measured with high accuracy. Using conservation of energy and momentum, the mass of the emitted neutrino can then be found. Sterile keV neutrinos would generate additional peaks in the spectrum, and the experiment uses high-resolution STJs to observe such a potential signal. During my internship, I contributed to different aspects of the BeEST experiment. I designed a UHV setup for X-ray spectroscopy analysis with STJs. It will allow high-resolution X-ray spectroscopy of different materials to examine the broadening in the BeEST spectra, the origins of which are unknown but may be related to material effects inside the STJ detectors. The design accounts for the specific demands of the BeEST project, but provides a degree of flexibility that makes it useful for other potential STJ projects. Additionally, to test the accuracy of the electronics and data analysis routine, an arbitrary waveform generator (AWG) will be used to generate simulated spectra that contain sterile neutrino signals of a specified mass and admixture. For this I programmed an AWG to output random STJ signals whose amplitudes follow the appropriate user-de ned probability distribution. The AWG signals can be input into the 32-channel preamplifier to simulate any spectrum of interest and test the uniformity of the amplifier chain, as well as the sensitivity of the data analysis.

42 ENGINEERING↗
DOE OSTI2023

Local adaptation of switchgrass drives trait relations to yield and differential responses to climate and soil environments

Abstract Switchgrass, a potential biofuel crop, is a genetically diverse species with phenotypic plasticity enabling it to grow in a range of environments. Two primary divergent ecotypes, uplands and lowlands, exhibit trait combinations representative of acquisitive and conservative growth allocation strategies, respectively. Whether these ecotypes respond differently to various types of environmental drivers remains unclear but is crucial to understanding how switchgrass varieties will respond to climate change. We grew two upland, two lowland, and two intermediate/hybrid cultivars of switchgrass at three sites along a latitudinal gradient in the central United States. Over a 4‐year period, we measured plant functional traits and biomass yields and evaluated genotype‐by‐environment (G × E) interaction effects by analyzing switchgrass responses to soil and climate variables. We found substantial evidence of G × E interactions on biomass yield, primarily due to deviations in the response of the southern lowland cultivar Alamo, which produced more biomass in hotter and drier environments relative to other cultivars. While lowland cultivars had the highest potential for yield, their yields were more variable year‐to‐year compared to other cultivars, suggesting greater sensitivity to environmental perturbations. Models comparing soil and climate principal components as explanatory variables revealed soil properties, especially nutrients, to be most effective at predicting switchgrass biomass yield. Also, positive correlations between biomass yield and conservative plant traits, such as high stem mass and tiller height, became stronger at lower latitudes where the climate is hotter and drier, regardless of ecotype. Lowland cultivars, however, showed a greater predisposition to exhibit these conservative traits. These results suggest switchgrass trait allocation trade‐offs that prioritize aboveground biomass production are more tightly associated in hot, dry environments and that lowland cultivars may exhibit a more specialized strategy relative to other cultivars. Altogether, this research provides essential knowledge for improving the viability of switchgrass as a biofuel crop.

09 BIOMASS FUELS↗
DOE OSTI2022

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) Version 4.1 Technical Reference Manual

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 4.1. The two most significant changes for HyRAM+ version 4.1 from HyRAM+ version 4.0 are direct incorporation of unconfined overpressure into the QRA calculations and modification of the models for cryogenic liquid flow through an orifice. In QRA mode, the user no longer needs to input peak overpressure and impulse values that were calculated separately; rather, the unconfined overpressure is estimated for the given system inputs, leak size, and occupant location. The orifice flow model now solves for the maximum mass flux through the orifice at constant entropy while conserving energy, which does not require a direct speed of sound calculation. This does not affect the mass flow for all-gaseous releases; the method results in the same speed of sound for choked flow. However, this method does result in a higher (and more realistic) mass flow rate for a given leak size for liquid releases than was previously calculated.

08 HYDROGEN↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

MOOSE Navier–Stokes module

The MOOSE Navier–Stokes module solves mass, momentum, energy, and passive scalar conservation equations in the context of fluid flow. The module supports solution of these equations in both free flow and porous medium contexts and for a range of fluid compressibility. The conservation equations can be discretized in space using continuous Galerkin finite elements or with cell centered finite volumes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

TDCOSMO 2025: Cosmological constraints from strong lensing time delays

We present cosmological constraints from eight strongly lensed quasars (hereafter, the TDCOSMO-2025 sample). Building on previous work, our analysis incorporated new deflector stellar velocity dispersions measured from spectra obtained with the James Webb Space Telescope (JWST), the Keck Telescopes, and the Very Large Telescope (VLT), utilizing improved methods. We used integrated JWST stellar kinematics for five lenses, VLT-MUSE for 2, and resolved kinematics from Keck and JWST for RX J1131−1231. We also considered two samples of non-time-delay lenses: 11 from the Sloan Lens ACS (SLACS) sample with Keck-KCWI resolved kinematics; and four from the Strong Lenses in the Legacy Survey (SL2S) sample. We improved our analysis of line-of-sight effects, the surface brightness profile of the lens galaxies, and orbital anisotropy, and corrected for projection effects in the dynamics. Our uncertainties are maximally conservative by accounting for the mass-sheet degeneracy in the deflectors’ mass density profiles. The analysis was blinded to prevent experimenter bias. Our primary result is based on the TDCOSMO-2025 sample, in combination with Ωm constraints from the Pantheon+ Type Ia supernovae (SN) dataset. In the flat Λ cold dark matter (CDM), we find H0 = 71.6+3.9−3.3 km s−1 Mpc−1. The SLACS and SL2S samples are in excellent agreement with the TDCOSMO-2025 sample, improving the precision on H0 in flat ΛCDM to 4.6%. Using the Dark Energy Survey SN Year-5 dataset (DES-SN5YR) or DESI-DR2 baryonic acoustic oscillations (BAO) likelihoods instead of Pantheon+ yields very similar results. We also present constraints in the open ΛCDM, wCDM, w0waCDM, and wϕCDM cosmologies. The TDCOSMO H0 inference is robust and consistent across all presented cosmological models, and our cosmological constraints in them agree with those from the BAO and SN.Key words: cosmological parameters / cosmology: observations / dark energy / distance scale⋆⋆ Brinson Fellow.⋆⋆⋆ NHFP Einstein Fellow.

Birrer, Simon [SUNY, Stony Brook] (ORCID:000000033↗
DOE OSTI2026

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗
DOE OSTI2024

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗
DOE OSTI2020

Measurement of the mass difference and the binding energy of the hypertriton and antihypertriton

According to the CPT theorem, which states that the combined operation of charge conjugation, parity transformation and time reversal must be conserved, particles and their antiparticles should have the same mass and lifetime but opposite charge and magnetic moment. In this work, we test CPT symmetry in a nucleus containing a strange quark, more specifically in the hypertriton. This hypernucleus is the lightest one yet discovered and consists of a proton, a neutron and a Λ hyperon. With data recorded by the STAR detector at the Relativistic Heavy Ion Collider, we measure the Λ hyperon binding energy B Λ for the hypertriton, and find that it differs from the widely used value and from predictions, where the hypertriton is treated as a weakly bound system. Our results place stringent constraints on the hyperon–nucleon interaction and have implications for understanding neutron star interiors, where strange matter may be present. A precise comparison of the masses of the hypertriton and the antihypertriton allows us to test CPT symmetry in a nucleus with strangeness, and we observe no deviation from the expected exact symmetry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗