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 901 records · Page 50

Simulation of proton radiation belt formation during the March 24, 1991 SSC

The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Hudson, M. K.↗

A Finite Element Study of Stress Reduction Techniques in REBCO HTS Conductor on a Round Cable (CORC) Cable

ReBCO high-temperature superconducting (HTS) tape is critical for achieving the high magnetic fields needed in next-generation particle accelerators. Enhancing the mechanical performance of ReBCO tape increases its critical current by reducing internal stress, especially in the superconducting layer. A finite element study examined how copper layer properties affect stress in ReBCO conductor on a round core (CORC) cables. The cable was modeled as a doubly supported beam under uniform compressive stresses."cable was modeled as a doubly supported beam under uniform load to simulate bending. A staged modeling approach from a single tape to a six-layer stack enabled validation and efficient parameter studies. Increasing the yield strength and Young s modulus of the copper layers reduced peak stress in the ReBCO layer. These results support development of improved tape stacks for high-field accelerator magnets

Mueller, Scott [Northern Illinois U.]↗

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

This is a new project which started in May 1996. The main objective of the experimental/numerical study is to improve the understanding of the physics of two-way coupling between the dispersed phase and turbulence in a prototypical turbulent shear flow - homogeneous shear, laden with small liquid droplets (in gas) or gaseous bubbles (in liquid). The method of direct numerical simulation (DNS) is used to solve the full three-dimensional, time-dependent Navier-Stokes equations including the terms describing the two-way coupling between the dispersed phase and the carrier flow. The results include the temporal evolution of the three-dimensional energy and dissipation spectra and the rate of energy transfer across the energy spectrum to understand the fundamental physics of turbulence modulation, especially the effects of varying the magnitude of gravitational acceleration. The mean-square displacement and diffusivity of the droplets (or bubbles) of a given size and the preferential accumulation of droplets in low vorticity regions and bubbles in high vorticity regions will be examined in detail for different magnitudes of gravitational acceleration. These numerical results which will be compared with their corresponding measured data will provide a data base from which a subgrid-scale (SGS) model can be developed and validated for use in large-eddy simulation (LES) of particle-laden shear flows. Two parallel sets of experiments will be conducted: bubbles in an immiscible liquid and droplets in air. In both experiments homogeneous shear will be imposed on the turbulent carrier flow. The instantaneous velocities of the fluid and polydispersed-size particles (droplets or bubbles) will be measured simultaneously using a two-component Phase-Doppler Particle Analyzer (PDPA). Also, the velocity statistics and energy spectra for the carrier flow will be measured.

Elghobashi, Said↗

On the Life Expectancy Ot High-power CW Magnetrons for Superconducting Accelrators

Modern CW or pulse Superconducting RF (SRF) accelerators require efficient RF sources controllable in phase and power with a reduced cost. Therefore, utilization of the high-power CW magnetrons as RF sources in SRF accelerator projects was proposed in a number of works, e.g., [1, 2]. But typically, the CW magnetrons are designed as RF sources for industrial heating, and the lifetime of the tubes is not the first priority as it is required for high-energy accelerators. The high-power industrial CW magnetrons use the cathodes made of pure tungsten. The emission properties of the tungsten cathodes are not deteriorated much by electron and ion bombardments, but the latter causes sputtering of the cathode in the magnetron crossed fields. The sputtered cathode material covers the magnetron interior. This leads to sparks and discharges that limit the life of the magnetrons. We considered an analysis of magnetron failure modes vs. output power [3]. We developed a model of ionization of the residual gas in the magnetrons interaction space and simulated the spattering of the cathode in 100 kW CW magnetrons to estimate the life expectancy. Basing on results we proposed ways to increase the CW magnetrons longevity for SRF accelerators.

43 PARTICLE ACCELERATORS↗

On the Life Expectance of High-power CW Magnetrons for SRF Accelerators

Modern CW or pulse Superconducting RF (SRF) accelerators require efficient RF sources controllable in phase and power with a reduced cost. Therefore, utilization of the high-power CW magnetrons as RF sources in SRF accelerator projects was proposed in a number of works, e.g., [1, 2]. But typically, the CW magnetrons are designed as RF sources for industrial heating, and the lifetime of the tubes is not the first priority as it is required for high-energy accelerators. The high-power industrial CW magnetrons use the cathodes made of pure tungsten. The emission properties of the tungsten cathodes are not deteriorated much by electron and ion bombardments, but the latter causes sputtering of the cathode in the magnetron crossed fields. The sputtered cathode material covers the magnetron interior. This leads to sparks and discharges that limit the life of the magnetrons. We considered an analysis of magnetron failure modes vs. output power [3]. We developed a model of ionization of the residual gas in the magnetrons interaction space and simulated the spattering of the cathode in 100 kW CW magnetrons to estimate the life expectancy. Basing on results we proposed ways to increase the CW magnetrons longevity for SRF accelerators.

43 PARTICLE ACCELERATORS↗

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Updated Simulation of the Booster Neutrino Beam

For accelerator neutrino experiments, an accurate prediction of the incoming neutrino flux is crucial for reducing uncertainties for all physics measurements. In this exciting period for the Short-Baseline Neutrino program at Fermilab, with far detector (ICARUS) already operating and the near detector (SBND) nearing operation, an updated flux model for the Booster Neutrino Beam (BNB) is presented. The BNB delivers 8 GeV protons to a beryllium target, subsequently producing neutrinos from the decay of the secondary beam of hadrons. A updated Monte-Carlo simulation of beam production in GEANT4 has been created, allowing predictions to be made for detectors with different baselines, offsets and sizes. This new simulation contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, with a view to storing all resulting particles - including neutral mesons - from the proton-Beryllium scatter to allow the study of exotic BSM scenarios. Results are presented, with comparisons to the original flux simulated for the MiniBooNE experiment.

Paton, Josephine↗

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗

US Department of Energy, Office of Science High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility

Oak Ridge National Laboratory's (ORNL's) Leadership Computing Facility (OLCF) continues to surpass its operational target goals: supporting users; delivering fast, reliable computational ecosystems; creating innovative solutions for high performance computing (HPC) needs; and managing risks, safety, and security associated with operating some of the most powerful computers in the world. The results can be seen in the cutting-edge science conducted by users and the praise from the research community. Calendar year (CY) 2019 was a big year as OLCF staff ran five world-class resources (the leadershipclass computers Titan and Summit, the large analysis cluster called Eos, and the massive parallel filesystems called Atlas and Alpine)) and also began power and cooling upgrades for a 2021 exascale system called Frontier. While continuing exceptional operation of Titan, Eos, and Rhea, the OLCF released the Summit supercomputer for production on January 1, 2019. Summit debuted as the most capable and efficient system in its class and has been recognized as the most powerful system in the world for its performance on both the high performance linpack (HPL) and conjugate gradient (HPCG) benchmark applications since June 2018 according to TOP500. Summit represents the culmination of a multiyear effort between the OLCF, IBM, NVIDIA, and Mellanox to deliver a system that is unmatched for modeling, simulation, data analysis, and learning. To hit the ground running with science-ready applications on day one, application teams worked closely with the OLCF through the Center for Accelerated Application Readiness (CAAR) program for years in advance of the Summit deployment. CY 2019 was filled with outstanding results and accomplishments: a very high rating from users on overall satisfaction for the sixth year in a row; a tremendous amount of core-hours delivered to researchers from two leadership-class systems; and success in delivering on the allocation split of roughly 60%, 30%, and 10% of core-hours offered for the Innovative and Novel Computational Impact on Theory and Experiment (INCITE), Advanced Scientific Computing Research Leadership Computing Challenge (ALCC), and Director's Discretionary (DD) programs, respectively (see Operational Performance section). These accomplishments, coupled with the high utilization rates (overall and capability usage), represent the fulfillment of the promise of both leadership-class machines: efficient facilitation of leadership-class computational applications. Table ES.1 presents a summary of the 2019 OLCF metric targets and the associated results. More information can be found in the Operational Performance section for each OLCF resource. The scientific accomplishments of OLCF users are a strong indication of long-term operational success, with publications this year in such notable journals and publications as Nature, Nature Physics, Nature Plants, Physical Review X, Journal of the American Physical Society, Cell, Nano Letters, and Trends in Biotechnology. Crucial domain-specific discoveries facilitated by resources at the OLCF are described in the High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility (OAR) Strategic Results section. For example, researchers used Summit to pinpoint and understand the production of proteins from genetic information, including mutations and the functional expression of disease (Section 8.2).

97 MATHEMATICS AND COMPUTING↗

PIP-II Beam Current Monitor Fault Case & Beam Position Monitor Linearity Studies in CST Studio

The use of non-invasive sensors & systems to measure particle beam characteristics is a crucial part of modern accelerator control systems. As such, simulations predicting the behavior of these sensors are essential for guiding beamline control system design. This poster details the results returned by two beam sensor models created using CST Studio software: the signal linearity of an elliptical beam position monitor (BPM), and the fault cases of an AC current transformer (ACCT) beam current monitor (BCM).

Rouzky, Adam↗

Advance prototype silver ion water bactericide system

An advance prototype unit was designed and fabricated to treat anticipated fuel cell water. The unit is a single canister that contains a membrane-type prefilter and a silver bromide contacting bed. A seven day baseline simulated mission test was performed; the performance was satisfactory and the effluent water was within all specifications for potability. After random vibrations another seven day simulated mission test was performed, and results indicate that simulated launch vibrations have no effects on the design and performance of the advanced prototype. Bench tests and accelerated breadboard tests were conducted to define the characteristics of an upgraded model of the advance prototype unit which would have 30 days of operating capability. A preliminary design of a silver ion generator for the shuttle orbiter was also prepared.

Jasionowski, W. J.↗

A comparison of the results of dynamic wind-tunnel tests with theoretical predictions for an aeromechanical gust-alleviation system for light airplanes

Dynamic wind tunnel tests have been conducted on a 1/6-scale model of a general aviation airplane equipped with an all-mechanical gust alleviation system which uses auxiliary aerodynamic surfaces to drive the flaps. The longitudinal short period motions were studied under simulated gust conditions in order to verify the mathematical model used in a previous study to predict the performance of the full scale system and determine the amount of normal acceleration alleviation which could be attained. The model responses were measured for different configurations with the system active and without the system active for comparison. The tests confirmed the general relationships between the experimental variables and the model responses predicted by the mathematical model, but there were significant differences in the magnitudes of the responses. The experimental results for the model were used to estimate a reduction of 30 percent in the rms normal acceleration response of a similar full scale airplane in atmospheric turbulence.

Stewart, E. C.↗

A unique approach to space thermal simulation

The Space Systems Lab of McDonnell Douglas developed several innovative techniques to perform the DELTA Star System thermal vacuum test. Test design constraints included an accelerated test schedule and an ambitious test plan focusing on cyclic absorbed orbital flux simulation. Test design and fabrication was completed in only 4 months, producing a 21 zone computer modeled refectorless lamp array system. Array zone control was based on discreet absorbed flux measurements, provided by adiabatic coupons conceived to eliminate the need for conventional radiometers. Test requirements, methodology, and test thermal model flight data correlation are presented.

Durant, D. Q.↗

Parallelizing Timed Petri Net simulations

The possibility of using parallel processing to accelerate the simulation of Timed Petri Nets (TPN's) was studied. It was recognized that complex system development tools often transform system descriptions into TPN's or TPN-like models, which are then simulated to obtain information about system behavior. Viewed this way, it was important that the parallelization of TPN's be as automatic as possible, to admit the possibility of the parallelization being embedded in the system design tool. Later years of the grant were devoted to examining the problem of joint performance and reliability analysis, to explore whether both types of analysis could be accomplished within a single framework. In this final report, the results of our studies are summarized. We believe that the problem of parallelizing TPN's automatically for MIMD architectures has been almost completely solved for a large and important class of problems. Our initial investigations into joint performance/reliability analysis are two-fold; it was shown that Monte Carlo simulation, with importance sampling, offers promise of joint analysis in the context of a single tool, and methods for the parallel simulation of general Continuous Time Markov Chains, a model framework within which joint performance/reliability models can be cast, were developed. However, very much more work is needed to determine the scope and generality of these approaches. The results obtained in our two studies, future directions for this type of work, and a list of publications are included.

Nicol, David M.↗

The Influence of Convection on Magnetotail Variability

This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model. In this model the plasma mantle continuously supplies particles to the magnetotail, the ion current periodically updates the magnetic field using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the model's near-Earth X-line. This variability results from the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. The characteristic time scale of variability in the magnetotail is on the order of 4 - 5 minutes. We also investigate how the magnetotail's topology responds to increased convection electric fields, and show examples of observations of variability in the magnetotail.

Peroomian, Vahe↗

Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms

This paper presents a final report on Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms. The focus of this research is to derive and implement: 1) Wavelet based methodologies for the compression, transmission, decoding, and visualization of three dimensional finite element geometry and simulation data in a network environment; 2) methodologies for interactive algorithm monitoring and tracking in computational mechanics; and 3) Methodologies for interactive algorithm steering for the acceleration of large scale finite element simulations. Also included in this report are appendices describing the derivation of wavelet based Particle Image Velocity algorithms and reduced order input-output models for nonlinear systems by utilizing wavelet approximations.

Kurdila, Andrew J.↗