Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear accelerators”

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 1,531 records · Page 85

Implicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids

An implicit code for computing inviscid and viscous incompressible flows on unstructured grids is described. The foundation of the code is a backward Euler time discretization for which the linear system is approximately solved at each time step with either a point implicit method or a preconditioned Generalized Minimal Residual (GMRES) technique. For the GMRES calculations, several techniques are investigated for forming the matrix-vector product. Convergence acceleration is achieved through a multigrid scheme that uses non-nested coarse grids that are generated using a technique described in the present paper. Convergence characteristics are investigated and results are compared with an exact solution for the inviscid flow over a four-element airfoil. Viscous results, which are compared with experimental data, include the turbulent flow over a NACA 4412 airfoil, a three-element airfoil for which Mach number effects are investigated, and three-dimensional flow over a wing with a partial-span flap.

Anderson, W. Kyle↗

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗

Accelerated convergence of structured banded systems using constrained corrections

An efficient iterative method for solving a structured banded system of equations is described. The method was developed for a full potential flow program and uses a basic interation step, a dynamic relation step, and a multigrid concept of constraining iterative corrections. The solution of a large linear system of equations is examined. Efficient iterative methods have become attractive for large problems. In the nonlinear cases, these iterations may be effectively merged to improve convergence rates.

Kneile, K.↗

Accelerating GNNs on GPU Sparse Tensor Cores through N:M Sparsity-Oriented Graph Reordering

Recent advancements in GPU hardware support have introduced the capability to leverage N:M sparse patterns for substantial performance gains. Graphs in Graph Neural Networks (GNNs) are typically sparse, but the sparsity is often irregular, not conforming to such sparse patterns. In this paper, we propose a novel graph reordering algorithm, the first of its kind, to reshape irregular graph data into the N:M structured sparse pattern at the tile level, allowing linear-algebra-based graph operations in GNNs to benefit from the N:M sparse hardware. The optimization is lossless, maintaining the accuracy of GNN. It can remove 98-100\% violations of the N:M sparse patterns at the vector level, and increase the proportion of conforming graphs in SuiteSparse collection from 5-9\% to 88.7-93.5\%. On A100 GPUs, the optimization accelerates Sparse Matrix Matrix (SpMM) by up to 43X (2.3X -- 7.5X on average) and speeds up the key graph operations in GNNs on real graphs by as much as 8.6X (3.5X on average).

artificial intelligence, graph neural networks↗

Non linear viscoelasticity applied for the study of durability of polymer matrix composites

A methodology is described for the durability analysis of polymer matrix composites, based on nonlinear viscoelasticity theory. The durability analysis is performed on the basis of a certain number of tests carried out on limited and, if possible, short time scale by the use of accelerating factors. The method was applied to thermomatrix composites under uniaxial and biaxial loadings, showing satisfactory agreement between the life-time predictions and the published data on real-time behavior.

Cardon, A.↗

Parametric reduced-order modeling for component-oriented treatment and localized nonlinear feature inclusion

Abstract We propose coupling a physics-based reduction framework with a suited response decomposition technique to derive a component-oriented reduction (COR) approach, which is suitable for assembly systems featuring localized nonlinearities. Dependencies on influencing parameters are injected into the reduced-order model (ROM), thus ensuring robustness and validity over a domain of parametric inputs, while capturing nonlinear effects. The implemented approach employs individual component modes to capture localized features while additionally relying on reduced modes of a global nature to approximate the system’s dynamics accurately. The global modes are derived from a linear monolithic system, defined as a result of a coordinate separation scheme, which permits the proposed COR-ROM to naturally couple the response between linear and nonlinear subdomains. The derived low-order representation utilizes a proper orthogonal decomposition projection and is additionally reinforced with the inclusion of a hyper-reduction technique to capture the underlying high-fidelity model response while providing accelerated computations. The resulting approach is exemplified in the synthetic case studies of a four-story shear frame with multiple nonlinear regions driven by hysteresis and a large-scale kingpin connection featuring plasticity.

Vlachas, Konstantinos (ORCID:000000029124219X)↗

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗

A skew quadrupole for the AGS to minimize the polarization losses of the polarized beams

The partial helices installed in the Alternating Gradient Synchrotron (AGS) 9 [1] eliminate both, the imperfection and the vertical intrinsic spin resonances thus yielding a 70% polarized proton beam with rigidity 79.4 [Tm] with 2x10 11 protons/bunch at the end of the AGS acceleration cycle. The initial beam polarization at the exit of the 200 MeV Linac is measured to be 80%. The 10% loss of the polarization at the end of the AGS acceleration cycle is due to the horizontal spin resonances which are caused by the betatron oscillations of the beam in the presence of the partial helices. To reduce the effect of these horizontal spin resonances on the beam polarization, the “jump Quads method” is applied in the AGS [2] which eliminates, almost all, the horizontal spin resonances and increases the final polarization of the proton beam to the value of 70%. A study [3] shows that these horizontal spin resonances can be totally eliminated by introducing into the AGS ring, skew quadrupoles which linearly couple the beam motions to excite new horizontal spin resonances which minimize the horizontal-spin-resonances caused by the partial-helices.

43 PARTICLE ACCELERATORS↗

Simulation of the thermoelectric effect in a multi-metallic superconducting cavity

Superconducting radio-frequency accelerating cavities made with different material layers, such as copper, Nb or Nb₃Sn, are susceptible to thermoelectric effects due to differences in Seebeck coefficients between the metals. A temperature gradient across the surfaces can drive thermoelectric currents, which may impact the cavity performance. A layered Cu/Nb/Nb3Sn single-cell cavity was tested with cryocoolers in 2022. Three heaters were mounted on the cavity surface at different locations and three single-axis cryogenic fluxgate magnetometers were attached close to the cavity equator. A linear increase in the magnetic field was measured while increasing the heaters' power. The cavity setup was analyzed with COMSOL and the results showed a trend similar to that observed in the experiment. This contribution details the approach chosen for the simulation and some of the challenges encountered.

Accelerator Physics↗

Flat Spectra of Energetic Particles in Interplanetary Shock Precursors

The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extending upstream of the shock for up to seven flux e-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy as ϵ -1 , consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivity κ on the particle flux F (nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ 3/2 that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗

Electric Drive Technologies Consortium (EDTC)/ Cost competitive, high-Performance, highly Reliable (CPR) Power Devices on 4H-SiC (Final Report)

4H-Silicon carbide (4H-SiC) is a wide bandgap semiconductor that offers superior material properties over silicon, including higher critical electric field, thermal conductivity, and electron saturation velocity. These advantages make 4H-SiC highly attractive for high-voltage, high-efficiency power electronics. However, realizing the full potential of SiC requires device technologies that are not only high-performing but also manufacturable and reliable under real-world operating conditions. This report summarizes the outcomes of a five-year R&D effort funded by the U.S. Department of Energy (DOE) under the Electric Drive Technologies Consortium (EDTC), focused on developing cost-competitive, high-performance, and highly reliable (CPR) power devices on 4H-SiC substrates. The program targeted scalable and manufacturable 1.2 kV-class SiC MOSFETs optimized for next-generation electric vehicles, renewable energy systems, and industrial power conversion. The project delivered transformative advancements in SiC power device performance and ruggedness. Particularly, Specific on-resistance (R on,sp ) was reduced by up to 37%, from ~4.0 m$\Omega \cdot$cm 2 in earlier designs to an industry-leading 2.40 m$\Omega \cdot$cm 2 , driven by optimized doping, refined JFET widths, and layout engineering. Breakdown voltages (BV) exceeded 1600 V, marking improvement over legacy baselines, and demonstrating the robustness of newly implemented junction profiles and edge terminations. Short-circuit withstand time (SCWT) saw a remarkable 4$\times$ increase, from ~2 $\mu$s to over 8 $\mu$s, achieved through the successful deployment of deep P-well structures (~1.8–2.0 $\mu$m) via channeling implantation. This innovative process breakthrough enabled precise junction formation without MeV-class implantation tools, reduced leakage under high field stress, and allowed even the shortest-channel devices (down to 0.3 $\mu$m) to achieve both high BV and excellent ruggedness—breaking the traditional trade-off between conduction efficiency and blocking capability. Several novel architectures pushed the performance envelope further. JBSFETs—featuring embedded Schottky portions—eliminated bipolar degradation and drastically reduced third-quadrant leakage, while Ladder MOSFETs introduced a clever orthogonal conduction path that achieved a 15.4% reduction in R on,sp over standard linear designs. Switching performance reached new benchmarks: short-channel devices showed a 31% reduction in total switching energy compared to 0.5 $\mu$m counterparts, while maintaining manageable gate drive requirements. Layout-optimized structures not only improved transconductance but also accelerated switching transitions, pointing to real-world benefits in converter-level efficiency. The devices also passed rigorous reliability validation. Stress-tested across TDDB, HTGB, HTRB, HVP, and burn-in, the devices screened under 30 V/10 hr and 43 V/1 s protocols consistently exhibited tighter lifetime distributions and long-term oxide robustness. These screening techniques proved effective in identifying latent defects and ensuring deployment-grade reliability. Meanwhile, advanced 3D TCAD simulations revealed and resolved electric field hotspots—particularly in HEXFET corners—where fields exceeding 4.8 MV/cm were mitigated through geometry-aware layout corrections. Overall, the results of this project demonstrate a manufacturable and scalable SiC power device platform that addresses key DOE performance targets for efficient, robust, and reliable 1.2kV 4H-SiC Power Devices. The developed technologies represent a meaningful step forward in the commercial readiness of high-voltage SiC solutions and provide a strong foundation for continued advancement in wide bandgap power electronics.

42 ENGINEERING↗

On possible Mn-53 heterogeneity in the early solar system

The effects of influence of shock wave propagation on the energy spectrum of accelerated particles that lead to different production rates of radionuclides, in particular, Mn-53, on small scales in the early solar system are shown. Search for evidence for extinct Mn-53 has stimulated investigations of Cr isotope anomalies in meteorites. The linear correlation between the magnitude of the Cr-53* excesses and the Mn/Cr ratio that unambiguously proves the in situ decay of Mn-53 was detected, really, in different mineral phases of some carbonaceous and enstatite chondrites, primitive achondrites, pallasites and iron meteorites. However, the data on the Cr-53* excess rarely defines a single linear array on a Mn-53-Cr-52 evolution diagram even for meteorites of the same chemical group. A clear isochron with Mn-53/Mn-55 = 4.4 plus or minus 1.0 x 10(exp -5) (in range of approximately 2.4 to approximately 9 x 10(exp -5)) is observed for CAI of the Allende C3-chondrite while the data for the Murchison C2- and Orgueil C1-chondrites fall much lower corresponding rather to Mn-53/Mn-55 less than 2 x 10(exp -5). In the case of iron meteorites it ranges from less than 5 x 10(exp -8) to less than 5 x 10(exp -5).

Lavrukhina, A. K.↗

Linear Proof Mass Actuator

This paper describes the mechanical design, analysis, fabrication, testing, and lessons learned by developing a uniquely designed spaceflight-like actuator. The linear proof mass actuator (LPMA) was designed to attach to both a large space structure and a ground test model without modification. Previous designs lacked the power to perform in a terrestrial environment while other designs failed to produce the desired accelerations or frequency range for spaceflight applications. Thus, the design for a unique actuator was conceived and developed at NASA Langley Research Center. The basic design consists of four large mechanical parts (mass, upper housing, lower housing, and center support) and numerous smaller supporting components including an accelerometer, encoder, and four drive motors. Fabrication personnel were included early in the design phase of the LPMA as part of an integrated manufacturing process to alleviate potential difficulties in machining an already challenging design. Operating testing of the LPMA demonstrated that the actuator is capable of various types of load functions.

Holloway, Sidney E., III↗

Linear Proof Mass Actuator

This paper describes the mechanical design, analysis, fabrication, testing, and lessons learned by developing a uniquely designed spaceflight-like actuator. The Linear Proof Mass Actuator (LPMA) was designed to attach to both a large space structure and a ground test model without modification. Previous designs lacked the power to perform in a terrestrial environment while other designs failed to produce the desired accelerations or frequency range for spaceflight applications. Thus, the design for a unique actuator was conceived and developed at NASA Langley Research Center. The basic design consists of four large mechanical parts (Mass, Upper Housing, Lower Housing, and Center Support) and numerous smaller supporting components including an accelerometer, encoder, and four drive motors. Fabrication personnel were included early in the design phase of the LPMA as part of an integrated manufacturing process to alleviate potential difficulties in machining an already challenging design. Operational testing of the LPMA demonstrated that the actuator is capable of various types of load functions.

Holloway, S. E., III↗

McMillan electron lens in a system with space charge

Space charge (SC) forces of a circulating beam in a ring have both linear (defocusing) and nonlinear components, due to a nonuniform beam distribution. The linear component of SC forces produces a betatron tune shift, which is the largest for a zero-amplitude particle, while the nonlinear component produces an amplitude-dependent betatron tune spread. These SC effects are responsible for several undesirable phenomena in accelerators: emittance growth, particle losses, beam halo, etc. In this paper, we investigate the possibility to mitigate the distributed SC forces by a thin McMillan lens, providing an axially-symmetric kick, which is qualitatively opposite to the accumulated effect of beam's own SC. Experimentally, the proposed concept can be tested in Fermilab's IOTA ring. A thin McMillan lens can be implemented by a short (70) insertion of an electron beam with a specifically chosen density distribution in transverse directions. In this article, to test if McMillan lenses can reduce the tune spread induced by SC, we make several simulations with a 6-D particle tracking code, Synergia. We choose such beam and lattice parameters that the SC tune spread is roughly 0.5 and the emittance growth due to the half-integer resonance is clearly observed without the SC compensation. Then, we focus on reducing the emittance growth by adjusting the bare betatron tunes using the ring quadrupoles, and reducing the tune spread by the McMillan lenses. The results of reducing a large tune spread (≈ 0.5), reported here, are not perfect, but substantial. There is still room for further investigation. The simulations performed so far in this work indicate that McMillan lenses can cope with an SC tune spread of ≈ 0.1 per lens.

43 PARTICLE ACCELERATORS↗

Medium-energy electrons and heavy ions in Jupiter's magnetosphere - Effects of lower hybrid wave-particle interactions

A theory of medium-energy (about keV) electrons and heavy ions in Jupiter's magnetosphere is presented. Lower hybrid waves are generated by the combined effects of a ring instability of neutral wind pickup ions and the modified two-stream instability associated with transport of cool Iogenic plasma. The quasi-linear energy diffusion coefficient for lower hybrid wave-particle interactions is evaluated, and several solutions to the diffusion equation are given. Calculations based on measured wave properties show that the noise substantially modifies the particle distribution functions. The effects are to accelerate superthermal ions and electrons to keV energies and to thermalize the pickup ions on time scales comparable to the particle residence time. The S(2+)/S(+) ratio at medium energies is a measure of the relative contribution from Iogenic thermal plasma and neutral wind ions, and this important quantity should be determined from future measurements. The theory also predicts a preferential acceleration of heavy ions with an accleration time that scales inversely with the root of the ion mass. Electrons accelerated by the process contribute to further reionization of the neutral wind by electron impact, thus providing a possible confirmation of Alfven's critical velocity effect in the Jovian magnetosphere.

Barbosa, D. D.↗