Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PARTICLE MODELS”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗
DOE OSTI2017

PIC simulations of wave-particle interactions with an initial electron velocity distribution from a kinetic ring current model

Whistler wave-particle interactions play an important role in the Earth inner magnetospheric dynamics and have been the subject of numerous investigations. By running a global kinetic ring current model (RAM-SCB) in a storm event occurred on Oct 23–24 2002, we obtain the ring current electron distribution at a selected location at MLT of 9 and L of 6 where the electron distribution is composed of a warm population in the form of a partial ring in the velocity space (with energy around 15 keV) in addition to a cool population with a Maxwellian-like distribution. The warm population is likely from the injected plasma sheet electrons during substorm injections that supply fresh source to the inner magnetosphere. These electron distributions are then used as input in an implicit particle-in-cell code (iPIC3D) to study whistler-wave generation and the subsequent wave-particle interactions. Here, we find that whistler waves are excited and propagate in the quasi-parallel direction along the background magnetic field. Several different wave modes are instantaneously generated with different growth rates and frequencies. The wave mode at the maximum growth rate has a frequency around 0.62ω ce , which corresponds to a parallel resonant energy of 2.5 keV. Linear theory analysis of wave growth is in excellent agreement with the simulation results. These waves grow initially due to the injected warm electrons and are later damped due to cyclotron absorption by electrons whose energy is close to the resonant energy and can effectively attenuate waves. The warm electron population overall experiences net energy loss and anisotropy drop while moving along the diffusion surfaces towards regions of lower phase space density, while the cool electron population undergoes heating when the waves grow, suggesting the cross-population interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Using Particle Image Velocimetry to Determine Turbulence Model Parameters

The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.

42 ENGINEERING↗
DOE OSTI2021

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗
NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗
DOE OSTI2017

Search for a massive invisible particle X{sup 0} in B{sup +} → e{sup +}X{sup 0} and B{sup +} → μ{sup +}X{sup 0} decays

We present a search for a non-Standard-Model invisible particle X{sup 0} in the mass range 0.1–1.8 GeV/c{sup 2} in B{sup +} → e{sup +}X{sup 0} and B{sup +} → μ{sup +}X{sup 0} decays. The results are obtained from a 711 fb{sup −1} data sample which corresponds to 772 × 10{sup 6}BB̄ pairs, collected at the Υ(4S) resonance with the Belle detector at the KEKB e{sup +}e{sup −} collider. One B meson is fully reconstructed in a hadronic mode to determine the momentum of the lepton of the signal decay in the rest frame of the recoiling partner B meson. We find no evidence of a signal and set upper limits on the order of 10{sup −6}.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

ELECTRON ACCELERATION IN PULSAR-WIND TERMINATION SHOCKS: AN APPLICATION TO THE CRAB NEBULA GAMMA-RAY FLARES

The γ -ray flares from the Crab Nebula observed by AGILE and Fermi -LAT reaching GeV energies and lasting several days challenge the standard models for particle acceleration in pulsar-wind nebulae because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Previous modeling has suggested that the synchrotron limit can be exceeded if the electrons experience electrostatic acceleration, but the resulting spectra do not agree very well with the data. As a result, there are still some important unanswered questions about the detailed particle acceleration and emission processes occurring during the flares. We revisit the problem using a new analytical approach based on an electron transport equation that includes terms describing electrostatic acceleration, stochastic wave-particle acceleration, shock acceleration, synchrotron losses, and particle escape. An exact solution is obtained for the electron distribution, which is used to compute the associated γ -ray synchrotron spectrum. We find that in our model the γ -ray flares are mainly powered by electrostatic acceleration, but the contributions from stochastic and shock acceleration play an important role in producing the observed spectral shapes. Our model can reproduce the spectra of all the Fermi -LAT and AGILE flares from the Crab Nebula, using magnetic field strengths in agreement with the multi-wavelength observational constraints. We also compute the spectrum and duration of the synchrotron afterglow created by the accelerated electrons, after they escape into the region on the downstream side of the pulsar-wind termination shock. The afterglow is expected to fade over a maximum period of about three weeks after the γ -ray flare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Simple physics of the partly pinned fluid systems

In this paper, we consider some aspects of the physics of the partly pinned (PP) systems obtained by freezing in place particles in equilibrium bulk fluid configurations in the normal (nonglassy) state. We first discuss the configurational overlap and the disconnected density correlation functions, both in the homogeneous and heterogeneous cases, using the tools of the theory of adsorption in disordered porous solids. The relevant Ornstein-Zernike equations are derived, and asymptotic results valid in the regime where the perturbation due to the pinning process is small are obtained. Second, we consider the homogeneous PP lattice gas as a means to make contact between pinning processes in particle and spin systems and show that it can be straightforwardly mapped onto a random field Ising model with a strongly asymmetric bimodal distribution of the field. Possible implications of these results for studies of the glass transition based on PP systems are also discussed.

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

Improved Creep-­Fatigue Models on Advanced Materials for SFR Applications (Final Report)

This report provides an update on the materials performance criteria and methodology relevant to sodium-cooled fast reactors (SFRs). The report is the second deliverable (Level 2) in FY11 (M2A11AN040303) under the work package A-11AN040303 “Materials Performance Criteria and Methodology” as part of the Advanced Structural Materials Program for the Advanced Reactor Concepts. The overall objective of the Materials Performance Criteria and Methodology work project is to evaluate the key requirements for the ASME Code qualification and the Nuclear Regulatory Commission (NRC) approval of advanced structural materials in support of the design and performance of sodium-cooled fast reactors. Advanced materials are a critical element in the development of fast reactor technologies. Enhanced materials performance not only improves safety margins and provides design flexibility, but also is essential for the economics of future advanced fast reactors. Qualification and licensing of advanced materials are prominent needs for the development and implementation of advanced fast reactor technologies. Nuclear structural component designs in the U.S. comply with the ASME Boiler and Pressure Vessel (B&PV) Code Section III (Rules for Construction of Nuclear Facility Components), and the NRC grants licensing. As the SFRs will operate at higher temperatures than the current light water reactors (LWRs), the design of elevated-temperature components must comply with ASME Section III Subsection NH (Class 1 Components in Elevated Temperature Service). A number of technical issues relevant to materials performance criteria and high temperature design methodology in the SFR were identified and presented in earlier reports. A viable approach to resolve these issues and the R&D priority were also recommended. The development of mechanistically based creep-fatigue interaction models for life prediction and reliable data extrapolation was chosen to be the central focus in near-term efforts. Our current focus is on the creep-fatigue damage issue in high-strength ferritic/martensitic steels such as mod.9Cr-1Mo (G91) and NF616 (G92) steels. The current ASME creep-fatigue design rule puts severe limits of fatigue and creep loads for G91 steel, the lead structural material for fast reactors. High-strength ferritic/martensitic steels behave fundamentally differently from austenitic stainless steels, for which the current ASME creep-fatigue design rules were developed. The unique deformation and damage characteristics in G91 steel, e.g. cyclic softening, degradation of creep and rupture strength during cyclic service, demands a new creep-fatigue design procedure that explicitly accounts for the material’s unique creep-fatigue behavior. G92 steel, a variant of G91 steel, is the lead candidate in the Advanced Alloy Development program. There is currently no creep-fatigue design rule available for this advanced alloy. To support the predictive model development and resolve the over-conservative issue with the ASME design rule for G91 steel, we recovered stress relaxation data from thirteen creep-fatigue tests conducted in the late 1980s and early 1990s at Oak Ridge National Laboratory, and conducted extensive data analysis in FY10. Based on this database and available literature data, we have developed a Cyclic Softening Model and a Stress Relaxation Model specifically applied to G91 steel. The Cyclic Softening Model describes the cyclic stress variation as a function of cycle number during creep-fatigue loading, and well captures the cyclic softening effect in G91 steel. The Stress Relaxation Model predicts the stress relaxation curve during the hold time of cyclic loading, and its dependence on the cyclic softening effect of G91 steel. Two creep-fatigue test systems have been set up to conduct model-oriented creep-fatigue experiments. Results from three creep-fatigue experiments on Heat 30176 of G91 steel are presented in the report: (1) a test performed at 550°C with total strain amplitude of 0.25% and a hold time of 60 s applied at peak strains in both tension and compression, (2) a test performed at 400°C with total strain amplitude of 0.5% and a tension and compression hold time of 60 s, and (3) a test performed at 20°C with total strain amplitude of 0.5% and a tension and compression hold time of 60 s. Three different test temperatures were chosen for two primary reasons: to determine the transition temperature of creep-fatigue interaction from a plasticity-dominated process to a creep-dominated process; to develop unified creep-fatigue models applicable for a wide range of temperatures. With the detailed information obtained in the creep-fatigue experiments, we have further refined the Cyclic Softening Model and Stress Relaxation Model. The Cyclic Softening Model and Stress Relaxation Model were incorporated into the Bilinear Creep-Fatigue Damage Model, forming a powerful framework for the life prediction by an improved creep-fatigue evaluation methodology. The cyclic softening and consecutive stress relaxation responses in G91 steel can be considered for each cycle of the entire cyclic process in the model calculation. The creep-fatigue damage in G91 steel was evaluated using the improved Bilinear Creep-Fatigue Damage Model without any safety margins. It was found that due to the significant cyclic softening in G91 steel, the unit creep damage in each cycle decreased continuously with increasing number of cycles. This continuous change in creep damage with cycling has significant implication in the evaluation of the total creep-fatigue damage. The total creep damage can be one order of magnitude lower when the unit creep damage at half-life was used in the evaluation in comparison with the calculations using the unit creep damage at cycle 10. An important conclusion was that the accurate assessment of the creep-fatigue damage in G91 steel must consider the history of cycling. This analysis provides a more realistic assessment of creep-fatigue damage and identifies critical factors that must be considered in the development of future advanced models for high-strength ferritic/martensitic alloys. An accelerated creep-fatigue testing methodology is proposed for the design of model-oriented experiments to guide and validate model developments. Future work will focus on creep-fatigue experiments of both G91 and G92 steels, and development of creep-fatigue models applicable to advanced ferritic/martensitic alloys. The Cyclic Softening Model and Stress Relaxation Model will be improved by incorporating microstructural parameters of dislocation density, subgrain structure, and second-phase particles. Collaborative modeling of combined Cyclic Softening Model, Stress Relaxation Model, and Interactive Damage Rate Model will be pursued for life prediction under creep-fatigue loading with extended hold time. Accelerated, well-controlled creep-fatigue tests will be conducted to provide high-quality experimental data to refine key material parameters and validate models.

Li, Meimei↗
DOE OSTI2013

SciDAC Center for Plasma Edge Simulation

This project with a total funding of $592,998 for six years has partially supported four postdoctoral researchers at the University of California, Irvine (UCI). The UCI team has formulated electrostatic and electromagnetic global gyrokinetic particle simulation models with kinetic electrons, implemented these models in the edge code XGC1, performed benchmark between GTC and XGC1, developed computational tools for gyrokinetic particle simulation in tokamak edge geometry, and initiated preparatory study of edge turbulence using GTC code. The research results has been published in 12 papers and presented at many international and national conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensor Placement Optimization Study for the Built Environment: Operational Use Cases

Systems of fixed-position radiation sensors can provide information that assists emergency responders following nuclear and radiological incidents. State, local, tribal, and territorial (SLTT) government agencies that implement systems of fixed-position sensors are faced with numerous decisions regarding sensor selection, quantity, and placement. To develop guidance on implementation of radiation detection systems, we simulated the release of radioactive material in an urban environment using a combination of three models: the Weather Research Forecasting (WRF), Quick Urban and Industrial Complex (QUIC), and Monte Carlo N-Particle (MCNP) models. We then evaluated the performance of several hypothetical sensor systems. The small number of simulations we conducted are not sufficient to generate definitive design guidance for radiation sensor systems, but we did identify trends that would be of interest to emergency planners. For a scenario that releases 1000 curies of Cs-137, radiation detectors were needed at 500-meter intervals to have a high likelihood of event detection and to estimate source location and plume detection. We also noted that optimal detector altitude varied with distance to the source. We recommend additional research in this area be conducted to support developing sensor placement guidelines that expand on a range of locations, isotopes, activity levels, and different weather conditions. Original simulation strategies included a range of environments, additional radioisotopes (Am241 and AmBe), and a larger selection of sensor types. These types of expansions would support SLTT guidance on sensor system recommendations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1990

Fundamentals of chemistry modeling applicable to a vectorized particle simulation

This paper describes the fundamentals of extending the vectorized particle simulation method derived by Baganoff and McDonald (1990), McDonald and Baganoff (1988), and McDonald (1989) for modeling chemically reacting flows. Details of reaction mechanics per reaction are presented, with particular attention given to the quantum nature of the vibrational mode. The models of reactive flows developed here were verified through a simulation of a superheated diatomic gas relaxing thermochemically to equilibrium in a reservoir.

Haas, Brian L.↗
NASA NTRS1991

On the observability of the gamma-ray line flux from dark matter annihilation

The limits on the possible cosmic gamma-ray line flux from the two-photon annihilation of dark matter in the Galaxy are discussed. These limits are derived using both particle physics and cosmological constraints on dark matter candidates which arise in supersymmetric extensions of the standard model of particle physics. Results are given in terms of allowed and prescribed areas in the flux-energy plane. Then these bounds are used to consider the observability of the line flux above continuum background fluxes using future high-resolution gamma-ray telescopes.

Rudaz, S.↗
NASA NTRS2000

Analysis of MIR-18 results for physical and biological dosimetry: radiation shielding effectiveness in LEO

We compare models of radiation transport and biological response to physical and biological dosimetry results from astronauts on the Mir space station. Transport models are shown to be in good agreement with physical measurements and indicate that the ratio of equivalent dose from the Galactic Cosmic Rays (GCR) to protons is about 3/2:1 and that this ratio will increase for exposures to internal organs. Two biological response models are used to compare to the Mir biodosimetry for chromosome aberration in lymphocyte cells; a track-structure model and the linear-quadratic model with linear energy transfer (LET) dependent weighting coefficients. These models are fit to in vitro data for aberration formation in human lymphocytes by photons and charged particles. Both models are found to be in reasonable agreement with data for aberrations in lymphocytes of Mir crew members: however there are differences between the use of LET dependent weighting factors and track structure models for assigning radiation quality factors. The major difference in the models is the increased effectiveness predicted by the track model for low charge and energy ions with LET near 10 keV/micrometers. The results of our calculations indicate that aluminum shielding, although providing important mitigation of the effects of trapped radiation, provides no protective effect from the galactic cosmic rays (GCR) in low-earth orbit (LEO) using either equivalent dose or the number of chromosome aberrations as a measure until about 100 g/cm 2 of material is used.

Flight Experiment↗
NASA NTRS2018

Current Activities in the Advanced Exploration Systems ISRU Project

Media Filter MOXIE support: completed flow performance tests and initial dust loading tests of the MOXIE HEPA filter in the Mars Flow Loop Designed and fabricated prototype full-scale scroll filter for testing of renewable filter performance Mars Flow Loop upgrades: installed more sensitive instrumentation, improved imaging, and increased run duration capability Electrostatic Precipitator (ESP) Plasma-physics based model predicts particle charge and trajectory Characterizing dust environment from fluidized-bed injection dust environment using laser sheet visualization and Fine Particle Analyzer Sabatier Design Study Defined the Sabatier design space, including reactor type, thermal management, gas recycling/separation Modeled Sabatier systems with one or two reactors with different types of thermal management Adding thermal management schemes and recycling and gas separation to model Catalyst Screening Performing compression testing on new and used catalyst pellets Performing vibration testing on different catalyst pellet types under different load conditions Preparing to test catalysts for performance and degradation

In-Situ Resource Utilization↗
DOE OSTI2024

Old-Aged groundwater contributes to mountain hillslope hydrologic dynamics

Understanding connectivity between the soil and deeper bedrock groundwater is needed to accurately predict a watershed’s response to perturbation, such as drought. Yet, the bedrock groundwater dynamics in mountainous environments are typically under-constrained and excluded from watershed hydrologic models. Here, we investigate the role of groundwater characterized with decadal and longer water ages on the hydrologic and mass-transport processes within a steep snow-dominated mountain hillslope in the Central Rocky Mountains (USA). We quantify subsurface and surface water mass-balance, groundwater flowpaths, and age distributions using the ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models, which are compared to hydrometric and environmental tracer observations. An ensemble of models with varied soil and hydrogeologic parameters reproduces observed groundwater levels and century-scale mean ages inferred from environmental tracers. The numerical models suggest soil water near the toe of the hillslope contains considerable (>60 % of the mass-flux) contributions from bedrock flowpaths characterized with water ages >10 years. Flowpath connectivity between the deeper bedrock and soil systems is present throughout the year, highlighting the potentially critical role of groundwater with old ages on processes such as evapotranspiration and streamflow generation. The coupled numerical model and groundwater age observations show the bedrock groundwater system influences the hillslope hydrodynamics and should be considered in mountain watershed conceptual and numerical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Fast-ignition transport studies: Realistic electron source, integrated particle-in-cell and hydrodynamic modeling, imposed magnetic fields

Transport modeling of idealized, cone-guided fast ignition targets indicates the severe challenge posed by fast-electron source divergence. The hybrid particle-in-cell (PIC) code Zuma is run in tandem with the radiation-hydrodynamics code Hydra to model fast-electron propagation, fuel heating, and thermonuclear burn. The fast electron source is based on a 3D explicit-PIC laser-plasma simulation with the PSC code. This shows a quasi two-temperature energy spectrum and a divergent angle spectrum (average velocity-space polar angle of 52 Degree-Sign ). Transport simulations with the PIC-based divergence do not ignite for >1 MJ of fast-electron energy, for a modest (70 {mu}m) standoff distance from fast-electron injection to the dense fuel. However, artificially collimating the source gives an ignition energy of 132 kJ. To mitigate the divergence, we consider imposed axial magnetic fields. Uniform fields {approx}50 MG are sufficient to recover the artificially collimated ignition energy. Experiments at the Omega laser facility have generated fields of this magnitude by imploding a capsule in seed fields of 50-100 kG. Such imploded fields will likely be more compressed in the transport region than in the laser absorption region. When fast electrons encounter increasing field strength, magnetic mirroring can reflect a substantial fraction of them and reduce coupling to the fuel. A hollow magnetic pipe, which peaks at a finite radius, is presented as one field configuration which circumvents mirroring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗