Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal diffusion”

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

Thermal stability of 3D interface Cu/Nb nanolaminates

Nanocrystalline alloys are promising structural materials yet lack thermal stability in many cases. Recent work shows that interface structure has an outsize effect on the thermal behavior of nanostructured alloys. Here, this work focuses on the role of controlled heterophase interface structure in the thermal evolution of model Cu/Nb nanolaminates. We introduce 3D interfaces containing nanoscale heterogeneities in all spatial dimensions between Cu and Nb, forming 3D Cu/Nb. TEM, nanoindentation, and DSC are used in tandem to establish thermal stability and to identify shifts in microstructure as a function of static annealing temperature. 3D interfaces are shown to survive annealing to 300 °C for 1 hr., while 3D Cu/Nb microstructure evolves to form low-density and voided regions correlating to the onset of layer pinch-off between 500 and 600 °C annealing temperatures. A diffusivity- and vacancy energetics-based mechanism is developed to explain void formation driven by 3D interface degradation at elevated temperature.

3D interfaces↗
DOE OSTI2023

Multi-Task Learning of Scanning Electron Microscopy and Synthetic Thermal Tomography Images for Detection of Defects in Additively Manufactured Metals

One of the key challenges in laser powder bed fusion (LPBF) additive manufacturing of metals is the appearance of microscopic pores in 3D-printed metallic structures. Quality control in LPBF can be accomplished with non-destructive imaging of the actual 3D-printed structures. Thermal tomography (TT) is a promising non-contact, non-destructive imaging method, which allows for the visualization of subsurface defects in arbitrary-sized metallic structures. However, because imaging is based on heat diffusion, TT images suffer from blurring, which increases with depth. We have been investigating the enhancement of TT imaging capability using machine learning. In this work, we introduce a novel multi-task learning (MTL) approach, which simultaneously performs the classification of synthetic TT images, and segmentation of experimental scanning electron microscopy (SEM) images. Synthetic TT images are obtained from computer simulations of metallic structures with subsurface elliptical-shaped defects, while experimental SEM images are obtained from imaging of LPBF-printed stainless-steel coupons. MTL network is implemented as a shared U-net encoder between the classification and the segmentation tasks. Results of this study show that the MTL network performs better in both the classification of synthetic TT images and the segmentation of SEM images tasks, as compared to the conventional approach when the individual tasks are performed independently of each other.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature dependence of heavy quark diffusion from (2+1)-flavor lattice QCD

We present a lattice determination of the heavy-quark diffusion coefficient in (2+1)-flavor QCD with almost physical quark masses. The momentum and spatial diffusion coefficients are extracted for a wide temperature range, from T = 163 MeV to 10 GeV. The results are in agreement with previous works from the HotQCD collaboration, and show fast thermalization of the heavy quark inside the QGP. Near the chiral crossover temperature T c ≃ 150 MeV, our results are close to the AdS/CFT estimation computed at strong coupling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stationary Oxygen Vacancy Construction toward a Superior-Performance Ultrahigh Nickel Single-Crystal Cathode

Oxygen vacancies exert a complex and profound influence on the layered cathodes, especially those with ultrahigh nickel content. They can facilitate lithium-ion transport and enhance electronic conductivity, while aggressive oxygen vacancy formation causes structural degradation and electrolyte decomposition. Herein, taking ultrahigh nickel single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92) as a model material, we propose a pinning strategy to harness the benefits of oxygen vacancies while mitigating their detrimental effects. Through a carefully controlled thermal process, both oxygen vacancies and pinning atoms are successfully introduced into the surface region. The resulting anchored oxygen vacancies, capitalizing on their inherent advantages, improve conductivity and lithium-ion diffusion. Simultaneously, the neighboring pinning atoms effectively increase the migration barrier and suppress the adverse effects of these vacancies, including electrolyte decomposition and structural degradation during long-term electrochemical cycling. Consequently, oxygen vacancy-anchored single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92-OV) demonstrates significantly improved high-voltage electrochemical performance, with 86.16% capacity retention after 200 cycles at 4.6 V and 1 C in a half-cell and 90.71% after 300 cycles at 4.5 V and 1 C in a full cell. Furthermore, this study not only provides valuable insights into the chemistry of oxygen vacancy but also introduces a viable strategy for leveraging oxygen vacancies to achieve stable high-voltage performance in ultrahigh nickel single-crystal cathodes.

defects in solids↗
DOE OSTI2026

Extreme Temperature Cryptography Based On Nitrogen-Incorporated Ultrananocrystalline Diamond

Physical entropy sources that remain stable under extreme temperatures are essential for cryptography in emerging technological frontiers in deep space exploration, geothermal energy harvesting, and nuclear energy. However, conventional semiconductor platforms fail to generate stable and reliable cryptographic keys above 200 degrees C due to performance degradation. Here, we report a diamond-based cryptographic primitive that exploits the defect-rich sp 2 -bonded grain boundary network in nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) film as a robust entropy source to generate cryptographic keys that remain operationally stable even after enduring extreme temperatures of 700 degrees C for 54 h while also surviving thermal cycling between room temperature and 700 degrees C for 48 h. The strength of the generated keys is assessed through several cryptographic metrics such as bit uniformity, entropy, hamming distances, and correlation coefficients, all of which are found to be near their respective ideal values. Moreover, the generated keys pass the NIST SP 800 and SP 800-90B tests and are also resilient to supply bias variations and a regression-based machine learning attack model based on the Fourier series. The robustness of the keys is attributed to the better thermal stability and chemical inertness of the n-UNCD film. This is supported by high-resolution energy-dispersive X-ray spectroscopy (EDS), which shows no significant lateral diffusion of metal atoms into the n-UNCD layer, and by Raman spectroscopy, which reveals no significant changes in the bonding configuration of the n-UNCD structure. Our findings highlight the remarkable potential of n-UNCD film for extreme environment cryptography by expanding the operational limits of conventional hardware security platforms.

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

DEVELOPMENT OF STABLE ANODE AND CATHODE MATERIALS FOR RECHARGEABLE BATTERIES

Batteries have been employed in a variety of applications, such as portable electronics, electric vehicles (EVs), and stationary energy storage to preserve energy from other renewable sources (like wind or solar energy). The ultimate goal is to develop high energy density, long life span, better safety, and low cost of the batteries. However, current commercialized batteries (like lead-acid, zinc-alkaline, lithium-ion batteries (LIBs)) could not fulfill all the demands of diversified applications. LIBs predominate the market because of their high energy density. To achieve a higher capacity of the cell, high-nickel layered (Ni > 90%) cathode materials are promising candidates since they compose high specific capacity and discharge voltage. In chapter 1, an introduction to cell energy density and the development of high-nickel layered cathode materials along with associated obstacles of poor cycling stability and thermal stability have been discussed. Associated works like doping or surface coating have also been mentioned in this section. In chapter 2, Al doping in high-nickel layered cathode materials to enhance the structural and thermal stability was introduced. Uniform incorporation of Al doping is achieved by mechanical fusion and calcination processes. The Al doping not only decreased the Li/Ni mixing ratio but also enhance the thermal resistance to oxygen evolution because of strong Al-O bonding, which leads to elevated electrochemical and thermal stability. In chapter 3, TiN is implemented as a Ti dopant for high-nickel layered cathode materials to improve the electrochemical performance in the LIBs. The Ti not only diffused within the bulk structure but also formed segregation on the surface, which improved the structural stability and led to better cycling performance. Although LIBs deliver high energy density, safety concerns of flammable organic electrolytes have not been resolved yet. Therefore, the aqueous rechargeable zinc-ion batteries (ZIB) have been praised for their safe, low-cost, eco-friendly stationary energy storage, which is considered as a complementary system to LIBs. In chapter 4, the mechanism, advantages, and challenges of ZIBs would be given and the strategies for solving the zinc metal anode issues have been discussed in this section. In chapter 5, a polymer coating method was reported to facilitate Zn deposition/ stripping by coordination with Zn2+ and prevented direct contact with aqueous electrolyte to block corrosion side-reactions. With this coating, a boost in a lifetime (up to 400 hours) and lower polarization under extremely high current conditions (10 mA cm-2) have been achieved in repeated Zn deposition/ stripping cycling.

25 ENERGY STORAGE↗
DOE OSTI2020

Applications of Hyporheic Zone Sampling for Assessing VOC Plume Engineered Remediation and Natural Attenuation Effectiveness - 20414

Determining the efficacy of engineered remedial actions and/or natural attenuation as an acceptable EPA Record of Decision (ROD) must address long-term temporal and spatial variability. Since final plume downslope destinations are aquifers and related surface streams and rivers, hyporheic zone sampling is essential. Comprehensive hyporheic sampling documents plume groundwater remediation prior to upwelling into stream benthic strata and natural attenuation within this dynamic zone itself. Our intensive sampling of the Savannah River Site (SRS) Chemical, Metals, and Pesticides (CMP) Pits VOC plume-fringe from 2005 to 2018 documented PCE natural attenuation and effectiveness of In Situ Thermal Treatment (ISTT). Most plume hyporheic upwelling was delimited to a 20 m inflow reach along Pen Branch. Our study utilized corroborative sampling methodologies: ambient hole-water grab samples: passive diffusion bags (PDB); and hole-core sediment samples to confirm station to station spatial and temporal trends. VOC hot spot detections were confirmed from both groundwater and soil core samples. VOC pathways through the Pen Branch hyporheic zone were assessed at hole-installation using sequential soil samples from augered cores and grab-samples of hole-bottom water. Longer term, composite water samples were collected from PDBs deployed for a minimum two-week equilibration period. Soil samples were collected from the 15 and 65 cm core levels in December 2017 and June 2018 at six stations. Additional core depths of 40 and 85 cm were sampled in June 2018 at the highest VOC station (5DB80). All samples were analyzed in an EPA certified lab by purge and trap GC-MS. Effectiveness of previous ISTT and soil vapor extraction (SVE) for this CMP Pits PCE plume, was clearly documented by significant downward trends in VOC at hyporheic zone stations prior to December 2015. Achieving final stage conversions of PCE daughter products to ethene and CO{sub 2} may be limited depending upon suitable interactions between hyporheic zone sediment morphology, hydrology, and redox conditions. Because declines in VOC daughter products at our hyporheic stations may also be caused by physical actions of dilution and dispersion rather than true biochemical dechlorination, full quantification of natural attenuation requires the use of compound-specific isotope ratio analysis (CSIA) and microbial 16S rRNA with qPCR. To assist these next-phase efforts, a natural attenuation index (NAI) was used to help produce quicker, more cost-effective decision making for locating CSIA and microbial sampling efforts. The NAI used natural logarithm transformed GC-MS results for station VOC degradation compounds versus VOC parent compounds. Wide differences were observed between stations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electron-phonon coupling and superconductivity in the doped topological crystalline insulator ( Pb 0.5 Sn 0.5 ) 1 – x In x Te

Here, we present a neutron-scattering study of phonons in single crystals of ( Pb 0.5 Sn 0.5 ) 1 – x In x Te with x = 0 (metallic, but nonsuperconducting) and x = 0.2 (nonmetallic normal state, but superconducting). We map the phonon dispersions (more completely for x = 0 ) and find general consistency with theoretical calculations, except for the transverse and longitudinal optical (LO) modes at the Brillouin-zone center. At low temperature, both modes are strongly damped but sit at a finite energy ( ≈ 4 meV in both samples), shifting to higher energy at room temperature. These modes are soft due to a proximate structural instability driven by the sensitivity of Pb-Te and Sn-Te p -orbital hybridization to off-center displacements of the metal atoms. The impact of the soft optical modes on the low-energy acoustic modes is inferred from the low thermal conductivity, especially at low temperature. Given that the strongest electron-phonon coupling is predicted for the LO mode, which should be similar for both studied compositions, it is intriguing that only the In-doped crystal is superconducting. In addition, we observe elastic diffuse (Huang) scattering that is qualitatively explained by the difference in Pb-Te and Sn-Te bond lengths within the lattice of randomly distributed Pb and Sn sites. We also confirm the presence of anomalous diffuse low-energy atomic vibrations that we speculatively attribute to local fluctuations of individual Pb atoms between off-center sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗
DOE OSTI2020

Determination of Stress-Free Lattice Spacing (D₀) for Residual Stress Relaxation Measurement in Ni-based Superalloys by Neutron Diffraction

The residual stress and its relaxation in critical engineering components are key to structural materials reliability. The determination of residual stress in bulk engineering structures like Ni-based superalloys can be conducted by measuring the lattice strain with penetrating neutron diffraction. The state-of-the-art engineering diffraction beamline VULCAN at the Oak Ridge National Laboratory (ORNL) Spallation Neutron Source (SNS) provides the critical capability to evaluate the residual strain/stress relaxation in engineering component thanks to the high flux and event data features. By the means of residual stress/strain calculation, it is critical to accurately determine the stress-free lattice spacing d₀, which can be altered by the change of chemistry, unrelieved stress, and measurement scheme. Here we reported d₀ determination for determining residual stress relaxation and distribution in Alloy 718 superalloys after different quenching treatments. Selection of locations for d₀ measurement was discussed by considering the neutron path, attenuation and sample alignment. A dynamic d₀ that resulted from atom diffusion and chemistry change was estimated as function of time for the in-situ relaxation characterization. It demonstrated how the dynamic d₀ values may influence on the strain calculation in different thermally treated samples. The results highlight the importance of dynamic d₀ for in-situ relaxation involving chemical changes and provide guidance to the dynamic d₀ measurement and simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗
DOE OSTI2023

Post-irradiation examination of the Sirius-1 nuclear thermal propulsion fuel test

Nuclear Thermal Propulsion (NTP) systems hold promise in reducing transit times for exploration of the solar system by both crewed and uncrewed missions. NTP systems currently under investigation include a once-through high-temperature gas-cooled fission reactor to provide thermal energy to heat the coolant, which also serves as the propellant. The fuel system consists of angular UN fuel particles dispersed in a matrix of W/Re, creating a ceramic and metallic composite or cermet, which has been irradiated in Idaho National Laboratory's (INL) Transient Reactor Test Facility (TREAT). This enabled evaluation of these materials under representative nuclear heating rates (~95 K/s) and peak surface temperatures (~2527 K). Further, these surface conditions were indicative of achievement of the target peak internal temperature of approximately 2850 K. These tests named Sirius-1 (UN–W/Re), have been irradiated and this paper will present post-irradiation examination results. The Sirius-1 test produced minor, stable cracks in the fuel specimen and spalling of surface material. Volatilized uranium ‘soot’ was found deposited on the inner wall of the irradiation capsule indicating loss of some fissile material from the fuel specimen. Spalling from the surfaces was also noted upon visual inspection. Uranium diffusion from the fuel particles resulted in the formation of U/Re phases and the production of a laminar microstructure on the edge of the fuel system. Overall, the specimen was stable and remained within a coolable geometry upon conclusion of multiple thermal cycles to prototypical operating temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗
DOE OSTI2024

High Pressure Melting Curve of Fe‐Si: Implication for the Thermal Properties in Mercury's Core

The motion of liquid iron (Fe) alloy materials in the outer core drives the dynamo, which generates Mercury's magnetic field. The assessment of core models requires laboratory measurements of the melting temperature of Fe alloys at high pressure. Here, we experimentally determined the melting curve of Fe9wt%Si and Fe17wt%Si up to 17 GPa using in situ and ex situ measurements of intermetallic fast diffusion that serves as the melting criterion in a large-volume press. Our determined melting slopes are comparable with previous studies up to about 17 GPa. However, when extrapolated, our melting slopes significantly deviate from previous studies at higher pressures. For Mercury's core with a model composition of Fe9wt%Si, the melting temperature-depth profile determined in our study is lower by ∼150–250 K when compared with theoretical calculations. Using the new melting curve of Fe9wt%Si and the electrical resistivity values from a previous study of Fe8.5wt%Si, we estimate that the electronic thermal conductivity of liquid Fe9wt%Si is 30 Wm −1 K −1 at the Mercury's CMB pressure of 5 GPa and 37 Wm −1 K −1 at an assumed ICB of 21 GPa, corresponding to heat flux values of 23 mWm −2 and 32 mWm −2 , respectively. These values provide new constraints on the core models.

58 GEOSCIENCES↗
DOE OSTI2020

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation is independent of the number of spatial points and the number of processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗