PROTON-PULSE INDUCED DYNAMIC STRAIN PREDICTION, MEASUREMENT, AND SIMULATION VALIDATION OF SNS TARGET WITH A NOVEL GAS INJECTOR
Explore the source record for details and available documents.
Jiang, Hao [ORNL] (ORCID:0000000232929338)↗
SEARCH · Search NASA
Results for “TARGET SIMULATION”
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.
2D kinetic-ion simulations of inverted corona fusion targets
Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Assessment of simulated and observed cavitation-induced erosion damage in Spallation Neutron Source target vessels
Cavitation-induced erosion damage in different Spallation Neutron Source (SNS) target designs are simulated using explicit finite element–based techniques and compared with observations of erosion in targets after operation. The efficacy of the previously developed method, called saturation time, was evaluated using erosion-damaged samples from new target designs. A new metric called maximum bubble size was implemented under the rationale that larger cavitation bubbles will collapse more intensely. The maximum cavitation bubble size over 1 ms of simulated time was calculated based on the Rayleigh–Plesset equation for each element integration point and presented as a contour map at the vessel surface for assessing with erosion observations. SNS targets are now operated with helium gas injection to reduce cavitation damage. A simulation method using a material model for the mixture of mercury and gas bubbles was recently developed and used to account for the effect of small gas bubbles on the structural response of the target vessel. Furthermore, this work compares the new method's results with observed cavitation damage. Maps of the calculated maximum bubble size for targets operated with and without gas injection were compared with photographs of erosion damage observed in SNS targets. The patterns in maximum bubble size maps correlated well with observations of erosion patterns in target vessels after service. Advantages and challenges of the maximum bubble size simulation technique are provided, and differences between results from the previous and the newly proposed metric are discussed.
Jiang, Hao↗
Quantification and evaluation of strain reduction from small-bubble gas injection in Spallation Neutron Source target vessels
Small-bubble gas injection has been routinely utilized in the operation of Spallation Neutron Source (SNS) mercury targets since 2017 to mitigate cavitation-induced erosion damage to target vessels. Strain measurements of target vessels collected in-situ during initial operation with gas injection were used to study the gas injection effect on the structural response of targets to proton pulses. A significant strain reduction owing to gas injection was found by comparing the strain measurement data during operation with and without gas injection. The research presented here focuses on quantifying strain reductions in SNS targets and evaluating the effect of small-bubble gas injection by comparing different bubbler types and target designs. The strain measurement results show the gas injection significantly reduced strain in SNS target vessels; strain values decreased by 30% to 80% for targets operating with gas injection. Stress and strain responses of SNS targets were simulated to numerically evaluate the gas injection effect. Based on the predicted stresses with and without gas injection, the fatigue lifetimes of SNS jet-flow design target were estimated using fe-safe fatigue analysis software. The simulations show these reductions should improve the fatigue life of target vessels and allow SNS targets to meet their fatigue design goal.
Fatigue life estimation↗
Multiphysics Simulations of Nanofibrous High Power Targets for the Inference of Failure Modes
Fermilab s High Power Targetry Research and Development (HPT R&D) group have been developing and studying an electrospun nanofiber target concept to support the need for robust targets in future fixed target facilities. These nanofiber mats have demonstrated resistance to radiation damage, and the free motion of the individual fibers is expected to mitigate the cyclic stresses induced by a pulsed, high-power beam. To evaluate the efficacy of this concept, nanofiber mat samples have been sent by the HPT R&D group to the HiRadMat facility at CERN for prototypic thermal shock testing; the outcomes of these experiments show that the survivability of a nanofiber target depends on its construction parameters, in particular the packing density of the fibers. Samples with higher packing densities have consistently been destroyed by exposure to the HiRadMat beam, with a hole at the beam center visible, and layers of nanofibers peeled away from the center hole, whereas samples with lower densities have survived with limited damage. The exact reason for the failure of the higher density targets was unclear at the time of the original experiments, but the results of our recent multiphysics simulations which recreate the experiments support the hypothesis that the expansion and pressurization of the air inside the target after being heated by the pulsed beam is the cause; in a high density nanofiber mat, the motion of air through the pores of the mat is much more restricted, and induces a larger pressure on the fibers, blowing the mat apart. In this talk, we ll share the results of these simulations and discuss how they support this hypothesis.
Asztalos, Will↗
Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF
The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.
Bremsstrahlung↗
Simulations of positron injector for Ce+BAF
A baseline concept for a continuous wave (CW) polarized positron injector was developed for the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. This concept is based on the generation of CW longitudinally polarized positrons by a high-current, polarized electron beam (1 mA, 130‑370 MeV, and 90% longitudinal polarization) that passes through a rotating, water-cooled, tungsten target. The simulation results for the Ce+BAF injector at the Low Energy Recirculator Facility (LERF) are presented, including positron beam generation, capture, energy selection, and acceleration to 123 MeV. The positron yield (or positron current) and longitudinal polarization are calculated considering the longitudinal and transverse CEBAF acceptances (<1% energy spread, <1 mm bunch length and normalized emittance of <100 mm mrad). The impact of target thickness, drive electron beam energy, and transverse size on positron yield within the required emittance limit is evaluated.
Accelerator Physics↗
High-Latitude, Low-Altitude SAI: Overview of G6-1.5K-HiLLA Simulations in E3SMv3
This report details a contribution of Energy Exascale Earth System Model version 3 (E3SMv3) simulations to a proposed model intercomparison project funded and organized by Reflective, a nonprofit group studying the possible global impacts of SAI. Using annual feedback control and seasonally-variable injection sites to achieve a desired global near-surface temperature target, the simulated SAI campaign indicates robust maintenance of the 2020-2039 climatic state over 60 years into the future.
54 ENVIRONMENTAL SCIENCES↗
Positron Injector for Ce+BAF
A baseline concept for a continuous wave (CW) polarized positron injector was developed for the CEBAF at Jefferson Lab. This concept is based on the generation of CW longitudinally polarized positrons by a high-current, polarized electron beam 1 mA, 130-370 MeV, and 90% longitudinal polarization that passes through a rotating, water-cooled, tungsten target. The simulations for the Ce+BAF injector at the Low Energy Recirculator Facility (LERF) are performed, including positron beam generation, capture, energy selection, and acceleration to 123 MeV. The positron yield (or positron current) and longitudinal polarization are calculated considering the longitudinal and transverse CEBAF acceptances (<1% energy spread, 1.2 mm bunch length and normalized emittance of 100 mm·mrad). The impact of target thickness, electron beam energy, and transverse size on positron yield within the required emittance limit is evaluated. Target thermal and structural FEM analyses (ANSYS Fluent) are performed to determine the maximum electron beam current and minimum transverse beam size on the target.
Benesch, Jay [Thomas Jefferson National Accelerato↗
PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas
We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST
ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.
17 WIND ENERGY↗
Foundation models for atomistic simulation of chemistry and materials
Conventional computational methods for modeling chemical and materials systems are limited by system size and timescale, forcing a trade-off between quantum-mechanical accuracy and the sampling needed for realistic observables. Large language and vision foundation models — pre-trained on massive datasets using transformer architectures — have revolutionized many fields. It is thus interesting to ask whether a foundation model — subject to suitable data, parameter scaling and training — could enable learned simulations of chemistry and materials. Here, in this study, we review the field of machine-learned interatomic potentials (MLIPs) and posit that scaling up large and diverse chemical and materials datasets and highly expressive architectures using advanced training strategies should result in models that are: more efficient, transferable, robust to out-of-distribution scenarios, and easier to fine-tune to a variety of downstream physical observables than models trained from scratch on small datasets corresponding to specific, targeted atomistic simulation tasks. We provide specific criteria for creating such large-scale MLIP foundation models, coordinated strategies for their development, evaluation and deployment, and highlight potential emergent capabilities that could transform predictive simulations in chemistry and materials science and accelerate discovery across multiple technological domains.
Yuan, Eric C.-Y. [University of California, Berkel↗
Structural and Dynamical Insights into the Formation Process of a Cross-Linked Polymer Network in Acrylic Adhesives During Thermal Curing
Many modern adhesives, sealants, and coatings rely on the controlled transition from a liquid to a solid state by forming a three-dimensional cross-linked polymer network, often referred to as curing. The curing process, which is initiated by the mixing of reactive components or an external trigger, defines the structure of a network and further controls the final mechanical properties of cured materials. However, the curing mechanism is not fully understood yet due to the lack of experimental tools capable of directly probing the structure and dynamics of a network over relevant time- and length scales. Here, in this paper, we report the curing process of a commercial two-component methyl methacrylate (MMA) adhesive using in operando X-ray photon correlation spectroscopy (XPCS), a method that closely simulates the target manufacturing environment of the adhesive. The results are then integrated with those obtained by rheology, differential scanning calorimetry (DSC), and transmission electron microscopy to establish the structure–dynamics–process–property relationship. The XPCS results identify four distinct stages in the curing process after the mixing, extrusion, and deposition of the acrylic adhesive: (i) At a cure time (or “aging time”, t age ) of less than 1 min, nanodomains of polymerized MMA are formed within a liquid monomeric MMA matrix. The average size is several nm and remains constant over t age , while the dynamics of the nanodomains are slowed down with t age due to an increase in the viscosity of the MMA matrix. (ii) After t age > 1 min, the size of the nanodomains increases with t age until the gel point (= 6.3 min after mixing as determined by rheology). The dynamics of the nanodomains also increase due to the heat generated by the exothermic reaction. (iii) At the gel point, the nanodomains begin to interconnect each other, resulting in a network structure with a characteristic length of about 100 nm. This characteristic network size does not change for the rest of the curing process up to t age = 500 min. The dynamics of the network structure, however, show a rapid slowing down with t age up to t age ≈ 12 min, corresponding to the onset of vitrification (as determined by rheology). (iv) At t age > 12 min, when the DSC and rheology data can no longer provide meaningful information, the XPCS data show a further slowing down of the network dynamics associated with vitrification. Our results provide rich and complex insights into the physics and material design of thermosets in commercially relevant processes, which are essential for future industrial applications.
36 MATERIALS SCIENCE↗
On the role of self-generated magnetic fields in focused ion beams driven by intense picosecond laser pulses
We investigate the physics of ion beam focusing in the interaction of intense picosecond laser pulses with concave foils, using two-dimensional particle-in-cell simulations. Curved targets have been observed to deliver brighter multi-MeV ion beams than flat targets, and they could be used to ignite pre-compressed ICF fusion capsules in a scheme called ion fast ignition. Focusing is achieved by shaping the target and relying on an ion acceleration mechanism in which the ions are accelerated normal to the target surface; subsequently, the plasma expansion grows magnetic fields that advect with the plasma and deflect ions away from the target symmetry axis, causing the focal spot to shift downstream from its geometric location. The drop in both electric and magnetic field strengths in the expanding plasma makes the ion beam envelope asymmetric with respect to the focusing plane.
Focused ion beam↗
First study of neutrino angle reconstruction using quasielasticlike interactions in MicroBooNE
We investigate the expected precision of the reconstructed neutrino direction using a 𝜈 𝜇 -argon quasielasticlike event topology with one muon and one proton in the final state and the reconstruction capabilities of the MicroBooNE liquid argon time projection chamber. This direction is of importance in the context of DUNE sub-GeV atmospheric oscillation studies. MicroBooNE allows for a data-driven quantification of this resolution by investigating the deviation of the reconstructed muon-proton system orientation with respect to the well-known direction of neutrinos originating from the Booster Neutrino Beam with an exposure of 1.3 ×10 21 protons on target. Using simulation studies, we derive the expected sub-GeV DUNE atmospheric-neutrino reconstructed simulated spectrum by developing a reweighting scheme as a function of the true neutrino energy. We further report flux-integrated single- and double-differential cross section measurements of charged-current 𝜈 𝜇 quasielasticlike scattering on argon as a function of the muon-proton system angle using the full MicroBooNE data sets. We also demonstrate the sensitivity of these results to nuclear effects and final state hadronic reinteraction modeling.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Structure refinement and anisotropic atomic displacement parameters of 1M Illite: Rietveld and pair distribution function analysis using synchrotron X-ray radiation
Illite, a widespread clay mineral, plays a pivotal role in geological processes, notably as an indicator in diagenetic and hydrothermal alteration environments, and possesses significant industrial relevance in applications including ceramics, construction and catalysis. However, challenges including its nanoscale crystallinity, structural disorder and frequent interstratification with other clay minerals have hindered detailed structural characterization using conventional X-ray diffraction (XRD) techniques. This study employs integrated synchrotron XRD and pair distribution function (PDF) analysis to elucidate the crystal structure of the 1M illite polytype, yielding the first determination of its anisotropic atomic displacement parameters (U aniso ). TheseU aniso parameters provide critical insights into atomic dynamics and static disorder within the structure, enabling a more refined understanding of structure–property relationships. This integrated approach, combining synchrotron XRD, Rietveld refinement and PDF analysis, yields a comprehensive structural characterization, capturing both average crystallographic and local atomic arrangements. Considering illite's widespread geological occurrence and industrial importance, this high-precision structural dataset, especially the determinedU aniso values, provides a crucial benchmark for future modeling and simulation efforts targeting accurate prediction of its physicochemical behavior.
Chemistry↗
Imaging Enhanced Simulation Validation for Directed Energy Deposition Additive Manufacturing (Final CRADA Report – NFE-22-09107)
Directed Energy Deposition (DED) is a welding-based metal Additive Manufacturing (AM) process that relies on the programmed rastering of an electric arc or laser induced weld pool to construct a component in a layerwise fashion. The induced complex thermal field and uneven thermal contraction depends on the printed geometry and scan pattern, and as such, accumulated residual stresses and distortions are complex and difficult to predict. Several prior works have resulted in tools to combat this issue; ANSYS has developed a thermoplastic simulation package targeting DED AM, and ORNL has developed an in-situ imaging sensor package ‘Stereo Correlated Optical and Pyrometric System’ (SCOPS) that can spatially monitor temperature and full field strain. In this CRADA, these tools are compared in order to validate the results and complimentarily address the weaknesses in each other.
36 MATERIALS SCIENCE↗
Implementation of disruptive designs for gas turbine components using direct energy deposition additive manufacturing
This research aims to develop a framework for establishing the correlation between in-situ monitoring data, process parameters, and microstructure evolution in blown-powder laser-directed energy deposition (DED) additive manufacturing (AM). To achieve this, a comprehensive manufacturing framework has been developed, spanning from in-situ data acquisition, melt-pool simulation, microstructure modeling, and statistical microstructure quantification. A machine learning-based surrogate model is constructed to predict melt pool geometry directly from in-situ coaxial camera data. The surrogate model is trained using outputs from a high-fidelity melt pool simulation, which provides accurate melt pool dimension data under varying process conditions. The predicted melt pool geometry is then used as input to a microstructure model to predict microstructural features. To rigorously compare and analyze microstructures, the project introduces statistical metrics that quantify differences based on key features such as morphology and texture. Microstructures are represented using advanced statistical descriptors including angular chord length distribution, two-point spatial statistics, orientation distribution function, and global spherical harmonic. These representations are used to compute four distinct “dissimilarity scores” that quantitatively capture differences in texture and morphology. This framework is demonstrated to enable automated calibration of simulation parameters by minimizing discrepancies between simulated and target microstructures. The technology developed in this project enables direct correlation between in-situ monitoring data and resulting microstructure, paving the way for adaptive microstructure control in metal AM. This capability strengthens the connection between process parameters and final material properties, facilitating more precise and reliable material design.
36 MATERIALS SCIENCE↗