Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 37 records · Page 2

High gradient tests of metallic mm-wave accelerating structures

This study explores the physics of vacuum rf breakdowns in high gradient mm-wave accelerating structures. We performed a series of experiments with 100 GHz and 200 GHz metallic accelerating structures, at the Facility for Advanced Accelerator Experimental Tests (FACET) at the SLAC National Accelerator Laboratory. This paper presents the experimental results of rf tests of 100 GHz travelling-wave accelerating structures, made of hard copper-silver alloy. The results are compared with pure hard copper structures. The rf fields were excited by the FACET ultra-relativistic electron beam. The accelerating structures have open geometries, 10 cm long, composed of two halves separated by a variable gap. The rf frequency of the fundamental accelerating mode depends on the gap size and can be changed from 90 GHz to 140 GHz. The measured frequency and pulse length are consistent with our simulations. When the beam travels off-axis, a deflecting field is induced in addition to the decelerating longitudinal field. We measured the deflecting forces by observing the displacement of the electron bunch and used this measurement to verify the expected accelerating gradient. We present the first quantitative measurement of rf breakdown rates in 100 GHz copper-silver accelerating structure, which was 10 –3 per pulse, with peak electric field of 0.42 GV/m, an accelerating gradient of 127 MV/m, at a pulse length of 2.3 ns. The goal of our studies is to understand the physics of gradient limitations in order to increase the energy reach of future accelerators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An integrated approach to understanding RF vacuum arcs

Although used in the design and costing of large projects such as linear colliders and fusion tokamaks, the theory of vacuum arcs and gradient limits is not well understood. Almost 120 years after the isolation of vacuum arcs, the exact mechanisms of the arcs and the damage they produce are still being debated. We describe our simple and general model of the vacuum arc that can incorporate all active mechanisms and aims to explain all relevant data. Our four stage model, is based on experiments done at 805 MHz with a variety of cavity geometries, magnetic fields, and experimental techniques as well as data from Atom Probe Tomography and failure analysis of microelectronics. The model considers the trigger, plasma formation, plasma evolution and surface damage phases of the RF arc. This paper also examines how known mechanisms can explain the observed sharp field dependence, fast breakdown times and observed surface damage. We update the model and discuss new features while also pointing out where new data would be useful in extending the model to a wider range of frequencies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Prospects for H-mode inhibition in negative triangularity tokamak reactor plasmas

Instability to high toroidal mode number ($\textit{n}$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity (NT) (δ) shaping, preventing access to strong H-mode regimes when δ is sufficiently negative. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$\textit{n}$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical δ required for avoiding 2nd stability to high-$\textit{n}$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium inverse aspect ratio, elongation and squareness, equilibria with NT are robustly prohibited from accessing the 2nd stability region, offering the prediction that NT reactors should maintain L-mode-like operation. In order to access high-$\textit{n}$ 2nd stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Finally, scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stoichiometry and thickness dependence of superconducting properties of niobium nitride thin films

The current technology used in linear particle accelerators is based on superconducting radio frequency (SRF) cavities fabricated from bulk niobium (Nb), which have smaller surface resistance and therefore dissipate less energy than traditional nonsuperconducting copper cavities. Using bulk Nb for the cavities has several advantages, which are discussed elsewhere; however, such SRF cavities have a material-dependent accelerating gradient limit. In order to overcome this fundamental limit, a multilayered coating has been proposed using layers of insulating and superconducting material applied to the interior surface of the cavity. The key to this multilayered model is to use superconducting thin films to exploit the potential field enhancement when these films are thinner than their London penetration depth. Such field enhancement has been demonstrated in MgB{sub 2} thin films; here, the authors consider films of another type-II superconductor, niobium nitride (NbN). The authors present their work correlating stoichiometry and superconducting properties in NbN thin films and discuss the thickness dependence of their superconducting properties, which is important for their potential use in the proposed multilayer structure. While there are some previous studies on the relationship between stoichiometry and critical temperature T{sub C}, the authors are the first to report on the correlation between stoichiometry and the lower critical field H{sub C1}.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Vacuum RF Breakdown of Accelerating Cavities in Multi-Tesla Magnetic Fields

Ionization cooling of intense muon beams requires the operation of high-gradient, normal-conducting RF structures within multi-Tesla magnetic fields. The application of strong magnetic fields has been shown to lead to an increase in vacuum RF breakdown. This phenomenon imposes operational (i.e. gradient) limitations on cavities in ionization cooling channels, and has a bearing on the design and operation of other RF structures as well, such as photocathodes and klystrons. We present recent results from Fermilab's MuCool Test Area (MTA), in which 201 and 805 MHz cavities were operated at high power both with and without the presence of multi-Tesla magnetic fields. We present an analysis of damage due to breakdown in these cavities, as well as measurements related to dark current and their relation to a conceptual model describing breakdown phenomena.

43 PARTICLE ACCELERATORS↗

Microphonics in the LCLS-II Superconducting Linac

The LCLS-II project has installed a new superconducting linac at SLAC that consists of 35 1.3 GHz cryomodules and 2 3.9 GHz cryomodules. The linac will provide a 4 GeV electron beam for generating soft and hard X-ray pulses. Cavity detuning induced by microphonics was a significant design challenge for the LCLS-II cryomodules. Cryomodules were produced that were within the detuning specification (10 Hz for 1.3 GHz cryomodules) on test stands. Here we present first measurements of the microphonics in the installed LCLS-II superconducting linac. Overall, the microphonics in the linac are manageable with 94% of cavities coming within the detune specification. Only two cavities are gradient limited due to microphonics. We identify a leaking cool down valve as the source of microphonics limiting those two cavities.

43 PARTICLE ACCELERATORS↗

Investigation of Heat Transfer From

The convective heat transfer from the surface of an ellipsoidal forebody of fineness ratio 3 and 20-inch maximum diameter was investigated in clear air for both stationary and rotating operation over a range of conditions including air speeds up to 240 knots, rotational speeds up to 1200 rpm, and angles of attack of 0 deg, 3 deg, and 6 deg. The results are presented in the form of heat-transfer coefficients and the correlation of Nusselt and Reynolds numbers. Both a uniform surface temperature and a uniform input heater density distribution were used. The experimental results agree well with theoretical predictions for uniform surface temperature distribution. Complete agreement was not obtained with uniform input heat density in the laminar-flow region because of conduction effects. No significant effects of rotation were obtained over the range of airstream and rotational speeds investigated. Operation at angle of attack had only minor effects on the local heat transfer. Transition from laminar to turbulent heat transfer occurred over a wide range of Reynolds numbers. The location of transition depended primarily on surface roughness and pressure and temperature gradients. Limited transient heating data indicate that the variation of surface temperature with time followed closely an exponential relation.

Lewis, James P.↗

Explicit upwind algorithm for the parabolized Navier-Stokes equations

A new explicit upwind algorithm based on Roe's flux-difference splitting (FDS) method has been developed for the three-dimensional Parabolized Navier-Stokes (PNS) equations. For three-dimensional flows, FDS's are determined separately for the two nonmarching directions and modified to account for the calculated shock angle in the crossflow plane. Second-order FDS is applied to the pressure and convection terms with the streamwise pressure gradient limited in the subsonic region to maintain a hyperbolic inviscid equation set. Second-order central differencing is obtained in the two-step algorithm for the shear and heat flux terms. The new algorithm is demonstrated for three laminar flow test cases: supersonic flow over a flat plate, hypersonic flow over a 15 deg ramp, and hypersonic flow past a 10 deg cone at a 24 deg angle of attack. The computed results agree well with experimental measurements.

Korte, John J.↗

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

Using Nanoparticle X-ray Spectroscopy to Probe the Formation of Reactive Chemical Gradients in Diffusion-Limited Aerosols

For aerosol particles that exist in highly viscous, diffusion-limited states, steep chemical gradients are expected to form during photochemical aging in the atmosphere. Under these conditions, species at the aerosol surface are more rapidly transformed than molecules residing in the particle interior. To examine the formation and evolution of chemical gradients at aerosol interfaces, the heterogeneous reaction of hydroxyl radicals (OH) on ~200 nm particles of pure squalane (a branched, liquid hydrocarbon) and octacosane (a linear, solid hydrocarbon) and binary mixtures of the two are used to understand how diffusion limitations and phase separation impact the particle reactivity. Aerosol mass spectrometry is used to measure the effective heterogeneous OH uptake coefficient (γ eff ) and oxidation kinetics in the bulk, which are compared with the elemental composition of the surface obtained using X-ray photoemission. When diffusion rates are fast relative to the reaction frequency, as is the case for squalane and low-viscosity squalane–octacosane mixtures, the reaction is efficient (γ eff ~ 0.3) and only limited by the arrival of OH to the interface. However, for cases, where the diffusion rates are slower than reaction rates, as in pure octacosane and higher-viscosity squalane–octacosane mixtures, the heterogeneous reaction occurs in a mixing-limited regime and is ~10× slower (γ eff ~ 0.03). This is in contrast to carbon and oxygen K edge X-ray absorption measurements that show that the octacosane interface is oxidized much more rapidly than that of pure squalane particles. The O/C ratio of the surface (estimated to be the top 6–8 nm of the interface) is measured to change with rate constants of (3.0 ± 0.9) × 10 –13 and (8.6 ± 1.2) × 10 –13 cm 3 molecule –1 s –1 for squalane and octacosane particles, respectively. The differences in surface oxidation rates are analyzed using a previously published reaction-diffusion model, which suggests that a 1–2 nm highly oxidized crust forms on octacosane particles, whereas in pure squalane, the reaction products are homogeneously mixed within the aerosol. This work illustrates how diffusion limitations can form particles with highly oxidized surfaces even at relatively low oxidant exposures, which is in turn expected to influence their microphysics in the atmosphere.

74 ATOMIC AND MOLECULAR PHYSICS↗

The Thick Level-Set model for dynamic fragmentation

The Thick Level-Set (TLS) model is implemented to simulate brittle media undergoing dynamic fragmentation. This non-local model is discretized by the finite element method with damage represented as a continuous field over the domain. A level-set function defines the extent and severity of damage, and a length scale is introduced to limit the damage gradient. Numerical studies in one dimension demonstrate that the proposed method reproduces the rate-dependent energy dissipation and fragment length observations from analytical, numerical, and experimental approaches. In conclusion, additional studies emphasize the importance of appropriate bulk constitutive models and sufficient spatial resolution of the length scale.

97 MATHEMATICS AND COMPUTING↗

Origins of quench in buffered chemical polished and low temperature baked SRF cavities

Electropolishing (EP) and buffered chemical polishing (BCP) are conventional surface preparation techniques for superconducting radiofrequency (SRF) cavities that remove damaged material from the cavity surface. One main issue with EP and BCP treated SRF cavities is high field Q-slope (HFQS), a drop in quality factor at high gradients that limits quench field. High gradient performance in EP cavities can be improved by applying a low temperature bake (LTB), but LTB does not consistently remove HFQS in BCP cavities. There is no consensus as to the why LTB is not effective on BCP prepared cavities, and the cause of HFQS in BCP cavities is not well understood. We examine the origins of quench in EP, BCP, EP+LTB, and BCP+LTB treated SRF cavities. We also show the effect of these treatments on the onset of HFQS, heating within the cavity up to quench, concentration of free hydrogen, and surface roughness.

43 PARTICLE ACCELERATORS↗

Iron transformation mediates phosphate retention across a permafrost thaw gradient

Phosphorus limits primary productivity in many (Sub-)Arctic ecosystems and may constrain biological carbon sequestration. Iron (III) oxides strongly bind phosphate in soils but can dissolve under flooded, reducing conditions induced by permafrost thaw and ground collapse. The ability for iron to regulate phosphate storage and solubility in thawing permafrost landscapes remains unclear. Here, iron-rich sediments containing iron oxides and organic-bound iron were incubated with or without added phosphate in soils along a permafrost thaw gradient to evaluate how iron-phosphate associations respond to thaw-induced redox shifts. Iron oxides partially dissolved and released sorbed phosphate when incubated in soils underlain by degraded permafrost. Iron complexed by organic matter remained stable but provided no phosphate binding capacity. Phosphate addition enhanced iron oxide dissolution and phosphorus concentrations in associated microbial biomass. Our study demonstrates that the capacity for iron oxides to immobilize and retain phosphate in permafrost peatlands decreases with permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗