Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1964

Dynamic Model Investigation of the Landing Characteristics of a Manned Spacecraft

Investigations were made to study the water-landing and certain grounds-surface landing characteristics of a Gemini spacecraft model. The water landing experiments were made by simulating paraglider and parachute letdowns with two 1/6- scale model configurations. Parameters included various combinations of attitude, horizontal speed, vertical speed, and landing skids extended and retracted. Investigations were made in calm water and in waves. The paraglider landings at horizontal speeds of 63 feet per second (19.8 m/sec) which resulted in a noseover or tumbling shortly after initial water contact. The maximum longitudinal acceleration of the model in calm water was about 14g units, and the maximum angular acceleration was 66 radians per second squared. In the parachute landings with the heat shield forward, the model skidded along the water surface on the heat shield. Parachute landings with the small end forward resulted in behavior similar to that of the paraglider landings. The ground-surface landings were made with a 1/3-scale model by simulating a parachute letdown with braking rockets, which were fired prior to touchdown to dissipate vertical velocity. In these landings, control of timing and aligning the rockets on the model was very critical, and violent behavior resulted when either rocket alignment or timing was in error. In the landings that were correctly controlled, the model either remained upright or slowly rolled over on its side.

EXPERIMENTS↗
DOE OSTI2026

A generative machine learning model for designing metal hydrides applied to hydrogen storage

Developing new metal hydrides is a critical step toward efficient hydrogen storage in carbon-neutral energy systems. However, existing materials databases, such as the Materials Project, contain a limited number of well-characterized hydrides, which constrains the discovery of optimal candidates. This work presents a framework that integrates causal discovery with a lightweight generative machine learning model to generate novel metal hydride candidates that may not exist in current databases. Using a dataset of 450 samples (270 training, 90 validation, and 90 testing), the model generates 1000 candidates. After ranking and filtering, six previously unreported chemical formulas and crystal structures are identified, four of which are validated by density functional theory simulations and show strong potential for future experimental investigation. Overall, the proposed framework provides a scalable and time-efficient approach for expanding hydrogen storage datasets and accelerating materials discovery.

generative model↗
DOE OSTI2023

Preliminary Design of a LAMP DTL

A preliminary Drift Tube Linac (DTL) layout was designed to create an algorithm for developing the conceptual design of the proposed LANSCE Accelerator Modernization Project (LAMP) final section: the proton linear accelerator from 3 MeV to 100 MeV. Preceding reports describe the proposed layouts of the LEBT, the RFQ, and the MEBT subsections of the linac. Initial estimates of the needed RF power for the DTL are also presented in the report. Present report describes the initial layout of the DTL and the first longitudinal beam dynamics results of the simulations n simplified models. The follow-up reports will include the transverse focusing scheme and beam dynamics in details. Present report in the last section will only describe the planned focusing scheme, based on the existing LANCE DTL.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Proton beam dynamics in bare IOTA with intense space-charge

We are commissioning a 2.5-MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MADX simulation codes. Our findings reveal agreement within a factor of 2 in core phase-space density predictions up to 100 synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗
NASA NTRS1989

Dynamics and control of the tether elevator/crawler system

This paper investigates the dynamics and acceleration levels of a new tethered system for micro- and variable-gravity applications. The system consists of two platforms tethered on opposite sides to the Space Station. A fourth platform, the elevator, is placed in between the Space Station and the upper platform. Variable-g levels on board the elevator are obtained by moving this facility along the upper tether, while microgravity experiments are carried out on board the Space Station. By controlling the length of the lower tether the position of the system center of mass can be maintained on board the Space Station despite variations of the system's distribution of mass. The paper illustrates the mathematical model, the environmental perturbations and the control techniques which have been adopted for the simulation and control of the system dynamics. Two sets of results from two different simulation runs are shown. The first set shows the system dynamics and the acceleration spectra on board the Space Station and the elevator during station-keeping. The second set of results demonstrates the capability of the elevator to attain a preselected g-level.

Lorenzini, E. C.↗
NASA NTRS1988

Acceleration levels on board the Space Station and a tethered elevator for micro and variable-gravity applications

This paper investigates the dynamics and acceleration levels of a new tethered system for micro and variable-gravity applications. The system consists of two platforms tethered on opposite sides to the Space Station. A fourth platform, the elevator, is placed in between the Space Station and the upper platform. Variable-g levels on board the elevator are obtained by moving this facility along the upper tether, while micro-g experiments are carried out on board the Space Station. By controlling the length of the lower tether the position of the system CM can be maintained on board the Space Station despite variations of the station's distribution of mass. The paper illustrates the mathematical model, the environmental perturbations and the control techniques which have been adopted for the simulation and control of the system dynamics. Two sets of results from two different simulation runs are shown. The first set shows the system dynamics and the acceleration spectra on board the Space Station and the elevator during station-keeping. The second set of results demonstrates the capability of the elevator to attain a preselected g-level.

Lorenzini, E. C.↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
NASA NTRS2000

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. III. Responses after labyrinthectomy

The horizontal angular vestibuloocular reflex (VOR) evoked by high-frequency, high-acceleration rotations was studied in four squirrel monkeys after unilateral labyrinthectomy. Spontaneous nystagmus was measured at the beginning and end of each testing session. During the period that animals were kept in darkness (4 days), the nystagmus at each of these times measured approximately 20 degrees /s. Within 18-24 h after return to the light, the nystagmus (measured in darkness) decreased to 2.8 +/- 1.5 degrees /s (mean +/- SD) when recorded at the beginning but was 20.3 +/- 3.9 degrees /s at the end of the testing session. The latency of the VOR measured from responses to steps of acceleration (3,000 degrees /s(2) reaching a velocity of 150 degrees /s) was 8.4 +/- 0.3 ms for responses to ipsilesional rotations and 7.7 +/- 0.4 ms for contralesional rotations. During the period that animals were kept in darkness after the labyrinthectomy, the gain of the VOR measured during the steps of acceleration was 0.67 +/- 0.12 for contralesional rotations and 0.39 +/- 0.04 for ipsilesional rotations. Within 18-24 h after return to light, the VOR gain for contralesional rotations increased to 0.87 +/- 0.08, whereas there was only a slight increase for ipsilesional rotations to 0.41 +/- 0. 06. A symmetrical increase in the gain measured at the plateau of head velocity was noted after the animals were returned to light. The VOR evoked by sinusoidal rotations of 2-15 Hz, +/-20 degrees /s, showed a better recovery of gain at lower (2-4 Hz) than at higher (6-15 Hz) frequencies. At 0.5 Hz, gain decreased symmetrically when the peak amplitude was increased from 20 to 100 degrees /s. At 10 Hz, gain was decreased for ipsilesional half-cycles and increased for contralesional half-cycles when velocity was raised from 20 to 50 degrees /s. A model incorporating linear and nonlinear pathways was used to simulate the data. Selective increases in the gain for the linear pathway accounted for the recovery in VOR gain for responses at the velocity plateau of the steps of acceleration and for the sinusoidal rotations at lower peak velocities. The increase in gain for contralesional responses to steps of acceleration and sinusoidal rotations at higher frequencies and velocities was due to an increase in the contribution of the nonlinear pathway. This pathway was driven into cutoff and therefore did not affect responses for rotations toward the lesioned side.

NASA Discipline Neuroscience↗
DOE OSTI2026

Assessment and Improvement Proposals for High-Pressure Rinsing Processes during Cleanroom Preparation of Superconducting Radio-Frequency Cavities for Particle Accelerators

The goal of this internship is to develop a fluid dynamics characterization of the high-pressure flow used for cleaning the internal surface of the SRF cavities. Initially, the High Pressure Rinsing (HPR) process is introduced, focusing on the problems of the SRF cavities and how this method can be used to increase their acceleration efficiency and performance. The characterization of the problem is then presented, including an introduction to the multiphase flow modeling and the analysis of the Fermilab facility for the HPR. To accomplish the goal, a preliminary estimation of the variables was performed using MATLAB. Then, multiple 2D simulations were developed with Ansys Fluent, and the influence of the main parameter was investigated. The simulation results are then presented and analysed, with a focus on the possible future developments of the HPR process fluid dynamic characterization.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Proton Beam Dynamics in Bare IOTA with Intense Space-Charge

We are commissioning a 2.5~MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MAD-X simulation codes. Our findings reveal agreement within 30\% in core phase-space density predictions up to 100~synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Machine Learning for Predicting Multipactor Susceptibility in Planar RF Structures

Multipactor discharge is a persistent challenge in high-power microwave (HPM) and accelerator systems, where secondary electron avalanches can cause heating, vacuum degradation, and failure. This work presents the first supervised machine learning (ML) framework for multipactor prediction, trained on high-fidelity 3D Particle-in-Cell (PIC) simulation data in planar geometries. The model maps operational, geometric, and material-dependent secondary electron yield (SEY) parameters to the time-averaged electron growth rate, enabling rapid reconstruction of susceptibility charts. Among the models evaluated, tree-based ensemble methods such as Random Forest and Extra Trees demonstrate superior generalization to unseen materials compared to neural networks such as multilayer perceptron (MLP). Performance metrics, including Intersection over Union (IoU), Structural Similarity Index Measure (SSIM), and Pearson correlation, show close agreement with simulation benchmarks. Principal Component Analysis attributes generalization limits to material feature-space disjointedness.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Audacity of huge: overcoming challenges of data scarcity and data quality for machine learning in computational materials discovery

Machine learning (ML)-accelerated discovery requires large amounts of high-fidelity data to reveal predictive structure–property relationships. For many properties of interest in materials discovery, the challenging nature and high cost of data generation has resulted in a data landscape that is both scarcely populated and of dubious quality. Data-driven techniques starting to overcome these limitations include the use of consensus across functionals in density functional theory, the development of new functionals or accelerated electronic structure theories, and the detection of where computationally demanding methods are most necessary. When properties cannot be reliably simulated, large experimental data sets can be used to train ML models. In the absence of manual curation, increasingly sophisticated natural language processing and automated image analysis are making it possible to learn structure–property relationships from the literature. Finally, models trained on these data sets will improve as they incorporate community feedback.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Study of solid molecular deuterium D 2 growth under gas pressure

The injection of high-speed cryogenic pellets made of frozen hydrogen-isotopes, represents to date the most effective method to fuel magnetically confined thermonuclear fusion plasmas. Additionally, the injection of very large pellets composed of cryogenic solid of some suitable impurity (typically a noble-gas such as H 2 , Ne, or H 2 /Ne, D 2 /Ne mixtures), shattered in relatively small fragments just before entering the plasma, seems to be the most promising method to reduce the damage risks for the plasma-facing components in case of a plasma disruption. This technology, known as "Shattered Pellet Injection" (SPI), allows to spread out the plasma energy and mitigate possible damage to the in-vessel components, as well as to densify the plasma to suppress the formation of runaway electrons, and/or dissipate their energy. Several techniques to produce and launch cryogenic pellets have been investigated in the past decades. "Pipe gun" injectors are reliable and relatively simple devices are still commonly used today. They make use of single- or two-stage pneumatic light-gas guns to accelerate the pellet at high speeds. In these injectors, the cryogenic pellets are formed “ in situ ” (i.e., inside the launching barrel), by de-sublimating them directly from the gas phase, i.e., at temperatures and pressures below those of the triple point. The simplest case is pure deuterium pellets ($\mathcal{T}$ < 18.7 K, $\mathcal{P}$ < 171.3 hPa). The production of good quality solid deuterium, capable of withstanding the mechanical stress during the acceleration of the pellets, is a key issue. To this end the phase transition of deuterium from gas to solid (and vice versa) is modeled with extensive molecular-dynamics (MD) simulations. Furthermore, the solid growth from the gas phase is simulated in an ample range of temperatures and pressures, to find the best compromise between growth velocity and mechanical properties of the resulting solid system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Measurement of exclusive 𝜋 + -argon interactions using ProtoDUNE-SP

We present the measurement of 𝜋 + -argon inelastic cross sections using the ProtoDUNE single-phase liquid argon time projection chamber in the incident 𝜋 + kinetic energy range of 500–800 MeV in multiple exclusive channels (absorption, charge exchange, and the remaining inelastic interactions). The results of this analysis are important inputs to simulations of liquid argon neutrino experiments such as the Deep Underground Neutrino Experiment and the Short Baseline Neutrino program at Fermi National Accelerator Laboratory. They will be employed to improve the modeling of final state interactions within neutrino event generators used by these experiments, as well as the modeling of 𝜋 + -argon secondary interactions within the liquid argon. This is the first measurement of 𝜋 + -argon absorption at this kinetic energy range as well as the first ever measurement of 𝜋 + -argon charge exchange.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Transformational Regional-Scale Earthquake Simulations with the DOE EarthQuake SIMulation Exascale Framework

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

58 GEOSCIENCES↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2022

Criticality in the Shallow-to-Deep Transition of Simulated Tropical Marine Convection

Idealized simulations of tropical, marine convection depict shallow, nonprecipitating cumuli located beneath the 0°C level transitioning into cumulonimbi that reach up to 12 km and higher. The timing of the transition was only weakly related to environmental stability, and 13 of the 15 simulations run with 5 different lapse-rate profiles had rain develop at nearly the same time after model start. The key quantity that apparently controlled deep convective formation was vertical acceleration inside cloudy updrafts between cloud base and the 0°C level. Below a critical value of updraft vertical acceleration, little rainfall occurred. Just as the domain-mean updraft acceleration reached the critical value, the first convection quickly grew to past 12 km altitude. Then, as acceleration increased above the critical value, rain rate averaged in the model domain increased quickly over about a 3-h-long period. The specific value of the critical updraft acceleration depended on how updrafts were defined and in what layer the acceleration was averaged; however, regardless of how criticality was defined, a robust relationship between domain-mean updraft vertical acceleration and rain rate occurred. Positive acceleration of updrafts below the 0°C level was present below 2.75 km and was largest in the 500 m above cloud base. However, the maximum difference between updraft and environmental temperatures occurred between 2 and 3 km. The domain-mean Archimedean buoyancy of updrafts relative to some reference state was a poor predictor for domain-mean rain rate. In conclusion, the exact value of the critical updraft acceleration likely depends on numerous other factors that were not investigated.

54 ENVIRONMENTAL SCIENCES↗