Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

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 919 records · Page 51

Direct detection of dark matter substructure in strong lens images with convolutional neural networks

Strong gravitational lensing is a potential way of uncovering the nature of dark matter, by finding perturbations to images that cannot be accounted for well by modeling the lens galaxy without additional structure, be it subhalos (smaller halos within the smooth lens) or line-of-sight (LOS) halos. Here, we introduce results attempting to infer the presence of substructure from images without requiring an intermediate step in which a smooth model has to be subtracted, using a simple convolutional neural network (CNN). We find that the network is only able to infer the presence of subhalos with greater than 75% accuracy when they have masses of greater than or equal to 5 × 10 9 M ⊙ if they lie within the main lens galaxy. Since less massive foreground LOS halos can have the same effect as higher-mass subhalos, the CNN can probe lower masses in the halo mass function. The accuracy does not improve significantly if we add a population of less massive subhalos. With the expectation of experiments such as Hubble Space Telescope and Euclid yielding thousands of high-quality strong lensing images in the next years, having a way of analyzing images quickly to identify candidates that merit further analysis to determine individual subhalo properties while preventing extensive resources being used for images that would yield null detections could be very useful. By understanding the sensitivity as a function of substructure mass, nondetections could be combined with the information from images with substructure to constrain the cold dark matter scenario, in particular if the sensitivity can be pushed to lower masses.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma

The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. Here, the analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP $β$ (NBI injection power), EP energy (NBI voltage) and the radial location of the EP density profile gradient (NBI radial deposition). The analysis indicates a transition from the soft (local plasma relaxation) to the hard MHD (global plasma relaxation) limit if the simulation EP $β$ $\unicode{x2A7E}$ $0.02$, leading to bursting MHD activity caused by radial AEs overlapping. MHD bursts cause an enhancement of the EP transport showing ballistic-like features as avalanche-like events. Simulations in the soft MHD limit show an increment of the EP density gradient as the EP $β$ increases. On the other hand, there is a gradient upper limit in the hard MHD limit, consistent with the critical-gradient behavior. AEs induce shear flows and zonal current leading to the deformation of the flux surfaces and the safety factor profile, respectively, particularly strong for the simulation in the hard MHD limit. Simulations in the hard MHD regime show a decrease of the AE frequency in the saturation phase; this is caused by the destabilization of a transitional mode between a $9/3-10/3$ TAE and a $9/3$ RSAE that may explain the AE frequency down-sweeping observed in some DIII-D discharges. Reducing the EP energy in the nonlinear simulations leads to a weakening of the plasma perturbation. On the other hand, increasing the EP energy causes the opposite effect. Nonlinear simulations of off-axis NBI profiles indicate a lower plasma perturbation as the EP density gradient is located further away from the magnetic axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Demonstrate Capability of NEAMS Tools to Generate Reactor Kinetics Parameters for Pebble-Bed HTGRs Transient Modeling

The system analysis of anticipated operating occurrences and design basis accident for pebblebed reactor systems requires knowledge of the neutron kinetic parameters. These parameters take into account various types of feedback from the pebble-bed core and are typically integrated into system analysis tools such as the System Analysis Module (SAM) in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tool package. These parameters are typically generated by a higher fidelity full-core coupled neutronics/thermal fluids analysis, such as using Griffin and Pronghorn, in 2D or 3D. The delayed neutron fractions (beta) and the neutron lifetimes (lambda) are typically generated by applying the adjoint solution of the neutron flux. Temperature reactivity coefficients are the most important feedback for all anticipated operating occurrences and design basis accident, then depending on the transient other effects might be important as rector control system worth, s-curves, and Xenon generation. This work package will establish a methodology and process for calculating the reactor kinetics parameters for high-temperature gas-cooled reactor and provide them in a suitable form for system level analyses of typical anticipated operating occurrences and design basis accident.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗

Splitting regions and shrinking islands from higher point constraints

We study constraints from higher-point amplitudes on 2 → 2 scattering in the context of effective field theory (EFT) using the perturbative numerical S-matrix bootstrap. Specifically, we investigate the class of weakly coupled EFTs with amplitudes that obey the hidden zero and split conditions that are known to hold both for Tr(Φ 3 ) theory and for certain string tree amplitudes, including at 4-point the beta function. Requiring the splitting condition for the 5-point amplitude not only fixes nearly all its contact terms, but it also imposes non-linear constraints among the 4-point EFT Wilson coefficients. When included in the bootstrap, the resulting allowed region consistent with positivity is no longer convex but is restricted to a smaller non-convex region — which has a sharp corner near the string beta function! Assuming the absence of an infinite spin tower at the mass gap, the allowed region bifurcates into a trivial region (with states only above a chosen cutoff) and an island that continues to shrink around the string as more constraints are included in the bootstrap. The numerics indicate that in the absence of single-mass infinite spin towers the string beta function is the unique 4-point amplitude compatible with hidden zero and the 5-point splitting constraints. The analysis provides a prototype example for how features of higher-point amplitudes constrain the bootstrap of 4-point amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

Amplicon Sequencing Assessment to Measure Microbial Community Response from Heavy Metal Contaminated Soils in Savannah River site, Tims Branch watershed

Long-term presence of heavy metal contamination in soils is correlated with changes in microbial community structures and can lead to the reduction of indigenous species. Tolerance to soil heavy metal toxicity varies among different microbial communities. It remains unclear whether long term effects caused by metal perturbations are associated with soil microbial community dysregulation in Tims Branch watershed. Research Questions: Does heavy metal contamination in Tims Branch soils decrease microbial species diversity significantly compared to non-contaminated soils? Does the species richness decrease in contaminated samples? between samples? Is diversity in low, mid and high contamination soils significantly different compared to control samples? Hypothesis: The relative abundance and diversity of species are significantly altered in soils which are contaminated with heavy metals in the Tims Branch watershed. Objectives: Use amplicon sequencing technology along with bioinformatics to compare bacterial communities between four sites in Tims Branch watershed. Evaluate the percent relative abundance and diversity across the site. Conclusions: All bacterial communities exposed to different heavy metal concentrations were dominated by four major groups (Proteobacteria, Acidobacteria, Chloroflexi and Verrucomicrobia). The dominant phyla Proteobacteria, Acidobacteria and Chloroflexi accounted for 61 % of the relative abundance. Soil bacterial α-diversity, expressed as observed species richness metric and Shannon's diversity index, was the highest in sample location four (high contamination). The lowest observed was at the control location. Species richness increased when concentrations of heavy metals were localized, inferring that the metals can alter microbial community structure. PERMANOVA analysis between contaminated groups tested significant using the Adonis (F (4,29) =1.9712; p < 0.001). There is a 20 percent variation explained between groups (R2 = 0.214). This analysis proves to be beneficial in detecting microbial communities altered by contaminated soils. Moreover, certain concentrations of heavy metals can alter microbial community structures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Measurement of the Eta Meson Radiative Decay Width via the Primakoff Effect

The ? meson is an interesting tool to study fundamental symmetries in Quantum Chromodynamics (QCD). In particular, its radiative decay width, ? p? Ñ ??q, is an important quantity that can be predicted in the framework of Chiral Perturbation Theory. A precision measurement of this quantity would provide critical inputs to understanding the mixing of the ? and ?1 mesons and extracting constants with wide-ranging applications in low-energy QCD. This decay width has been measured in the past using two different experimental techniques. The more popular technique utilized e`e´ collisions to produce ? mesons through electromagnetic interactions. Today, the Particle Data Group (PDG) averages the results of five such experiments to obtain their currently-accepted value of the decay width as: 0.515?0.018 keV. However the first measurement of this quantity was obtained from a fixed-target experiment that measured the cross section for photoproduction of ? mesons on a nuclear target via the Primakoff effect. Their result of 0.324?0.046 keV shows strong tension with the average of the collider measurements, motivating a new, high precision measurement using the Primakoff method. For this purpose, the PrimEx-eta experiment was conducted in Hall D of the Thomas Jefferson National Accelerator Facility (Jefferson Lab or JLab). The data is currently being analyzed to measure the differential cross section for the photoproduction of ? mesons on a liquid, 4He target. Preliminary results obtained from the analysis of the first phase of the PrimEx-eta experiment show reasonable agreement with the currently-accepted PDG value of the radiative decay width. However, as will be discussed, there are many challenges to this precision measurement which must be studied before any results can be finalized and compared with previous measurements. In parallel to the ? decay width measurement, the PrimEx-eta experiment measured the total cross section for the fundamental, Quantum Electrodynamics (QED) process of Compton scattering from the atomic electrons inside the target. The results obtained from this measurement are in strong agreement with the next-to-leading order QED calculations, and the total combined uncertainties are below 3% for incident photon energies between 7-10 GeV. In addition to providing the first precision measurement of the total Compton scattering cross section within this beam energy range, this measurement verifies the capability of the PrimEx-eta experimental setup to perform absolute cross section measurements at forward angles, and serves as a reference process for the calibration of systematic uncertainties.

Smith, Andrew↗

Small-𝑥 gluon GPD constrained from deeply virtual 𝐽/𝜓 production and gluon PDF through universal-moment parametrization

We phenomenologically constrain the small-𝑥 and small-𝜉 gluon generalized parton distributions (GPDs) with the deeply virtual 𝐽/𝜓 production (DV⁢𝐽/𝜓⁢P) in the framework of GPDs through universal moment parameterization (GUMP). We use a hybrid cross-section formula combining collinear factorization to the next-to-leading order (NLO) accuracy of the strong coupling 𝛼 𝑠 , with corrections from nonrelativistic QCD to account for the power corrections due to the heavy 𝐽/𝜓 mass. We reach reasonable fit to the measured differential cross sections of DV⁢𝐽/𝜓⁢P by H1 at the Hadron-Electron Ring Accelerator as well as forward gluon parton distribution functions (PDFs) from JAM22 global analysis. We find that both NLO and nonrelativistic corrections are significant for heavy vector meson productions. Of course, the gluon GPD we obtain still contains considerable freedom in need of inputs from other constraints, particularly in the distribution-amplitude-like region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Weak structure functions in $\nu_l-N$ and $\nu_l-A$ scattering with nonperturbative and higher order perturbative QCD effects

We study the effect of various perturbative and nonperturbative QCD corrections on the free nucleon structure functions [FiNWI(x,Q2);i=1–3] and their implications in the determination of nuclear structure functions. The evaluation of the nucleon structure functions has been performed by using the MMHT 2014 parton distribution functions (PDFs) parametrization, and the target mass correction (TMC) and higher twist (HT) effects are incorporated following the works of Schienbein et al. and Dasgupta et al., respectively. These nucleon structure functions are taken as input in the determination of nuclear structure functions. The numerical calculations for the νl/ν¯l-A deep inelastic scattering (DIS) process have been performed by incorporating the nuclear medium effects like Fermi motion, binding energy, nucleon correlations, mesonic contributions, shadowing, and antishadowing in several nuclear targets such as carbon, polystyrene scintillator, iron, and lead, which are being used in MINERνA, and in argon nuclei, which is relevant for the ArgoNeuT and DUNE experiments. The differential scattering cross sections d2σAWIdxdy and (dσAWIdx/dσCHWIdx) have also been studied in the kinematic region of the MINERνA experiment. The theoretical results are compared with the recent experimental data of MINERνA and the earlier data of the NuTeV, CCFR, CDHSW, and CHORUS Collaborations. Moreover, a comparative analysis of the present results for the ratio (dσAWIdx/dσCHWIdx), and the results from the Monte Carlo (MC) generator GENIE and other phenomenological models of Bodek and Yang, and Cloet et al., has been performed in the context of the MINERνA experiment. The predictions have also been made for the ν¯l-A cross section relevant for the MINERνA experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗

Giant nonlinear optical responses from photon-avalanching nanoparticles

Avalanche phenomena use steeply nonlinear dynamics to generate disproportionately large responses from small perturbations, and are found in a multitude of events and materials. Photon avalanching enables technologies such as optical phase-conjugate imaging, infrared quantum counting and efficient upconverted lasing. However, the photon-avalanching mechanism underlying these optical applications has been observed only in bulk materials and aggregates, limiting its utility and impact. In this study we report the realization of photon avalanching at room temperature in single nanostructures-small, Tm 3+ -doped upconverting nanocrystals-and demonstrate their use in super-resolution imaging in near-infrared spectral windows of maximal biological transparency. Avalanching nanoparticles (ANPs) can be pumped by continuous-wave lasers, and exhibit all of the defining features of photon avalanching, including clear excitation-power thresholds, exceptionally long rise time at threshold, and a dominant excited-state absorption that is more than 10,000 times larger than ground-state absorption. Beyond the avalanching threshold, ANP emission scales nonlinearly with the 26th power of the pump intensity, owing to induced positive optical feedback in each nanocrystal. This enables the experimental realization of photon-avalanche single-beam super-resolution imaging with sub-70-nanometre spatial resolution, achieved by using only simple scanning confocal microscopy and without any computational analysis. Pairing their steep nonlinearity with existing super-resolution techniques and computational methods, ANPs enable imaging with higher resolution and at excitation intensities about 100 times lower than other probes. The low photon-avalanching threshold and excellent photostability of ANPs also suggest their utility in a diverse array of applications, including sub-wavelength imaging and optical and environmental sensing.

36 MATERIALS SCIENCE↗

Octet and decuplet baryon σ terms and mass decompositions

We present a comprehensive analysis of the SU(3) octet and decuplet baryon masses and σ terms using high-precision lattice QCD data and chiral SU(3) effective theory with finite range regularization. The effects of various systematic uncertainties, including from the scale setting of the lattice data and the regularization prescriptions, are quantified. We find the pion-nucleon and strange nucleon σ terms to be σ πN = 44(3)(3) MeV and σ Ns = 50(6)(1) MeV, respectively. Furthermore, the results provide constraints on the energy-momentum tensor mass decompositions of the SU(3) octet and decuplet baryons, where we find that the trace anomaly and quark-gluon energies decrease for strange baryons due to their larger strange σ terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Crucial roles of eastward propagating environments in the summer MCS initiation over the U.S. Great Plains

Mesoscale convective systems (MCSs) frequently occur over the U.S. Great Plains during summer. An analysis using self-organizing map is conducted to identify four types of summer MCS initiation environments during 2004-2017. The first two types feature favorable large-scale environments at both upper and low levels, while Type-3 has favorable lower-level and surface conditions but unfavorable upper-level circulation, and Type-4 features the most unfavorable large-scale environments for MCS initiation. Despite the unfavorable large-scale environment, the convection-centered environments in Type-4 are favorable for MCS initiation and similar to the first two types, suggesting a role of sub-synoptic disturbances as an MCS precursor. All four types of the MCS initiation delineate a clear eastward propagating feature in many fields, such as upper-level potential vorticity/geopotential height, surface pressure and surface equivalent potential temperature, upstream up to 25°-longitude west of and ~36 hours before the MCS initiation. The propagating environments and local, non-propagating low-level moisture are found to be important in MCS initiation at the foothill of the Rocky Mountains, but over the central Great Plains, it is the coupling of dynamical and moisture anomalies associated with propagating waves that results in the MCS initiation. By tracking MCSs and mid-tropospheric perturbations (MPs), a type of sub-synoptic disturbances with Rocky Mountains origin, ~30% of MPs is associated with MCS initiation, mostly in Type-4. Although MPs are related to a small fraction of MCS initiation, MCSs that are associated with MPs tend to produce more rainfall in a larger area with a stronger convective intensity, suggesting MPs to be a source of intense MCSs in summer.

Song, Fengfei↗

Simulations of Water Flow in Relation to a Steady Criticality in an Unsaturated Alluvial Repository

Oak Ridge National Laboratory (ORNL) is examining the conditions under which a criticality event could be ruled out for large packages placed in an unsaturated alluvial repository. This work specifically examines the water level in a stylized waste package that is assumed (1) to fail on the top, and (2) to empty as water is driven away—not by boiling—but by fluid flow phenomena driven by thermal gradients. This analysis uses the PFLOTRAN code package. Specifically, this year’s work includes examination of the heat rate at which water is driven away for various percolation rates using the same assumptions listed in the 2022 annual report. The STANDARDS (formerly UNF-ST&DARDS) package has predicted dual-purpose canister (DPC) neutron multiplication factors (k-effectives) ranging above 1.1 (unitless) with failed components. The DPC designs were not created with the intent to use for disposal purposes. Emplacing such packages may be worth consideration if it can be shown that criticality is either extremely unlikely or of very low consequence to the repository’s safety case. This work aims to help establish hydrologic boundary conditions for a highly hypothetical, highly stylized subsurface criticality event. The approach to critical is assumed to be slow, producing steady heating rates under a few kilowatts. For a DPC containing capillary media, the simulated power required for dry-out was proportional to the percolation rate but was unaffected by the size of a large breech on top of the DPC. Results for single cases with no capillary media in the DPC suggest that higher infiltration rates may be required in those cases to cause flooding. Bounding mechanisms and assumptions must be identified and characterized further by performing additional perturbations to various input parameters to demonstrate robustness of the results. This might be accomplished by examining responses to power pulses to determine the system transfer function required to enable a reduced order feedback model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗