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 1,495 records · Page 83

Development of a Trajectory-Centric CFD-RBD Framework for Advanced Multidisciplinary/Multiphysics Simulation

The desire to model increasingly complex unsteady flow phenomena drives coupling of physics-based disciplinary analysis tools, such as coupled aerodynamics-rigid body dynamics simulations. This paper documents the creation of a framework linking the six-degree-of-freedom trajectory propagator POST2 with NASA’s FUN3D computational fluid dynamics flow solver. Cross-code verification between the framework and a CFD-centric 6DOF code is performed using the Army-Navy Finner projectile experiencing unsteady accelerating flow. Free-flight simulations of an entry vehicle ballistic range test are validated against physical and computational experiments.

Zachary J Ernst↗

Fluid-Structure Interaction Simulations Of Supersonic Parachute Inflation: Model Sensitivities

Fluid-structure interaction simulations have the potential to accelerate the certification process of parachute systems for interplanetary robotic exploration missions while helping to mitigate risks by covering a wider range of scenarios than can be flight tested. The Launch, Ascent, and Vehicle Aerodynamics team is working toward actualizing this potential by further developing the capability to perform such fluid-structure interaction simulations and validating it with the best supersonic parachute flight test data available. We perform several simulations of the Advanced Supersonic Parachute Inflation Research Experiments’ third flight test (SR03) where we independently vary several parameters to determine their sensitivity on the predicted drag. We investigate sensitivity of the drag to the added mass of the confluence fitting, to the effective porosity, to how the gore seams are modeled, to viscous effects, and to deceleration. Modeling the added mass of the confluence fitting has virtually no impact on the drag, but is worthwhile to include to better match the overall system’s mass, which is important for the simulations with deceleration. Parachute drag is only marginally sensitive to porosity and viscous effects. Maximum parachute drag is most sensitive to how the gore seams are modeled, causing a nearly 20% drop when they are modeled as stiffer radials, and to deceleration, causing an additional 10% reduction. Results show improved agreement with flight test measurements when the FSI model is closest to the as-built configuration and physics experienced in flight: with the confluence fitting, viscous effects, the gore seams modeled as stiff radials, and including deceleration.

ESM↗

Modeling laser-wakefield accelerators using the time-averaged ponderomotive approximation in a Lorentz boosted frame

Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.

Laser Wakefield Acceleration↗

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗

Advances in the simulation of space-charge effects

The nonlinear space-charge effects are an important topic in high intensity accelerator beam dynamics and have been extensively studied using macroparticle tracking simulations. Here, we report on recent advances in the simulation of space-charge effects using a symplectic multiparticle space-charge model. The transverse space-charge limit was explored for a periodic focusing and defocusing lattice. The artificial numerical emittance growth in the macroparticle space-charge simulation was analyzed using a one-dimensional model.

47 OTHER INSTRUMENTATION↗

Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.

43 PARTICLE ACCELERATORS↗

Pandemic drugs at pandemic speed: infrastructure for accelerating COVID-19 drug discovery with hybrid machine learning- and physics-based simulations on high-performance computers

The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods, in this case, developed for linear accelerators, and physics-based methods. The two in silico methods, each have their own advantages and limitations which, interestingly, complement each other. Here, we present an innovative infrastructural development that combines both approaches to accelerate drug discovery. The scale of the potential resulting workflow is such that it is dependent on supercomputing to achieve extremely high throughput. We have demonstrated the viability of this workflow for the study of inhibitors for four COVID-19 target proteins and our ability to perform the required large-scale calculations to identify lead antiviral compounds through repurposing on a variety of supercomputers.

97 MATHEMATICS AND COMPUTING↗

Microscopic Processes On Radiation from Accelerated Particles in Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The jitter'' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

Supervised paralled-in-time algorithm for long-time Lagrangian simulations of stochastic dynamics: Application to hydrodynamics

Lagrangian particle methods based on detailed atomic and molecular and nanoscale systems. However, the maximum time step is limited by the smallest oscillation period of the fastest atomic motion, rendering long-time simulations very expensive. To resolve this bottleneck, we propose a supervised parallel-in-time algorithm for stochastic dynamics to accelerate long-time Lagrangian particle simulations. Our method is inspired by bottom-up coarse-graining projections that yield mean-field hydrodynamic behavior in the continuum limit.

Multiscale modeling, Lagrangian method, Parallel-i↗

Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries

Lithium plating is one of the main concerns for cell durability, as cyclable lithium and cell capacity are decreased during the process. Thermal and mechanical loads can influence cell behavior, including lithium plating driving forces, as suggested by the Arrhenius equation and the Bruggeman relation, respectively. Hence, a 1D and a 3D model are employed to study the thermal and mechanical effects on lithium plating in a lithium-ion battery with a high areal capacity (4 mAh/cm 2 ). The 1D model assumes uniform mechanical and temperature profile, while the 3D model is capable of capturing the nonuniformity of stress and temperature in the cell. Further, the effects of mechanical deformation and different temperatures on the process of lithium plating are investigated in this study. With the help of the 1D charging simulations, it is observed that lower temperatures increase the risk of lithium plating and the presence of mechanical deformations can accelerate its initiation. Furthermore, 3D charging simulations indicated that folds and boundaries of jellyrolls are more susceptible to lithium plating due to the higher stress concentration and lower temperature, respectively.

25 ENERGY STORAGE↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a Green Parabolic Trough Collector

This report describes in detail a research project that assesses the potential of a parabolic trough solar collector with a lumber construction to meet ambitious performance, cost, and durability goals. Instead of exploring lumber as a general question for solar collectors, the project focuses on the design optimization, laboratory testing, outdoor testing, and economic modeling of a particular scale of parabolic trough, namely, the 6m aperture and 12m long module that has become the standard in the utility-scale power plant market. As will be appreciated in the report, it is also a scale that matches well the capabilities of lumber as a structural material and the commonly available cross sections. The design and its analysis extend far beyond traditional considerations for wooden structures in the built environment, and includes a computationally efficient coupled structural-optical simulation and the empirical characterization of joint stiffness at small deformations. Laboratory testing is used to characterize cyclic joint behavior, accelerated weathering, and to validate the simulation used as the basis for building a full-scale prototype. The construction and outdoor testing of this prototype is described. A total installed cost model of the collector is presented and the impact of the cost is seen in the context of a full plant techno-economic model. The results of this model in combination with the experience with the first prototype suggest that with improvements to the design and assembly process, this collector type has the potential to lower the levelized cost of heat for industrial solar steam plants used for process heating or cogeneration.

14 SOLAR ENERGY↗

On the high-energy gamma-ray signature of cosmic-ray sources

Monte Carlo simulations of the gamma-ray emission from hypothetical cosmic-ray sources are performed. Sources which might correspond to acceleration by supernova shocks in 'average' interstellar conditions and deep within giant molecular clouds are considered. The consequences of dropping the common assumption that the cosmic-ray spectrum at the sources is the same as that observed at earth are examined. Spectral effects which can be related to the depth of the material shroud and the population of accelerated particles are explored using these simulations and are described. The results are compared with the COS B catalog of gamma-ray sources, and the implications for the underlying particle populations and source mechanisms are discussed.

Ormes, J. F.↗

Coupler Studies for PBG Fiber Accelerators

Photonic band gap (PBG) fiber with hollow core defects are being designed and fabricated for use as laser driven accelerators because they can provide gradients of several GeV/m for picosecond pulse lengths. We expect to produce fiber down to {lambda} = 1.5-2.0 {micro}m wavelengths but still lack a viable means for efficient coupling of laser power into such structures due to the very different character of the TM-like modes from those used in the telecom field and the fact that the defect must function as both a longitudinal waveguide for the accelerating field and a transport channel for the particles. We discuss the status of our work in pursuing both end and side coupling. For both options, the symmetry of these crystals leads to significant differences with the telecom field. Side coupling provides more options and appears to be preferred. Our goals are to test gradients, mode content and coupling efficiencies on the NLCTA at SLAC. While there are many potential types of fiber based on very different fabrication methods and materials we will concentrate on 2D axisymmetric glass with hexagonal symmetry but will discuss several different geometries including 2D and 3D planar structures. Since all of these can be fabricated using modern techniques with a variety of dielectric materials they are expected to have desirable optical and radiation hardness properties. Thus, we expect a new generation of very high gradient accelerators that extends the Livingston-Panofsky chart of exponential growth in energy vs. time at greatly reduced costs. For illustration, Fig.1 shows a simulation of our first engineered fiber with an accelerating mode expected near 7.3 {micro}m that is now ready to test on the NLCTA. In this example, one sees the uniform longitudinal accelerating field in the central defect as first shown by Lin3 together with a hexagonal array of surrounding hot spots. Contrary to what one expects from the telecom field, Ng et al. have shown4 that the ideal end-coupling scheme for this structure appears as shown in Fig. 2 with a six-fold array of laser spots focused inside the end of the fiber. While convenient for an on-axis particle beam, this is inconvenient for the laser drive field as well as the tolerances it places on the end cleave of the fiber. The importance of the crystal symmetry is clearly shown so that one might expect side coupling to reflect a similar pattern which we find that it does unless the hexagonal symmetry is perturbed sufficiently. This can be done in several ways and will be discussed further.

36 MATERIALS SCIENCE↗

The mechanisms of electron heating and acceleration during magnetic reconnection

The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell simulations with non-zero guide fields so that electrons remain magnetized. In this regime, electric fields parallel to B accelerate particles directly, while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection, while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20% of the magnitude of the reconnecting component), the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases, the heating by the gradient B drift is negligible in magnitude. It produces net cooling because the conservation of the magnetic moment and the drop of B during reconnection produce a decrease in the perpendicular electron energy. Heating by the curvature drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. This suggests that acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The International Linear Collider: A Global Project

The International Linear Collider (ILC) is now under consideration as the next global project in particle physics. In this report, we review of all aspects of the ILC program: the physics motivation, the accelerator design, the run plan, the proposed detectors, the experimental measurements on the Higgs boson, the top quark, the couplings of the W and Z bosons, and searches for new particles. We review the important role that polarized beams play in the ILC program. The first stage of the ILC is planned to be a Higgs factory at 250 GeV in the centre of mass. Energy upgrades can naturally be implemented based on the concept of a linear collider. We discuss in detail the ILC program of Higgs boson measurements and the expected precision in the determination of Higgs couplings. We compare the ILC capabilities to those of the HL-LHC and to those of other proposed e + e - Higgs factories. We emphasize throughout that the readiness of the accelerator and the estimates of ILC performance are based on detailed simulations backed by extensive RandD and, for the accelerator technology, operational experience.

43 PARTICLE ACCELERATORS↗

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

An Updated Simulation of the Booster Neutrino Beam

For accelerator neutrino experiments, an accurate prediction of the incoming neutrino flux is crucial for reducing uncertainties for all physics measurements. In this exciting period for the Short-Baseline Neutrino program at Fermilab, with far detector (ICARUS) already operating and the near detector (SBND) nearing operation, an updated flux model for the Booster Neutrino Beam (BNB) is presented. The BNB delivers 8 GeV protons to a beryllium target, subsequently producing neutrinos from the decay of the secondary beam of hadrons. A updated Monte-Carlo simulation of beam production in GEANT4 has been created, allowing predictions to be made for detectors with different baselines, offsets and sizes. This new simulation contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, with a view to storing all resulting particles - including neutral mesons - from the proton-Beryllium scatter to allow the study of exotic BSM scenarios. Results are presented, with comparisons to the original flux simulated for the MiniBooNE experiment.

Paton, Josephine↗