Search NASA⌕ Search

SEARCH · Search NASA

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

Oxidation kinetics of SPS-densified U 3 Si 2 fuels—Microstructure impact

Here, U 3 Si 2 is a potential candidate for accident tolerant fuels because of its high uranium density and excellent thermal conductivity in comparison to UO 2 . However, U 3 Si 2 suffers from oxidation, steam corrosion, and subsequent disintegration/pulverization. The detailed investigation of kinetics that incorporates fundamental treatment of oxidation of U 3 Si 2 is scarcely reported, and the oxidation mechanisms have not been fully elucidated. In this paper, the oxidation behavior of microcrystalline (mc - ) and nanocrystalline (nc - ) U 3 Si 2 have been systematically investigated using a thermogravimetric analysis (TGA) apparatus through a series of isothermal and non-isothermal kinetic studies. The isothermal kinetic study with a model-fitting approach indicates oxidation activation energy of 85 kJ/mol for dense mc - U 3 Si 2 and 96.4 kJ/mol for nc - U 3 Si 2 pellets, while the isoconversional approach leads to an activation energy in the range of 70–85 kJ/mol for mc - U 3 Si 2 and 75–86 kJ/mol for nc - U 3 Si 2 with three most common model-free methods, including Kissinger–Akahira–Sunose, Flynn–Wall–Ozawa, and Friedman methods. The derivation of oxidation activation energies using both isothermal and isoconversional methods highlights the approach to evaluate the oxidation resistance of nuclear materials using TGA quantitatively and makes it possible to compare among various nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Basis for the ICRP’s updated biokinetic model for systemic astatine

The International Commission on Radiological Protection (ICRP) recently updated its biokinetic models for workers in a series of reports called the OIR (occupational intakes of radionuclides) series. A new biokinetic model for astatine (At), the heaviest member of the halogen family, was adopted in OIR Part 5 (ICRP in press). Occupational intakes of radionuclides: Part 5). Furthermore, this paper provides an overview of available biokinetic data for At; describes the basis for the ICRP's updated model for At; and tabulates dose coefficients for intravenous injection of each of the two longest lived and most important At isotopes, 211 At and 210 At. At-211 (T 1/2 = 7.214 h) is a promising radionuclide for use in targeted α-particle therapy due to several favourable properties including its half-life and the absence of progeny that could deliver significant radiation doses outside the region of α-particle therapy. At-210 (T 1/2 = 8.1 h) is an impurity generated in the production of 211 At in a cyclotron and represents a potential radiation hazard via its long-lived progeny 210 Po (T 1/2 = 138 days). Tissue dose coefficients for injected 210 At and 211 At based on the updated model are shown to differ considerably from values based on the ICRP's previous model for At, particularly for the thyroid, stomach wall, salivary glands, lungs, spleen, and kidneys.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

MCNP calculations concerning the UW 137 Cs release event (Final Report)

This document supersedes earlier reports (LA-UR-19-32358 and LA-UR-20-22216) that discussed MCNP-calculated conversion factors (count rate per unit activity per unit area) for collimated 2”x2” NaI detectors (Ludlum model 44-10 and Thermo model SPA3) in contact with duct and plenum walls. Specifically, for duct and plenum dimensions and construction details of interest to the University of Washington ¹³⁷Cs release event. This document also corrects some of the data presented in the earlier LA-UR reports.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Physics-based model of irradiation creep for ferritic materials under fusion energy operation conditions

Irradiation creep is known to be an important process for structural materials in nuclear environments, potentially leading to creep failure at temperatures where thermal creep is generally negligible. While there is a great deal of data for irradiation creep in steels and zirconium alloys in light water reactor conditions, much less is known for first wall materials under fusion energy conditions. Lacking suitable fusion neutron sources for detailed experimentation, modeling, and simulation can help bridge the dose-rate and spectral-effects gap and produce quantifiable expectations for creep deformation of first wall materials under standard fusion conditions. In this paper, we develop a comprehensive model for irradiation creep created from merging a crystal plasticity representation of the dislocation microstructure and a defect evolution simulator that accounts for the entire cluster dimensionality space. Both approaches are linked by way of a climb velocity that captures dislocation-biased defect absorption and a dislocation strengthening term that reflects the accumulation of defect clusters in the system. Here, we carry out our study in Fe under first wall fusion reactor conditions, characterized by a fusion neutron spectrum with average recoil energies of 20 keV and a damage dose rate of ≈3×10 –7 dpa/s at temperatures between 300 and 800 K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding the Applicability of Gas/Particle Partitioning Theory to SOA using Novel Chemically-Speciated Partitioning Measurements (Final Report)

The main objective of this research grant project was to investigate the physical-chemical processes and properties of particles and partitioning gases that influence gas/particle (G/P) partitioning and timescales for equilibration, to directly determine the conditions under which equilibrium G/P partitioning theory can be applied to secondary organic aerosol (SOA). SOA constitutes a major component of submicron particles in the atmosphere, which have important implications for air quality, human and ecosystem health, radiation balance, and climate; however SOA formation, evolution and losses are not well understood. Use of novel on-line measurements of chemically-speciated partitioning measurements of both gas-phase and particle-phase components was a central component to advancing these scientific goals. The key approach was to conduct environmental chamber studies with systems of increasing complexity to determine G/P partitioning rates and equilibrium states for chemically-speciated SOA components. The project builds on our previous work on gas-wall partitioning that must be understood in order to study G/P partitioning in chambers, and to help extrapolate the results to be atmospherically relevant. The goal was to employ measurements and modeling to systematically determine when G/P partitioning theory accurately describes the distribution of oxidized organic compounds in the gas and aerosol phases, and when and why it does not (e.g., aerosol phase-state kinetic limitations, non-ideal mixing, phase separation, oligomer formation,), for chemical and physical conditions typically found throughout the continental boundary layer.

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

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2024

Irradiation creep model for α-Fe: Dislocation climb in the presence of stress fields (Final Report)

Irradiation creep is known to be an important process for structural materials in nuclear environments, potentially leading to creep failure at temperatures where thermal creep is generally negligible. While there is a great deal of data for irradiation creep in steels and zirconium alloys in light water reactor conditions, much less is known for first wall materials under fusion energy conditions. Lacking suitable fusion neutron sources for detailed experimentation, modeling, and simulation can help bridge the dose-rate and spectral-effects gap and produce quantifiable expectations for creep deformation of first wall materials under standard fusion conditions. In this work, we present a comprehensive model for irradiation creep created from merging a crystal plasticity representation of the dislocation microstructure and a defect evolution simulator that accounts for the entire cluster dimensionality space. Both approaches are linked by way of a climb velocity that captures dislocation-biased defect absorption and a dislocation strengthening term that reflects the accumulation of defect clusters in the system. We carry out our study in Fe under first wall fusion reactor conditions, characterized by a fusion neutron spectrum at temperatures between 300 and 800°K.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗
DOE OSTI2026

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental apparatus and methodology to test and quantify thermal performance of micro and macro-encapsulated phase change materials in building envelope applications

Thermal energy storage (TES) is used as a viable technology to shift peak electricity demand caused by the space cooling requirements in buildings. Passive TES is implemented in building envelope via micro and macroencapsulation methods. This study describes the use of a state-of-the-art laboratory to test different PCM inclusions in simplified building walls. A microencapsulated PCM and two macroencapsulated PCMs are tested in a controlled environment to gather data for validation purposes of PCM modelling algorithms in building energy modelling programs. Data indicates that the chamber environment and the heating and cooling system can conduct full-cycle tests of wall panels with PCM inclusions. This study also generates data from parallel tests on 4 wall panels which can be used in building energy modelling programs to validate the PCM modelling algorithms. The cyclic tests also capture the thermal effects of PCMs and complex PCM behaviors like sub-cooling in PCM hydrate-salts.

25 ENERGY STORAGE↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Preliminary neutronic and thermal analysis of Inertial Fusion Energy systems with a focus on thick liquid wall concepts

This paper provides an assessment of the neutronic behavior of Inertial Fusion Energy systems using the Serpent 2 Monte Carlo code. The study explores different types of reactor concepts, including dry-wall and thick liquid wall systems. In particular, it characterizes the tritium breeding and energy multiplication performance, as well as the effect of long-term irradiation on structural materials. For this purpose, the displacement per atom (DPA) and the gas production rate in structural materials are calculated. A simplified geometry with varying geometrical parameters is used as a representation for the different kinds of devices. Additionally, the case of the HYLIFE-II concept is used to assess modeling inaccuracies arising from certain approximations. A comparison between the simplified geometry and a realistic one is performed, thanks to Serpent’s capability to run transport calculations on unstructured meshes. The study then addresses the impact of these discrepancies on the temperature distribution in the first wall using the GeN-Foam code for conjugate heat transfer analysis. The results highlight the necessity to consider the heterogeneous nature of the liquid wall, as simplified models tend to underestimate radiation damage metrics and heat loads.

Conjugate heat transfer↗
DOE OSTI2025

Localized material compression to correct distortion in wire arc additive manufacturing

Wire Arc Additive Manufacturing (WAAM) is an advanced manufacturing technology which utilizes welding systems to generate three dimensional geometries in a layer-by-layer fashion. Distortion or warping of a print substrate and WAAM components due to thermally induced residual stresses is an ongoing challenge limiting the widespread adoption of WAAM technologies for producing components. In this manuscript, a novel approach is described to address thermal distortion in deposited components by applying lateral compressions along the length of the deposited material. To demonstrate this method, a series of single-track walls were printed and compressed at evenly spaced intervals using a modified hydraulic cutter tool. The jaws of the tool were modified to compress material rather than to shear it. A mathematical model was developed to relate the curvature of the deposited material to the volume of compression required to eliminate this distortion. Validation of this model was performed using 3D scan data to compare the change in wall curvature induced by compression to the volume of the applied compressions. Substrate deflection was also compared against a control wall, and implementation of wall compression reduced maximum deflections by 93% across a series of four depositions and subsequent compressions. Wall cross sections were also analyzed to determine the impact of compression on material hardness and grain structure. The results demonstrate that successively placed lateral compressions can effectively control and potentially eliminate bending distortion in printed parts. This methodology can be further developed to form a robust model for correction of thermally-induced distortion in WAAM components.

Additive manufacturing↗
DOE OSTI2025

A mixture parameterized biologically based dosimetry model to predict body burdens of polycyclic aromatic hydrocarbons in developmental zebrafish toxicity assays

Polycyclic aromatic hydrocarbons (PAHs) are a group of environmental toxicants found ubiquitously as complex mixtures in human-impacted environments. Developmental zebrafish exposures have been used widely to study PAH toxicity, but most studies report nominal exposure concentrations. Nominal exposure concentrations can be unreliable dose metrics due to differences in toxicant bioavailability resulting from disparate exposure methodologies and chemical properties. Toxicokinetic modeling can predict toxicant tissue doses to facilitate comparison between exposures of different chemicals, methodologies, and biological models. We parameterize a biologically based dosimetry model for developmental zebrafish toxicity assays for 9 PAHs. The model was optimized with measurements from media, tissue, and plastic plate walls throughout a static developmental exposure to a mixture of 10 PAHs of high abundance within the Portland Harbor Superfund Site. Plate binding, volatilization, zebrafish permeability, and tissue—media partitioning coefficients vary widely between PAHs. Model predictions accounted for 83% and 54% of 48 hpf body burdens within a factor of 2 resulting from exposures to mixtures and individual PAHs, respectively. Accounting for solubility significantly improves model performance. Competition for active sites in metabolizing enzymes may change biotransformation kinetics between individual PAH and mixture exposures. Area under the curve estimations of concentrations in zebrafish resulted in altered hazard rankings from nominal exposure concentrations. Future work will be oriented to generalizing the model to other PAHs. This PAH dosimetry model improves the interpretability of developmental zebrafish toxicity assays by providing time-resolved body burdens from nominal exposure concentrations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗