SEARCH · Engineering Papers

Results for “phase space dynamics”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗

DOE OSTI2026

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

DOE OSTI2026

Third integer resonant extraction transit time simulation studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab]↗

DOE OSTI2020

Coordinated Multiscale Modeling and Synthesis of Novel Nanostructured Composite Membranes for Solar Fuels Generation (Final Report)

Generation of carbon-based liquid fuels from reduction of CO 2 using solar energy requires the development of improved membranes that provide good ionic conductivity and mechanical properties while minimizing the crossover of gases and reduction products. Using closely coordinated multiscale modeling, synthesis and characterization studies we propose to investigate and develop novel nanostructured composite membranes with applications to conversion of carbon dioxide to storable chemical fuels. The envisioned proton exchange membranes (PEMs) will be based on a support of self- assembling polymer-modified nanoparticles with controllable porosity and tortuosity. The use of a nanoparticle support will provide control of transport and mechanical properties of the membranes. Tuning the copolymer modification of nanoparticles (polymer architectures, brush grafting density) will allow to optimize the ion (proton) transport via the Grotthuss mechanism through narrow hydrated channels while minimizing crossover of gases and reaction products due to the high density of the grafted brush and use of a non-ionic conducting glassy polymers to fill spaces between polymer-grafted nanoparticles. In Phase I synthesis of the copolymers and their attachment to the nanoparticle support will be guided by multiscale simulations that include atomistic molecular dynamics and continuum-level transport simulations. The former will include reactive and non-reactive atomistic simulations that will provide key insight into nanoscale polymer morphology and ion/molecular transport mechanisms need to optimize the copolymer structure. The latter will provide insight into the role of nanoparticle self-assembly on the global transport of ionic and molecular species. The proton conduction and CO 2 /methanol permeability of the membranes will be experimentally characterized as well as simulated and will be correlated with morphology and polymer structure. Prototypes of optimal membranes will be fabricated and characterized at the end of the project. Successful development of improved nanostructured composite PEMs and demonstration of the simulation-guided materials-by-design paradigm will result in our ability to design and synthesize membranes for a wide variety of solar fuels generation and related electrochemical applications.

14 SOLAR ENERGY↗

DOE OSTI2025

Third Integer Resonant Extraction Transit Time Simulation Studies

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

DOE OSTI2025

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Taming the Pseudoelastic Response of Nitinol Using Ion Implantation

Implantation of Ni 50.5 Ti 49.5 wire with 30 MeV Ni 6+ ions at doses (< 0.1 DPA) typically smaller than employed in the literature is shown to systematically alter the pseudoelastic response, with extrema in Berkovich nanoindentation load (+50%), hysteresis (–60%), and recoverable displacement (–19%) occurring at ~ 3.6 μm below the implantation surface. These extraordinary values are attributed to ~ 10 to 20 nm amorphous clusters that constrain the stress-induced B2-B19' phase transformation. This is substantiated by phase field simulations of crystalline-amorphous composites and molecular dynamics simulations of crystalline-vacancy cluster composites showing the spatial refinement of martensite caused by nm-scale defects. Finally, the results suggest that ion implantation may potentially expand the processing and performance space for NiTi, by creating amorphous defects at smaller length scales than dislocation substructures produced by conventional deformation processing.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Capture Cavities for the CW Polarized Positron Source Ce+BAF

The initial design of the capture cavities for a continuous wave (CW) polarized positron beam for the Continuous Electron Beam Accelerator Facility (CEBAF) up-grade at Jefferson Lab is presented. A chain of standing wave multi-cell copper cavities inside a solenoid channel are selected to capture positrons in CW mode. The cavity shunt impedance is surveyed by tuning the cavity geometry while considering accommodating large phase space distribution positron beams with large beam pipe radius while ensuring a large enough passband mode separation. The RF field wall loss power and maximum wall loss power density are considered in cavity and waveguide design. A range of design parameters are given for larger system optimization when the capture cavities are considered together with thermal calculation and beam dynamics in next phase of work.

Wang, S.↗

DOE OSTI2024

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Dynamic Manipulation in Piezoresponse Force Microscopy: Creating Nonequilibrium Phases with Large Electromechanical Response

We discuss how domain walls and topological defects in ferroelectric materials have emerged as a powerful tool for functional electronic devices including memory and logic. Similarly, wall interactions and dynamics underpin a broad range of mesoscale phenomena ranging from giant electromechanical responses to memory effects. Exploring the functionalities of individual domain walls, their interactions, and controlled modifications of the domain structures is crucial for applications and fundamental physical studies. However, the dynamic nature of these features severely limits studies of their local physics since application of local biases or pressures in piezoresponse force microscopy induce wall displacement as a primary response. Here, we introduce an approach for the control and modification of domain structures based on automated experimentation, whereby real-space image-based feedback is used to control the tip bias during ferroelectric switching, allowing for modification routes conditioned on domain states under the tip. This automated experiment approach is demonstrated for the exploration of domain wall dynamics and creation of metastable phases with large electromechanical response.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Evaluation of the anisotropic grain boundaries and surfaces of α -U via molecular dynamics

Alloys based on uranium-zirconium are gaining renewed interest as fuels for the Versatile Test Reactor and a number of microreactor designs. Implementing metallic fuel in reactors creates the need for robust descriptive and predictive fuel performance modeling. Here, the current state of metallic fuel performance modeling relies on empirical equations derived from historical experiments, which may be unreliable when applied outside of their temperature, power, and composition phase space. One area where such data is lacking is the irradiation behavior of alpha -U, specifically tearing and porosity formation at the early stages of irradiation. While grain boundaries likely play a key role in this fuel behavior, relatively little is known about grain boundaries in alpha-U. Thus, we evaluate the grain boundary, surface energy, and work of adhesion of alpha -U utilizing molecular dynamics. Symmetric tilt grain boundaries (STGBs) are analyzed with the tilt plane oriented along each major crystallographic axis, for a total of eighty unique grain boundaries. The effect of temperature, tilt plane, and misorientation angle on interfacial energies are analyzed. The interfacial energies typically increase with temperature and there is significant variance as a function of misorientation angle, irrespective of the tilt plane. At 500 K, the average surface energy (1.23 J/m 2 ) is approximately 1.5 times the grain boundary energy (0.79 J/m 2 ), and the work of adhesion is approximately twice the grain boundary energy (1.68 J/m 2 ). Orientations for the likely formation of twins and likely failure planes are identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Self-assembly of a triply periodic continuous mesophase with Fddd symmetry in simple one-component liquids

Triply periodic continuous morphologies (networks) arising as a result of the microphase separation in block copolymer melts have so far never been observed self-assembled in systems of particles with spherically symmetric interaction. Here, we report a molecular dynamics simulation where two simple one-component liquids form upon cooling an equilibrium network with the Fddd space group symmetry. This complexity reduction in the liquid network formation in terms of the particle geometry and the number of components evidences the generic nature of this class of phase transition suggesting opportunities for producing these structures in a variety of new systems.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multi-Physics Topology Optimization for Dual Flow Heat Exchangers

Current heat exchanger (HX) designs are subjected to strict constraints pertaining to their manufacturability through conventional lengthy and expensive processes consisting of subtractive, forming and brazing processes. Therefore, existing HXs are typically restricted to the classical and rather simple geometries, such as plate-fin or shell-tube geometries, with lower efficiency and larger size. To increase HX efficiencies, novel designs are needed. Recent advances in topology optimization (TO) and additive manufacturing (AM) offer significant flexibility in design and fabrication of complex geometries, with potentially superior thermo-mechanical performance and high-power density. On the other hand, there is no commercially available software tool that addresses the combined coupling of structural, fluid dynamics, and thermal with multiple materials (hot and cold streams, solid, and void space) in TO of HXs. The goal of this project is to develop a multi-physics multi-phase TO approach to enhance the performance of HXs. The optimized design will be additively manufactured and tested in a physical environment to quantify its performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Status of the Spallation Neutron Source beam test facility and progress of beam dynamics studies

The Spallation Neutron Source (SNS) Beam Test Facility (BTF) supports the study of beam dynamics in the front end of a high power LINAC. The BTF combines a replica of the SNS front end, including nearly-identical ion source, RFQ and MEBT, with extensive phase space diagnostics and a FODO transport line. Diagnostic capabilities include direct measurement of 6D phase space distribution and detection of halo distributions to a sensitivity of greater than one part-per-million. The goal of on-going BTF studies is to demonstrate accurate particle-in-cell modeling of halo growth and evolution by leveraging unprecedented accuracy in the description of the initial beam distribution. This work is motivated by operational experience at the SNS, which currently operates with beam loss that cannot be described by any model. This paper summarizes progress in the BTF beam study program as well as diagnostics development and recent upgrades to the beamline configuration.

Ruisard, Kiersten↗

DOE OSTI2022

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A machine learning study on spinodal clumping in heavy ion collisions

Possible observables of baryon number clustering due to the instabilities occurring at a first order QCD phase transition are discussed. The dynamical formation of baryon clusters at a QCD phase transition can be described by numerical fluid dynamics, augmented with a gradient term and an equation of state with a mechanically unstable region. It is shown that the dynamical description of this phase transition, in nuclear collisions, will lead to the formation of dense baryon clusters at the phase boundary. State-of-the-art machine learning methods find that the coordinate space clumping leaves characteristic imprints on the spatial net density distribution in almost every event. On the other hand the momentum distributions do not show any clear event-by-event features. Lastly, it is shown that the 'third order' cumulant, the skewness, shows a peak at the beam energy where the system, created in the heavy ion collision, reaches the deconfinement phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Gold-gold luminosity increase in RHIC for a beam energy scan with colliding beam energies extending below the nominal injection energy

The Beam Energy Scan phase II (BES-II), performed in the Relativistic Heavy Ion Collider (RHIC) from 2019 to 2021, explored the phase transition between quark-gluon plasma and hadronic gas. BES-II exceeded the goal of a fourfold increase in the average luminosity over that achieved during Beam Energy Scan phase I (BES-I), at five gold beam energies: 9.8, 7.3, 5.75, 4.59, and 3.85 GeV / nucleon. This was accomplished by addressing several beam dynamics effects, including intrabeam scattering, beam-beam, space charge, beam instability, and field errors induced by superconducting magnet persistent currents. Some of these effects are especially detrimental at low energies. BES-II achievements are presented, and the measures taken to improve RHIC performance are described. These measures span the whole RHIC complex, including ion beam sources, injectors, beam lifetime improvements in RHIC, and operation with the world’s first bunched beam Low Energy RHIC electron Cooler (LEReC).

43 PARTICLE ACCELERATORS↗