Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accelerator 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 919 records · Page 51

( Fermi Fermi)/LAT study of gamma-ray emission in the direction of the Monoceros Loop supernova remnant

We present an analysis of the gamma-ray measurements by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope in the region of the supernova remnant (SNR) Monoceros Loop (G205.5+0.5). The brightest gamma-ray peak is spatially correlated with the Rosette Nebula, which is a molecular cloud complex adjacent to the southeast edge of the SNR. After subtraction of this emission by spatial modeling, the gamma-ray emission from the SNR emerges, which is extended and fit by a Gaussian spatial template. The gamma-ray spectra are significantly better reproduced by a curved shape than a simple power law. The luminosities between 0.2 and 300 GeV are ∼4×10{sup 34} erg s{sup −1} for the SNR and ∼3×10{sup 34} erg s{sup −1} for the Rosette Nebula, respectively. We argue that the gamma-rays likely originate from the interactions of particles accelerated in the SNR. The decay of neutral pions produced in nucleon–nucleon interactions of accelerated hadrons with interstellar gas provides a reasonable explanation for the gamma-ray emission of both the Rosette Nebula and the Monoceros SNR.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure of Energetic Particle Mediated Shocks Revisited

The structure of collisionless shock waves is often modified by the presence of energetic particles that are not equilibrated with the thermal plasma (such as pickup ions [PUIs] and solar energetic particles [SEPs]). This is relevant to the inner and outer heliosphere and the Very Local Interstellar Medium (VLISM), where observations of shock waves (e.g., in the inner heliosphere) show that both the magnetic field and thermal gas pressure are less than the energetic particle component pressures. Voyager 2 observations revealed that the heliospheric termination shock (HTS) is very broad and mediated by energetic particles. PUIs and SEPs contribute both a collisionless heat flux and a higher-order viscosity. We show that the incorporation of both effects can completely determine the structure of collisionless shocks mediated by energetic ions. Since the reduced form of the PUI-mediated plasma model is structurally identical to the classical cosmic ray two-fluid model, we note that the presence of viscosity, at least formally, eliminates the need for a gas sub-shock in the classical two-fluid model, including in that regime where three are possible. By considering parameters upstream of the HTS, we show that the thermal gas remains relatively cold and the shock is mediated by PUIs. We determine the structure of the weak interstellar shock observed by Voyager 1 . We consider the inclusion of the thermal heat flux and viscosity to address the most general form of an energetic particle-thermal plasma two-fluid model.

79 ASTRONOMY AND ASTROPHYSICS↗

PARTICLE ACCELERATION AND HEATING BY TURBULENT RECONNECTION

Turbulent flows in the solar wind, large-scale current sheets, multiple current sheets, and shock waves lead to the formation of environments in which a dense network of current sheets is established and sustains “turbulent reconnection.” We constructed a 2D grid on which a number of randomly chosen grid points are acting as scatterers (i.e., magnetic clouds or current sheets). Our goal is to examine how test particles respond inside this large-scale collection of scatterers. We study the energy gain of individual particles, the evolution of their energy distribution, and their escape time distribution. We have developed a new method to estimate the transport coefficients from the dynamics of the interaction of the particles with the scatterers. Replacing the “magnetic clouds” with current sheets, we have proven that the energization processes can be more efficient depending on the strength of the effective electric fields inside the current sheets and their statistical properties. Using the estimated transport coefficients and solving the Fokker–Planck (FP) equation, we can recover the energy distribution of the particles only for the stochastic Fermi process. We have shown that the evolution of the particles inside a turbulent reconnecting volume is not a solution of the FP equation, since the interaction of the particles with the current sheets is “anomalous,” in contrast to the case of the second-order Fermi process.

79 ASTRONOMY AND ASTROPHYSICS↗

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗

The Composition of GRB Jets and the ICMART Model

Models of gamma ray bursts (GRBs) are drawn from observations of light curves, spectra, and spectral evolution. The ICMART (Internal Collision-induced MAgnetic Reconnection & Turbulence) model and some of its features are presented. Increasing evidence points towards Poynting-flux-dominated jets in at least some (even a good fraction of) GRBs. The main emission component (Band) is of a synchrotron emission origin, produced by electrons accelerated in the emission region. The data seem to require that magnetic reconnection in the moderately-high sigma regime is the mechanism to accelerate particles. Extensive numerical simulations are needed to verify physical details of such a model, and some encouraging results have been obtained.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating functional verification of an integrated circuit

Illustrative embodiments include a method, system, and computer program product for accelerating functional verification in simulation testing of an integrated circuit (IC). Using a processor and a memory, a serial operation is replaced with a direct register access operation, wherein the serial operation is configured to perform bit shifting operation using a register in a simulation of the IC. The serial operation is blocked from manipulating the register in the simulation of the IC. Using the register in the simulation of the IC, the direct register access operation is performed in place of the serial operation.

42 ENGINEERING↗

LIA Performance Limitations

Location: SLAC National Accelerator Laboratory Presenter: Evan R Scott Date: 29 August 2023 Purpose: Part of SDRD for collaboration with SLAC

SDRD, SLAC, linear induction accelerator, cathode,↗

Thermal Aging Degradation of Cast Stainless Steels in LWR Systems

Cast austenitic stainless steels (CASSs) have been extensively used for the large components of light water reactor (LWR) power plants such as primary coolant piping and pump casing. Since a large number of massive CASS components are installed in every modern nuclear power plant, any significant degradation of CASS components will raise serious concerns on the integrity of the entire power plant. Upon the presumption that the replacement cost of such massive components is prohibitive, therefore, thorough scientific understanding and systematic monitoring the thermal aging degradation of CASS components should be essential requirements for assessing the plant safety. In the past few decades many researches have concluded that the formation of Cr-rich α-phase by Spinoidal decomposition of δ-ferrite phase, which is formed in high temperature region during cooling, is the primary mechanism for the thermal embrittlement. Cracking mechanism in the thermally-embrittled duplex stainless steels consists of cleavage initiation at ferrite and crack propagation via separation of ferrite-austenite interphase. It has been known, however, that such thermal aging embrittlement is observed only in accelerated aging conditions and the cast stainless steel grades used in nuclear plants have performed well without experiencing real embrittlement. It has been speculated that long-term thermal aging for ≥ 60 years could lead to further decomposition of designed phases and increase of precipitation. This will result in significantly increased susceptibility to embrittlement or other degradation due to thermal aging-radiation synergistic effect, irradiation-induced degradation, stress-corrosion cracking, and general corrosion. Currently, therefore, the thermal embrittlement of CASS components is one of the most serious concerns for extended operation of nuclear power plants. A scientific research is underway to address the issues related to the thermal aging degradation of CASS components. The research aims to provide knowledge-based conclusive prediction for the integrity of the CASS components of LWR power plants during extended operation up to and beyond 60 years. This presentation intends to provide an introductory overview on the thermal aging phenomena in LWR relevant conditions and discussion on the measurements of mechanical property degradation after accelerated and LWR relevant thermal aging treatments. First, the thermal aging effect on toughness is discussed in terms of the causes of embrittlement and influential parameters. Second, the talk will include a comprehensive discussion on the thermal aging effects in different CASS grades (i.e., CF3, CF3M, CF8, and CF8M), which will focus on the influence of alloying elements such as molybdenum and carbon. Third, approximate analysis on thermal reaction in CASS was carried out using Arrhenius equation to decide the effective temperature range for accelerated aging experiments to simulate 60 and 80 years of services in limited time frame. The results will provide valuable insights in the degree of thermal aging. Fourth, calculation on equilibrium precipitation was performed for model CASS alloys using the CALPHAD program, and the results will be used to describe the precipitation behaviors in association with the degradation in mechanical properties. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order S{sub N} Equations in Near-Void Regions

The common definition of the diffusion coefficient as the inverse of three times the transport cross section is not compatible with voids. Morel introduced a non-local tensor diffusion coefficient that remains finite in voids. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coefficient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction. It was demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors have been applied to modeling the TREAT reactor confirming the findings of Larsen and Trahan. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S{sub N} equations discretized with the discontinuous finite element method (DFEM). This NDA method uses a scalar diffusion coefficient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it suffers from very large and potentially unbounded diffusion coefficients in the low order problem. However, it was noted that the choice of the diffusion coefficient does not influence consistency of the method at convergence and hence the diffusion coefficient is essentially a free parameter. The choice of the diffusion coefficient does, however, affect the convergence behavior of the nonlinear diffusion iterations. Within this work we use Morel's non-local diffusion coefficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant enhancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coefficients, the periodical horizontal interface problem is used. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

GPU Acceleration of Kernel Density Estimators in Monte Carlo Neutron Transport Simulations

Kernel Density Estimators (KDEs) have recently been explored for use in Monte Carlo radiation transport simulations as an alternative to histogram tallies for capturing spatially-resolved quantities such as neutron flux and reaction rates. Histogram tallies suffer from large uncertainties when detailed spatial resolution of the quantity of interest is required; increasing the resolution of the underlying mesh increases the uncertainty in each histogram bin. KDEs obtain estimates of underlying densities on user-defined points, with the uncertainty at those points being independent of the spatial resolution desired since collisions and particle-track segments are allowed to contribute to the score at multiple tally points. Thus, KDEs show potential for obtaining estimates of a smoother spatially-resolved with reduced variance when compared to a histogram. However, this ability that enables KDEs to obtain lower uncertainties per particle history is more computationally expensive than using a traditional histogram bins on a structured mesh. While KDEs are capable of producing lower tally uncertainties compared to histogram binning for a given number of histories, their increased computational cost may make them less efficient. Even so, it is possible to reduce this additional cost by exporting the KDE tally routines onto Graphics Processing Units (GPUs). GPUs are commonly found on high-performance computing machines as a means of increasing performance without adding additional compute nodes. While it is difficult for existing Monte Carlo radiation transport codes to take advantage of GPUs without rewriting large portions of the code, it is possible to leverage GPUs through heterogeneous computing. This work uses a modified version of OpenMC and exports the KDE tally routines onto a GPU to obtain problem-dependent speedups ranging between 1.6 and 5.0 for the problem studied in this paper. While this methodology is applied to KDEs, the same ideas can be applied to other computationally expensive vectorizable portions of Monte Carlo algorithms. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Implementation of Krylov acceleration method for 3D S{sub N} transport code AETIUS with unstructured tetrahedral mesh - 25557

We have been developing an in-house deterministic neutron and gamma transport solver, AETIUS (An Easy modeling Transport code using Unstructured tetrahedral mesh, Shared memory parallel). To achieve fast results, BiCGStab (one of the Krylov acceleration methods) was implemented as a within-group iteration solver. In this paper, two test problems (simple double cube problem and helical tube in the outer cube problem) were used to verify the acceleration performance. The number of iterations and the calculation time were demonstrated and compared with the source iteration. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Effect of resistivity gradient on laser-driven electron transport and ion acceleration

The effect of resistivity gradient on laser-driven electron transport and ion acceleration is investigated using collisional particle-in-cell simulation. The study is motivated by recent proton acceleration experiments [Gizzi et al., Phys. Rev. ST Accel. Beams 14, 011301 (2011)], which showed significant effect of the resistivity gradient in layered targets on the proton angular spread. This effect is reproduced in the present simulations. It is found that resistivity-gradient generation of magnetic fields and inhibition of electron transport is significantly enhanced when the feedback interaction between the magnetic field and the fast-electron current is included. Filamentation of the laser-generated hot electron jets inside the target, considered as the origin of the nonuniform proton patterns observed in the experiments, is clearly suppressed by the resistive magnetic field. As a result, the electrostatic sheath field at the target back surface acquires a relatively smooth profile, which contributes to the superior quality of the proton beams accelerated off layered targets in the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of Plasma Waveguides to High Energy Accelerators

This grant supported basic experimental, theoretical and computer simulation research into developing a compact, high pulse repetition rate laser accelerator using the direct laser acceleration mechanism in plasma-based slow wave structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relativistic electron acceleration by oblique whistler waves

Test-particle simulations of electrons interacting with finite-amplitude, obliquely propagating whistler waves are carried out in order to investigate the acceleration of relativistic electrons by these waves. According to the present findings, an efficient acceleration of relativistic electrons requires a narrow range of oblique propagation angles, close to the whistler resonance cone angle, when the wave amplitude is held constant at relatively low value. For a constant wave propagation angle, it is found that a range of oblique whistler wave amplitudes permits the acceleration of relativistic electrons to O(MeV) energies. An initial distribution of test electrons is shown to form a power-law distribution when plotted in energy space. It is also found that the acceleration is largely uniform in electron pitch-angle space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerically optimized structures for dielectric asymmetric dual-grating laser accelerators

Optical scale dielectric structures are promising candidates to realize future compact, low cost particle accelerators, since they can sustain high acceleration gradients in the range of GeV/m. Here, we present numerical simulation results for a dielectric asymmetric dual-grating accelerator. It was found that the asymmetric dual-grating structures can efficiently modify the laser field to synchronize it with relativistic electrons, therefore increasing the average acceleration gradient by ∼10% in comparison to symmetric structures. The optimum pillar height which was determined by simulation agrees well with that estimated analytically. The effect of the initial kinetic energy of injected electrons on the acceleration gradient is also discussed. Finally, the required laser parameters were calculated analytically and a suitable laser is proposed as energy source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗