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

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Static and dynamic magnetic properties of the spiin-$\frac{5}{2}$triangle lattice antiferromagnet Na 3 Fe(PO 4 ) 2 studied by 31 P NMR

We report 31 P nuclear magnetic resonance (NMR) measurements have been carried out to investigate the magnetic properties and spin dynamics of Fe 3+ (S = 5/2) spins in the two-dimensional triangular lattice (TL) compound Na 3 Fe(PO 4 ) 2 . The temperature (T) dependence of nuclear spin-lattice relaxation rates (1/T 1 ) shows a clear peak around Néel temperature, T N = 10.9 K, corresponding to an antiferromagnetic (AFM) transition. From the temperature dependence of NMR shift (K) above T N , an exchange coupling between Fe 3+ spins was estimated to be J/k B ≃ 1.9 K using the spin-5/2 Heisenberg isotropic-TL model. The temperature dependence of 1/T 1 T divided by the magnetic susceptibility (χ), 1/T 1 T χ , above T N proves the AFM nature of spin fluctuations below ~50 K in the paramagnetic state. In the magnetically ordered state below T N , the characteristic rectangular shape of the NMR spectra is observed, indicative of a commensurate AFM state in its ground state. The strong temperature dependence of 1/T 1 in the AFM state is well explained by the two-magnon (Raman) process of the spin waves in a 3D antiferromagnet with a spin-anisotropy energy gap of 5.7 K. The temperature dependence of sublattice magnetization is also well reproduced by the spin waves. Those results indicate that the magnetically ordered state of Na 3 Fe(PO 4 ) 2 is a conventional 3D AFM state, and no obvious spin frustration effects were detected in its ground state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
DOE OSTI2020

Adaptive Machine Learning for Bragg Coherent Diffraction Imaging (BCDI) of 3D Electron Density Maps with Application to La 2-x Ba x CuO 4 (LBCO) High Temperature Superconductor Studies [PowerPoint]

Understanding mesoscale heterogeneity is important for evaluating fatigue and failure of structural materials and for manufacturing processes. Data-driven tools can facilitate in guiding effort investment (measurements and computations) at various stages of the materials development. Our 3D reconstruction approach is to find a set of coefficients which define a bounding surface of the electron density. Data for training the 3D CNN was generated by sampling coefficients from uniform distributions. Test vs. prediction values are shown for the 28 even-valued coefficients for 1000 test structures, with the best and worst performers highlighted. The adaptive part of this work utilized a model independent extremum seeking (ES). The predictions of the CNN were used as the starting point for the ES algorithm. The convergence of the ES algorithm for 3 different structures is illustrated, and the robustness of the adaptive ML approach was demonstrated on experimentally measured 3D crystal from high energy diffraction microscopy. Reconstruction of a non-uniform density volume with different levels of Poisson noise is also shown. Preliminary attempts at reconstructing measured La 2-x Ba x CuO 4 diffraction patterns have so far failed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

New heat flux model for non-axisymmetric divertor infrared structures

A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Heuristic predictions of RMP configurations for ELM suppression in ITER burning plasmas and their impact on divertor performance

Abstract A subspace of resonant magnetic perturbation (RMP) configurations for edge localized mode (ELM) suppression is predicted for H-mode burning plasmas at 15 MA current and 5.3 T magnetic field in ITER. Perturbations to the core plasma can be reduced by a factor of 2 for equivalent edge stability proxies, while the perturbed plasma boundary geometry remains mostly resilient. The striated domain of perturbed field lines connecting from the main plasma (normalized poloidal flux < 1 ) to the divertor targets is found to be significantly larger than the expected heat load width in the absence of RMPs. This facilitates heat load spreading with peak values at an acceptable level below 10 MW m − 2 on the outer target already at moderate gas fueling and low Ne seeding for additional radiative dissipation of the 100 MW of power into the scrape-off layer (SOL). On the inner target, however, re-attachment is predicted away from the equilibrium strike point due to increased upstream heat flux, higher downstream temperature and less efficient impurity radiation.

3D plasma boundary modeling↗
DOE OSTI2025

Adhesion of Self-Complementary, Sinusoidal Surfaces Fabricated Using Two-Photon Polymerization

Microscale, pick-and-place assembly is a non-lithographic assembly method poised to impact diverse fields including flexible electronics, microfluidics and robotics. However, a major technological challenge is the need to deterministically control adhesion between parts. Here, switchable adhesion involving 3D-printed, self-complementary surfaces is demonstrated. Mechanical properties of metasurfaces pressed against flat, rigid substrates are modeled using finite element methods. A series of flat slabs and metastructured slabs with 2D sinusoidal surfaces are printed using two-photon polymerization (2PP) of a shape-memory resin. The surface frequency of featured slabs was varied between $3.\bar3$ mm −1 and $26.\bar6$ mm −1 with similar amplitudes. Adhesion between printed metasurfaces and glass and between printed, self-complementary metasurfaces is studied above and below the cured resin’s glass transition temperature (∼45 °C). Simple heating of adhering surfaces to above 60 °C lowers adhesion, and compression of surfaces while above the glass transition temperature followed by cooling to room temperature elevates adhesion. The nominal adhesive strength between printed, self-complementary surfaces, as determined by the maximum observable pull-off stress, exceeds 3 MPa. Further tailoring complementary surfaces for adhesion control may facilitate microscale disassembly for recovery of components or precious metals.

amorphous materials↗
DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Coaxial Ceramic Direct Ink Writing on Heterogenous and Rough Surfaces: Investigation of Core–Shell Interactions

In this work, coaxial conductor–ceramic direct ink writing enables the printing of sensitive or encapsulated materials onto heterogeneous and rough substrates. While encasing the core fluid within a stiff ceramic shell, continuity may be maintained, even while printing onto conventionally challenging substrates. Here, we report the development of a coaxial ceramic direct ink writing suite and explore coflow interrelationships based on microfluidic principles. A coaxial nozzle is designed to facilitate the coextrusion of an alumina shell, whereas indium–tin-oxide inks constitute the core. In this manner, a core–shell ceramic element may be printed onto rough substrates for future high-temperature applications. Colloidal inks are engineered to provide the required rheological and sintering performance. Moreover, flow simulations in conjunction with microfluidic coflow principles are used to explore the coaxial printing processing space, thus controlling the core–shell architectures. Physical modeling is further used to analyze core deformations and eccentricity. As a result, simulations are validated experimentally, and the analyses are used to deposit coaxial ceramic features onto heterogeneous, high-temperature ceramic substrates.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the 3D printability of one-part moisture-curable polyurethanes via direct ink writing (DIW)

Direct ink writing (DIW) is an ambient temperature additive manufacturing (AM) method that accommodates many possible ink materials. Here, in this work, we demonstrate using a moisture-curable commercially available polyurethane (PU) sealant as an alternative ink for DIW. We discussed the fundamentals of PU chemistry and determined the best 3D printing parameters. Studies on rheological, thermogravimetric, spectroscopic characterization, and initial finite element analyses (FEA) showed properties expected from a performance sealant with high elongation and low modulus of a 3D-printed object. This affirms the flexibility of the DIW technique as an accessible AM method amenable for future materials development from commercial model formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2024

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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↗