Search NASA⌕ Search

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.

At least 163 records · Page 9

Implementation of INCL nuclear model in GENIE Generator

The Liège Intranuclear Cascade (INCL) model is a nuclear-physics model that simulates hadron (baryon, anti-baryon and meson) reactions on nuclei, for incident energies ranging from a few tens of MeV to 10-20 GeV. The INCL model has been well validated by hadron scattering data. In my work, I implement an interface in GENIE to use the INCL nuclear model in the simulations of both the initial state of the target nucleus and the Final State Interaction in neutrino-nucleus interaction. It has a consistent treatment of nuclear models in both neutrino interaction and hadron rescattering. A full event record including neutrino vertex and each vertex of hadron rescattering has been accomplished. Several processes, e.g. cluster production, Delta transportation and de-excitation will also be included as benefits of the implementation of the INCL model in GENIE. I will show some initial simulation results showcasing the new GENIE features and discuss plans for making them available for use in experimental analyses.

Liu, Liang [Fermilab] (ORCID:000000026753925X)↗

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗

Constraining hydrodynamic models of inertial confinement fusion implosions using capsule surrogate experiments

We conduct capsule surrogate experiments at the National Ignition Facility to calibrate radiation hydrodynamic simulations to infer hydrodynamic conditions that are not observable in indirect drive ignition implosions. We tune the simulations by applying laser power and cross beam energy transfer (CBET) saturation multipliers to match the observables from capsule surrogate experiments. Shock timing, velocity, and symmetry are measured in liquid D 2 filled Keyhole capsule surrogate experiments and implosion trajectory, stagnation time, and shape time history are measured in in-flight 2D backlit x-ray radiography experiments (“2DConA”) of D 2 gas filled capsule implosions. Calibrated simulations suggest that the N210808 ignition implosion (fusion target gain = 0.7) had a shell mass remaining at stagnation of less than the nominal %5 (3.8%) and resulted in less confinement. For N221204, the shell was made 5.75 μm thicker to trade implosion velocity for increased confinement and resulted in a target gain = 1.5 with a shell mass remaining of 5.7%. Furthermore, a single adjusted model can reproduce all shock timing data as changes are made to shell thickness (79–85 μm) and laser wavelength separation (1.8–4.0 Å). However, for the 2DConA implosions, a 5% variation in the peak power laser multipliers and a 30% variation in late-time CBET between experiments are needed to match the observed stagnation times, in-flight $P_2$ shape, and hot-spot $P_2$ shape. While progress is being made to improve the models in simulations using focused experiments, capsule surrogate experiments will continue to be needed to optimize future ignition designs.

Lasers↗

Design of a high power 8 MeV electron beam injector for the LERF

We present the beamline design for a high power 8 MeV electron injector to be installed at the Low Energy Recirculator Facility (LERF). The injector is designed to provide beam power of up to 80 kW for testing positron source target candidates. We discuss the results of beam dynamics simulations conducted using the General Particle Tracer (GPT) software, where electromagnetic beamline elements were optimized using multi-objective genetic algorithm tools. Additionally, we present preliminary studies on secondary particle generation from both rotating solid and free surface liquid metal jet targets described elsewhere in this workshop.

Lizarraga-Rubio, Victor [Thomas Jefferson National↗

Design of a high power 8 MeV electron beam injector for the LERF

Here, we present the beamline design for a high power 8 MeV electron injector to be installed at the Low Energy Recirculator Facility (LERF). The injector is designed to provide beam power of up to 80 kW for testing positron source target candidates. We discuss the results of beam dynamics simulations conducted using the General Particle Tracer (GPT) software, where electromagnetic beamline elements were optimized using multi-objective genetic algorithm tools. Additionally, we present preliminary studies on secondary particle generation from both rotating solid and free surface liquid metal jet targets described elsewhere in this workshop.

Lizarraga-Rubio, Victor Manuel [Thomas Jefferson N↗

Swelling and Fission Gas Release of U-10Mo and U-17Mo Following Neutron Irradiation at 250 – 500°C

The performance of metallic fuel alloys, U-10Mo and U-17Mo, was examined using the MiniFuel test system in the High Flux Isotope Reactor. Approximately 0.8 mm thick disks were irradiated at target temperatures of 250°C, 350°C, 450°C, and 500°C up to three different fission densities, culminating at a maximum fission density of 6.8 × 10 20 cm −3 . Fission rates decayed from 3 – 6 × 10 13 cm −3 s −1 to 2 – 4 × 10 13 cm −3 s −1 over the course of the longest, eight cycle, irradiation as the 235 U was consumed and the 239 Pu concentration accumulated. Actual irradiation temperature on the last day of irradiation was measured via dilatometry using the SiC passive thermometry recovered from the MiniFuel subcapsules and compared favorably with the thermal calculations using as-built geometry and test conditions. Furthermore, average simulated temperatures were within 50°C of the target temperatures, except for the 500°C irradiation, for which the temperature variation was 62°C, 32°C, and 90°C in the in two, four, and eight cycle irradiations, respectively. Fission gas release (FGR) measurements showed no release above recoil for any U-17Mo fuels or for the U-10Mo fuel irradiated at 250°C. Finally, significant (40%–80%) FGR was found for the medium- to highest-burnup U-10Mo samples irradiated at target temperatures 350°C–500°C. Significant FGR correlated with sample thickness swelling, which was as high as 13%–35% for high-release samples and below 7% for all other (low–gas release) samples.

FGR↗

Radioisotope production at the Spallation Neutron Source: Design concept of experimental target station

Completion of the Proton Power Upgrade Project for the Spallation Neutron Source (SNS) accelerator at Oak Ridge National Laboratory opens an opportunity to utilize reserve beam power of more than 100 kW for applications beyond neutron production. One of these applications is the production of critical radionuclides. To demonstrate the feasibility of using the reserve beam power to produce radioisotope at SNS, a design concept of a small-scale experimental target station in the Linac Dump area has been developed. This experimental facility will provide isotope yield benchmarking data using protons in the GeV range. It will also enable additional research and development in isotope handling and radiochemical separation. The target station consists of a target module enclosed in a vessel and concrete shielding. Particle transport calculations and thermo-mechanical simulations are used to determine beam parameters, decay time, isotope yield, shielding dimensions, and target design parameters. Calculations verified that the irradiated capsule can be handled manually using hands-off tools and transported to a hot cell in a shielded container for post-irradiation characterizations.

Lee, Yong Joong [ORNL] (ORCID:0000000298381723)↗

SPADES (Scalable Parallel Discrete Events Simulation) [SWR-24-99]

SPADES (Solver for PArallel Discrete Event Simulation) is an open-source parallel discrete event simulation (PDES) package built on the AMReX library. Targeted at solving discrete event systems in parallel, this software package aims to be performance portable and scalable on heterogeneous computing architectures, e.g., graphic processing units (GPU). SPADES implements optimistic synchronization with rollback through an implementation of the Time Warp algorithm. An alternative conservative synchronization approach is also implemented using the Lower Bound on Incoming Time Stamp. In our implementation, logical processes are represented as cells in a grid and event messages are represented as particles. SPADES supports various parallel decomposition strategies, including the use of the Message Passing Interface (MPI) and OpenMP threading. All major GPU architectures (e.g., Intel, AMD, NVIDIA) are supported through the use of performance portability functionalities implemented in AMReX. The SPADES software is released in NREL Software Record SWR-24-99 “SPADES (Scalable Parallel Discrete Events Simulation)”.

Henry de Frahan, Marc [National Renewable Energy L↗

AuriDESI: mock catalogues for the DESI Milky Way Survey

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

dynamics↗

Automated shaker placement and regularized input estimation for MIMO testing.

Multi-input, multi-output (MIMO) testing is used in component qualification to reproduce operational responses in the laboratory. It is often preferred to single-input and base-shake testing because of the potential for equivalent or better tests using smaller actuators and shorter test suites. Given a target response, two key steps in MIMO test design are selecting actuator locations and solving for input loads. Actuator locations are often manually selected using expert judgment. If an automatic method is used, locations are usually determined by simulating the vibration control problem and minimizing a combination of the input energy and control residuals. To select a configuration, the relative importance of input energy and residuals must be specified. Specifying relative weights is, in general, a manual and subjective process. This paper develops an objective function that compares actuator configurations based on control accuracy and required input energy without any manual parameter tuning. The objective function uses an optimally selected tradeoff parameter for each candidate configuration. To choose actuator locations using the new objective function, a pivoting algorithm for integer programming problems is developed. Starting with an initial configuration (such as the one generated by a greedy algorithm), the pivoting algorithm guarantees an objective function decrease in each iteration until convergence is reached. In a simulation featuring a structure excited by a diffuse acoustic field, electrodynamic shaker locations and regularized inputs are solved for without any analyst-specified parameters. Simulations are performed in MIMO configurations where the number of target responses is less than, equal to, and greater than the number of actuators.

Multi-input multi-output↗

CryoTEN: efficiently enhancing cryo-EM density maps using transformers

Abstract Motivation Cryogenic electron microscopy (cryo-EM) is a core experimental technique used to determine the structure of macromolecules such as proteins. However, the effectiveness of cryo-EM is often hindered by the noise and missing density values in cryo-EM density maps caused by experimental conditions such as low contrast and conformational heterogeneity. Although various global and local map-sharpening techniques are widely employed to improve cryo-EM density maps, it is still challenging to efficiently improve their quality for building better protein structures from them. Results In this study, we introduce CryoTEN—a 3D UNETR++ style transformer to improve cryo-EM maps effectively. CryoTEN is trained using a diverse set of 1295 cryo-EM maps as inputs and their corresponding simulated maps generated from known protein structures as targets. An independent test set containing 150 maps is used to evaluate CryoTEN, and the results demonstrate that it can robustly enhance the quality of cryo-EM density maps. In addition, automatic de novo protein structure modeling shows that protein structures built from the density maps processed by CryoTEN have substantially better quality than those built from the original maps. Compared to the existing state-of-the-art deep learning methods for enhancing cryo-EM density maps, CryoTEN ranks second in improving the quality of density maps, while running >10 times faster and requiring much less GPU memory than them. Availability and implementation The source code and data are freely available at https://github.com/jianlin-cheng/cryoten.

Biochemistry & Molecular Biology↗

Characterizing electron-collecting CdTe for use in a 77 ns burst-rate imager

The Keck-PAD (pixel array detector) was developed at Cornell as a burst-rate imager capable of recording images from successive electron bunches (153 ns period) from the Advanced Photon Source (APS). Both Si and hole-collecting Schottky CdTe have been successfully bonded to this ASIC (application-specific integrated circuit) and used with this frame rate. The facility upgrades at the APS will lower the bunch period to 77 ns, which will require modifications to the Keck-PAD electronics to image properly at this reduced period. In addition, operation at high X-ray energies will require a different sensor material having a shorter charge collection time. For the target energy of 40 keV for this project, simulations have shown that electron-collecting CdTe should allow >90% charge collection within 35 ns. This collection time will be sufficient to sample the signal from one frame and prepare for the next. 750 µm-thick electron-collecting Schottky CdTe has been obtained from Acrorad and bonded to two different charge-integrating ASICs developed at Cornell, the Keck-PAD and the CU-APS-PAD. Carrier mobility has been investigated using the detector response to single X-ray bunches at the Cornell High Energy Synchrotron Source and to a pulsed optical laser. The tests indicate that the collection time will meet the requirements for 77 ns imaging.

36 MATERIALS SCIENCE↗

Transverse BBU Suppression with Feedback for Energy-Recovery Linacs

The Electron–Ion Collider (EIC) achieves its design luminosity by cooling the ion beam with a high-current electron beam generated in an energy-recovery linac (ERL). The baseline ERL lattice employs a BNL five-cell cavity that is frequency-scaled to 197 MHz, 591 MHz, and 1.773 GHz, raising concerns about multibunch beam-breakup (BBU) instabilities. Threshold currents for each frequency option are established with two independent BBU tracking codes, providing cross-validated operating margins that guide cavity selection. To further increase the current limit, we incorporate a digital transverse feedback (FB) system that targets the dominant higher-order modes (HOMs). Simulations show the FB raises the BBU threshold by roughly an order of magnitude: power-spectral-density analysis of the beam centroid at the linac exit confirms strong suppression of the aliased HOM peak, although some spectral growth appears at secondary frequencies. This combined study quantifies baseline BBU limits, demonstrates effective active mitigation, and charts a practical path toward robust, high-current operation of the EIC ERL.

Accelerator Physics↗

MATHUSLA: An External Long-Lived Particle Detector to Maximize the Discovery Potential of the HL-LHC

We present the current status of the MATHUSLA (MAssive Timing Hodoscope for Ultra-Stable neutraL pArticles) long-lived particle (LLP) detector at the HL-LHC, covering the design, fabrication and installation at CERN Point 5. MATHUSLA40 is a 40 m-scale detector with an air-filled decay volume that is instrumented with scintillator tracking detectors, to be located near CMS. Its large size, close proximity to the CMS interaction point and about 100 m of rock shielding from LHC backgrounds allows it to detect LLP production rates and lifetimes that are one to two orders of magnitude beyond the ultimate reach of the LHC main detectors. This provides unique sensitivity to many LLP signals that are highly theoretically motivated, due to their connection to the hierarchy problem, the nature of dark matter, and baryogenesis. Data taking is projected to commence with the start of HL-LHC operations. We summarize the new 40m design for the detector that was recently presented in the MATHUSLA Conceptual Design Report, alongside new realistic background and signal simulations that demonstrate high efficiency for the main target LLP signals in a background-free HL-LHC search. We argue that MATHUSLA's uniquely robust expansion of the HL-LHC physics reach is a crucial ingredient in CERN's mission to search for new physics and characterize the Higgs boson with precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Super-X and conventional divertor configurations in MAST-U ohmic L-mode; a comparison facilitated by interpretative modelling

Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,\textrm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16 $\,\pm\,7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e = 5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream n e and T e profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,\textrm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗

Production Target Studies at the M4 Diagnostic Absorber

High power production targetry is becoming essential to the future of HEP, but our current physics simulations are lacking validation in the lower energy ranges that many proposed secondary beam experiments will use. We propose a study of target materials at the M4 diagnostic absorber to monitor target temperatures during irradiation by an 8 GeV proton beam and provide some experimental validation for our simulations at lower energies as well as potentially perform post-irradiation examination of samples to understand changes to material properties.

Bloomer, Madeleine [UC, Davis] (ORCID:000900080329↗