SEARCH · Engineering Papers

Results for “Physical quantities”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Snowmass2021 theory frontier white paper: Astrophysical and cosmological probes of dark matter

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Gastight rotating cylinder electrode: Toward decoupling mass transport and intrinsic kinetics in electrocatalysis

Abstract Decoupling and understanding the various mass, charge, and heat transport phenomena involved in the electrocatalytic transformation of small molecules (i.e., CO 2 , CO, H 2 , N 2 , NH 3 , O 2 , and CH 4 ) is challenging but it can be readily achieved using dimensionless quantities (i.e., Reynolds, Sherwood, Schmidt, Damköhler, Nusselt, Prandtl, and Peclet Numbers) to simplify the characterization of systems with multiple interacting physical phenomena. Herein we report the development of a gastight rotating cylinder electrode cell with well‐defined mass transport characteristics that can be applied to experimentally decouple mass transfer effects from intrinsic kinetics in electrocatalytic systems. The gastight rotating cylinder electrode cell enables the dimensionless analysis of electrocatalytic systems and should enable the rigorous research and development of electrocatalytic technologies.

Jang, Joonbaek↗

DOE OSTI2024

Sustainable production of radionuclidically pure antimony-119

Background Radiopharmaceutical therapy (RPT) uses radionuclides that decay via one of three therapeutically relevant decay modes (alpha, beta, and internal conversion (IC) / Auger electron (AE) emission) to deliver short range, highly damaging radiation inside of diseased cells, maintaining localized dose distribution and sparing healthy cells. Antimony-119 ( 119 Sb, t 1/2 = 38.19 h, EC = 100%) is one such IC/AE emitting radionuclide, previously limited to in silico computational investigation due to barriers in production, chemical separation, and chelation. A theranostic (therapeutic/diagnostic) pair can be formed with 119 Sb’s radioisotopic imaging analogue 117 Sb (t 1/2 = 2.80 h, E γ = 158.6 keV, I γ = 85.9%, β + = 262.4 keV, I β+ = 1.81%). Results Within, we report techniques for sustainable and cost-effective production of pre-clinical quality and quantity, radionuclidically pure 119 Sb and 117 Sb, novel low energy photon measurement techniques for 119 Sb activity determination, and physical yields for various tin target isotopic enrichments and thicknesses using (p, n) and (d, n) nuclear reactions. Additionally, we present a two-column separation providing a radioantimony yield of 73.1% ± 6.9% (N = 3) and tin separation factor of (6.8 ± 5.5) x 10 5 (N = 3). Apparent molar activity measurements for deuteron produced 117 Sb using the chelator TREN-CAM were measured at 42.4 ± 25 MBq 117 Sb/µmol (1.14 ± 0.68 mCi/µmol), and we recovered enriched 119 Sn target material at a recycling efficiency of 80.2% ± 5.5% (N = 6) with losses of 11.6 mg ± 0.8 mg (N = 6) per production. Conclusion We report significant steps in overcoming barriers in 119 Sb production, chemical isolation and purification, enriched target material recycling, and chelation, helping promote accessibility and application of this promising therapeutic radionuclide. We describe a method for 119 Sb activity measurement using its low energy gamma (23.87 keV), negating the need for attenuation correction. Finally, we report the largest yet-measured 119 Sb production yields using proton and deuteron irradiation of natural and enriched targets and radioisotopic purity > 99.8% at end of purification.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2022

Quantum annealing for jet clustering with thrust

Quantum computing holds the promise of substantially speeding up computationally expensive tasks, such as solving optimization problems over a large number of elements. In high-energy collider physics, quantum-assisted algorithms might accelerate the clustering of particles into jets. In this study, we benchmark quantum annealing strategies for jet clustering based on optimizing a quantity called “thrust” in electron-positron collision events. Here, we find that quantum annealing yields similar performance to exact classical approaches and classical heuristics, after tuning the annealing parameters. Without tuning, comparable performance can be obtained through a hybrid quantum/classical approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Biochar alters hydraulic conductivity and impacts nutrient leaching in two agricultural soils

Biochar is purported to provide agricultural benefits when added to the soil, through changes in saturated hydraulic conductivity (K sat ) and increased nutrient retention through chemical or physical means. Despite increased interest and investigation, there remains uncertainty regarding the ability of biochar to deliver these agronomic benefits due to differences in biochar feedstock, production method, production temperature, and soil texture. In this project, a suite of experiments was carried out using biochars of diverse feedstocks and production temperatures, in order to determine the biochar parameters which may optimize agricultural benefits. Sorption experiments were performed with seven distinct biochars to determine sorption efficiencies for ammonium and nitrate. Only one biochar effectively retained nitrate, while all biochars bound ammonium. The three biochars with the highest binding capacities (produced from almond shell at 500 and 800 °C (AS500 and AS800) and softwood at 500 °C (SW500)) were chosen for column experiments. Biochars were amended to a sandy loam and a silt loam at 0 % and 2 % (w/w), and K sat was measured. Biochars reduced K sat in both soils by 64 %–80 %, with the exception of AS800, which increased K sat by 98 % in the silt loam. Breakthrough curves for nitrate and ammonium, as well as leachate nutrient concentration, were also measured in the sandy loam columns. All biochars significantly decreased the quantity of ammonium in the leachate, by 22 % to 78 %, and slowed its movement through the soil profile. Ammonium retention was linked to high cation exchange capacity and a high oxygen-to-carbon ratio, indicating that the primary control of ammonium retention in biochar-amended soils is the chemical affinity between biochar surfaces and ammonium. Biochars had little to no effect on the timing of nitrate release, and only SW500 decreased total quantity, by 27 % to 36 %. The ability of biochar to retain nitrate may be linked to high micropore specific surface area, suggesting a physical entrapment rather than a chemical binding. Together, this work sheds new light on the combined chemical and physical means by which biochar may alter soils to impact nutrient leaching and hydraulic conductivity for agricultural production.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Investigating Influential Parameters for High-Purity Germanium Crystal Growth

This paper focuses on the research and development of high-purity germanium (HPGe) crystals for detector fabrication, specifically targeting applications in rare-event physics searches. The primary objective was to produce large-scale germanium crystals weighing >1 kg with a controlled diameter of ∼10 cm and an impurity range of approximately 1010/cm 3. Ensuring structural integrity and excellent crystalline quality requires a thorough assessment of dislocation density, a critical aspect of the crystal development process. Dislocation density measurements play a crucial role in maximizing the sensitivity of HPGe detectors, and our findings confirmed that the dislocation density fell within acceptable ranges for detector fabrication. Additionally, this paper examines the segregation coefficient of various contaminants during the crystal development process. Comprehensive analysis of impurity segregation is essential for reducing contaminant quantities in the crystal lattice and customizing purification processes. This, in turn, minimizes undesired background noise, enhancing signal-to-noise ratios for rare-event physics searches and overall detector performance. The investigation included the segregation coefficients of three major acceptors and one donor in crystals grown at the University of South Dakota, providing valuable insights for optimizing crystal purity and detector efficiency.

Bhattarai, Sanjay↗

DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗

DOE OSTI2024

The impact of collisionality on the runaway electron avalanche during a tokamak disruption

The exponential growth (avalanching) of runaway electrons (REs) during a tokamak disruption continues to be a large uncertainty in RE modeling. The present work investigates the impact of tokamak geometry on the efficiency of the avalanche mechanism across a broad range of disruption scenarios. It is found that the parameter ν *,crit ⁠, describing the collisionality at the critical energy to run away, delineates how toroidal geometry impacts RE formation. In particular, utilizing a reduced but self-consistent description of plasma power balance, it is shown that for a high-density deuterium-dominated plasma, ν *,crit is robustly less than one, resulting in a substantial decrease in the efficiency of the RE avalanche compared to predictions from slab geometry. In contrast, for plasmas containing a substantial quantity of neon or argon, ν *,crit ≳1, no reduction of the avalanche is observed due to toroidal geometry. This sharp contrast in the impact of low- vs high-Z material results primarily from the relatively strong radiative cooling from high-Z impurities, enabling the plasma to be radiatively pinned at low temperatures and thus large electric fields, even for modest quantities of high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗

DOE OSTI2023

Entropy and Seebeck signals meet on the edges

Here, we explore the electronic entropy per particle s and Seebeck coefficient S in zigzag graphene ribbons. Pristine and edge-doped ribbons are considered using tight-binding models to inspect the role of edge states in the observed thermal transport properties. As a band gap opens when the ribbons are doped at one or both edges, due to asymmetric edge potentials, we find that s and S signals are closely related to each other: both develop sharp dip-peak line shapes as the chemical potential lies in the gap, while the ratio s/S exhibits a near-constant value equal to the elementary charge e at low temperatures. This constant ratio suggests that S can be seen as the transport differential entropy per charge, as suggested by some authors. Our calculations also indicate that measurement of s and S may be useful as a spectroscopic probe of different electronic energy scales involved in such quantities in gapped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Hybrid Approaches for Data Reduction of Spatiotemporal Scientific Applications

Scientists conduct large-scale simulations to compute derived quantities from primary data. Thus, it is crucial that data compression techniques maintain bounded errors on these derived quantities or quantities of interest (QOI). For many spatiotemporal applications, these QOIs are binary in nature and represent presence or absence of a physical phenomenon. In this work, we propose to use a hybrid approah for differential compression for such applications. We use a neural network (NN) approach to determine regions-of-interest (ROIs) where the binary QOIs are going to be prevalent. This is then used with traditional approaches that compress at a lower level (and higher accuracy) for these ROIs as compared to other regions.

Li, Xiao↗

DOE OSTI2025

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Theory Manual for the Fuel Cycle Analysis Code REBUS

REBUS is a system of programs designed for the fuel cycle analysis of fast reactors and accelerator-driven systems. The first version of REBUS was developed in 1970 for calculating the equilibrium operating conditions of liquid metal fast breeder reactors. During the past five decades, various extensions and refinements have been made in the simulation of material motions in fuel cycle facilities and the solution of the neutronics and isotopic transmutation equations that govern the reactor behavior. The latest version complies with the standard code practices and interface dataset specifications of the committee on computer code coordination (CCCC). This report provides a detailed description of the physical processes and their mathematical models and solution methods of the current production version of REBUS. The models of the facilities and physical processes involved in fast reactor fuel cycle are first discussed. Then the mathematical formulation of these physical models is presented, followed by detailed discussions of solution methods of the resulting set of equations. Finally the output edit quantities are discussed along with their definitions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗