Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative capture”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Watershed zonation through hillslope clustering for tractably quantifying above- and below-ground watershed heterogeneity and functions

Abstract. In this study, we develop a watershed zonation approach for characterizing watershed organization and functions in a tractable manner by integrating multiple spatial data layers. We hypothesize that (1) a hillslope is an appropriate unit for capturing the watershed-scale heterogeneity of key bedrock-through-canopy properties and for quantifying the co-variability of these properties representing coupled ecohydrological and biogeochemical interactions, (2) remote sensing data layers and clustering methods can be used to identify watershed hillslope zones having the unique distributions of these properties relative to neighboring parcels, and (3) property suites associated with the identified zones can be used to understand zone-based functions, such as response to early snowmelt or drought and solute exports to the river. We demonstrate this concept using unsupervised clustering methods that synthesize airborne remote sensing data (lidar, hyperspectral, and electromagnetic surveys) along with satellite and streamflow data collected in the East River Watershed, Crested Butte, Colorado, USA. Results show that (1) we can define the scale of hillslopes at which the hillslope-averaged metrics can capture the majority of the overall variability in key properties (such as elevation, net potential annual radiation, and peak snow-water equivalent – SWE), (2) elevation and aspect are independent controls on plant and snow signatures, (3) near-surface bedrock electrical resistivity (top 20 m) and geological structures are significantly correlated with surface topography and plan species distribution, and (4) K-means, hierarchical clustering, and Gaussian mixture clustering methods generate similar zonation patterns across the watershed. Using independently collected data, we show that the identified zones provide information about zone-based watershed functions, including foresummer drought sensitivity and river nitrogen exports. The approach is expected to be applicable to other sites and generally useful for guiding the selection of hillslope-experiment locations and informing model parameterization.

58 GEOSCIENCES↗
DOE OSTI2022

Capture probe

A system and method for sampling a sample material includes a device for directing sample into a capture probe. The device for supplying sample material to the probe can be a device for radiating energy to the surface to eject sample from the sample material. A probe includes an outer probe housing having an open end. A liquid supply conduit has an outlet positioned to deliver liquid to the open end. An exhaust conduit removes liquid from the open end of the housing. The liquid supply conduit can be connectable to a liquid supply for delivering liquid at a first volumetric flow rate to the open end of the housing. A liquid exhaust system can be in fluid connection with the liquid exhaust conduit for removing liquid from the liquid exhaust conduit at a second volumetric flow rate such that gas with sample is withdrawn with the liquid.

Kertesz, Vilmos↗
DOE OSTI2024

Capture probe

A system and method for sampling a sample material includes a device for directing sample into a capture probe. The device for supplying sample material to the probe can be a device for radiating energy to the sample material to eject sample from the sample material. A probe includes an outer probe housing having an open end. A liquid supply conduit has an outlet positioned to deliver liquid to the open end. An exhaust conduit removes liquid from the open end of the housing. The liquid supply conduit can be connectable to a liquid supply for delivering liquid at a first volumetric flow rate to the open end of the housing. A liquid exhaust system can be in fluid connection with the liquid exhaust conduit for removing liquid from the liquid exhaust conduit at a second volumetric flow rate such that gas with sample is withdrawn with the liquid.

Kertesz, Vilmos↗
DOE OSTI2024

Desert-Urban SysTem IntegratEd AtmospherIc Monsoon (DUSTIEAIM) in the Southwestern United States

Phoenix, Arizona, is the fifth-largest city in the United States with about 1.6 million residents. It is also considered to be the hottest major city in the country, having set a record of 31 consecutive days with temperatures above 110 F (43 C) in summer 2023. Phoenix depends on the Colorado River for over a third of its water supply, but there are questions about how much water will eventually be available to serve the rapidly growing metropolitan area, especially in such extreme heat. Current earth system models struggle to accurately capture precipitation in the Southwestern United States, especially in relation to the summertime North American Monsoon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗
DOE OSTI2021

Capture efficiency and bias from the defect dynamics near grain boundaries in BCC Fe using mesoscale simulations

The capture efficiency describes the capability of a sink, such as a grain boundary (GB), dislocation, and void, to absorb point defects (PDs). The bias defines the difference in capture efficiency between the absorption of a vacancy and dumbbell at a sink. Complete kinetic information on PDs, including diffusion barriers and diffusion orientations, as well as accurate saddle points, are needed to determine the capture efficiency and bias at a sink accurately, which is computationally demanding. Here, in the present study, the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method was used to investigate the defect dynamics of PDs near different types of grain boundaries (GBs) (with both $\langle$100$\rangle$ and $\langle$110$\rangle$ families) accurately in body-centered cubic (BCC) iron (Fe). The capture efficiency, sink strength, and bias factor of different types of GBs were determined in Fe, which, different from traditional rate theory estimation, showed a distinct capture efficiency, sink strength, and bias in different GBs. The results demonstrate a strong positive correlation between the capture efficiency and the GB strain width, instead of the GB misorientation, GB energy, or GB-PD binding energy, which have been investigated previously. This work provides valuable insight into the radiation-induced microstructural evolution of GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Noninvasive liquid level sensing with laser generated ultrasonic waves

This article proposes a noninvasive liquid level sensing technique using laser-generated ultrasound waves for nuclear power plant applications. Liquid level sensors play an important role of managing the coolant system safely and stably in the plant structure. Current sensing techniques are mostly intrusive, performing inside the fluidic structure, which is disadvantageous in terms of the regular maintenance of the plant system. Furthermore, typical intrusive sensors do not perform stably under varying environmental conditions such as temperature and radiation. In this study, sensing units are attached to the outer surface of a liquid vessel to capture guided ultrasound waves in a nonintrusive manner. The signal intensity of the guided wave dissipates when the signal interacts with the internal liquid media. The sensing mechanism is mathematically expressed as an index value to correlate the liquid level with the sensor signal. For the acoustic wave generation, laser-generated ultrasound was adopted instead of using typical contact type transducers. Following the simulation validation of the proposed concept, the performance of the developed sensor was confirmed through experimental results under elevated liquid temperature conditions. Finally, the nonlinear multivariable regression exhibited the best-fit to the datasets measured under the variable liquid level and temperature conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗
DOE OSTI2026

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗
DOE OSTI2022

Integration of the Kromek D3S Detector and Spot Robot For Secondary Inspections

Inspecting vehicles and containers for the presence of nuclear material is a challenging task for border control and security. When performed manually by inspectors, this task also has an associated risk of exposing the inspectors to unknown radiation. With the advent of agile, easy-to-program, quadruped robots like the Boston Dynamics Spot, automation of secondary inspection can improve the efficiency of the inspection process and alleviates the radiation risks to inspectors. In this project, Brookhaven National Laboratory and the University of Massachussetts at Lowell explored how to automate a simple secondary inspection mission. The Spot robot comes with its own software development kit (SDK) that allows clients/users to write custom code in the Python programming language to control the robot. Spot also has a payload computer called Spot-CORE, which runs the Ubuntu Linux operating system and allows users to integrate external sensors, such as a radiation detector, with Spot. In this study, the Kromek D3S detector has been integrated with Spot via the Spot-CORE, allowing Spot to capture gamma spectra and neutron counts for a specified acquisition period. Two custom routines, search and confirmation, have been developed and executed in this specified order. The search routine directs Spot to go around the nearest obstacle, e.g., vehicle and container, in a preset distance and step to collect gamma and neutron gross counts with the D3S detector. The radiation data and the robot location corresponding to each step are stored and fed to the confirmation routine at the end of the search. The confirmation routine then navigates Spot to the locations of the highest gamma or neutron counts to perform a long, e.g., one minute, measurement, and gives the operators the signature gamma spectra and neutron counts at the hotspots. This paper presents a detailed description of this automated system along with results of the preliminary tests in identifying the location and signature of a 137Cs radiation source in a vehicle.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Performance and flexibility improvements of Staged Pressurized Oxy-Combustion

Staged Pressurized Oxy-Combustion (SPOC) is a low-carbon coal combustion power technology being developed by Washington University in St. Louis (WUSTL). Oxy-combustion plants enable straightforward capture of carbon dioxide (CO 2 ) by removing most of the nitrogen in the combustion air prior to use, thereby burning fuel in near-pure oxygen instead of air, producing a flue gas containing primarily CO 2 and water. CO 2 capture at amounts > 90% is possible, often using cryogenic air separation. Oxy-combustion typically relies on flue gas recycle (FGR) to reduce the peak temperature and radiation that would otherwise occur in a fuel/oxygen only flame. SPOC reduces the peak temperatures of combustion by utilizing two or more pressurized boiler modules connected in series to produce fuel staging; hence, only a portion of the fuel is combusted in any given furnace module. This means that the thermal energy released at each stage can be captured and removed from the gases prior to subsequent stages, when more fuel is introduced. This allows the SPOC process to operate with minimal FGR, avoiding the associated efficiency losses and additional costs. Also, the process operates at an elevated gas-side pressure, reducing boiler size, enhancing heat transfer to achieve a compact boiler configuration as compared to an atmospheric-pressure boiler design, and allowing for recovery of the latent heat of the water from the flue gas at a temperature useful to the steam cycle. The resultant net efficiency of the system is over 3 percentage points greater than traditional atmospheric-pressure oxy-combustion, and 7 percentage points greater than the post combustion variant, representing a step-change improvement over first-generation capture technologies. To further develop the concept, WUSTL and the Electric Power Research Institute, Inc., organized a project with American Air Liquide, Inc., Doosan Babcock Limited, and the U.S. Department of Energy to investigate a practicable and workable boiler design. The team has identified the potential for enhanced process flexibility for controlling power generation over a wider load range than is normally available to conventional coal-fired power plants due to the staged nature of the heat release. With increasing intermittent renewable generator contribution, on-demand generators need to be highly flexible to participate in the future energy market, requiring extensive operation at reduced load. Conventional coal-fired steam generators typically face challenges in maintaining temperature control of the reheat steam and main steam at reduced loads. This results in inefficient operation, both in terms of the boiler efficiency and steam turbine heat rate. The results of this project show the SPOC process is capable of exceptional turndown, both on a stage basis and with the ability to bypass entire stages. Oxygen-supply flexibility was also investigated, as this is also a key consideration for the overall flexibility of the SPOC process given the operating constraints of conventional air separation units. A boiler design concept assessment was conducted and was focused on delivering compact and constructible design. The assessment checked appropriate tube operating metal temperatures at full load and at lower operating loads, balanced against the needs of efficient coal combustion, and the resultant slagging and ash environments. Combustion testing in the 100-kWth pressurized combustion test rig at WUSTL was carried out to validate the combustion, heat flux profiles and burnout at multiple loads. Combustion parameters investigated were flame stability, fuel burnout, ash composition, radiative heat flux, and temperature profiles. The results of these tests formed the basis of a full-scale boiler design that will encompass improvements in both efficiency and flexibility over conventional oxy-combustion processes. The air separation unit flexibility was investigated, and associated cost implications were addressed. Detailed economic assessment results for a 550 MWe net power block are also provided, allowing for a comparison against the baseline NETL oxy-combustion and post combustion capture cases.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kilauea 2018

All instrumented basaltic caldera collapses have generated M w > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. Here, we present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kilauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward-moving caldera block and magma column, and post-collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near-field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Positron Emission Intensity in the Decay of 86 gY for Use in Dosimetry Studies

The β + -emitting radionuclide 86 gY (t 1/2 = 14.7 h) forms a matched-pair with the β – -emitting therapeutic radionuclide 90Y (t 1/2 = 2.7 d) for theranostic application in medicine. This approach demands a precise knowledge of the positron emission probability of the PET nuclide which was until recently rather uncertain for 86 gY. In this work, an 86 gY source of high radionuclidic purity was prepared and a direct measurement of the positron emission intensity per 100 decay of the parent (hereafter “positron emission intensity”) was performed using high-resolution HPGe detector γ-ray spectroscopy. The electron capture intensity was also determined as an additional check by measuring the K α and K β X-rays of energies 14.1 and 15.8 keV, respectively, using a low energy HPGe detector. From those measurements, normalized values of 27.2 ± 2.0% for β + -emission and 72.8 ± 2.0% for EC were obtained. These results are in excellent agreement with values recently reported in the literature based on a detailed decay scheme study.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Evaluating Mesoscale Convective Systems Over the US in Conventional and Multiscale Modeling Framework Configurations of E3SMv1

Organized mesoscale convective systems (MCSs) contribute a significant amount of precipitation in the Central and Eastern US during spring and summer, which impacts the availability of freshwater and flooding events. However, current global Earth system models cannot capture MCSs well and misrepresent the statistics of precipitation in the region. In this study, we investigate the representation of MCSs in three configurations of the Energy Exascale Earth System Model (E3SMv1) by tracking individual storms based on outgoing longwave radiation using a new application of TempestExtremes. Our results indicate that conventional parameterizations of convection, implemented in both low (LR; ~150 km) and high (HR; ~25 km) resolution configurations, fail to capture almost all MCS-like events, in-part because they underestimate high-level cloud ice associated with deep convection. On the other hand, the multiscale modeling framework (MMF; cloud-resolving models embedded in each grid-column of ~150 km resolution E3SMv1) configuration represents MCSs and their annual cycle better. Nevertheless, relative to observations, the E3SMv1-MMF spatial distribution of MCSs and associated precipitation is shifted eastward, and the diurnal timing is lagged. A comparison between the large-scale environment in E3SMv1-MMF and ERA5 reanalysis suggests that the biases during the summer in E3SMv1-MMF are associated with biases in low-level humidity and meridional moisture transport within the low-level jet. The fact that conventional parameterizations of convection, even with high-resolution, cannot capture MCSs over the US suggests that methods with explicit representation of kilometer-scale convective organization, such as the MMF, may be necessary for improving the simulation of these convective systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A laser powder bed fusion system for operando synchrotron x-ray imaging and correlative diagnostic experiments at the Stanford Synchrotron Radiation Lightsource

Laser powder bed fusion (LPBF) is a highly dynamic multi-physics process used for the additive manufacturing (AM) of metal components. Improving process understanding and validating predictive computational models require high-fidelity diagnostics capable of capturing data in challenging environments. Synchrotron x-ray techniques play a vital role in the validation process as they are the only in situ diagnostic capable of imaging sub-surface melt pool dynamics and microstructure evolution during LPBF-AM. In this article, a laboratory scale system designed to mimic LPBF process conditions while operating at a synchrotron facility is described. The system is implemented with process accurate atmospheric conditions, including an air knife for active vapor plume removal. Significantly, the chamber also incorporates a diagnostic sensor suite that monitors emitted optical, acoustic, and electronic signals during laser processing with coincident x-ray imaging. The addition of the sensor suite enables validation of these industrially compatible single point sensors by detecting pore formation and spatter events and directly correlating the events with changes in the detected signal. Experiments in the Ti–6Al–4V alloy performed at the Stanford Synchrotron Radiation Lightsource using the system are detailed with sufficient sampling rates to probe melt pool dynamics. X-ray imaging captures melt pool dynamics at frame rates of 20 kHz with a 2 µm pixel resolution, and the coincident diagnostic sensor data are recorded at 470 kHz. Furthermore, this work shows that the current system enables the in situ detection of defects during the LPBF process and permits direct correlation of diagnostic signatures at the exact time of defect formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multigroup Scattering in Monte Carlo Radiation Transport Codes [Slides]

Legendre truncations to multigroup scattering distributions are not amenable to Monte Carlo sampling due to negative values. We have implemented two moment-preserving methods in MGMC that capture the shape of the truncation, are non-negative over [-1, 1], and can be efficiently sampled on CPUs and GPUs. MGMC can now simulate neutrons with anisotropic scattering mechanics and MGMC shows good agreement with LANL production codes PARTISN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗