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 163 records · Page 9

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗

Tailoring resonant magnetic perturbation to optimize fast-ion confinement during ELM control in KSTAR

Abstract 3D resonant magnetic perturbation (RMP) is one promising way to control edge localized modes that can cause excessive material erosion of tokamak first walls. However, RMP can lead to undesired degradation of plasma confinement, including fast-particle losses, which can impact the performance and safety of the reactor. This work investigates the optimization of the poloidal spectrum of the 3D field to optimize fast ion confinement during edge localized mode (ELM) suppression. In the initial step, the validity of the modeling framework is tested against experimental data. Simulations successfully replicate an increase in poloidal limiter temperature with different poloidal spectra. Then, the simulation shows improvement of fast ion confinement with a reduction of core resonant response, while edge resonant magnetic fields are maintained above the threshold to sustain the ELM suppression. Reduction of the core resonant fields keeps the Kolmogorov–Arnold–Moser surface and reduces the fast particle losses due to the stochastic magnetic field lines. The results highlight the potential of edge localization of the resonant fields to enhance the performance of fusion reactors, but further investigation is needed to improve the validation of this approach.

3D field↗

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗

Prediction of transport in the JET DTE2 discharges with TGLF and NEO models using the TGYRO transport code

Abstract The JET Deuterium-Tritium-Experiment Campaign 2 (DTE2) has demonstrated the highest-ever fusion energy production. To forecast the transport dynamics within these discharges, the TGLF and NEO models within the TGYRO transport code were employed. A critical development in this study is the new quasilinear transport model, TGLF-SAT2, specifically designed to resolve discrepancies identified in JET deuterium discharges. This model accurately describes the saturated three-dimensional (3D) fluctuation spectrum, aligning closely with a database of nonlinear CGYRO turbulence simulations, thereby enhancing the predictive accuracy of TGYRO simulations. In validating against the JET DTE2 discharges across two primary operating scenarios, TGYRO effectively predicted the temperature profiles within a broad radial window ( ρ ∼ 0.2–0.85), though with minor ion temperature discrepancies near the core. However, a consistent underprediction of electron density profiles by 20% across the simulation domain was noted, indicating areas for future refinement. To achieve a self-consistent steady-state solution based on the JET DTE2 discharges, an integrated modeling workflow TGYRO-STEP within the OMFIT framework was introduced. This workflow iterates among the core transport, the pedestal pressure and the MHD equilibrium, ultimately yielding a converged solution that significantly reduces dependence on experimental boundary conditions for temperature and density profiles. The integrated simulation results show negligible differences in electron density and temperature profiles compared to standalone TGYRO modeling, while the ion temperature profile is lower due to the updated boundary condition in TGYRO-STEP. The application of the TGYRO-STEP workflow to JET DTE2 discharges serves as a crucial test to validate its robustness and highlights its limitations, providing valuable insights for its potential future application in ITER and Fusion Power Plant deuterium and tritium prediction modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vector-level feedforward control of LPBF melt pool area using a physics-based thermal model

Laser powder bed fusion (LPBF) is an additive manufacturing technique that has gained popularity thanks to its ability to produce geometrically complex, fully dense metal parts. However, these parts are prone to internal defects and geometric inaccuracies, stemming in part from variations in the melt pool. Here, this paper proposes a novel vector-level feedforward control framework for regulating melt pool area in LPBF. By decoupling part-scale thermal behavior from small-scale melt pool physics, the controller provides a scale-agnostic prediction of melt pool area and efficient optimization over it. This is done by operating on two coupled lightweight models: a finite-difference thermal model that efficiently captures vector-level temperature fields and a reduced-order, analytical melt pool model. Each model is calibrated separately with minimal single-track and 2D experiments, and the framework is validated on a complex 3D geometry in both Inconel 718 and 316L stainless steel. Results showed that feedforward vector-level laser power scheduling reduced geometric inaccuracy in key dimensions by 62%, overall porosity by 16.5%, and photodiode root-mean-squared deviation by 38.5% on average. Overall, this modular, data-efficient approach demonstrates that proactively compensating for known thermal effects can significantly improve part quality while remaining computationally efficient and readily extensible to other materials and machines.

Additive manufacturing↗

Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2

Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Low-cost fin-tube heat exchanger design for building thermal energy storage using phase change material

Heat transfer in phase change materials (PCMs) is complex because the melting and freezing fronts change as functions of stored or released heat. In prior attempts to optimize heat exchangers (HXs) in one or two dimensions, complex geometry has often been used to maximize the melt and freeze front area. This complex geometry is difficult and hence expensive to construct. This paper proposes a multiple-scale 3D finite element modeling approach to design fin-tube HXs for low-cost latent thermal energy storage applications. Further, the optimal fin and tube designs were determined at three scales (unit-scale, medium-scale, and large-scale) by modeling the melt and freeze front in three dimensions and using measured bulk thermal properties. The finite element model was validated by comparing it with the experimental data for a referenced design of a similar type. The results indicate that commercially available organic PCMs with low conductivity (<0.3 W/m·K) can have charge and discharge times appropriate for building thermal energy storage (i.e., 4–5 h) with fin-tube HX designs at costs <$26/kWh, even when the temperature difference (5.56°C) between the heat transfer fluid and the PCM phase change temperature is small. However, as the HX increases in length, the temperature reduction along the tube limits some larger-scale designs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Predicting 3D moisture sorption behavior of materials from 1D investigations

Abstract Moisture in materials can be a source of future outgassing and exacerbate unwanted changes in physical and chemical properties. Here, we investigate the effect of sample size and shape on the moisture transport phenomena through a combined experimental and modeling approach. Several different materials varying in size and shape were investigated over a wide range of relative humidities (0–90%) and temperatures ( $$30{-}70^{\,\circ } \hbox {C}$$ 30 - 70 ∘ C ) using gravimetric type dynamic vapor sorption (DVS). A dynamic triple-mode sorption model, developed previously, was employed to describe the experimental results with good success; the model includes absorption, adsorption, pooling (clustering) of species, and molecular diffusion. Here we show that the full triple-mode sorption model is robust enough to predict the dynamic uptake and outgassing of 3-dimensional (3D) samples using parameters derived from quasi-1D samples. This successful demonstration on three different materials (filled polydimethylsiloxane (PDMS), unfilled PDMS, and ceramic inorganic composite) illustrates that the model is robust at describing the scale-independent physics and chemistry of moisture sorption and diffusion materials. This work demonstrates that while sorption mechanisms manifest in testing of all sample sizes, some of these mechanisms were so subtle that they were overlooked in our initial modeling and assessment, illustrating the importance of multi-scale experiments in the development of robust predictive capabilities. Our study also outlines the challenges and viable solutions for global optimization of a multi-parameter model. The ability to quantify moisture sorption and diffusion, independent of scale, using 1D lab-scale experiments enables prediction of long-term bulk materials behavior in real applications.

36 MATERIALS SCIENCE↗

3D printed carbon fiber reinforced carbon as an energy efficient alternative to graphite for EFAS tooling

As Electric Field Assisted Sintering (EFAS) gains more industrial acceptance and use, it becomes more important to develop more efficient means to implement this technology. To this aim, 3D printed continuous carbon fiber reinforced carbon (CCC) was manufactured and fabricated into tooling for EFAS systems as an alternative to traditional graphite tooling. The impact of fiber orientation on the thermal and electrical properties of the CCC was characterized. Sample material was sintered in Tokai G535 graphite tooling, under common processing conditions and compared with CCC tooling. There was nearly 50 % energy savings compared to graphite while maintaining equivalent sample density and microstructure plus keeping ram temperatures 39 % cooler. This is due to spatial control of generated heat and thermal diffusivity within the molds, by means of fiber orientation anisotropy. Finite element modeling of the tooling design supported the experimental results as well as displays the effect of optimization of this 3D printed CCC material.

36 MATERIALS SCIENCE↗

Antiferromagnetic transitions of Dirac fermions in three dimensions

We use determinant quantum Monte Carlo simulations to study the role of electron-electron interactions on three-dimensional (3D) Dirac fermions based on the π-flux model on a cubic lattice. We show that the Hubbard interaction drives the 3D Dirac semimetal to an antiferromagnetic (AF) insulator only above a finite critical interaction strength and the long-range AF order persists up to a finite temperature. We evaluate the critical interaction strength and temperatures using finite-size scaling of the spin structure factor. The critical behaviors are consistent with the (3+1)-dimensional Gross-Neveu universality class for the quantum critical point and 3D Heisenberg universality class for the thermal phase transitions. We further investigate correlation effects in the birefringent Dirac fermion system. It is found that the critical interaction strength U c is decreased by reducing the velocity of the Dirac cone, quantifying the effect of velocity on the critical interaction strength in 3D Dirac fermion systems. Our findings unambiguously uncover correlation effects in 3D Dirac fermions and may be observed using ultracold atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗

Advanced modeling of field enhanced thermionic emission

Shaped emitters are of interest to a broad range of applications in vacuum electronic devices. In particular, thermionic energy converters (TECs) take advantage of shaped emitters to increase the local surface field, thereby extracting more current for a given cathode temperature and applied voltage. However, modeling these devices is challenging; Warp [J.-L. Vay, D. P. Grote, R. H. Cohen, and A. Friedman, Comput. Sci. Discov. 5, 014019 (2012)] is a fully 3D particle-in-cell code capable of handling a wide range of physics problems and is well suited to modeling TECs. Additionally, recent improvements to Warp have enabled the accurate modeling of emitters with arbitrary curved surfaces. Specifically, the inclusion of subgrid resolution for computing the electrostatic potential and the ability to apply mesh refinement for specific areas of interest allow for a more accurate solution to the fields on these surfaces. These improvements coupled with Warp’s ability to handle variable particle weights make it an ideal candidate for simulating these complex devices. In this paper, the authors study the applicability of different subgrid configurations for simulating shaped emission surfaces and field convergence for different mesh-refinement techniques. They then implement a custom weighting algorithm that allows for uniform sampling of emission surfaces with a large variation in the surface electric field. They then use this algorithm to study emission for curved emitters in both the field-enhancement regime and the space-charge regime.

Edelen, Jonathan P. (ORCID:0000000215180652)↗

Tardigrade: Summary (LL20-AM-sensor)

Tardigrade is an effort to create small, passive remote sensor for material detection and identification. The microbolometer based sensor design uses custom bandpass filter bands combined with tuned thermal expansion materials in a ruggedized lens housing to provide passive mechanical correction for optical changes across a wide temperature range. This requires developing computational tools to model and predict thermal performance of meta-materials composed of multiple structured base materials and void-space, along with matching and developing 3D print technology to demonstrate the housing. Current efforts in the project are in 3 areas: material design algorithm development (transitioning into full 3D material design space), multi-material additive manufacturing process exploration, and sensor/optics design and engineering.

42 ENGINEERING↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Atomic structure and Mott nature of the insulating charge density wave phase of 1T-TaS 2

Here, using x-ray pair distribution function (PDF) analysis and computer modeling, we explore structure models for the complex charge density wave (CDW) phases of layered 1T-TaS 2 that both well capture their atomic-level features and are amenable to electronic structure calculations. The models give the most probable position of constituent atoms in terms of 3D repetitive unit cells comprising a minimum number of Ta–S layers. Structure modeling results confirm the emergence of star-of-David (SD) like clusters of Ta atoms in the high-temperature incommensurate (IC) CDW phase and show that, contrary to the suggestions of recent studies, the low-temperature commensurate (C) CDW phase expands upon cooling thus reducing lattice strain. The C-CDW phase is also found to preserve the stacking sequence of Ta–S layers found in the room temperature, nearly commensurate (NC) CDW phase to a large extent. DFT based on the PDF refined model shows that bulk C-CDW 1T-TaS 2 also preserves the insulating state of individual layers of SD clusters, favoring the Mott physics description of the metal-to-insulator (NC-CDW to C-CDW) phase transition in 1T-TaS 2 . Our work highlights the importance of using precise crystal structure models in determining the nature of electronic phases in complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

RELAP5-3D Solutions to Exercise 1 of the OECD-NEA HTTF Benchmark

The high-temperature gas-cooled reactor (HTGR) is an advanced reactor concept that has received considerable attention over the last 60 years. HTGRs boast high coolant outlet temperatures, passive safety, and large margins to fuel failure. Recent years have seen increasing interest in HTGRs, whether large-scale HTGRs for electricity or process heat applications or as microreactors serving remote communities. HTGRs have considerable potential to fulfill these applications, but modeling and simulation tools typically used for reactor safety analysis have not been validated for HTGR modeling. The High-Temperature Test Facility (HTTF) was constructed at Oregon State University (OSU) to provide validation data for HTGR modeling. Recently, Idaho National Laboratory (INL), OSU, Argonne National Laboratory, Canadian Nuclear Laboratories, and the University of Tennessee Knoxville have collaborated to develop an HTGR thermal hydraulics benchmark based on HTTF experiments. This summary presents results from a portion of that benchmark using RELAP5-3D

22 GENERAL STUDIES OF NUCLEAR REACTORS↗