Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “second phase particles”

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 199 records · Page 11

Smart Pixels: In-pixel AI for on-sensor data filtering

We present a smart pixel prototype readout integrated circuit (ROIC) designed in CMOS 28 nm bulk process, with in-pixel implementation of an artificial intelligence (AI) / machine learning (ML) based data filtering algorithm designed as proof-of-principle for a Phase III upgrade at the Large Hadron Collider (LHC) pixel detector. The first version of the ROIC consists of two matrices of 256 smart pixels, each 25$\times$25 µm\textsuperscript{2} in size. Each pixel consists of a charge-sensitive preamplifier with leakage current compensation and three auto-zero comparators for a 2-bit flash-type ADC. The frontend is capable of synchronously digitizing the sensor charge within 25 ns. Measurement results show an equivalent noise charge (ENC) of $\sim$30e\textsuperscript{-} and a total dispersion of $\sim$100e\textsuperscript{-} The second version of the ROIC uses a fully connected two-layer neural network (NN) to process information from a cluster of 256 pixels to determine if the pattern corresponds to highly desirable high-momentum particle tracks for selection and readout. The digital NN is embedded in-between analog signal processing regions of the 256 pixels without increasing the pixel size and is implemented as fully combinatorial digital logic to minimize power consumption and eliminate clock distribution, and is active only in the presence of an input signal. The total power consumption of the neural network is $\sim$ 300 $\mu$W. The NN performs momentum classification based on the generated cluster patterns and even with a modest momentum threshold, it is capable of 54.4\% – 75.4\% total data rejection, opening the possibility of using the pixel information at 40MHz for the trigger. The total power consumption of analog and digital functions per pixel is $\sim$ 6 $\mu$W per pixel, which corresponds to $\sim$ 1 W/cm\textsuperscript{2} staying within the experimental constraints.

Parpillon, Benjamin↗

New proxies for second-order cumulants of conserved charges in heavy-ion collisions within the EPOS4 framework

Proxies for cumulants of baryon number 𝐵, electric charge 𝑄, and strangeness 𝑆 are usually measured in heavy-ion collisions via moments of net-number distribution of given hadronic species. Since these cumulants of conserved charges are expected to be sensitive to the existence of a critical point in the phase diagram of nuclear matter, it is crucial to ensure that the proxies used as substitutes are as close to them as possible. Hence, we use the EPOS 4 framework to generate Au + Au collisions at several collision energies of the BNL Relativistic Heavy Ion Collider beam energy scan. We compute second-order net cumulants of 𝜋, 𝐾, and 𝑝, for which experimental data have been published as well as the corresponding conserved charge cumulants. We then compare them with proxies, defined in previous lattice QCD and hadron resonance gas model studies, which are shown to reproduce more accurately their associated conserved charge cumulants. We investigate the impact of hadronic rescatterings occurring in the late evolution of the system on these quantities, as well as the amount of signal actually originating from the bulk medium which endures a phase transition.

Physics↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗

Nb SRF cavities applied to searches for beyond the standard model physics

The SRF technology developed for accelerators can be successfully applied to new applications, including quantum computing, dark matter searches and beyond the standard model physics. The ultra-high quality factor of SRF cavities can al-low to achieve unprecedented sensitivity in fields outside of the usual accelerator applications, for examples fundamental and beyond the standard model physics. Applications of SRF cavities for gravitational waves searches are also being investigated. The SQMS Physics and Sensing thrust is working on this effort as it strives to combine SRF cavities with QIS, with a focus on BSM and fundamental physics. Using Nb SRF cavities, we run a proof of principle search for the dark photon dark matter candidate at fixed frequency demonstrating the deepest sensitivity to wavelike dark mater [1]. The Dark SRF experiment uses two 1.3GHz SRF cavities to conduct a light shining through wall (LSW) experiment to attempt producing and detecting dark photons in the laboratory. The results of this first superconducting LSW search were recently published [2] and we are now working on the second phase of the experiment, which will be run in a dilution refrigerator using 2.6GHz SRF cavities. In addition to the limits we have achieved, we will also present the experiments we are currently developing, such as the heterodyne axion DM search [3, 4], the SERAPH series of experiments, and gravitational waves searches based on SRF cavities as proposed in [5-8].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

KIPT ADS Target Cooling Analyses by MCNP/Fluent

This report presents multi-physics analyses of the target cooling system of the accelerator driven system (ADS) of the Kharkiv Institute of Physics and Technology (KIPT) using MCNP and Fluent computer programs. MCNP has been used to transport electrons, gammas, and neutrons, and to calculate the energy deposition in the target using a mesh-tally. The mesh-tally data have been imported in Fluent by a C subroutine that has been compiled and linked to Fluent as a user defined function. The KIPT ADS is in operation and Fluent used the CAD model from the manufacturing process of the target. The Fluent results for the reference case match very well the literature results obtained by STAR-CCM+ during the design phase. Other cases that differ from the reference one have been analyzed; in these cases, it is assumed a malfunction of the electron accelerator control or the primary water coolant pump. The target cooling system operates normally for all the analyzed cases except when the inlet water mass flow rate is decreased. If the water coolant mass flow decreases to 25%, the water maximum temperature exceeds the boiling point. However, the transient analysis showed that the target cooling system can operate for 180 seconds when the inlet water mass flow rate is decreased down by 50%.

42 ENGINEERING↗

Structure and dynamical properties of two-dimensional dusty plasmas on one-dimensional periodic substrates

In the presence of the substrate, various structural and dynamical properties of two-dimensional dusty plasma (2DDP) were investigated using Langevin dynamical simulations. This paper reviews a series of results of the structural and dynamical properties of 2DDP modified by one-dimensional periodic substrates (1DPSs) as follows. First, when the depth of the 1DPS increases gradually, it is found that the static structural order of 2DDP along each potential well of the 1DPS increases first, and then decreases gradually. When the width of the 1DPS increases gradually, the particle arrangement in each potential well of the 1DPS changes from one straight row gradually to two rows, including the stable zigzag structure. Second, when there are two rows in each potential well of the 1DPS, the phonon spectra are split into two branches, corresponding to the breathing and sloshing modes, respectively. In addition, due to the pure repulsion between dust particles, from the obtained phonon spectra, the sloshing wave propagates backward at small wave numbers. Third, the calculated mean-square displacement (MSD) shown that, at the intermediate timescale between the initial ballistic and final diffusive motions, the particle motion tends to be more subdiffusive while the depth of the 1DPS increases. While the width of the 1DPS increases gradually, the long-time diffusive motion first increases, then decreases, and finally increases again, exhibiting the oscillation-like diffusion, due to the stable zigzag structure. Finally, when an external direct-current driving force is applied on all particles of 2DDP, three different depinning phases are discovered, which are the pinned, disordered plastic flow, and moving ordered states, respectively, as the driving force increases from zero. In addition, the continuous/discontinuous property of the phase transition between these different depinning states is investigated, showing that the transition property is modulated by the depth of the applied 1DPS

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗

The pNab experiment and the quest for ever better neutron beam polarization

The Nab and pNAB collaborations are conducting a program of studies of free neutron beta decay, with the primary goal of testing the unitarity of the Cabbibo-Kobayashi-Maskawa matrix that describes quark mixing due to the weak interaction. For this purpose, a large, novel electromagnetic spectrometer, the Nab spectrometer, has been designed, built, and placed in use to determine the correlation coefficients in unpolarized neutron beta decay: "a", the neutrino-electron correlation coefficient; and "b", the Fierz term. The subject of this paper is pNAB, the second phase of the program, that will deploy the same spectrometer with a polarized neutron beam to determine "A", the beta asymmetry; and "B", the neutrino asymmetry coefficients. A focus of this paper will be on the strategies to provide a high and precisely known neutron beam polarization.

Baessler, Stefan↗

Apex hydrogen bonds in dendron assemblies modulate close-packed mesocrystal structures

The close-packed mesocrystal structures from soft-matter assemblies have recently received attention due to their structural similarity to atomic crystals, displaying various sphere-packing Frank–Kasper (FK) and quasicrystal structures. Herein, diverse mesocrystal structures are explored in second-generation dendrons (G2-X) designed with identical wedges, in which the terminal functionalities X = CONH 2 and CH 2 NH 2 represent two levels of the strong and weak hydrogen-bonding apexes, respectively. The cohesive interactions at the core apex, referred to as the core interactions, are effectively modulated by forming heterogeneous hydrogen bonds between these two functional units. For the dendron assemblies compositionally close to each pure component of G2-CONH 2 and G2-CH 2 NH 2 , their own FK A15 and C14 phases dominate other phases, respectively. We show the existence of the wide-range FK σ including the dodecagonal quasicrystal (DDQC) phases from the dendron mixtures between G2-CONH 2 and G2-CH 2 NH 2 , providing an experimental phase sequence of A15–σ–DDQC–C14 as the core interactions are alleviated. Intriguingly, the temperature dependence of particle sizes shows that the high plateau values of particle sizes are maintained equivalently until each threshold temperature (T th ), followed by a prompt decrease above the T th . A decrease in T th by alleviating the core interactions and its composition dependence suggest that the more size-dispersed particles, the more susceptibility to chain exchange with increasing temperature. Here our results on the formation of supramolecular dendron assemblies provide a guide to understand the core-interaction-dependent mesocrystal structures toward the fundamental principle underlying the temperature dependence of their particle sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the Nature of Neutrinos with the Deep Underground Neutrino Experiment (DUNE)

The Deep Underground Neutrino Experiment (DUNE) is a next-generation experimental program designed to study the behaviour of neutrino oscillation. DUNE will utilize a neutrino beam originating at Fermilab, near Chicago, and will leverage a detector at Fermilab (Near Detector) and a detector 1300 km away in South Dakota (Far Detector), south of Saskatchewan. In the first phase of DUNE, its Far Detector will comprise of two 10,000 ton (fiducial) liquid argon (LAr) time-projection chamber (TPC) modules – powerful tracking calorimeter detectors – placed nearly a mile underground. With this large, sensitive, underground detector, DUNE aims to collect a high statistics and pure sample of neutrinos at the Far Detector. This setup also offers the potential to study non-beam physical processes via e.g. neutrinos produced in the atmosphere, supernova neutrino bursts, and/or solar neutrinos, etc. A second phase will aim to add more detector mass and expand the program. The Near Detector will consist of a LAr TPC module as well: critical to constraining systematic uncertainties in the oscillation analysis. However, this LAr TPC will have a novel design using a pixel-based readout instead of the traditional wire-based readout. This and the segmentation of the LAr TPC into multiple units are crucial in mitigating the high multiplicity of neutrino interactions expected in any readout window given its proximity to the beam. The Near Detector will feature additional components and capability beyond the LAr TPC, allowing one to deeply characterize the neutrino flux. Due to the complexity of this experimental program, several smaller-scale prototype detectors have been operating to test, validate, and improve both the technical designs and software for processing and analyzing events. By operating in charged particle test beams or neutrino beams, several of the prototypes are also capable of producing valuable results. Canadian institutions are involved in the realization of the DUNE through efforts with both the Near and Far Detectors and prototypes. DUNE is anticipated to begin operating near the end of this decade/the beginning of the next. This talk will focus on the overall DUNE program, for example its ultimate plans, status, and the efforts with prototypes.

Howard, Bruce [York U., Canada; Fermilab]↗

Statistical Properties of Undulator Radiation

Two experiments were carried out to study the statistical properties of undulator radiation in the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. The first experiment studied the turn-to-turn fluctuations in the power of the radiation generated by an electron bunch. The magnitude of these fluctuations depends on the 6D phase-space distribution of the electron bunch. In IOTA, we demonstrated that this effect can be used to measure some electron bunch parameters, small transverse emittances in particular. In the second experiment, a single electron was stored in the ring, emitting a photon only once per several hundred turns. In this regime, any classical interference-related collective effects were eliminated, and the quantum fluctuations could be studied in detail to search for possible deviations from the expected Poissonian photon statistics. In addition, the photocount arrival times were used to track the longitudinal motion of a single electron and to compare it with simulations. This allowed us to determine several dynamical parameters of the storage ring such as the rf cavity phase jitter and the dependence of the synchrotron motion period on amplitude.

43 PARTICLE ACCELERATORS↗

Improved calibration of rf cavities for relativistic electron beams: Effects of secondary corrections and experimental verification

In the aspect of longitudinal beam bunching, the bunching strength can be controlled by the rf cavity phase and voltage. However, these machine parameters are different from those that interact with the beam itself. In order to gain control of the beam-cavity interaction, cavity calibration must be performed. Furthermore, it relies on fitting the beam energy gain versus cavity phase to a calibration function. Under the conventional assumption of relativistic beam conditions, the calibration function is a first harmonic sinusoidal function (a sinusoidal function with a period of 2π). However, this expression is insufficient for a high-voltage bunching cavity. Due to beam acceleration inside the cavity, an energy bias and a second harmonic function should be included to modify the conventional calibration function, even for a relativistic electron beam. In this paper, we will derive this modification and provide a comparison to both the Coherent Electron Cooling Experiment and the IMPACT T simulation, respectively.

43 PARTICLE ACCELERATORS↗

Constraining the particle-scale diversity of black carbon light absorption using a unified framework

Abstract. Atmospheric black carbon (BC), the strongest absorber of visible solar radiation in the atmosphere, manifests across a wide spectrum of morphologies and compositional heterogeneity. Phenomenologically, the distribution of BC among diverse particles of varied composition gives rise to enhancement of its light absorption capabilities by over twofold in comparison to that of nascent or unmixed homogeneous BC. This situation has challenged the modeling community to consider the full complexity and diversity of BC on a per-particle basis for accurate estimation of its light absorption. The conventionally adopted core–shell approximation, although computationally inexpensive, is inadequate not only in estimating but also capturing absorption trends for ambient BC. Here we develop a unified framework that encompasses the complex diversity in BC morphology and composition using a single metric, the phase shift parameter (ρBC), which quantifies how much phase shift the incoming light waves encounter across a particle compared to that in its absence. We systematically investigate variations in ρBC across the multi-space distribution of BC morphology, mixing state, mass, and composition as reported by field and laboratory observations. We find that ρBC>1 leads to decreased absorption by BC, which explains the weaker absorption enhancements observed in certain regional BC compared to laboratory results of similar mixing state. We formulate universal scaling laws centered on ρBC and provide physics-based insights regarding core–shell approximation overestimating BC light absorption. We conclude by packaging our framework in an open-source Python application to facilitate community-level use in future BC-related research. The package has two main functionalities. The first functionality is for forward problems, wherein experimentally measured BC mixing state and assumed BC morphology are input, and the aerosol absorption properties are output. The second functionality is for inverse problems, wherein experimentally measured BC mixing state and absorption are input, and the morphology of BC is returned. Further, if absorption is measured at multiple wavelengths, the package facilitates the estimation of the imaginary refractive index of coating materials by combining the forward and inverse procedures. Our framework thus provides a computationally inexpensive source for calculation of absorption by BC and can be used to constrain light absorption throughout the atmospheric lifetime of BC.

54 ENVIRONMENTAL SCIENCES↗

Interplanetary Protons versus Interacting Protons in the 2017 September 10 Solar Eruptive Event

We analyze the relativistic proton emission from the Sun during the eruptive event on 2017 September 10, which caused a ground-level enhancement (GLE 72) registered by the worldwide network of neutron monitors. Using the neutron monitor data and interplanetary transport modeling both along and across interplanetary magnetic field (IMF) lines, we deduce parameters of the proton injection into the interplanetary medium. The inferred injection profile of the interplanetary protons is compared with the profile of the >100 MeV γ -ray emission observed by the Fermi Large Area Telescope, attributed to pion production from the interaction of >300 MeV protons at the Sun. GLE 72 started with a prompt component that arrived along the IMF lines. This was followed by a more prolonged enhancement caused by protons arriving at the Earth across the IMF lines from the southwest. The interplanetary proton event is modeled using two sources—one source at the root of the Earth-connected IMF line and another source situated near the solar western limb. The maximum phase of the second injection of interplanetary protons coincides with the maximum phase of the prolonged >100 MeV γ -ray emission that originated from a small area at the solar western limb, below the current sheet trailing the associated coronal mass ejection (CME). A possible common source of interacting protons and interplanetary protons is discussed in terms of proton acceleration at the CME bow shock versus coronal (re-)acceleration in the wake of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗

Flow regime and Reynolds number variation effects on the mixing behavior of parallel flows

The hydraulic single-phase mixing of three parallel rectangular channels is experimentally investigated at various Reynolds numbers (Re) and flow regime combinations. Particle Image Velocimetry results for seven mixing cases are presented and discussed with varying Re combinations ranging from 1,824 to 20,844. While all cases result in the same Re ratio of ~0.69 between the inner and outer flows, two cases represent multi-regime mixing with the inner-outer regime pair of laminar-transitional and transitional-turbulent, while the other 5 cases are all characteristic of turbulent mixing with varying levels of turbulence. The outer channels initially share characteristics with a backward facing step. The center channel is found to initially behave like a slot jet, but then sees a significant increase in velocity decay. This inner flow velocity decay increased dramatically in the laminar-transitional mixing case, whose centerline velocity decay was ~6 times larger than the decay in the turbulent mixing cases. Second order statistics revealed a consistent mixing layer thickness of ~0.1 hydraulic diameters for all the cases but showed more intense shearing in the multi-regime mixing cases. The combined point and thereby the mixing layer length is determined using centerline velocity decay profiles, which show a much more aggressive mixing in multi-regime flows. Multi-regime mixing demonstrated superior characteristics relative to turbulent mixing due to a more dramatic velocity decay in the inner flow and a shorter mixing length. The contributions of this work include communicating the benefits of multi-regime mixing and providing detailed characterization efforts that can serve future efforts for validating computational models. Here this research also lays the groundwork for future studies aimed at achieving high levels of mixing without a severe penalty in pressure drop.

42 ENGINEERING↗

PI Loop Resonance Control for Dark Photon Experiment at 2 K Using a 2.6 GHz SRF Cavity

Two 2.6 GHz SRF cavities are being used for a dark photon search at the vertical test stand (VTS) in FNAL, for the second phase of the Dark SRF experiment. During testing at 2 K the cavities experience frequency detuning caused by microphonics and slow frequency drifts. The experiment requires that the two cavities have the same frequency within the cavity's bandwidth. These two cavities are equipped with frequency tuners consisting of three piezo actuators. The piezo actuators are used for fine-fast frequency tuning. A proportional-integral (PI) loop utilizing the three piezos on the emitter was used to stabilize the cavity frequency and match the receiver cavity frequency. The results from this implementation will be discussed. The integration time was also calculated via simulation.

43 PARTICLE ACCELERATORS↗

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗