Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D temperature modeling”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Magneto-electroluminescence response in 2D and 3D hybrid organic–inorganic perovskite light emitting diodes

We have studied the magneto-electroluminescence (MEL) response in light emitting diodes based on 2D-(PEA) 2 PbI 4 and 3D-MAPbI 3 hybrid organic–inorganic perovskites at cryogenic temperatures. We found that the MEL is negative, i.e., the EL decreases with the applied field strength, B. In addition, the MEL(B) response has a Lorentzian line shape whose width depends on the perovskite used. We interpret the MEL(B) response using the “Δg mechanism” in which the spin of the injected electron–hole (e–h) pairs oscillates between singlet and triplet configurations due to different precession frequencies of the electron and hole constituents that originate from the difference, Δg, in the electron and hole gyromagnetic constants, g. In this model, the MEL(B) linewidth is inversely proportional to the spin lifetime and Δg. The model used is validated by directly measuring the spin lifetime of photogenerated e–h pairs using the circularly polarized pump–probe transmission technique with 100 fs resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pressure tuning of Kitaev spin liquid candidate Na 3 Co 2 SbO 6

The search for Kitaev’s quantum spin liquid in real materials has recently expanded with the prediction that honeycomb lattices of divalent, high-spin cobalt ions could host the dominant bond-dependent exchange interactions required to stabilize the elusive entangled quantum state. The layered honeycomb Na 3 Co 2 SbO 6 has been singled out as a leading candidate provided that the trigonal crystal field acting on Co 3d orbitals, which enhances non-Kitaev exchange interactions between $J$ eff = $\frac{1}{2}$ spin-orbital pseudospins, is reduced. Here we show that applied pressure leads to anisotropic compression of the layered structure, significantly reducing the trigonal distortion of CoO 6 octahedra. Ferromagnetic correlations between pseudospins are enhanced in the spin-polarized (3 Tesla) phase up to about 60 GPa. Higher pressures drive a high-spin to low-spin transition destroying the $J$ eff = $\frac{1}{2}$ moments required to map the spin Hamiltonian into Kitaev’s model. The spin transition strongly suppresses the low-temperature magnetic susceptibility and appears to stabilize a paramagnetic phase driven by frustration. The possible emergence of frustrated magnetism of localized $S$ = $\frac{1}{2}$ moments opens the door for exploration of novel magnetic quantum states in compressed honeycomb lattices of divalent cobaltates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Three-dimensional diagnostics and measurements of inertial confinement fusion plasmas

Recent inertial confinement fusion measurements have highlighted the importance of 3D asymmetry effects on implosion performance. One prominent example is the bulk drift velocity of the deuterium–tritium plasma undergoing fusion (“hotspot”), v HS . Upgrades to the National Ignition Facility neutron time-of-flight diagnostics now provide v HS to better than 1 part in 10 4 and enable cross correlations with other measurements. Here we present the impact of v HS on the neutron yield, downscatter ratio, apparent ion temperature, electron temperature, and 2D x-ray emission. The necessary improvements to diagnostic suites to take these measurements are also detailed. The benefits of using cross-diagnostic analysis to test hotspot models and theory are discussed, and cross-shot trends are shown.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Resolving atomic transport through iron oxide under irradiation using isotopic tracers

Structural materials in nuclear environments are exposed to complex and harsh conditions, including high temperatures, corrosive solutions, and irradiation. Material stability and likewise degradation are fundamentally controlled by atomic transport mechanisms. Here, we investigate the influence of radiation-induced, non-equilibrium point defect populations on mass transport within model oxide systems using isotopic tracers and atom probe tomography. An isotopically labelled hematite (a-Fe2O3) film was synthesized with 16O- and 18O-rich layers (>30nm and 10 nm, respectively) sequentially deposited via molecular beam epitaxy. One region of the film was proton irradiated to ~0.1 dpa at 450°C, and regions outside the beam spot provided a non-irradiated control of thermal annealing only. 3D APT was used to directly visualize the elemental and O-isotopic redistribution following both thermal and combined thermal and irradiation conditions. Complementary scanning transmission electron microscopy (STEM) was also employed to study the microstructural evolution. Error function fits of the 18O-concentration profiles enabled quantification of diffusion coefficients.

isotopic tracers, diffusion, proton irradiation, h↗
DOE OSTI2023

Long‐Range Electron Transport Rates Depend on Wire Dimensions in Cytochrome Nanowires

Abstract The ability to redirect electron transport to new reactions in living systems opens possibilities to store energy, generate new products, or probe physiological processes. Recent work by Huang et al. showed that 3D crystals of small tetraheme cytochromes (STC) can transport electrons over nanoscopic to mesoscopic distances by an electron hopping mechanism, making them promising materials for nanowires. However, fluctuations at room temperature may distort the nanostructure, hindering efficient electron transport. Classical molecular dynamics simulations of these fluctuations at the nano‐ and mesoscopic scales allowed us to develop a graph network representation to estimate maximum electron flow that can be driven through STC wires. In longer nanowires, transient structural fluctuations at protein‐protein interfaces tended to obstruct efficient electron transfer, but these blockages are ameliorated in thicker crystals where alternative electron transfer pathways become more efficient. The model implies that more flexible proteinprotein interfaces limit the required minimum diameter to carry currents commensurate with conventional electronics.

Chemistry↗
DOE OSTI2025

Boosting H I -Galaxy Cross-Clustering Signal through Higher-Order Cross-Correlations

After reionization, neutral hydrogen (${\rm H\, \small {I}}$) traces the large-scale structure (LSS) of the Universe, enabling ${\rm H\, \small {I}}$ intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study ${\rm H\, \small {I}}$ clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of k-nearest neighbor (kNN) optical galaxies and the ${\rm H\, \small {I}}$ brightness temperature field smoothed at relevant scales (the kNN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate ${\rm H\, \small {I}}$ data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with k ∥ < k min,∥ . Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the ${\rm H\, \small {I}}$-galaxy cross-correlation signal in our simulations, using the kNN-field framework, is detectable at >30σ across r = [3, 12] h –1 Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp k ∥ filter can spuriously boost detection to 8σ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that kNN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

4D structural changes and pore network model of biomass during pyrolysis

Biochar is an engineered carbon-rich substance used for soil improvement, environmental management, and other diverse applications. To date, the understanding of how biomass affects biochar microstructure has been limited due to the complexity of analysis involved in tracing the changes in the physical structure of biomass as it undergoes thermochemical conversion. In this study, we used synchrotron x-ray micro-tomography to visualize changes in the internal structure of biochar from diverse feedstock (miscanthus straw pellets, wheat straw pellets, oilseed rape straw pellets, and rice husk) during pyrolysis by collecting a sequence of 3D scans at 50 °C intervals during progressive heating from 50 °C to 800 °C. The results show a strong dependence of biochar porosity on feedstock as well as pyrolysis temperature, with observed porosity in the range of 7.41–60.56%. Our results show that the porosity, total surface area, pore volume, and equivalent diameter of the largest pore increases with increasing pyrolysis temperature up to about 550 °C. The most dramatic development of pore structure occurred in the temperature range of 350–450 °C. This understanding is pivotal for optimizing biochar’s properties for specific applications in soil improvement, environmental management, and beyond. By elucidating the nuanced variations in biochar’s physical characteristics across different production temperatures and feedstocks, this research advances the practical application of biochar, offering significant benefits in agricultural, environmental, and engineering contexts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Electronic structure, dimer physics, orbital-selective behavior, and magnetic tendencies in the bilayer nickelate superconductor $\mathrm{La_3Ni_2O_7}$under pressure

Motivated by the recently reported high-temperature superconductivity in the bilayer La 3 Ni 2 O 7 (LNO) under pressure, we comprehensively study this system using ab initio techniques. The Ni 3d orbitals have a large bandwidth at ambient pressure, increasing by ~22% at 29.5 GPa. Without electronic interactions, the Ni d 3z 2 -r 2 orbitals form a bonding-antibonding molecular orbital state via the O p z inducing a “dimer” lattice in the LNO bilayers. The Fermi surface consists of two-electron sheets with mixed e g orbitals and a hole pocket defined by the d 3z 2 -r 2 orbital, suggesting a Ni two-orbital minimum model. Different from the infinite-layer nickelate, we obtained a large interorbital hopping between d 3z 2 -r 2 and d x 2 -y 2 states in LNO, caused by the ligand “bridge” of in-plane O p x or p y orbitals connecting those two orbitals, inducing d-p σ-bonding characteristics. The competition between the intraorbital and interorbital hoppings leads to an interesting dominant spin stripe (π,0) order because of bond ferromagnetic tendencies via the recently discussed “half-empty” mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2022

Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events

Abstract We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Pressure tuning of Kitaev spin liquid candidate Na$_3$Co$_2$SbO$_6$goo

The search for Kitaev's quantum spin liquid (KQSL) state in real materials has recently expanded with the prediction that honeycomb lattices of divalent, high-spin cobalt ions could host the dominant bond-dependent exchange interactions required to stabilize the elusive entangled quantum state. The layered honeycomb Na$_3$Co$_2$SbO$_6$ has been singled out as a leading candidate provided that the trigonal crystal field acting on Co $3d$ orbitals, which enhances non-Kitaev exchange interactions between $J_{\rm eff}=\frac{1}{2}$ spin-orbital pseudospins, is reduced. We find that applied pressure leads to anisotropic compression of the layered structure, significantly reducing the trigonal distortion of CoO$_6$ octahedra. A strong enhancement of ferromagnetic correlations between pseudospins is observed in the spin-polarized (3 Tesla) phase up to about 60 GPa. Higher pressures drive a spin transition into a low-spin state destroying the $J_{\rm eff}=\frac{1}{2}$ local moments required to map the spin Hamiltonian into Kitaev's model. The spin transition strongly suppresses the low-temperature magnetic susceptibility and appears to stabilize a paramagnetic phase driven by frustration. Although applied pressure fails to realize a KQSL state, the possible emergence of frustrated magnetism of localized, low-spin $S=\frac{1}{2}$ moments opens the door for exploration of novel magnetic quantum states in compressed honeycomb lattices of divalent cobaltates.

FOS: Physical sciences↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2024

Unraveling Reactivity Differences: Room-Temperature Ring-Opening Metathesis Polymerization (ROMP) versus Frontal ROMP

Here, in this study, we explore the distinct reactivity patterns between frontal ring-opening metathesis polymerization (FROMP) and room-temperature solventless ring-opening metathesis polymerization (ROMP). Despite their shared mechanism, we find that FROMP is less sensitive to inhibitor concentration than room-temperature ROMP. By increasing the initiator-to-monomer ratio for a fixed inhibitor/initiator quantity, we find reduction in the ROMP background reactivity at room temperature (i.e., increased resin pot life). At elevated temperatures where inhibitor dissociation prevails, accelerated frontal polymerization rates are observed because of the concentrated presence of the initiator. Surprisingly, the strategy of employing higher initiator loading enhances both pot life and front speeds, which leads to FROMP rates exceeding prior reported values by over 5 times. This counterintuitive behavior is attributed to an increase in the proximity of the inhibitor to the initiator within the bulk resin and to whether the temperature favors coordination or dissociation of the inhibitor. A rapid method was developed for assessing resin pot life, and a straightforward model for active initiator behavior was established. Modified resin systems enabled direct ink writing of robust thermoset structures at rates much faster than previously possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Two-Dimensional and Three-Dimensional Tetrel-Arsenide Frameworks Templated by Li and Cs Cations

Two new ternary tetrel pnictide phases, Li 1.08 Ge 4.92 As 7 and Cs 0.11 Zn 0.05 Si 0.95 As, have been synthesized. Both compounds crystallize in new structure types, which demonstrate the flexibility of Ge(Si)–As layers when accommodating electropositive cations of various sizes. Li 1.08 Ge 4.92 As 7 displays a Ge-As layer motif that has not been observed in ternary Ge-As systems, while Cs 0.11 Zn 0.05 Si 0.95 As exhibits a 3D framework composed of fragments of the 2D parent-phase SiAs. For Li 1.08 Ge 4.92 As 7 , a combination of computational modeling and solid-state 7 Li MAS NMR was used to establish the correct crystal structure with one mixed occupied Ge/Li site. The compositions of the phases were confirmed through energy-dispersive X-ray spectroscopy. Transport properties and UV–Vis spectroscopy reveal both reported compounds to be semiconductors. Cs 0.11 Zn 0.05 Si 0.95 As shows high thermal stability with congruent melting at 1300 K and an exceptionally low thermal conductivity of 0.5 W m –1 K –1 at 300 K. Temperature-dependent Li-ion hopping dynamics in Li 1.08 Ge 4.92 As 7 was characterized using 7 Li solid-state NMR.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extensive 3D mapping of dislocation structures in bulk aluminum

Abstract Thermomechanical processing such as annealing is one of the main methods to tailor the mechanical properties of materials, however, much is unknown about the reorganization of dislocation structures deep inside macroscopic crystals that give rise to those changes. Here, we demonstrate the self-organization of dislocation structures upon high-temperature annealing in a mm-sized single crystal of aluminum. We map a large embedded 3D volume ( $$100\times 300\times 300$$ 100 × 300 × 300 $$\upmu$$ μ m $$^3$$ 3 ) of dislocation structures using dark field X-ray microscopy (DFXM), a diffraction-based imaging technique. Over the wide field of view, DFXM’s high angular resolution allows us to identify subgrains, separated by dislocation boundaries, which we identify and characterize down to the single-dislocation level using computer-vision methods. We demonstrate how even after long annealing times at high temperatures, the remaining low density of dislocations still pack into well-defined, straight dislocation boundaries (DBs) that lie on specific crystallographic planes. In contrast to conventional grain growth models, our results show that the dihedral angles at the triple junctions are not the predicted 120 $$^{\circ }$$ ∘ , suggesting additional complexities in the boundary stabilization mechanisms. Mapping the local misorientation and lattice strain around these boundaries shows that the observed strain is shear, imparting an average misorientation around the DB of $$\approx$$ ≈ 0.003 to 0.006 $$^{\circ }$$ ∘ .

36 MATERIALS SCIENCE↗