Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
DOE OSTI2024

Towards a quantitative understanding of bonding in supersonic single particle impacts: A three-dimensional FIB-SEM exploration

Particle bonding is crucial to coating quality in cold spray, but it has been a challenge to accurately quantify bonding even in single particle impacts. This paper uses FIB-SEM to explicitly map the particle-substrate interface for Cu-on-Cu single microparticle impacts in a full 3D rendering that spans a wide range of impact velocities. This approach permits a detailed quantification of the total bonding area and all of its associated components. In addition to revealing why prior 2D characterization efforts have missed important details about impact bonding, these data quantitatively reveal the evolution of bonding from its onset at the “critical velocity” V cr (where bonding is generally poor, ∼6 %) to its peak at around 1.3‧V cr (where almost 90 % of the particle bonds). Further increase in the velocity to 1.5‧V cr and beyond finds the onset of hydrodynamic penetration and a decrease in bonding. These data then support the development of a simple analytical model based on oxide rarefication and extrusion of bare metal through gaps in the oxide layer as driving the development of bonding. As a result, the model reproduces the experiments and provides guidance on optimization of bonding as a function of material and process parameters.

Analytical methods↗
DOE OSTI2025

Observation of kinetic mix enhancement in thin-shell OMEGA implosions

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D 3 He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fall line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D 3 He-p mix yield of 1.7 × 10 9 as compared to the experimental value of 9.3⁢ (±2.1) × 10 8 . This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D 3 He-p yields of 5×10 5 . This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Synchrotron self-Compton in a radiative-adiabatic fireball scenario: modelling the multiwavelength observations in some Fermi /LAT bursts

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide a unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium, we present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of $1\lt p\lt 2$ and $2\le p$. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\ge 100$ MeV in a sample of nine bursts from the second Fermi/LAT GRB catalogue that exhibited temporal and spectral indices with $\gtrsim 1.5$ and $\approx 2$, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modelled the multiwavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the mechanism responsible for explaining the energetic GeV photons.

acceleration of particles↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. IV. Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at z ≤ 2

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift z = 2 and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at z = 4 and 3, and before the last major merger, focusing on the formation of well-defined rotationally supported disks, the mass–metallicity relation, the specific star formation rate, the gas metallicity gradients, and the nonaxisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semianalytic models at z ~ 2. We see that almost all the hydro codes develop rotationally supported structures at low redshifts. Most agree within 0.5 dex with the observed mass–metallicity relation at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the intercode differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to z = 2 and the auxiliary data will be made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Physics-based modeling and data analytics [Slides]

This presentation contains a summary of ongoing work within the physics-based modeling and data analytics work package within the Nuclear Materials Discovery and Qualification initiative (NMDQi). Topics include work on MOOSE-based crystal plasticity, molecular dynamics modeling of recombination in metals and alloys, the MOOSE Stochastic Tools Module, and machine learning and atomistic modeling to predict thermo-kinetic properties of nuclear structural materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

SPT-3G D1: A Measurement of Secondary Cosmic Microwave Background Anisotropy Power

We report new measurements of millimeter-wave temperature power spectra in the angular multipole range $1700 \le \ell \le 11,000$ (wavelengths $13^\prime \gtrsim λ\gtrsim 2^\prime$). We use two years of data in three observing bands centered near 95, 150, and 220 GHz from the SPT-3G receiver on the South Pole Telescope that cover a 1646 deg$^2$ region of the Southern sky. Using the measured power spectra, we present constraints on the thermal and kinematic Sunyaev-Zel'dovich (SZ) effects, radio galaxies, and cosmic infrared background (CIB). We find that inferred SZ powers are dependent on the detailed modeling of the thermal SZ-CIB correlation, and to a lesser extent on the assumed angular dependence of the SZ spectra. We report constraints for simulation-based model templates as well as fits where the angular dependencies of the SZ and CIB power spectra are allowed to vary. In the latter case at $\ell=3000$, we find thermal SZ power at 143 GHz of $D_{3000}^{\rm tSZ} = 4.91\pm0.37\, μ{\rm K}^2$ and kinematic SZ power of $D_{3000}^{\rm kSZ} =1.75\pm0.86\, μ{\rm K}^2$. We use the measured kinematic SZ power to estimate the duration of reionization, noting that the reionization inferences are sensitive to the model choices and assumed level of homogeneous kinematic SZ power from the late-time universe. We find a 95% limit on the duration from an ionization fraction of 25% to 75% of $Δ^{50} z_{\rm re} <\,3.8$ based on a semi-analytic model, or a limit on the duration from an ionization fraction of 5% to 95% of $Δ^{90} z_{\rm re} <\,6.1$ based on the AMBER simulations.

Chaubal, P. [Melbourne U.]↗
DOE OSTI2024

Ne viii in the Warm-hot Circumgalactic Medium of FIRE Simulations and in Observations

Abstract The properties of warm-hot gas around ∼ L * galaxies can be studied with absorption lines from highly ionized metals. We predict Ne viii column densities from cosmological zoom-in simulations of halos with masses in ∼10 12 and ∼10 13 M ☉ from the Feedback in Realistic Environments (FIRE) project. Ne viii traces the volume-filling, virial-temperature gas in ∼10 12 M ☉ halos. In ∼10 13 M ☉ halos the Ne viii gas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially at M ⋆ ≳ 10 10.5 M ☉ . Median Ne viii columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne viii profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne viii columns. The circumgalactic Ne viii observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

Wijers, Nastasha A. (ORCID:0000000163747185)↗
DOE OSTI2025

Global models predict clouds at the wrong time of day: Does it matter for radiation and climate?

Accurate prediction of future climate change hinges upon the ability of Earth system models (ESMs) to simulate clouds and their radiative effects. Even if an ESM can simulate the correct clouds, a systematic error in the amount of sunlight reflected by clouds (and, thus, cloud radiative effect) can exist if the clouds are simulated at the incorrect time of day. In this work, we develop an analytical model connecting diurnal cloud biases to emergent mean state radiative biases. With the use of satellite observations, we demonstrate that there are errors in the time of day that clouds are occurring in ESMs that would cause bias in shortwave cloud radiative effect (SWCRE) that is greater than 45% of the total SWCRE bias, but such errors in the cloud diurnal cycle are masked by other compensating errors, indicating that these ESMs are getting the right answer for the wrong reasons.

Geosciences↗
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 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 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↗
DOE OSTI2024

Weak lensing combined with the kinetic Sunyaev–Zel’dovich effect: a study of baryonic feedback

ABSTRACT Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev–Zel’dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, ‘baryonification’. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model, and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters, with the joint WL + kSZ analysis constraining $S_8=0.823^{+0.019}_{-0.020}$. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

79 ASTRONOMY AND ASTROPHYSICS↗