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 487 records · Page 27

Assessment of Artemis-1 Pogo Flight Instrumentation System

As the space industry continues to strive for more efficient launch vehicles, it must rely on increasingly accurate predictive models. Verification of models typically requires physical testing. Flight data measurements offer the most real and therefore the most accurate data for model correlation. As NASA prepares for the inaugural launch of Space Launch System (SLS), Artemis-1, they must rely heavily on predictive system models to ensure flight safety. NASA has implemented a Development Flight Instrumentation (DFI) system in hopes of recovering useful flight data to aid in model correlation. Historically, some of the most important flight measurements are those that monitor the potentially destructive dynamic interaction of the structural and propellant modes – a phenomenon known as pogo. Pogo is a dynamic instability that can occur on a launch vehicle during any phase of ascent. During this investigation an end-to-end assessment of the Artemis-I Pogo-related DFI was performed to identify any obstacles inherent in the current instrumentation system which may prevent successful measurement of the data necessary to validate the current predictive models of the Main Propulsion System (MPS). Input drive signals were approximated and applied to a system-level SLS state-space model to derive predicted pressure and acceleration responses. These predicted responses were then fed through a simulation of the data acquisition process in order to recover predicted measurements. Finally, a mock Flight Data Analysis (FDA) was performed to assess the ability of these measurements to meet the Flight Test Objectives (FTO).

Adam Johnson↗

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗

A simple, flexible technique for RF cavity wake-field calculations

It is typical in the accelerator field to model machine components, especially RF cavities, as parallel RLC resonators. In the interest of simulating and diagnosing beam instabilities, knowledge of the time-domain voltage waveform over an equivalent resonator by a bunch current often proves useful. This waveform may be found by convolving the bunch current with the RLC resonator impulse response. While analytical and quasi-analytical expressions are available in this regime for common distributions such as the Gaussian, analogous results for less standard distributions are difficult to obtain using direct methods, which opens the door for the development of a more generalized technique. In this paper, a formulation is created that allows for the simple computation of the time-domain voltage waveform of and RLC resonator. The formulation uses the Cauchy Residue Theorem to extract the convolution result from the Fourier Domain, and it only requires that the current distribution Fourier Transform be holomorphic and known at one specific evaluation point. This greatly simplifies the computation of the time domain voltage for a large amount of bunch distributions both common and uncommon.

43 PARTICLE ACCELERATORS↗

Developing Multiphysics, Integrated, High-Fidelity, Massively Parallel Computational Capabilities for Fusion Applications Using MOOSE

As the need for fusion as a clean, sustainable, and abundant energy source grows internationally, so does the need for multiphysics, computational tools to model, study, and predict the complex interactions between plasma, materials, and engineering processes. These tools have a crucial role to play in solving scientific and engineering challenges and accelerating fusion energy deployment. To address these needs, modeling capabilities should enable massively parallel, multiphysics, fully integrated high-fidelity simulations of fusion systems. Additional attributes, such as being open source and modular while maintaining high software quality assurance standards will maximize impact by ensuring accessibility for all and wide acceptance, rapid expansion and development, as well as reliability, efficiency, and robustness. In this paper, we describe how the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, which has a track record of success in the fission space thanks to the attributes listed above, can be leveraged in the fusion energy field. We highlight key successes of the MOOSE application in the fission space and describe how MOOSE has been and is being applied to fusion applications in the United States---e.g., Tritium Migration Analysis Program, version 8 (TMAP8), MOOSE Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX)---and the United Kingdom---e.g., AURORA, Achlys, Apollo. These efforts aim to establish a suite of tools that can be further extended to accelerate fusion energy deployment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

OpenFerro v0.1.0

OpenFerro is a Python package for on-lattice atomistic dynamics simulation of ferroic materials. OpenFerro is based on JAX, a high-performance linear algebra package supporting auto-differentiation and GPU acceleration. OpenFerro is designed to minimize the effort required to build on-lattice Hamiltonian models, and to perform molecular dynamics (MD) and Landau-Lifshitz-Gilbert simulations. Unlike existing codes, OpenFerro provides a unified interface to model different types of local order parameters.

Xie, Pinchen [Lawrence Berkeley National Laborator↗

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗

Evaluation of Portable Acceleration Solutions for LArTPC Simulation Using Wire-Cell Toolkit

The Liquid Argon Time Projection Chamber (LArTPC) technology plays an essential role in many current and future neutrino experiments. Accurate and fast simulation is critical to developing efficient analysis algorithms and precise physics model projections. The speed of simulation becomes more important as Deep Learning algorithms are getting more widely used in LArTPC analysis and their training requires a large simulated dataset. Heterogeneous computing is an efficient way to delegate computationally intensive tasks to specialized hardware. However, as the landscape of compute accelerators quickly evolves, it becomes increasingly difficult to manually adapt the code to the latest hardware or software environments. A solution which is portable to multiple hardware architectures without substantially compromising performance would thus be very beneficial, especially for long-term projects such as the LArTPC simulations. In search of a portable, scalable and maintainable software solution for LArTPC simulations, we have started to explore high-level portable programming frameworks that support several hardware backends. In this paper, we present our experience porting the LArTPC simulation code in the Wire-Cell Toolkit to NVIDIA GPUs, first with the CUDA programming model and then with a portable library called Kokkos. Preliminary performance results on NVIDIA V100 GPUs and multi-core CPUs are presented, followed by a discussion of the factors affiecting the performance and plans for future improvements.

Yu, Haiwang↗

Frostwing Co-Operation in Aircraft Icing Research

e aerodynamic effects of Cold Soaked Fuel Frost have become increasingly significant as airworthiness authorities have been asked to allow it during aircraft take-off. The Federal Aviation Administration and the Finnish Transport Safety Agency signed a Research Agreement in aircraft icing research in 2015 and started a research co-operation in frost formation studies, computational fluid dynamics for ground de/anti-icing fluids, and de/anti-icing fluids aerodynamic characteristics. The main effort has been so far on the formation and aerodynamic effects of CSFF. To investigate the effects, a generic high-lift common research wind tunnel model and DLR-F15 airfoil, representing the wing of a modern jet aircraft, was built including a wing tank cooling system. Real frost was generated on the wing in a wind tunnel test section and the frost thickness was measured with an Elcometer gauge. Frost surface geometry was measured with laser scanning and photogrammetry. The aerodynamic effect of the frost was studied in a simulated aircraft take-off sequence, in which the speed was accelerated to a typical rotation speed and the wing model was then rotated to an angle of attack used at initial climb. Time histories of the lift coefficient were measured with a force balance. The experiments showed that depending on the ambient temperature the frost may evaporate/melt during the take-off sequence. Lift losses after rotation with CSFF contamination at ambient temperatures of 4˚ to 7˚C above freezing point were measured to be 4 to 5 % for roughness values, k/c, below 10(exp -3). For comparison, lift loss tests with typical anti-icing fluids were to roughly equal lift losses. This paper gives an overview of the performed activities.

Koivisto, Pekka↗

Earthquake response of head-mounted equipment in advanced nuclear reactors

The seismic response of safety-related equipment mounted on the head of an advanced reactor, including pumps, control rod drive mechanisms, and reactor monitoring devices, will affect the design and layout of many advanced reactors. High earthquake-induced accelerations in such equipment may challenge their seismic qualification and trigger the need for additional support framing on the reactor head. Base isolation is a design solution that can drastically reduce seismic demands on equipment. This article describes a set of earthquake-simulator experiments conducted on a scale-model of a base-isolated reactor vessel including four representations of head-mounted equipment, with frequencies spanning from 4.5 to 27 Hz. Dynamic responses of the head-mounted equipment, including displacements, accelerations, and strains, were measured in the experiments for three support conditions: conventional, and seismically isolated using single concave Friction Pendulum (SFP) bearings and triple Friction Pendulum (TFP) bearings. Seismic isolation was effective at reducing equipment responses (accelerations, displacements, and strains) with respect to those in the conventionally supported vessel across a range of seismic inputs. Companion numerical studies highlight the accuracy to be expected in the calculation of different response quantities for lightly damped equipment. The importance of characterizing damping in head-mounted, safety-related equipment through physical experiments to support design and risk assessment is made clear through the numerical simulations.

Engineering↗

Evaluation of Mid-Size Male Hybrid III Models for use in Spaceflight Occupant Protection Analysis

In an effort to improve occupant safety during dynamic phases of spaceflight, the National Aeronautics and Space Administration (NASA) has worked to develop occupant protection standards for future crewed spacecraft. One key aspect of these standards is the identification of injury mechanisms through anthropometric test devices (ATDs). Within this analysis, both physical and computational ATD evaluations are required to reasonably encompass the vast range of loading conditions any spaceflight crew may encounter. In this study the accuracy of publically available mid-size male HIII ATD finite element (FE) models are evaluated within applicable loading conditions against extensive sled testing performed on their physical counterparts. Methods: A series of sled tests were performed at the Wright Patterson Air force Base (WPAFB) employing variations of magnitude, duration, and impact direction to encompass the dynamic loading range for expected spaceflight. FE simulations were developed to the specifications of the test setup and driven using measured acceleration profiles. Both fast and detailed FE models of the mid-size male HIII were ran to quantify differences in their accuracy and thus assess the applicability of each within this field. Results: Preliminary results identify the dependence of model accuracy on loading direction, magnitude, and rate. Additionally the accuracy of individual response metrics are shown to vary across each model within evaluated test conditions. Causes for model inaccuracy are identified based on the observed relationships. Discussion: Computational modeling provides an essential component to ATD injury metric evaluation used to ensure the safety of future spaceflight occupants. The assessment of current ATD models lays the groundwork for how these models can be used appropriately in the future. Identification of limitations and possible paths for improvement aid in the development of these effective analysis tools.

Putnam, Jacob B.↗

Evaluation of Mid-Size Male Hybrid III Models for use in Spaceflight Occupant Protection Analysis

Introduction: In an effort to improve occupant safety during dynamic phases of spaceflight, the National Aeronautics and Space Administration (NASA) has worked to develop occupant protection standards for future crewed spacecraft. One key aspect of these standards is the identification of injury mechanisms through anthropometric test devices (ATDs). Within this analysis, both physical and computational ATD evaluations are required to reasonably encompass the vast range of loading conditions any spaceflight crew may encounter. In this study the accuracy of publically available mid-size male HIII ATD finite element (FE) models are evaluated within applicable loading conditions against extensive sled testing performed on their physical counterparts. Methods: A series of sled tests were performed at the Wright Patterson Air force Base (WPAFB) employing variations of magnitude, duration, and impact direction to encompass the dynamic loading range for expected spaceflight. FE simulations were developed to the specifications of the test setup and driven using measured acceleration profiles. Both fast and detailed FE models of the mid-size male HIII were ran to quantify differences in their accuracy and thus assess the applicability of each within this field. Results: Preliminary results identify the dependence of model accuracy on loading direction, magnitude, and rate. Additionally the accuracy of individual response metrics are shown to vary across each model within evaluated test conditions. Causes for model inaccuracy are identified based on the observed relationships. Discussion: Computational modeling provides an essential component to ATD injury metric evaluation used to ensure the safety of future spaceflight occupants. The assessment of current ATD models lays the groundwork for how these models can be used appropriately in the future. Identification of limitations and possible paths for improvement aid in the development of these effective analysis tools.

Putnam, J.↗

The Role of Three-dimensional Effects on Ion Injection and Acceleration in Perpendicular Shocks

Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.

Orusa, Luca [Princeton Univ., NJ (United States); ↗

Transforming Energy Through Computational Excellence: Bringing Low Mach Number Reactive Flow Simulations at the Exascale

PeleLMeX's unique capabilities are allowing for reactive flow modeling at unprecedented scales and a reasonable time and cost. The code is currently being extended to tackle more practical, design-oriented simulations by implementing Large Eddy Simulation and data-driven chemical models, providing a fast but accurate tool for engineers considering the emergence of GPU-accelerated platforms. These extensions are critical for enabling the physical insight required to design the next generation of combustion devices as a key component of a renewable energy future.

MATHEMATICS AND COMPUTING↗

Generation and propagation of an electromagnetic pulse in the Trigger experiment and its possible role in electron acceleration

Instruments onboard the Trigger payload detected a large-amplitude, low-frequency, electric field pulse which was observed with a time delay consistent only with an electromagnetic wave. A model for this perturbation is constructed, and the associated field-aligned current is calculated as a function of altitude. This experiment may simulate the acceleration mechanism which results in the formation of auroral arcs, and possibly even other events in cosmic plasmas.

Kelley, M. C.↗

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Consequences of Alfven Waves and Parallel Potential Drops in the Auroral Zone

The goal of this research is to examine the causes of field-aligned plasma acceleration in the auroral zone using satellite data and numerical simulations. A primary question to be addressed is what causes the field-aligned acceleration of electrons (leading to precipitation) and ions (leading to upwelling ions) in the auroral zone. Data from the Fast Auroral SnapshoT (FAST) and Polar satellites is used when the two satellites are in approximate magnetic conjunction and are in the auroral region. FAST is at relatively low altitudes and samples plasma in the midst of the auroral acceleration region while Polar is at much higher altitudes and can measure plasmas and waves propagating towards the Earth. Polar can determine the sources of energy streaming earthward from the magnetotail, either in the form of field-aligned currents, electromagnetic waves or kinetic particle energy, that ultimately leads to the acceleration of plasma in the auroral zone. After identifying and examining several events, numerical simulations are run that bridges the spatial region between the two satellites. The code is a one-dimensional, long system length particle in cell simulation that has been developed to model the auroral region. A main goal of this research project is to include Alfven waves in the simulation to examine how these waves can accelerate plasma in the auroral zone.

Schriver, David↗

Adaptive sequential estimation with applications

A sequential estimator is presented and demonstrated which successfully tracks the system state and model errors in the presence of significant and unpredictable system or environmental variations. This adaptive estimation concept is shown to lead to a new and attractive approach to parameter identification problems. Simulation results are presented in an orbit determination problem, where the estimator tracks the orbit and unmodeled accelerations due to errors in the geopotential model. Some results in a re-entry trajectory estimation problem are also summarized.

Jazwinski, A. H.↗