Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modeling”

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 703 records · Page 39

Extended State Observer-Based Robust Model Predictive Velocity Control for Permanent Magnet Synchronous Motor

This article proposes an extended state observer based robust model predictive velocity control to decrease system prediction error under parameter uncertainties for permanent magnet synchronous motor (PMSM). We develop a new PMSM model that consists of velocity and acceleration to lump the system information and an external disturbance into a disturbance. The extended state observer (ESO) is designed to estimate the velocity, acceleration, and disturbance. By estimating the state variables and disturbance using the ESO, the model predictive control (MPC) finds the optimal control input by predicting future system behavior. Additionally, the direct current controller is designed so that the direct current converges to zero. Because the proposed method is not designed based on the cascade structure from the viewpoint of velocity control, the optimization control for the velocity and currents can be defined. Thus, the proposed method is robust against external disturbances and parameter uncertainties owing to feedback linearization, state feedback, and ESO-based MPC using the acceleration PMSM model. The proposed control algorithm was experimentally verified and it showed improved velocity tracking performance compared with ESO-based MPC using the conventional PMSM model.

42 ENGINEERING↗

A method to determine the kinematics of the lower limbs of a subject pedaling a bicycle using encoders and accelerometers

The goal of this research was to determine kinematic parameters of the lower limbs of a subject pedaling a bicycle. An existing measurement system was used as the basis to develop the model to determine position and acceleration of the limbs. The system consists of an ergometer instrumented to provide position of the pedal (foot), accelerometers to be attached to the lower limbs to measure accelerations, a recorder used for filtering, and a computer instrumented with an A/D board and a decoder board. The system is designed to read and record data from accelerometers and encoders. Software has been developed for data collection, analysis and presentation. Based on the measurement system, a two dimensional analytical model has been developed to determine configuration (position, orientation) and kinematics (velocities, accelerations). The model has been implemented in software and verified by simulation. An error analysis to determine the system's accuracy shows that the expected error is well within the specifications of practical applications. When the physical hardware is completed, NASA researchers hope to use the system developed to determine forces exerted by muscles and forces at articulations. This data will be useful in the development of countermeasures to minimize bone loss experienced by astronauts in microgravity conditions.

Liu, Shih-Ching↗

Lunar Scout Two spacecraft gravity experiment

Measurement of the gravity field of the Moon has a high science priority because of its implications for the internal structure and thermal history of the Moon, and it has a high priority for future exploration activities because of the influence of lunar gravity on spacecraft navigation and orbit maintenance. The current state of knowledge in the lunar gravity field (and the uncertainty in the knowledge) is based primarily on data accumulated from the Lunar Orbiter and Apollo programs. Data are sparse and emphasize the equatorial band (+/- 30 deg) on the near side of the Moon. There are no tracking data on the far side and only the Lunar Orbiter 5 provides a small amount of high inclination data. A host of gravity models developed from different combinations of tracking data have large discrepancies in their predictions of spacecraft motion and orbit lifetimes. There are also large disagreements in the Mercator projections of the gravity acceleration from each model, especially on the far side, where the contours tend to have no obvious relationship with the local topography. The science and engineering requirements for global gravity field mapping will be satisfied with continuous radio metric tracking of Lunar Scout 1 in a low polar orbit using the Deep Space Network and Lunar Scout 2 in a high elliptical orbit. The gravity field of the Moon will be mapped during the Scout Program using a two spacecraft concept. In the two spacecraft concept, one spacecraft is placed in a high altitude eccentric orbit while the second spacecraft is in a low altitude polar orbit. The gravity experiment requires a radio frequency that will permit two-way Doppler tracking between the spacecraft and the Deep Space Network (DSN). Both spacecraft carry NASA standard transponder systems for data transmission to Earth as well as for tracking and orbit determination. Data sufficient to produce a gravity field map could be acquired within one month with this system.

Cheng, Andrew F.↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Here, we use the public code EBHLIGHT to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

79 ASTRONOMY AND ASTROPHYSICS↗

Solar Energetic Particle Acceleration at a Spherical Shock with the Shock Normal Angle $θ$ $B_n$ Evolving in Space and Time

We present a 2D kinematic model to study the acceleration of solar energetic particles (SEPs) at a shock driven by a coronal mass ejection. The shock is assumed to be spherical about an origin that is offset from the center of the Sun. This leads to a spatial and temporal evolution of the angle between the magnetic field and the shock-normal direction (θ Bn ) as it propagates through the Parker spiral magnetic field from the lower corona to 1 au. We find that the high-energy SEP intensity varies significantly along the shock front due to the evolution of θ Bn . Generally, the west flank of the shock preferentially accelerates particles to high energies compared to the east flank and shock nose. This can be understood in terms of the rate of acceleration, which is higher at the west flank. Double power-law energy spectra are reproduced in our model as a consequence of the local acceleration and transport effects. These results will help us to better understand the evolution of SEP acceleration and provide new insights into large SEP events observed by multiple spacecraft, especially those close to the Sun, such as Parker Solar Probe and Solar Orbiter.

79 ASTRONOMY AND ASTROPHYSICS↗

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗

Modeling and Impact of Solar Energetic Particles in the Heliosphere and Geospace

Understanding the radiation environment due to solar energetic particles in the heliosphere and the Earth’s magnetosphere is a challenging and practically important task. Exposure to energetic particles often leads to malfunctions and unexpected failures of electronics onboard spacecraft. The most vulnerable are exploratory missions when outside of the Earth’s magnetosphere. Geomagnetic field deflects Solar Energetic Particles (SEPs) moving through geospace though some of these particles propagate to LEO and have a high penetrating capability, thus producing significant radiation hazards for human spaceflight. Solar energetic particles also have an essential effect on the composition and dynamics of the Earth’s atmosphere. Precipitating SEPs enhance the atmospheric concentration of NOx and HOx, which play a crucial role in the ozone balance in the middle atmosphere by destroying odd oxygen through catalytic reactions. Numerical modeling of the radiation environment due to SEPs in the inner heliosphere and geospace is a multifold problem. That includes simulating 1) solar wind dynamics and the interplanetary magnetic field, 2) global modeling of the Earth’s magnetosphere, and 3) modeling transport and acceleration of SEPs in the inner heliosphere and geospace. The lecture will outline the key physical processes that control the behavior of SEPs in the heliosphere and geospace, along with the contemporary numerical methods used for their modeling. It will primarily focus on describing the SEP population in geospace across different altitudes, ranging from Low Earth Orbit (LEO) through Medium Earth Orbit (MEO) and Geostationary Orbit (GEO), up to the magnetopause, while considering the realistic geomagnetic field. Additionally, the lecture will cover how the SEP population in geospace varies in response to geomagnetic activity. This includes the temporal trapping of SEPs in geospace and the reduction of the rigidity cutoff during geomagnetic storms.

solar energetic particles↗

The origin of the nonthermal radiation of quasi-stellar objects and galactic nuclei

The supermassive-star model for the origin of the nonthermal luminosity of quasars and nuclei of Seyfert galaxies is analyzed. In the general model, a moving magnetic field is the accelerating mechanism for ultrarelativistic particles and gravitational collapse is the energy source. Several specific models are considered for the magnetic field, including a frozen-in field that is twisted about the axis of a rotating supermassive star, a tightly wound magnetic spiral, and the rotating field of a giant pulsar. In the general model, an accelerating efficiency close to unity is assumed. However, it is found that observational and theoretical evidence indicate an accelerating efficiency of less than 0.01 for the specific magnetic-field models. Mass-loss rates of quasars and Seyfert nuclei are calculated on the basis of the lower efficiency value and found to be reasonable.

Opher, R.↗

MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8) for advanced open-source tritium transport and fuel cycle modeling

Tritium management is critical for the safety, sustainability, and economics of fusion energy systems, and advanced and reliable modeling tools help accelerate the development of tritium technologies. This paper presents the Tritium Migration Analysis Program, Version 8 (TMAP8), an open-source, MOOSE-based application developed to provide state-of-the-art tritium transport and fuel cycle modeling capabilities. TMAP8 aims to expand the capabilities of previous versions (i.e., TMAP4 and TMAP7) by leveraging modern computational techniques, ensuring high software quality assurance standards (key to building trust), and enabling multispecies, multiscale, and multiphysics simulations for integrated tritium transport modeling in complex geometries. This paper outlines TMAP8’s scope and rigorous development practices, emphasizing its transparency, accessibility, modularity, and reliability. We present the current suite of verification and validation cases based on those from TMAP4, demonstrating TMAP8’s accuracy and reliability against analytical solutions and experimental data. Additionally, the paper showcases TMAP8’s integrated fuel cycle modeling capabilities, highlighting its applicability at various scales and levels. The TMAP8 code and documentation are openly available, promoting collaborative development and widespread adoption within the fusion community. Future work will soon expand TMAP8’s verification and validation suite to include those from TMAP7 and other recent experimental studies for validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Spectral Single and Double Power-law Formation by Sequential Particle Acceleration in Flare Magnetic Islands

Spectral single and double power laws are common in high-energy phenomena, such as solar flares and solar energetic particles, including ground level enhancement events. It is not clear what determines the energy breaks and spectral indexes of these power laws. Here, we will describe a first-principles model of pitch-angle and energy distribution function evolution, which produces power laws and provides physical interpretation for such spectral features (Guidoni et al. 2022, ApJ). In this model, a prescribed fraction of particles sequentially “hops” between shrinking magnetic islands (accelerators) formed by flare reconnection. Each accelerator increases particles’ energies by a modest amount, but particles must visit only a few accelerators to increase their energies by orders of magnitude. Data from global magnetohydrodynamic simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. We will also describe the fully analytic method for forming and interpreting power laws, which requires only a few constrained physical parameters of the acceleration region and is independent of the acceleration model, as well as preliminary results extending the analytical model to the formation of double power laws.

S E Guidoni↗

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coherent ion acceleration at a density gradient

The paper describes a model, in which ions are accelerated against the direction of the electric field at a density gradient. A one-dimensional, inhomogeneous, steady-state is constructed, in which cold electrons stream relative to cold ions. The wavenumber and amplitude of the fluctuating potential are calculated, then used to compute the trajectories of representative ions trapped in the potential troughs. Ion acceleration to energies an order of magnitude greater than the applied potential is demonstrated.

Hsieh, S. L.↗

Beam breakup in a solid state powered linear induction accelerator

It has been proposed to use solid-state pulsed-power modules to power the new Scorpius accelerator. Since this highcurrent accelerator would have 360 gaps (more than five times the number of cells as present radiography LIAs), it is natural to be concerned about the beam-breakup (BBU) instability. Earlier calculations of this effect used a half-scale model of the accelerator to avoid the limitation to 200 gaps in the code. We have now updated the code to accommodate the extraordinary number of gaps in this proposed LIA. The results of these early simulations are presented here. They indicate that BBU can only be suppressed to the same level as in designs powered by conventional pulsed power, if the calculated cell impedance can truly be achieved in practice. The magnetic focusing field required for this is of the same order as the field required for a conventionally powered LIA only if the transverse impedance can indeed be constrained to the values obtained in some simulations of the cell.

42 ENGINEERING↗

Electron acceleration by magnetosonic waves in solar flares

Results of a simulation of electron acceleration by magnetosonic waves in a solar coronal loop are presented. The dependence of the energization rate on the strength and shape of the wave spectrum is discussed. Comparisons are made with the predictions of quasilinear theory and of a model describing stochastic acceleration by waves as isotropic diffusion in momentum space.

Zweibel, Ellen G.↗

Stellar winds driven by Alfven waves

Models of stellar winds were considered in which the dynamic expansion of a corona is driven by Alfven waves propagating outward along radial magnetic field lines. In the presence of Alfven waves, a coronal expansion can exist for a broad range of reference conditions which would, in the absence of waves, lead to static configurations. Wind models in which the acceleration mechanism is due to Alfven waves alone and exhibit lower mass fluxes and higher energies per particle are compared to wind models in which the acceleration is due to thermal processes. For example, winds driven by Alfven waves exhibit streaming velocities at infinity which may vary between the escape velocity at the coronal base and the geometrical mean of the escape velocity and the speed of light. Upper and lower limits were derived for the allowed energy fluxes and mass fluxes associated with these winds.

Belcher, J. W.↗

Mesoscale density variability in the mesosphere and thermosphere: Effects of vertical flow accelerations

A mechanistic one dimensional numerical (iteration) model was developed which can be used to simulate specific types of mesoscale atmospheric density (and pressure) variability in the mesosphere and the thermosphere, namely those due to waves and those due to vertical flow accelerations. The model was developed with the idea that it could be used as a supplement to the TGCMs (thermospheric general circulation models) since such models have a very limited ability to model phenomena on small spatial scales. The simplest case to consider was the integration upward through a time averaged, height independent, horizontally divergent flow field. Vertical winds were initialized at the lower boundary using the Ekman pumping theory over flat terrain. The results of the computations are summarized.

Revelle, D. O.↗

In-flight simulation with pilot-center of gravity offset and velocity mismatch

Similarity transformations which preserve modal characteristics and pilot's acceleration cues in in-flight simulation are presented. The model transformation for lateral acceleration matching is developed. A velocity-mismatch example, based on a VRA simulation of the Space Shuttle, illustrates that acceleration matching is achieved at the expense of mismatching in cues which are secondary to the simulated piloting task, while primarily cues are preserved. The approach is applicable for both implicit and explicit model-following, and it can easily be extended to the longitudinal case.

Stengel, R. F.↗