SEARCH · Engineering Papers

Results for “particle methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Evaluation of Candidate Theranostics for 227 Th/ 89 Zr Paired Radioimmunotherapy of Lymphoma

227 Th is a promising radioisotope for targeted α-particle therapy. It produces 5 α-particles through its decay, with the clinically approved 223 Ra as its first daughter. There is an ample supply of 227 Th, allowing for clinical use; however, the chemical challenges of chelating this large tetravalent f-block cation are considerable. Using the CD20-targeting antibody ofatumumab, we evaluated chelation of 227 Th 4+ for α-particle–emitting and radiotheranostic applications. Methods: We compared 4 bifunctional chelators for thorium radiopharmaceutical preparation: S-2-(4-Isothiocyanatobenzyl)-1,4,7,10-tetraazacyclododecane tetraacetic acid (p-SCN-Bn-DOTA), 2-(4-isothicyanatobenzyl)-1,2,7,10,13-hexaazacyclooctadecane-1,4,7,10,13,16-hexaacetic acid (p-SCN-Bn-HEHA), p-isothiacyanatophenyl-1-hydroxy-2-oxopiperidine-desferrioxamine (DFOcyclo*-p-Phe-NCS), and macrocyclic 1,2-HOPO N-hydroxysuccinimide (L804-NHS). Immunoconstructs were evaluated for yield, purity, and stability in vitro and in vivo. Tumor targeting of the lead 227 Th-labeled compound in vivo was performed in CD20-expressing models and compared with a companion 89 Zr-labeled PET agent. Results: 227 Th-labeled ofatumumab-chelator constructs were synthesized to a radiochemical purity of more than 95%, excepting HEHA. 227 Th-HEHA-ofatumumab showed moderate in vitro stability. 227 Th-DFOcyclo*-ofatumumab presented excellent 227 Th labeling efficiency; however, high liver and spleen uptake was revealed in vivo, indicative of aggregation. 227 Th-DOTA-ofatumumab labeled poorly, yielding no more than 5%, with low specific activity (0.08 GBq/g) and modest long-term in vitro stability (<80%). 227 Th-L804-ofatumumab coordinated 227 Th rapidly and efficiently at high yields, purity, and specific activity (8 GBq/g) and demonstrated extended stability. In vivo tumor targeting confirmed the utility of this chelator, and the diagnostic analog, 89 Zr-L804-ofatumumab, showed organ distribution matching that of 227 Th to delineate SU-DHL-6 tumors. Conclusion: Commercially available and novel chelators for 227 Th showed a range of performances. The L804 chelator can be used with potent radiotheranostic capabilities for 89 Zr/ 227 Th quantitative imaging and α-particle therapy.

227Th↗

DOE OSTI2023

Ultrafast radiographic imaging and tracking: An overview of instruments, methods, data, and applications

Ultrafast radiographic imaging and tracking (U-RadIT) use state-of-the-art ionizing particle and light sources to experimentally study sub-nanosecond transients or dynamic processes in physics, chemistry, biology, geology, materials science and other fields. These processes are fundamental to modern technologies and applications, such as nuclear fusion energy, advanced manufacturing, communication, and green transportation, which often involve one mole or more atoms and elementary particles, and thus are challenging to compute by using the first principles of quantum physics or other forward models. One of the central problems in U-RadIT is to optimize information yield through, e.g. high-luminosity X-ray and particle sources, efficient imaging and tracking detectors, novel methods to collect data, and large-bandwidth online and offline data processing, regulated by the underlying physics, statistics, and computing power. We review and highlight recent progress in: (a.) Detectors such as high-speed complementary metal-oxide semiconductor (CMOS) cameras, hybrid pixelated array detectors integrated with Timepix4 and other application-specific integrated circuits (ASICs), and digital photon detectors; (b.) U-RadIT modalities such as dynamic phase contrast imaging, dynamic diffractive imaging, and four-dimensional (4D) particle tracking; (c.) U-RadIT data and algorithms such as neural networks and machine learning, and (d.) Applications in ultrafast dynamic material science using XFELs, synchrotrons and laser-driven sources. Hardware-centric approaches to U-RadIT optimization are constrained by detector material properties, low signal-to-noise ratio, high cost and long development cycles of critical hardware components such as ASICs. Interpretation of experimental data, including comparisons with forward models, is frequently hindered by sparse measurements, model and measurement uncertainties, and noise. Alternatively, U-RadIT make increasing use of data science and machine learning algorithms, including experimental implementations of compressed sensing. Machine learning and artificial intelligence approaches, refined by physics and materials information, may also contribute significantly to data interpretation, uncertainty quantification and U-RadIT optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Destructive Analysis of TRISO Particles: Crush/Burn/Leach Followed by Davies-Gray Titration and IDMS

The accurate accounting of nuclear materials is a cornerstone of international nuclear safeguards. One emerging challenge in this domain is the fabrication of TRIstructural ISOtropic (TRISO) particle fuels. Although these innovative fuel forms are critical for advanced reactor applications, their robust refractory ceramics and coating compositions present significant obstacles to destructive analysis (DA) methods. Ensuring full and quantitative recovery from these particles is essential for accurate mass accountancy. The current study was initiated to address these challenges, first by validating a previously established destructive method developed by Oak Ridge National Laboratory (ORNL) for the quantitative recovery of uranium from TRISO particles and then following that process with uranium content determination through isotope dilution mass spectrometry (IDMS) and Davies-Gray titration. This study expands on the scope of a digestive method that was developed under the Advanced Gas Reactor Fuel Development and Qualification program and is currently implemented in both the Coated Particle Fuel Development Laboratory and Irradiated Fuels Examination Laboratory at ORNL. The success of the previous Advanced Gas Reactor work relied on developing a DA method to evaluate the fabrication process and reactor experiments. The methodology described in this report was designed to rigorously investigate the efficacy of the crush/burn/leach sample preparation of TRISO particles; it aims to quantify uranium recovery while also assessing the effects of TRISO constituents (e.g., silicon and zirconium) on analytical precision and accuracy. By comparing the results from the titration method and IDMS, we sought to determine whether existing analytical procedures accepted by the International Atomic Energy Agency (IAEA) could be effectively translated to TRISO fuel forms. The team employed an approach that involved processing replicate TRISO samples, optimizing the milling (i.e., crushing) step and performing serial leaches. The elemental composition of the analytical samples was examined to prepare for interference studies in the second year of this project. The integration of gamma spectrometry to verify residual uranium activity further strengthened the validation. Statistical methods were applied to the collected data to evaluate the uncertainties arising from sampling, sample preparation and uranium quantification. These uncertainties were then compared to the IAEA’s international target values (ITVs). Additional data collected in upcoming project work will strengthen the uncertainty estimates. Ultimately, it is hoped that this project will contribute materially to the body of work related to characterization of TRISO based fuels for the purpose of material accountancy and its applications to international safeguards.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Melting LHC detectors: a novel search for stopped long-lived particles

Particles at the TeV scale with lifetimes of a year or longer could have been abundantly produced at the LHC yet escaped detection because of backgrounds, and could still be trapped within detector materials. With gluinos in split-supersymmetry as a working example, we show that these trapped particles can be recovered from detector materials once prepared in liquid form, for example, by melting silicon detectors, extracting liquid argon from the electromagnetic calorimeter, or constructing a large water pool near ATLAS or CMS. These liquid samples can then be processed using iterative centrifugation followed by mass spectrometry, enabling single-particle sensitivity in macroscopic samples. This method can potentially discover gluinos up to 3 TeV in mass at the HL-LHC. It can also improve upon existing limits for other long-lived particles. For example, it can discover the stop up to 2 TeV, and will also be sensitive to integer-charged particles in the TeV range.

Gonski, Julia L. [SLAC]↗

DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗

DOE OSTI2024

Lie–Poisson Neural Networks (LPNets): Data-based computing of Hamiltonian systems with symmetries

An accurate data-based prediction of the long-term evolution of Hamiltonian systems requires a network that preserves the appropriate structure under each time step. Every Hamiltonian system contains two essential ingredients: the Poisson bracket and the Hamiltonian. Hamiltonian systems with symmetries, whose paradigm examples are the Lie–Poisson systems, have been shown to describe a broad category of physical phenomena, from satellite motion to underwater vehicles, fluids, geophysical applications, complex fluids, and plasma physics. The Poisson bracket in these systems comes from the symmetries, while the Hamiltonian comes from the underlying physics. We view the symmetry of the system as primary, hence the Lie–Poisson bracket is known exactly, whereas the Hamiltonian is regarded as coming from physics and is considered not known, or known approximately. Using this approach, we develop a network based on transformations that exactly preserve the Poisson bracket and the special functions of the Lie–Poisson systems (Casimirs) to machine precision. We present two flavors of such systems: one, where the parameters of transformations are computed from data using a dense neural network (LPNets), and another, where the composition of transformations is used as building blocks (G-LPNets). We also show how to adapt these methods to a larger class of Poisson brackets. We apply the resulting methods to several examples, such as rigid body (satellite) motion, underwater vehicles, a particle in a magnetic field, and others. The methods developed in this paper are important for the construction of accurate data-based methods for simulating the long-term dynamics of physical systems.

Data-based modeling↗

DOE OSTI2020

Pileup mitigation at CMS in 13 TeV data

With increasing instantaneous luminosity at the LHC come additional reconstruction challenges. At high luminosity, many collisions occur simultaneously within one proton-proton bunch crossing. The isolation of an interesting collision from the additional “pileup” collisions is needed for effective physics performance. In the CMS Collaboration, several techniques capable of mitigating the impact of these pileup collisions have been developed. Such methods include charged-hadron subtraction, pileup jet identification, isospin-based neutral particle “δβ” correction, and, most recently, pileup per particle identification. This paper surveys the performance of these techniques for jet and missing transverse momentum reconstruction, as well as muon isolation. The analysis makes use of data corresponding to 35.9 fb−1 collected with the CMS experiment in 2016 at a center-of-mass energy of 13 TeV. The performance of each algorithm is discussed for up to 70 simultaneous collisions per bunch crossing. Significant improvements are found in the identification of pileup jets, the jet energy, mass, and angular resolution, missing transverse momentum resolution, and muon isolation when using pileup per particle identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Stability of the particle-hole Pfaffian state and the 5 2 -fractional quantum Hall effect

Here we present a method for the exact construction of the fully particle-hole symmetric Pfaffian (PH-Pfaffian) ground state and its charged excitations on a sphere. We adopt the Moore-Read state, but with a nonholomorphic pairing component as in previous studies, and project it to the lowest Landau level. We study the energetics as well as other properties of these states and find that in a pure system interacting with the Coulomb forces the PH-Pfaffian cannot compete with either the Moore-Read state or its particle-hole conjugate, the anti-Pfaffian state, as an explanation for the 5/2 effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Ammonium-coordinated exchanger (ACE) functionalized silica sorbents for recovering/removing aqueous anionic contaminants

There are limited studies of functionalized silica anion exchange sorbents used for critical/heavy metal recovery/removal relative to polymeric and other inorganic materials. This work features ammonium-coordinated exchanger (ACE) anion exchange particle sorbents prepared by either acid-washing epoxy-crosslinked polyethylenimine (PEI) hydrogen bonded within/to a silica particle sorbent (two-step method) or reacting a di-chlorinated crosslinker, α,α-dichloro-p-xylene (DPX), with PEI within silica (single-step method). Energy dispersive X-ray spectroscopy (EDS) and infrared spectroscopy confirmed the presence of -NH 2 + ···Cl - and -NH 3 + ···Cl - groups, which removed oxyanionic species –arsenate, selenate, chromate, sulfate, phosphate, and nitrate– plus bromide from ideal solutions, authentic acid mine drainage (AMD), and authentic flue gas desulfurization (FGD) wastewater. Affinity of the anions for ACE varied across single- and mixed-element solutions. However, affinity for CrO 4 2- was among the highest in both cases. Total anion uptake by the optimized ACE, PEI-E3-HCl_1.1, reached 1.2 mmol anion/g-sorb. (0.56 mmol CrO 4 2- /g), or ∼2.3 mmol negative charge/g-sorb. This was close to the 0.52 mmol CrO 4 2- /g of a commercial anion exchange resin. Near-consistent removal of 20–80 % of each anion from FGD during an eight-cycle adsorption-desorption (1 M NaCl) test predicted good ACE viability for testing under practical conditions at larger scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Extended range X-ray pair distribution functions

Here we describe a dual detector system for high-energy X-ray, simultaneous, small and wide-angle X-ray scattering (SAXS and WAXS), designed to extract extended-range pair distribution functions (ER-PDF) for disordered materials. The hardware and software provides continuous reciprocal space coverage over atomic to nanometer length-scales. Details of the varying resolution, splicing of data and normalization on an absolute scale are outlined. In addition, the combination of SAXS and WAXS theory is considered with a view to enabling a direct Fourier transformation of the structure factor spanning multiple length-scales into real space. Important distinctions between the ER-PDF and the pair distance distribution function (PDDF) representations are demonstrated. It is shown that when the SAXS intensity in the structure factor, S(Q), is similar to the WAXS intensity, the contributions to the ER-PDF are minimal. However, when the SAXS S(Q) intensities are substantially stronger than the WAXS, the ER-PDF can provide important structural information on the local, intermediate and nanometer length-scales. Notably, the ER-PDF method provides direct information on particle sizes and their density distributions, overcoming the limitations of PDDF analysis for densely packed systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Highly Durable and Active Pt/Sb-Doped SnO 2 Oxygen Reduction Reaction Electrocatalysts Produced by Atomic Layer Deposition

Platinum supported on mixed-metal oxides (MMOs) are a class of active and durable cathode catalysts for proton exchange membrane fuel cell (PEMFC) due to a combination of the high oxidative stability of the supports and strong-metal-support interactions (SMSI) that enable them to exceed the activity of Pt/C. Herein, we solve a significant remaining challenge with Pt/MMO systems, namely the relatively low surface area and porosity. This is achieved by dispersing nearly-uniform Pt clusters using atomic layer deposition (ALD) on highly conductive (6.2 S/cm) and stable antimony doped tin dioxide (ATO) support. ALD-Pt/ATO exhibited significantly higher electrochemically active surface area (ECSA) (74 m 2 /g) and oxygen reduction reaction (ORR) catalytic activity (102 mA/mg Pt at 0.9 V vs. RHE) compared to Pt/ATO synthesized using ethylene glycol (ECSA=31 m 2 /g Pt , mass activity=52 mA/mg Pt at 0.9 V vs. RHE) and formic acid reduction methods (ECSA=28 m 2 /g Pt , mass activity=46 mA/mg Pt at 0.9 V vs. RHE). Further characterization showed that wet chemical methods resulted in poorer Pt particle dispersion, poor control over Pt particle size distribution and chemical degradation of the support (during Pt deposition). Given the near-ideal Pt particle size distribution of the ALD-Pt/ATO, particle size growth and loss of ECSA was found to be minimal over the course of rigorous potential cycling. Thus, after 10,000 potential cycles between 1V and 1.5V vs. RHE, ALD-Pt/ATO and other Pt/ATOs were found to retain 100% of their initial ECSA compared to 57.6% retention for Pt/C. Upon testing in a H 2 /air PEMFC, following 1,000 potential cycles, the change in ALD-Pt/ATO performance was negligible while Pt/C exhibited a 68.2% loss of initial peak power density. Thus, ALD-Pt/ATO is an active and highly durable ORR electrocatalyst in PEMFCs under start-up-shut down conditions.

25 ENERGY STORAGE↗

DOE OSTI2024

Tutorial: The equations of electron emission and their evaluation

Electron emission and transport through and over potential barriers is an essential process requiring modeling and simulation to meet the design needs and characterization of an exceedingly broad range of technologically important devices and processes. The simulation and description of thermal, field, and photoemission, and the related concerns of space–charge affected electron flow, often make use of specialized formulations developed in the early days of quantum mechanics. Advancements in the utilization of electron sources and particularly the simulation of devices and applications using advanced particle-in-cell and trajectory methods for beam optics codes create a strong need for a pedagogical account of the emission models to ensure correct numerical evaluation of their equations. This Tutorial starts from simple phenomenological accounts and progressively builds to comprehensive models emphasizing straightforward and often rapid calculation. It recommends formulations to supplant the canonical Richardson–Laue–Dushman (thermal), Fowler–Nordheim (field), Fowler–DuBridge (photo), and Baroody (secondary) equations and provides a useful formulation of space–charge affected flow commonly described by the Child–Langmuir relation that takes into account cathode dependence on surface field.

Physics↗

DOE OSTI2024

Cooper minima in high-Z atoms: effects of correlation and relativity on np photoionization

Abstract Photoionization dipole transition matrix elements pass through a zero or attain a minimum that leaves imprints on photoionization parameters like the cross-section, angular distribution asymmetry parameter, phase shift, and photoionization time delay. This minimum is commonly known as the ‘Cooper minimum’ (CM). The CM, in general, is strongly affected by relativistic and correlation effects. Previous works investigated CM in the 6pand 5psubshell photoionization up toZ= 100 using the single-particle Dirac-Slater (DS) method. The present work extends the earlier work toZup to 120 using more accurate methods; Dirac–Hartree–Fock (DHF) which includes the relativistic effects and exchange correlations, and the relativistic random phase approximation (RRPA) which includes both initial and final state electron-electron correlations along with relativistic effects. In addition to the study of photoionization from the 6pand 5psubshells, the 4psubshell has also been investigated in the present work. To demonstrate the prominent effects in the high-Z atoms, Rn (Z= 86), Ra (Z= 88), No (Z= 102), Cn (Z = 112), Og (Z= 118), and Ubn (Z= 120) are investigated.

Optics↗

DOE OSTI2025

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↗

DOE OSTI2022

Variational Adiabatic Gauge Transformation on Real Quantum Hardware for Effective Low-Energy Hamiltonians and Accurate Diagonalization

Effective low-energy theories represent powerful theoretical tools to reduce the complexity in modeling interacting quantum many-particle systems. However, common theoretical methods rely on perturbation theory, which limits their applicability to weak interactions. Here we introduce the Variational Adiabatic Gauge Transformation (VAGT), a nonperturbative hybrid quantum algorithm that can use nowadays quantum computers to learn the variational parameters of the unitary circuit that brings the Hamiltonian to either its block-diagonal or full-diagonal form. If a Hamiltonian can be diagonalized via a shallow quantum circuit, then VAGT can learn the optimal parameters using a polynomial number of runs. The accuracy of VAGT is tested through numerical simulations, as well as simulations on Rigetti and IonQ quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗