Search NASA⌕ Search

SEARCH · Search NASA

Results for “High gradient acceleration”

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 685 records · Page 38

A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation

In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad–Shafranov equation, the solution of which determines the poloidal component of the magnetic field. When the geometry of the confinement region is known, the problem becomes an interior Dirichlet boundary value problem. We propose a high order solver based on the Hybridizable Discontinuous Galerkin method. The resulting algorithm (1) provides high order of convergence for the flux function and its gradient, (2) incorporates a novel method for handling piecewise smooth geometries by extension from polygonal meshes, (3) can handle geometries with non-smooth boundaries and x-points, and (4) deals with the semi-linearity through an accelerated two-grid fixed-point iteration. The effectiveness of the algorithm is verified with computations for cases where analytic solutions are known on configurations similar to those of actual devices (ITER with single null and double null divertor, NSTX, ASDEX upgrade, and Field Reversed Configurations).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HD-Bind: Encoding of Molecular Structure with Low Precision, Hyperdimensional Binary Representations

Publicly available collections of drug-like molecules have grown to comprise tens of billions of compounds due to advances in combinatorial chemistry. Traditional methods for identifying "hit" molecules from a large collection of potential drug-like candidates have relied on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have the major drawback that they require exceptional computing capabilities for even relatively small collections of molecules. Hyperdimensional Computing (HDC) is a recently-proposed learning paradigm that represents data with high-dimension binary vectors; this allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas. We consider existing HDC approaches for molecular property classification and introduce two novel encodings of a commonly-used molecular representation, the extended connectivity fingerprint (ECFP). We show that HDC-based inference methods are as much as 91 times more efficient than traditional machine learning methods, and achieve an acceleration of nearly nine orders of magnitude compared to molecular docking. Our results show that HDC accelerated methods retain competitive accuracy on a number of well-studied tasks such as molecular property predictions using the MoleculeNet dataset, and bind/no-bind activity classification using the DUD-E and LIT-PCBA datasets. Our work thus motivates further investigation into molecular representation learning to develop ultraefficient pre-screening tools.

Jones, William↗

A Parallel Particle Swarm Optimization Algorithm Accelerated by Asynchronous Evaluations

A parallel Particle Swarm Optimization (PSO) algorithm is presented. Particle swarm optimization is a fairly recent addition to the family of non-gradient based, probabilistic search algorithms that is based on a simplified social model and is closely tied to swarming theory. Although PSO algorithms present several attractive properties to the designer, they are plagued by high computational cost as measured by elapsed time. One approach to reduce the elapsed time is to make use of coarse-grained parallelization to evaluate the design points. Previous parallel PSO algorithms were mostly implemented in a synchronous manner, where all design points within a design iteration are evaluated before the next iteration is started. This approach leads to poor parallel speedup in cases where a heterogeneous parallel environment is used and/or where the analysis time depends on the design point being analyzed. This paper introduces an asynchronous parallel PSO algorithm that greatly improves the parallel e ciency. The asynchronous algorithm is benchmarked on a cluster assembled of Apple Macintosh G5 desktop computers, using the multi-disciplinary optimization of a typical transport aircraft wing as an example.

Venter, Gerhard↗

The RHIC and RHIC pre-injectors controls systems: status and plans

For the past twelve years experiments at the Relativistic Heavy Ion Collider (RHIC) have recorded data from collisions of heavy ions and polarized protons, leading to important discoveries in nuclear physics and the spin dynamics of quarks and gluons. BNL is the site of one of the first and still operating alternating gradient synchrotrons, the AGS, which first operated in 1960. The accelerator controls systems for these instruments span multiple generations of technologies. In this report we will describe the current status of the Collider-Accelerator Department controls systems, which are used to control seven different accelerator facilities and multiple science programs (high energy nuclear physics, high energy polarized proton physics, NASA programs, isotope production, and multiple accelerator research and development projects). We will describe the status of current projects, such as the just completed Electron Beam Ion Source (EBIS), our R&D programs in superconducting RF and an Energy Recovery LINAC (ERL), innovations in feedback systems and bunched beam stochastic cooling at RHIC, and plans for future controls system developments.

43 PARTICLE ACCELERATORS↗

Operation Regime Analysis of Conduction Cooled Cavities through Multi-Physics Simulation

Euclid Techlabs in collaboration with Fermilab IARC (Batavia, IL) is developing industrial superconducting 10MeV electron linac [1, 2]. Cryocoolers are to be used for cooling instead of liquid helium bath to simplify the linac infrastructure [3]. The cavity linked to commercially available cryo-cooler cold head [4, 5] through highly conductive aluminium (AL) strips. However, this solution raises a problem of contact thermal resistance. This paper shows some results of Comsol multyphysics simulations of the cavity cooling by AL strips. Some insight was obtained on the acceptable range of contact resistance. Operation regimes were obtained at different accelerating gradients and cavity temperatures. The results of simula-tion are presented and discussed.

43 PARTICLE ACCELERATORS↗

SU-E-J-237: Image Feature Based DRR and Portal Image Registration

Purpose: Two-dimensional (2D) matching of the kV X-ray and digitally reconstructed radiography (DRR) images is an important setup technique for image-guided radiotherapy (IGRT). In our clinics, mutual information based methods are used for this purpose on commercial linear accelerators, but with often needs for manual corrections. This work proved the feasibility that feature based image transform can be used to register kV and DRR images. Methods: The scale invariant feature transform (SIFT) method was implemented to detect the matching image details (or key points) between the kV and DRR images. These key points represent high image intensity gradients, and thus the scale invariant features. Due to the poor image contrast from our kV image, direct application of the SIFT method yielded many detection errors. To assist the finding of key points, the center coordinates of the kV and DRR images were read from the DICOM header, and the two groups of key points with similar relative positions to their corresponding centers were paired up. Using these points, a rigid transform (with scaling, horizontal and vertical shifts) was estimated. We also artificially introduced vertical and horizontal shifts to test the accuracy of our registration method on anterior-posterior (AP) and lateral pelvic images. Results: The results provided a satisfactory overlay of the transformed kV onto the DRR image. The introduced vs. detected shifts were fit into a linear regression. In the AP image experiments, linear regression analysis showed a slope of 1.15 and 0.98 with an R2 of 0.89 and 0.99 for the horizontal and vertical shifts, respectively. The results are 1.2 and 1.3 with R2 of 0.72 and 0.82 for the lateral image shifts. Conclusion: This work provided an alternative technique for kV to DRR alignment. Further improvements in the estimation accuracy and image contrast tolerance are underway.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Ion probe beam experiments and kinetic modeling in a dense plasma focus Z-pinch

The Z-pinch phase of a dense plasma focus (DPF) emits multiple-MeV ions in a ∼cm length. The mechanisms through which these physically simple devices generate such high energy beams in a relatively short distance are not fully understood. We are exploring the origins of these large gradients using measurements of an ion probe beam injected into a DPF during the pinch phase and the first kinetic simulations of a DPF Z-pinch. To probe the accelerating fields in our table top experiment, we inject a 4 MeV deuteron beam along the z-axis and then sample the beam energy distribution after it passes through the pinch region. Using this technique, we have directly measured for the first time the acceleration of an injected ion beam. Our particle-in-cell simulations have been benchmarked on both a kJ-scale DPF and a MJ-scale DPF. They have reproduced experimentally measured neutron yields as well as ion beams and EM oscillations which fluid simulations do not exhibit. Direct comparisons between the experiment and simulations enhance our understanding of these plasmas and provide predictive design capability for accelerator and neutron source applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of gravity on sheared and nonsheared turbulent nonpremixed flames

The present numerical study is concerned with the fundamental physics of the multiway interaction between turbulence, chemical reaction, and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous, three-dimensional governing equations. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration. The effects of buoyancy include the generation of local shear, baroclinic production or destruction of vorticity, and countergradient heat and mass transport. Increased vorticity and small-scale turbulence promote further mixing and reaction. However, if the strain-rates become too high, local flame extinction can occur. Our objective is to gain an understanding of the complex interactions between the physical phenomena involved, with particular attention to the effects of buoyancy on the turbulence structure, flame behavior, and factors influencing flame extinction.

Elghobashi, Said↗

An Estimation of the Flying Qualities of the Kaiser Fleetwing All-Wing Airplane from Tests of a 1/7-Scale Model, TED No. NACA 2340

An investigation of a 1/7-scale powered model of the Kaiser Fleetwing all-wing airplane was made in the Langley full-scale tunnel to provide data for an estimation of the flying qualities of the airplane. The analysis of the stability and control characteristics of the airplane has been made as closely as possible in accordance with the requirements of the Bureau of Aeronautics, Navy Department's specifications, and a summary of the more significant conclusions is presented as follows. With the normal center of gravity located at 20 percent of the mean aerodynamic chord, the airplane will have adequate static longitudinal stability, elevator fixed, for all flight conditions except for low-power operation at low speeds where the stability will be about neutral. There will not be sufficient down-elevator deflection available for trim above speeds of about 130 miles per hour. It is probable that the reduction in the up-elevator deflections required for trim will be accompanied by reduced elevator hinge moments for low-power operation at low flight speeds. The static directional stability for this airplane will be low for all rudder-fixed or rudder-free flight conditions. The maximum rudder deflection of 30 deg will trim only about 15 deg yaw for most flight conditions and only 10 deg yaw for the condition with low power at low speeds. Also, at low powers and low speeds, it is estimated that the rudders will not trim the total adverse yaw resulting from an abrupt aileron roll using maximum aileron deflection. The airplane will meet the requirements for stability and control for asymmetric power operation with one outboard engine inoperative. The airplane would have no tendency for directional divergence but would probably be spirally unstable, with rudders fixed. The static lateral stability of the airplane will probably be about neutral for the high-speed flight conditions and will be only slightly increased for the low-power operation in low-speed flight. The airplane will not roll against the ailerons in a side-slip maneuver. Although the airplane would probably meet the minimum requirements of pb/2V of 0.07 at all speeds, there will be a loss in rolling ability of the airplane at high aileron deflections and at low flight speeds. It is estimated that the wing stall will be a gradual movement forward from the trailing edge and will be accompanied by no sudden pitching or rolling accelerations. Some stall warning may be indicated by reduction in the elevator and aileron force gradients and by the shaking of the controls caused by unsteady flow over the surfaces near the stall.

Brewer, Gerald W.↗

Status of High Power Tests of Normal Conducting Single-Cell Structures

We report the results of ongoing high power tests of single-cell standing wave structures. These tests are part of an experimental and theoretical study of rf breakdown in normal conducting structures at 11.4 GHz. The goal of this study is to determine the maximum gradient possibilities for normal-conducting rf powered particle beam accelerators. The test setup consists of reusable mode launchers and short test structures powered by SLACs XL-4 klystron. The mode launchers and structures were manufactured at SLAC and KEK and tested at the SLAC klystron test laboratory.

43 PARTICLE ACCELERATORS↗

4π Non-Coplanar Liver SBRT: A Novel Delivery Technique

Purpose: To improve the quality of liver stereotactic body radiation therapy (SBRT) treatments, a novel 4π framework was developed with accompanying algorithms to optimize non-coplanar beam orientations and fluences. The dose optimization is performed on a patient-specific deliverable beam geometry solution space, parameterized with patient and linear accelerator gantry orientations. Methods and Materials: Beams causing collision between the gantry and the couch or patient were eliminated by simulating all beam orientations using a precise computer assisted design model of the linear accelerator and a human subject. Integrated beam orientation and fluence map optimizations were performed on remaining beams using a greedy column generation method. Testing of the new method was performed on 10 liver SBRT cases previously treated with 50 to 60 Gy in 5 fractions using volumetric modulated arc therapy (VMAT). For each patient, both 14 and 22 non-coplanar fields were selected and optimized to meet the objective of ≥95% of the planning target volume (PTV) covered by 100% of the prescription dose. Doses to organs at risk, normal liver volumes receiving <15 Gy, integral dose, and 50% dose spillage volumes were compared against the delivered clinical VMAT plans. Results: Compared with the VMAT plans, the 4π plans yielded reduced 50% dose spillage volume and integral dose by 22% (range 10%-40%) and 19% (range 13%-26%), respectively. The mean normal liver volume receiving <15 Gy was increased by 51 cc (range 21-107 cc) with a 31% reduction of the mean normal liver dose. Mean doses to the left kidney and right kidney and maximum doses to the stomach and spinal cord were on average reduced by 70%, 51%, 67%, and 64% (P≤.05). Conclusions: This novel 4π non-coplanar radiation delivery technique significantly improved dose gradient, reduced high dose spillage, and improved organ at risk sparing compared with state of the art VMAT plans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Progress in the Development of the Nb 3 Sn MQXFB Quadrupole for the HiLumi Upgrade of the LHC

The high-luminosity upgrade of the Large Hadron Collider (HL-LHC) requires new high field and large-aperture quadrupole magnets for the low-beta inner triplets (MQXF). With a nominal operating gradient of 132.2 T/m in a 150 mm aperture and a conductor peak field of 11.3 T, the new quadrupole magnets are based on Nb 3 Sn superconducting technology. After a series of short models constructed in close collaboration by LARP (LHC Accelerator Research Program) and CERN, the development program is entering in the series production phase with CERN on one side and the US Accelerator Upgrade Project (US-AUP) on the other side assembling and testing full-length magnets. This paper describes the status of the development activities at CERN, in particular on the cold powering test of the first MQFXB prototype and on the construction of the second full scale prototype. Critical operations such as reaction heat treatment, coil impregnation and magnet assembly are discussed. Finally, the plan towards the series production is described.

43 PARTICLE ACCELERATORS↗

A modeling study of ocean thermal energy conversion resource and potential environmental effects around Kailua-Kona, Hawaii

Ocean Thermal Energy Conversion (OTEC) offers a promising renewable energy solution through a heat exchange process using the temperature difference between warm surface seawater and cold deep seawater. Because accurate resource characterization is critical for the optimal design and implementation of OTEC systems, a high-resolution numerical model is employed to better characterize the OTEC resource at Kona, Hawaii. Our model provides detailed spatial and temporal variability of the thermal gradient, which is essential for assessing the viability and efficiency of OTEC systems. The model results reveal distinct patterns and dynamics not captured by existing observations or models (e.g., lower-resolution information). These findings highlight the importance of using high-resolution models for accurate predictions of thermal gradient variability, ultimately supporting more efficient and sustainable OTEC deployment. Additionally, the study investigates the impacts of mixed water discharge from OTEC plants that can cause shock to organisms living in the surface water and potentially destabilize the water column. Understanding these effects is vital for minimizing any potential negative environmental consequences and ensuring the long-term viability of OTEC operations. Further, our model improves OTEC resource characterization, which can lead to optimal design and deployment of OTEC systems. The analysis of OTEC water discharge impacts can accelerate the development of OTEC technologies, overcoming permitting/consenting challenges. These findings contribute to the broader adoption of high-resolution modeling in ocean energy resource characterization, particularly for OTEC applications.

30 DIRECT ENERGY CONVERSION↗

Flux Compression Magnetic Nozzle

In pulsed fusion propulsion schemes in which the fusion energy creates a radially expanding plasma, a magnetic nozzle is required to redirect the radially diverging flow of the expanding fusion plasma into a rearward axial flow, thereby producing a forward axial impulse to the vehicle. In a highly electrically conducting plasma, the presence of a magnetic field B in the plasma creates a pressure B(exp 2)/2(mu) in the plasma, the magnetic pressure. A gradient in the magnetic pressure can be used to decelerate the plasma traveling in the direction of increasing magnetic field, or to accelerate a plasma from rest in the direction of decreasing magnetic pressure. In principle, ignoring dissipative processes, it is possible to design magnetic configurations to produce an 'elastic' deflection of a plasma beam. In particular, it is conceivable that, by an appropriate arrangement of a set of coils, a good approximation to a parabolic 'magnetic mirror' may be formed, such that a beam of charged particles emanating from the focal point of the parabolic mirror would be reflected by the mirror to travel axially away from the mirror. The degree to which this may be accomplished depends on the degree of control one has over the flux surface of the magnetic field, which changes as a result of its interaction with a moving plasma.

Thio, Y. C. Francis↗

Multi-variance replica exchange SGMCMC for inverse and forward problems via Bayesian PINN

Physics-informed neural network (PINN) has been successfully applied in solving a variety of nonlinear non-convex forward and inverse problems. However, the training is challenging because of the non-convex loss functions and the multiple optima in the Bayesian inverse problem. In this work, we propose a multi-variance replica exchange stochastic gradient Langevin dynamics method to tackle the challenge of the multiple local optima in the optimization and the challenge of the multiple modal posterior distribution in the inverse problem. Replica exchange methods are capable of escaping from the local traps and accelerating the convergence; two chains with different temperatures are designed where the low temperature chain aims for the local convergence, and the target of the high temperature chain is to travel globally and explore the whole loss function entropy landscape. However, it may not be efficient to solve mathematical inversion problems by using the vanilla replica method directly since the method doubles the computational cost in evaluating the forward solvers (likelihood functions) in the two chains. To address this issue, we propose to make different assumptions on the energy function estimation and this facilities one to use solvers of different fidelities in the likelihood function evaluation. More precisely, one can use a solver with low fidelity in the high temperature chain while using a solver with high fidelity in the low temperature chain. Our proposed method significantly lowers the computational cost in the high temperature chain, meanwhile preserving the accuracy and converging very fast. Here we give an unbiased estimate of the swapping rate and give an estimation of the discretization error of the scheme. To verify our idea, we design and solve four inverse problems which have multiple modes. The proposed method is also employed to train the Bayesian PINN to solve the forward and inverse problems; faster and more accurate convergence has been observed when compared to the stochastic gradient Langevin dynamics (SGLD) method and vanilla replica exchange methods.

97 MATHEMATICS AND COMPUTING↗

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗