Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

How sheath properties change with gas pressure: modeling and simulation

Abstract Particle-in-cell simulations are used to study how neutral pressure influences plasma properties at the sheath edge. The high rate of ion–neutral collisions at pressures above several mTorr are found to cause a decrease in the ion velocity at the sheath edge (collisional Bohm criterion), a decrease in the edge-to-center density ratio ( h l factor), and an increase in the sheath width and sheath potential drop. A comparison with existing analytic models generally indicates favorable agreement, but with some distinctions. One is that models for the h l factor need to be made consistent with the collisional Bohm criterion. With this and similar corrections, a comprehensive fluid-based model of the plasma boundary transition is constructed that compares well with the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Domain Dynamics under Ultrafast Electric-Field Pulses

Exploring the dynamic responses of a material is of importance to both understanding its fundamental physics at high frequencies and potential device applications. Here we develop a phase-field model for predicting the dynamics of ferroelectric materials and study the dynamic responses of ferroelectric domains and domain walls subjected to an ultrafast electric field pulse. We discover a transition of domain evolution mechanisms from pure domain growth at a relatively low field to combined nucleation and growth of domains at a high field. We derive analytical models for the two regimes which allow us to extract the effective mass and damping coefficient of ferroelectric domain walls. The exhibition of two regimes for the ferroelectric domain dynamics at low and high electric fields is expected to be a general phenomenon that would appear for ferroic domains under other ultrafast stimuli. Finally, the present work also offers a general framework for studying domain dynamics and obtaining fundamental properties of domain walls and thus for manipulating the dynamic functionalities of ferroelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary-/dual- pore network: Part II - numerical model

To increase the market share of electric vehicles, it is desirable to reduce the battery charge times, which are significantly limited by poor electrolyte transport. A high rate charging is achievable by using expensive and low energy density cells with thin electrodes. For higher energy density cells, new electrolytes with improved conductivity and diffusivity and/or electrodes with advanced architecture are required to boost the electrolyte transport, leading to a more uniform utilization of active materials. In our previous work, an analytical model was developed to investigate the effect of secondary pore network (SPN) on electrolyte transport and the configuration of SPN was optimized by enforcing equal characteristic diffusion times in through-plane and in-plane directions. Here, to evaluate the effect of SPN on the fast-charging capability of lithium-ion batteries, a 2D physics-based electrochemical model is developed with SPN in either one or both electrodes. Additionally, the effect of SPN on cell energy density and lithium plating is investigated for cells with different loadings and electrode porosities. Combining SPN with elevated charging temperatures, the model predicts that the volumetric discharge energy density of a 3 mA h/cm2 cell can reach 270 Wh/L after a 6C constant-current charging.

25 ENERGY STORAGE↗
DOE OSTI2023

Fermi-LAT Gamma-Ray Emission Discovered from the Composite Supernova Remnant B0453-685 in the Large Magellanic Cloud

Abstract We report the second extragalactic pulsar wind nebula (PWN) to be detected in the megaelectronvolt–gigaelectronvolt band by the Fermi-LAT, located within the Large Magellanic Cloud. The only other known PWN to emit in the Fermi band outside of the Milky Way is N157B, which lies to the west of the newly detected gamma-ray emission at an angular distance of 4°. Faint, pointlike gamma-ray emission is discovered at the location of the composite supernova remnant (SNR) B0453-685 with a ∼4 σ significance with energies ranging from 300 MeV–2 TeV. We present the Fermi-LAT data analysis of the new gamma-ray source, coupled with a detailed multiwavelength investigation to understand the nature of the observed emission. Combining the observed characteristics of the SNR and the physical implications from broadband modeling, we argue it is unlikely that the SNR is responsible for the gamma-ray emission. While the gamma-ray emission is too faint for a pulsation search, we try to distinguish between any pulsar and PWN component of SNR B0453-685 that could be responsible for the observed gamma-ray emission using semi-analytic models. We determine the most likely scenario is that the old PWN ( τ ∼ 14,000 yr) within B0453-685 has been impacted by the return of the SNR reverse shock with a possible substantial pulsar component below 5 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2022

Modeling interconnections of safety and financial performance of nuclear power plants, part 3: Spatiotemporal probabilistic physics-of-failure analysis and its connection to safety and financial performance

Here, this paper is a byproduct of a line of research by the authors to analyze interrelationships of safety and financial performance of nuclear power plants (NPPs). The result of this line of research is summarized in three parts: Part 1 covers a categorical review of relevant literature and the theoretical bases that support the methodological developments in Part 2. Part 2 introduces an Integrated Enterprise Risk Management (I-ERM) methodological framework to quantify the interconnections of safety and financial performance with a focus on operation and maintenance (O&M) of NPPs. Part 2 has also demonstrated the applicability and values of the I-ERM methodology through an NPP case study. This paper is Part 3, where detailed development and implementation of one of the I-ERM modules, i.e., probabilistic physics-of-failure (PPoF) analysis, and its connection with safety and financial performance is reported. In this article, the physical failure modeling for hardware components is advanced by incorporating finite element analysis (FEA) into PPoF analysis and coupling the FEA-based PPoF with the maintenance performance through a renewal process model. This article covers two scientific contributions: (i) first-of-its-kind incorporation of FEA into the PPoF model of thermal fatigue for NPP components; and (ii) advancing the interface between the PPoF analysis and the renewal process model in order to deal with spatiotemporal FEA outputs and to efficiently estimate the physical transition rates even when the PPoF outputs are dominated by success data. Through the incorporation of FEA, the resolution of the PPoF analysis is enhanced as spatiotemporal conditions such as stress and temperature can be considered explicitly instead of relying on simplified assumptions or analytical models with reduced spatiotemporal dimensions. To demonstrate an application of the FEA-based PPoF analysis and its coupling with maintenance through the renewal process model, a case study is conducted using excess letdown elbow piping in the chemical and volume control system of a Pressurized Water Reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Effects of radiation heat loss on planar and spherical hydrofluorocarbon/air flames

Hydrofluorocarbons (HFC), which are mildly flammable and pose potential fire risks, have received greater attention as a viable low global warming potential alternative to traditional refrigerant and fire-suppressant compounds. Therefore, there is a demand to accurately quantify their flammability and re-activity to establish proper safety metrics. Here, this study investigates the effects of radiation heat loss on slowly-propagating HFC/air laminar flames. Planar 1-D simulations of R-32/air and R-1234yf/air flames show significant reductions in laminar flame speed due to radiative heat losses from the flame zone. Simulations of spherically expanding flames (SEF) revealed that the radiation-induced flow needs to be considered when interpreting data from experiments. To this end, a Spherical-flame RADiation-Induced Flow (SRADIF) model was developed to estimate the burned gas inward flow velocities in constant-pressure SEFs, utilizing the optically thin limit assumption to model radiation heat loss. The model was validated against results from detailed numerical simulations of SEFs, from which radiation-induced inward flow was derived using a new formulation considering both the radiation heat loss and convective flow effects. Results show that SRADIF accurately predicts the inward flow velocity for R-32/air mixtures over a range of conditions and performs significantly better compared to existing analytical models. However, the model was unable to accurately predict flow velocities for R-1234yf/air flames and the reason for this is discussed. (224 words)

42 ENGINEERING↗
DOE OSTI2025

Quantifying and Optimizing the Energy Benefits of Mass Timber Construction

The International Mass Timber Alliance (IMTA) is a global organization of industry leaders, engineers, scientists, and associations dedicated to advancing mass timber construction. Its mission is to generate and disseminate scientific data supporting the development of standardized construction and energy efficient practices that promote the adoption of mass timber worldwide. IMTA collaborated with Oak Ridge National Laboratory (ORNL) to leverage ORNL’s expertise in building envelope modeling and testing to evaluate how mass timber construction can reduce peak heating and cooling demand, lower overall energy use, and improve resilience during power outages. A previous study of 80 mass timber buildings in Finland found measured energy use up to 50% lower than predicted by simulation. This project aimed to validate and extend those findings for U.S. buildings through analytical modeling, laboratory testing, and full-scale building evaluations. The research focused on the thermal performance of low-embodied-energy wall assemblies, such as cross-laminated timber (CLT) panels and log walls, with particular attention to the effects of thermal inertia on indoor comfort and energy performance. While mass timber’s structural and fire-resistance properties are well documented, its whole-building thermal behavior has received limited attention. Field data, simulation results, and resilience testing from this study will inform future modeling practices, design guidelines, and construction practices by quantifying the unique thermal and demand-flexibility benefits of mass timber construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities

The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effects of pillar design on the thermal performance of vacuum-insulated glazing

Vacuum-insulated glazing (VIG) with a low-emittance coating has a great market potential as an effective transparent insulator. The thermal insulating performance of VIG is determined by its design, including material selection and configuration of different components. Thermal conductance of the vacuum gap as a transport bottleneck is one of the primary factors controlling the thermal transport across VIG. In particular, because support pillars provide the main thermal transport channels across the vacuum gap, increasing the pillar thermal resistance is a key strategy for creating effective thermal insulation while maintaining the vacuum space. Here, the effects of various pillar design parameters, such as thermal conductivity, geometry, and arrangement, on the VIG thermal performance were comprehensively investigated via the finite element method. In addition, analytical models for thermal transport were examined and thermal conductance across the VIG unit was experimentally measured for validation. The pillar design parameters, especially the height, shape, spacing, and arrangement of the pillars, showed significant effects on the thermal performance of VIG. This research also shows that the smaller contact area for horizontal pillars can effectively decrease the heat loss by more than 30%. Because current VIG analytical equations of thermal performance are only applicable to cylindrical pillars, an analytical equation that can better describe the thermal performance of rectangular parallelepiped pillars is presented, along with a discussion about the mechanism of thermal transfer for different pillar shapes. Thermo-mechanical analyses based on 3D FEM simulations can provide valuable insights into the effect of various design parameters on the overall performance the VIG, allowing for the development of an optimal VIG design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Recovering non-Maxwellian particle velocity distribution functions from collective Thomson-scattered spectra

Collective optical Thomson scattering (TS) is a diagnostic commonly used to characterize plasma parameters. These parameters are typically extracted by a fitting algorithm that minimizes the difference between a measured scattered spectrum and an analytic spectrum calculated from the velocity distribution function (VDF) of the plasma. However, most existing TS analysis algorithms assume that the VDFs are Maxwellian, and applying an algorithm that makes this assumption does not accurately extract the plasma parameters of a non-Maxwellian plasma due to the effect of non-Maxwellian deviations on the TS spectra. We present new open-source numerical tools for forward modeling analytic spectra from arbitrary VDFs and show that these tools are able to more accurately extract plasma parameters from synthetic TS spectra generated by non-Maxwellian VDFs compared to standard TS algorithms. Estimated posterior probability distributions of fits to synthetic spectra for a variety of example non-Maxwellian VDFs are used to determine uncertainties in the extracted plasma parameters and show that correlations between parameters can significantly affect the accuracy of fits in plasmas with non-Maxwellian VDFs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Additive Manufactured Compact Microwave Absorbers

A high-performance, compact microwave absorber was created using Fused Deposition Modeling (FDM) 3D printing. Both a narrowband and a broadband absorber were created. The narrowband absorber was designed at 4.9 GHz, mid-band in WR-187 waveguide. The broadband absorber tried to achieve the best attenuation across the entire 3.95 to 5.85 GHz band. Two types of carbon loaded polylactic acid (PLA) plastic and one type of unloaded PLA were 3D printed with variable percentages of air to achieve different values of effective dielectric constant and loss tangent. The absorber comprised five or six rectangular pieces of these plastic materials. The thickness and fill factor values for each piece were optimized to minimize reflection through fast analytic modeling in MATLAB®. The results were then verified by HFSS® simulation as well. The stack progressed from the lowest loss and lowest dielectric constant to the highest at the shorting end. The final narrowband load had simulated return loss of 87 dB at 4.9 GHz with an analytic solution in MATLAB. The measured return loss of the 3D printed attenuator was 73 dB at 4.929 GHz. The total length of the absorber was 2.44 inches. A commercial absorber for WR-187 with return loss of 40 dB has length of 13 inches. The experiment proves that an effective and compact microwave absorber can be created using 3D printing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Architecture-Aware Models of AI Engines for High-Performance Matrix Matrix Multiplication

The AI Engine (AIE) architecture, available in systems from mobile SoCs to server-class FPGAs, aims to efficiently execute AI/ML tasks through a two-dimensional array of compute tiles. Previous work on AIEs has explored different approaches to mapping computation across spatial arrays, but the compute kernel running on each tile has not been the focus. Additionally, the AIE-ML architecture introduces memory tiles and omits programmable logic, requiring new approaches to staging and moving data throughout the array. In this work we update analytical models developed for CPUs to produce the design of high performance kernels while introducing new model considerations such as memory structure, throughput, and latency as required by the AIE hardware. We evaluate our models by developing AIE-ML kernels for matrix multiplication in low-precision data types showing performance up to 95% of compute peak for the kernel when data resides in local memory and above 90% of compute peak when data resides in main memory.

Binder, Elliott D. [Carnegie Mellon University, Pi↗