Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Breakup”

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 217 records · Page 12

Transverse BBU Suppression with Feedback for Energy-Recovery Linacs

The Electron–Ion Collider (EIC) achieves its design luminosity by cooling the ion beam with a high-current electron beam generated in an energy-recovery linac (ERL). The baseline ERL lattice employs a BNL five-cell cavity that is frequency-scaled to 197 MHz, 591 MHz, and 1.773 GHz, raising concerns about multibunch beam-breakup (BBU) instabilities. Threshold currents for each frequency option are established with two independent BBU tracking codes, providing cross-validated operating margins that guide cavity selection. To further increase the current limit, we incorporate a digital transverse feedback (FB) system that targets the dominant higher-order modes (HOMs). Simulations show the FB raises the BBU threshold by roughly an order of magnitude: power-spectral-density analysis of the beam centroid at the linac exit confirms strong suppression of the aliased HOM peak, although some spectral growth appears at secondary frequencies. This combined study quantifies baseline BBU limits, demonstrates effective active mitigation, and charts a practical path toward robust, high-current operation of the EIC ERL.

Accelerator Physics↗

THE LHeC ERL – OPTICS AND PERFORMANCE OPTIMIZATION

Unprecedently high luminosity of 1034 cm-2s-1, prom-ised by the LHeC accelerator complex poses several beam dynamics and lattice design challenges. Here we present beam dynamics driven approach to accelerator design, which in particular, addresses emittance dilution due to quantum excitations and beam breakup instability in a large scale, multi-pass Energy Recovery Linac (ERL). The use of ERL accelerator technology to provide im-proved beam quality and higher brightness continues to be the subject of active community interest and active accelerator development of future Electron Ion Colliders (EIC). Here, we employ current state of though for ERLs aiming at the energy frontier EIC. The main thrust of these studies was to enhance the collider performance, while limiting overall power consumption through explor-ing interplay between emittance preservation and effi-ciencies promised by the ERL technology.

Bogacz, Slawomir↗

Physics Goals of DWA Experiments at FACET-II

The dielectric wakefield acceleration (DWA) program at FACET produced a multitude of new physics results that range from GeV/m acceleration to the discovery of high field-induced conductivity in THz waves, and beyond, to a demonstration of positron-driven wakes. Here we review the rich program now developing in the DWA experiments at FACET-II. With increases in beam quality, a key feature of this program is extended interaction lengths, near 0.5 m, permitting GeV-class acceleration. Detailed physics studies in this context include beam breakup and its control through the exploitation of DWA structure symmetry. The next step in understanding DWA limits requires the exploration of new materials with low loss tangent, large bandgap, and improved thermal characteristics. Advanced structures with photonic features for mode confinement and exclusion of the field from the dielectric, as well as quasi-optical handling of coherent Cerenkov signals is discussed. Use of DWA for laser-based injection and advanced temporal diagnostics is examined.

43 PARTICLE ACCELERATORS↗

Multibunch Studies for LCLS-II High Energy Upgrade

The Linac Coherent Light Source (LCLS) X-ray free-electron laser at SLAC is being upgraded to LCLS-II with a superconducting linac and 1 MHz bunch repetition rate. The proposed high-energy upgrade (LCLS-II-HE) will increase the beam energy from 4 to 8 GeV, extending the reach of accessible X-ray photon energies. With the increased repetition rate and longer linac of LCLS-II-HE, multi-bunch effects are of greater concern. We use recently introduced capabilities in the beam transport code ELEGANT to study dipole and monopole beam breakup effects for LCLS-II HE beam parameters. The results indicate that resonant dipole kicks have steady-state settle times on the order of 500 bunches or less and appear manageable. We also consider a statistical variation of the cavity frequencies and transverse offsets of cavities and quadrupoles. Resonant emittance growth driven by monopole kicks is found to be disrupted by frequency variation between cavities.

43 PARTICLE ACCELERATORS↗

LINAC OPTICS OPTIMIZATION WITH MULTI-OBJECTIVE OPTIMIZATION

The beamline design of recirculating linacs requires special attention to avoid beam instabilities due to RF wakefields. A proposed high-energy, multi-pass energy recovery demonstration at CEBAF uses a low beam current. Stronger focusing at lower energies is necessary to avoid beam breakup(BBU) instabilities, even with this small beam current. The CEBAF linac optics optimization balances over-focusing at higher energies and beta excursions at lower energies. Using proper mathematical expressions, linac optics optimization can be achieved with evolutionary algorithms. Here, we present the optimization process of North Linac optics using multi-objective optimization.

Neththikumara, I.↗

A Preliminary Feasibility Study on Multi-Cavity Cryomodule Integration for the Electron Ion Collider Energy Recover Linac Cooler

The Electron-Ion Collider (EIC) is a state-of-the-art accelerator designed for collisions between highly polarized electrons and ions. To achieve optimal luminosity, the ion beam is cooled using an electron beam sourced from an energy recovery linac (ERL). In the current ERL baseline design, one BNL type RF cavity is used per cryomodule, leading to spatial and cost challenges. This study examines the feasibility of using more compact PERLE (Powerful Energy Recovery Linac for Experiments) type cavities, which can house multiple cavities within a single cryomodule, by evaluating their Beam Breakup (BBU) instability performance. Higher Order Modes (HOM) parameters were obtained through frequency scaling, assuming constant quality factor Q_L, R/Q, and linearly scaling HOM frequencies. To predict the BBU threshold current and ensure accuracy, two different BBU tracking simulations are used for cross-verification. Although the reduced footprint of PERLE-type cavities is advantageous, maintaining sufficient damping of HOMs remains crucial. Finally, we compare the HOM damping efficiency of both cavity types and suggest a pathway forward.

Deitrick, K.↗

Transverse Stability in an Alternating Symmetry Planar Dielectric Wakefield Structure

Dielectric Wakefield Acceleration (DWA) is a promising technique for realizing the next generation of linear colliders. It provides access to significantly higher accelerating gradients than traditional radio-frequency cavities. One impediment to realizing a DWA-powered accelerator is the issue of the transverse stability of the beams within the dielectric structure due to short-range wakefields. These short-range wakefields have a tendency to induce a phenomenon known as single-bunch beam breakup, which acts as its name implies and destroys the relevant beam. We attempt to solve this issue by leveraging the quadrupole mode excited in a planar dielectric structure and then alternating the orientation of said structure to turn an unstable system into a stable one. We examine this issue computationally to determine the limits of stability and based on those simulations describe a future experimental realization of this strategy.

43 PARTICLE ACCELERATORS↗

Investigation of the background in coherent $J/\psi$ production at the EIC

Understanding fundamental properties of nucleons and nuclei are among the most important scientific goals of the next-generation machine, the Electron-Ion Collider (EIC). With the unprecedented versatility provided by the EIC, it will provide answers to many standing puzzles and open questions in modern nuclear physics. One of the golden measurements at the EIC is coherent vector meson production in electron-nucleus (eA) scattering in order to obtain the spatial gluon density distribution in heavy nuclei. This requires the experiment to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus mostly breaks up into fragments of particles in the far-forward direction close to the hadron beam rapidity. In this report, we systematically study the rejection of incoherent J/\psi J / ψ production by vetoing products from the nuclear breakup - protons, neutrons, and photons, which is modeled with the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design.

Chang, Wan↗

Cathode to Target Simulations for Scorpius: I. Simulation Codes and Models

The new Scorpius linear induction electron accelerator is under development for multi-pulse flash radiography of large, explosively-driven hydrodynamic experiments. Beam physics from the cathode to the target was examined with computer simulations, including envelope, ray-trace, and particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. The results of these simulations is documented in a series of reports. In this report the computer codes and physical models used for these simulations are described. The conclusion of this study is that Scorpius will produce and accelerate a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory if the same engineering standards and construction details are upheld.

43 PARTICLE ACCELERATORS↗

Stagger Tuning Summary [Slides]

Beam Breakup (BBU) is dangerous for high-current linear induction accelerators (LIA) used for flash radiography. Cavity mode RF is coupled to the beam. For frequencies lower than the beam-pipe cutoff, the cavities only communicate via RF on the beam. This cumulative mode BBU grows exponentially with a number of e-foldings proportional to (Number of cavities) X (Beam current) X (Coupling factor) / (Magnetic focusing field).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of a Turbulent Liquid Spray Atomization Model for Diesel Engine Simulations (Final Technical Report)

This project addresses the systematic lack of predictive capabilities by spray models within engine CFD codes. We develop a new modeling approach to predict the breakup of diesel sprays based on recent literature showing that liquid turbulence plays a fundamental role in spray atomization. A new body of quantitative validation data is also developed as a critical element of the project, leveraging the joint capabilities of Georgia Tech’s high-pressure continuous-flow spray chamber and Argonne National Lab’s near-nozzle x-ray diagnostics at the Advanced Photon Source. This project contributes spatially-resolved measurements of drop size distribution within well-characterized diesel injectors, Spray A and D, from the Engine Combustion Network (ECN) to the engine combustion community for the first time. Utilizing this new body of measurements, we validate and demonstrate a new spray model for diesel sprays, termed the KH-Faeth model, that predicts global and local spray characteristic more accurately than the widely adopted and employed KH model. Predicted drop size distributions are seen to predict measured drops sizes both quantitatively and predictively, with accurate response in droplet size distributions over a wide range of ambient density, injection pressure, and injector nozzle size (Spray A and D) without model tuning. The KH-Faeth model can reduce error in the predicted centerline droplet size profile by up to 80% for ECN Spray D simulations when compared to use of the widely employed KH model.

33 ADVANCED PROPULSION SYSTEMS↗

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1

The fuel characterization models of SAS4A/SASSYS-1 (SAS) have recently been extended to include a new U-Pu-Zr metal fuel model, MFUEL. MFUEL is equipped with mechanistic physics based models to predict the pre-transient characterization and transient response of metal fuel, with emphasis on fuel melting, cladding failure, and the metal fuel’s impact on core reactivity. The MFUEL model will be available in the full version of SAS4A/SASSYS-1 5.7, which is scheduled to be released in June 2023. Fast reactor fuel pins that operated in EBR-II and FFTF with low smear density U-Zr and U-Pu-Zr metal fuels and irradiation resistant ferritic-martensitic cladding showed significant advantages in achieving high burnups and assuring inherent safety characteristics during anticipated transients, design basis events, and beyond design basis events. For safety analysis, a fuel performance model must be able to predict (1) Fuel pin mechanics and compositional and dimensional changes, (2) Clad failure, and (3) Fuel pin thermal resistance. Achieving these high level goals accurately is strongly related to the model performance of individual physical processes taking place within a fuel pin during its lifetime. Metal fuels typically operate above the mid-point temperature of melting during normal, as well as off-normal, conditions. At these elevated temperatures, the availability of thermal activation provides a driving force for various diffusional processes leading to complex phase transformations, micro-structure evolution, significant amounts of fuel swelling, interconnected porosity formation, excessive amounts of fission gas release, and fuel clad chemical interactions. Clad failure in fast reactors primarily occurs as a result of creep rupture augmented by clad wastage formation. The reaction is driven by thermal creep induced dislocation motion, grain boundary cavity nucleation, growth and breakup of grain boundaries. The high level complexity and limited available data requires introducing physics-based modeling approaches to gain extrapolation ability and sensitivity with respect to various conditions. The objective of this report is to perform validation of MFUEL using the experimental data for (1) Normal operation EBR-II fuel behavior, (2) Normal operation PHENIX fuel behavior, (3) HT9 Pressure tube ramp-and-hold creep rupture tests, (4) Whole Pin Furnace (WPF) creep strain, creep rupture and eutectic tests, (5) Fuel Behavior Test Apparatus (FBTA) eutectic tests, and (6) TREAT M5-7 OverPower tests up to clad failure. Section-2 includes a brief description of the MFUEL models. A detailed description of the MFUEL physics-based, semi-empirical models will be presented in the SAS V 5.7 theory manual. Section-3, 4, and 5 describes the validation effort for the pre-transient irradiation, furnace transients, and TREAT M-Series transients, respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Primer on Nuclear "Recoil" Data

For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Review of the LANL n-003_Li_006.endf evaluation of Jan 2020 Candidate for ENDF/B-VIII.1

This evaluation is a significant general improvement of descriptions of nuclear reactions of neutrons on Li6. This will be of great value in fusion applications such as the planned IFMIF. This new evaluation does not use the many pseudo-excited states of Li6* to describe the breakup of Li6 in the n+H2+He4 exit channel, as was done in ENDF/B-VII.0 and ENDF/B-VIII.0, but puts all the (uncorrelated) distributions into an MT=32 channel with MF=6 data, alongside the single specific excited state Li6_e1 at 3.562 MeV excitation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗