Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72

Discrimination of Small Earthquakes and Buried Single-Fired Chemical Explosions at Local Distances (<150 km) in the Western United States from Comparison of Local Magnitude ( M L ) and Coda Duration Magnitude ( M C )

Seismologists distinguish underground nuclear explosions from more commonly occurring earthquakes using moment tensor inversion, high-frequency P/S amplitude ratios, m b :M s comparisons, and P–pP differential travel times. These methods are generally successful for large seismic events (M >3-4) well recorded at regional-to-teleseismic distances (>150 km); however, it is unclear whether they can be modified to work for small events (M <3) well recorded only at local distances (<150 km). Here, we evaluate a recently proposed, local-distance seismic source discriminant—the difference between local magnitude (M L ) and coda duration magnitude (M C )—using seismograms of earthquakes and buried, single-fired chemical explosions recorded in three regions of the western United States. The quantity M L –M C was previously found to be sensitive to source depth, effectively discriminating mine blasts, induced earthquakes, and very shallow tectonic earthquakes from deeper crustal earthquakes. In this study, we report the first evaluation of M L –M C as a depth discriminant using data from buried, single-fired explosions that, unlike the seismic sources studied earlier, are good analogs for underground nuclear explosions. We find that even when using generic, uncalibrated methods of assigning magnitudes, M L –M C separates single-fired explosions and earthquakes. The area under the receiver operating characteristic curve is 0.92 for 19 explosions and 14 earthquakes in Washington, 0.90 for 22 explosions and 90 earthquakes in Wyoming, and 0.99 for three explosions and 149 earthquakes in Nevada. Finally, M L :M C comparisons have the potential to enhance discrimination based on high-frequency P/S amplitudes ratios—which perform less well at local than regional distances—because the two metrics have complementary sensitivities.

58 GEOSCIENCES↗

Unified differentiable digital twin for the IOTA/FAST facility

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or produce insights about accelerator state. One notable example of additional information are gradients of physical observables with respect to design parameters produced by differentiable simulations. The IOTA/FAST facility has recently begun a program to implement and experimentally validate a unified start-to-end differentiable digital twin to serve as a virtual accelerator test stand, allowing for rapid prototyping of new software and experiments with minimal beam time costs. In this contribution we will discuss our plans and progress. Specifically, we will cover the selection and benchmarking of both physics and ML codes, the development of generic interfaces between device models and surrogate or physics-based sections, and the export of the parameters through either a deterministic event loop or a fully asynchronous EPICS soft input/output controller. We will also discuss challenges in model calibration and uncertainty quantification, as well as future plans to support larger proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗

mosartwmpy: A Python implementation of the MOSART-WM coupled hydrologic routing and water management model

mosartwmpy is a Python implementation of the Model for Scale Adaptive River Transport with Water Management (MOSART-WM). This new version retains the functionality of the legacy model (written in FORTRAN) while providing new features to enhance user experience and extensibility. MOSART is a large-scale river-routing model used to study riverine dynamics of water, energy, and biogeochemistry cycles across local, regional, and global scales (Li et al., 2013). The WM component introduced by Voisin et al. (2013) represents river regulation through reservoir storage and release operations, diversions from reservoir releases, and allocation to sectoral water demands. Each reservoir release is independently calibrated using long-term mean monthly inflow into the reservoir, long-term mean monthly demand associated with this reservoir, and reservoir goals (flood control, irrigation, recreation, etc.). Generic monthly pre-release rules and storage targets are set up for individual reservoirs; however, those releases are updated annually for inter-annual variability (dry or wet year) and daily for environmental constraints such as flow minimum release and minimum/maximum storage levels. The WM model allows an evaluation of the impact of water management over multiple river basins at once (global, continental scales) and with consistent representation of human operations over the full domain.

97 MATHEMATICS AND COMPUTING↗

Analcime-wairakite formation during experimental cement-bentonite alteration at 200–300 °C

Engineered barrier system materials in nuclear waste repositories may undergo hydrothermal alteration in response to groundwater saturation and heating events over their long operational timescales. Hydrothermal interactions between engineered materials (e.g., bentonite buffers, cements, waste canister materials) and the host rock environment may drive alteration processes that affect the advantageous properties of some barrier materials. However, such alteration will also promote the formation of zeolites, hydrothermally formed minerals that may themselves act to isolate radionuclides. Understanding the environmental conditions that control radionuclide-sorbing properties of zeolites, such as Si/Al ratio or Na content, will be valuable for assessing the changing properties of engineered barrier materials in the case of an in situ heating event. Here, we present experimental work characterizing the formation of zeolites during hydrothermal interactions between generic clay and cement barrier materials, with a focus on the analcime-wairakite zeolite series, which has known radionuclide sorption and exchange properties. We present the results of hydrothermal experiments combining uncured ordinary Portland cement powder with steel and Wyoming bentonite in the presence of Opalinus clay and a synthesized Opalinus clay groundwater to simulate water-saturated conditions in an argillaceous rock repository. The experiments were conducted isothermally at 200 or 300 °C for 8 to 24 weeks. In contrast to analogous studies that did not include cement reactants, we observed the formation of analcime-group minerals in all experiments. The addition of cement resulted in lower Si/Al ratios in the zeolites compared to similar studies that did not include cement. The presence of cement as a reactant was interpreted to promote analcime formation at 200 °C. At 300 °C, we observed higher calcium and silica concentrations in the aqueous solutions as well as increased wairakite formation and decreased analcime formation compared to the experiments at 200 °C. These results show a fully realized analcime-wairakite solid solution that falls between Si/Al = 2 (ideal) and a trend of analcime minerals that have increasing Si/Al ratios with increasing Na/(Na+Ca). Our results show predictable relationships between the Si/Al ratio and the analcime-wairakite content under hydrothermal conditions and illustrate that repository material interactions may promote the formation of zeolites in the analcime-wairakite solid-solution series during heating events in the subsurface.

58 GEOSCIENCES↗

STrEAMlining EFT Matching

This paper presents STrEAM (SuperTrace Evaluation Automated for Matching), a Mathematica package that calculates all functional supertraces which arise when matching a generic UV model onto a relativistic Effective Field Theory (EFT) at one loop and to arbitrary order in the heavy mass expansion. STrEAM implements the covariant derivative expansion to automate the most tedious step of the streamlined functional matching prescription presented in Ref. [1]. The code and an example notebook are available at this link.

43 PARTICLE ACCELERATORS↗

EFT asymptotics: the growth of operator degeneracy

We establish formulae for the asymptotic growth (with respect to the scaling dimension) of the number of operators in effective field theory, or equivalently the number of S-matrix elements, in arbitrary spacetime dimensions and with generic field content. This we achieve by generalising a theorem due to Meinardus and applying it to Hilbert series---partition functions for the degeneracy of (subsets of) operators. Although our formulae are asymptotic, numerical experiments reveal remarkable agreement with exact results at very low orders in the EFT expansion, including for complicated phenomenological theories such as the standard model EFT. Our methods also reveal phase transition-like behaviour in Hilbert series. We discuss prospects for tightening the bounds and providing rigorous errors to the growth of operator degeneracy, and of extending the analytic study and utility of Hilbert series to EFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bootstrapping Mixed MN Correlators in 3D

The recent emergence of the modern conformal bootstrap method for the study of conformal field theories (CFTs) has enabled the revisiting of old problems in classical critical phenomena described by three-dimensional CFTs. The study of such CFTs with O(m)^n \rtimes S_n O ( m ) n ⋊ S n global symmetry, also known as models, is pursued in this work. Systems of mixed correlators involving scalar operators in two different representations of the global symmetry group are considered. Isolated allowed regions are found in parameter space for various values of m m and n n . These ``islands’’ can be separated into two qualitative groups: those close to the unitarity bound and those further away. As a by-product of our analysis generic tensor structures required to bootstrap any G^n \rtimes S_n G n ⋊ S n theory with G G arbitrary are worked out.

Kousvos, Stefanos Robert↗

Hilbert space fragmentation in a 2D quantum spin system with subsystem symmetries

We consider a 2D quantum spin model with ring-exchange interaction that has subsystem symmetries associated to conserved magnetization along rows and columns of a square lattice, which implies the conservation of the global dipole moment. The model is not integrable, but violates the eigenstate thermalization hypothesis through an extensive Hilbert space fragmentation, including an exponential number of inert subsectors with trivial dynamics, arising from kinetic constraints. While subsystem symmetries are quite restrictive for the dynamics, we show that they alone cannot account for such a number of inert states, even with infinite-range interactions. We present a procedure for constructing shielding structures that can separate and disentangle dynamically active regions from each other. Notably, subsystem symmetries allow the thickness of the shields to be dependent only on the interaction range rather than on the size of the active regions, unlike in the case of generic dipole-conserving systems.

Khudorozhkov, Alexey↗

Natural boundaries for scattering amplitudes

Singularities, such as poles and branch points, play a crucial role in investigating the analytic properties of scattering amplitudes that inform new computational techniques. In this note, we point out that scattering amplitudes can also have another class of singularities called natural boundaries of analyticity. They create a barrier beyond which analytic continuation cannot be performed. More concretely, we use unitarity to show that 2 \to 2 2 → 2 scattering amplitudes in theories with a mass gap can have a natural boundary on the second sheet of the lightest threshold cut. There, an infinite number of ladder-type Landau singularities densely accumulates on the real axis in the center-of-mass energy plane. We argue that natural boundaries are generic features of higher-multiplicity scattering amplitudes in gapped theories.

97 MATHEMATICS AND COMPUTING↗

Skyrmion crystals in the triangular Kondo lattice model

We present a systematic study of the formation of skyrmion crystals in a triangular Kondo Lattice model for generic electron filling fractions. Our results indicate that the four-sublattice chiral antiferromagnetic ordering that was reported more than one decade ago can be understood as the dense limit of a sequence of skyrmion crystals whose lattice parameter is dictated by the Fermi wave-vector. This observation has important implications for the ongoing search of skyrmion crystals in metallic materials with localized magnetic moments.

Wang, Zhentao↗

Probing fractional statistics in quantum simulators of spin liquid Hamiltonians

Recent advances in programmable quantum devices brought to the fore the intriguing possibility of using them to realise and investigate topological quantum spin liquid phases. This new and exciting direction brings about important research questions on how to probe and determine the presence of such exotic, highly entangled phases in a noisy quantum environment. One of the most promising tools is investigating the behaviour of the topological excitations, and in particular their fractional statistics. In this work we put forward a generic route to achieve this, and we illustrate it in the specific case of \mathbb{Z}_2 ℤ 2 topological spin liquids implemented with the aid of combinatorial gauge symmetry. We design a convenient architecture to study signatures of fractional statistics via quasiparticle interferometry, and we assess its robustness to diagonal and off-diagonal disorder, as well as to dephasing - effects that are generally pervasive in current quantum programmable devices. Interestingly, when turned on its head, our scheme provides a remarkably clear test of the “quantumness” of these devices, with robust signatures that crucially hinge on quantum coherence and quantum interference effects, and cannot be mimicked by classical stochastic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral anomalies and broken symmetries in maximally chaotic quantum maps

Spectral statistics such as the level spacing statistics and spectral form factor (SFF) are widely expected to accurately identify “ergodicity,” including the presence of underlying macroscopic symmetries, in generic quantum systems ranging from quantized chaotic maps to interacting many-body systems. By studying various quantizations of maximally chaotic maps that break a discrete classical symmetry upon quantization, we demonstrate that this approach can be misleading and fail to detect macroscopic symmetries. Notably, the same classical map can exhibit signatures of different random matrix symmetry classes in short-range spectral statistics depending on the quantization. While the long-range spectral statistics encoded in the early time ramp of the SFF are more robust and correctly identify macroscopic symmetries in several common quantizations, we also demonstrate analytically and numerically that the presence of Berry-like phases in the quantization leads to spectral anomalies, which break this correspondence. Finally, we provide numerical evidence that long-range spectral rigidity remains directly correlated with ergodicity in the quantum dynamical sense of visiting a complete orthonormal basis.

Shou, Laura [Univ. of Maryland, College Park, MD (↗

Anomalies of coset non-invertible symmetries

Anomalies of global symmetries provide important information on the quantum dynamics. We show the dynamical constraints can be organized into three classes: genuine anomalies, fractional topological responses, and integer responses that can be realized in symmetry-protected topological (SPT) phases. Coset symmetry can be present in many physical systems including quantum spin liquids, and the coset symmetry can be a non-invertible symmetry. We introduce twists in coset symmetries, which modify the fusion rules and the generalized Frobenius-Schur indicators. We call such coset symmetries twisted coset symmetries, and they are labeled by the quadruple (G, K, ω D+1 , α D ) in D spacetime dimensions where G is a group and K ⊂ G is a discrete subgroup, ω D+1 is a (D + 1)-cocycle for group G, and α D is a D-cochain for group K. We present several examples with twisted coset symmetries using lattice models and field theory, including both gapped and gapless systems (such as gapless symmetry-protected topological phases). We investigate the anomalies of general twisted coset symmetry, which presents obstructions to realizing the coset symmetry in (gapped) symmetry-protected topological phases. We show that finite coset symmetry G/K becomes anomalous when G cannot be expressed as the bicrossed product G = H ⋈ K , and such anomalous coset symmetry leads to symmetry-enforced gaplessness in generic spacetime dimensions. We illustrate examples of anomalous coset symmetries with A 5 / \Z 2 symmetry, with realizations in lattice models.

Hsin, Po-Shen [King’s College London (United Kingd↗

NuHepMC: A standardized event record format for neutrino event generators

Simulations of neutrino interactions are playing an increasingly important role in the pursuit of high-priority measurements for the field of particle physics. A significant technical barrier for efficient development of these simulations is the lack of a standard data format for representing individual neutrino scattering events. We propose and define such a universal format, named NuHepMC, as a common standard for the output of neutrino event generators. The NuHepMC format uses data structures and concepts from the HepMC3 event record library adopted by other subfields of high-energy physics. These are supplemented with an original set of conventions for generically representing neutrino interaction physics within the HepMC3 infrastructure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Innovations in optimization and control of accelerators using methods of differential geometry and genetic algorithms (Final Report)

Online tuning of particle accelerators is necessary in order to achieve optimal machine performance. However, it is also a major challenge due to the large parameter space which must be searched and the fact that many of the desir- able objectives compete with one another, and so will not reach their optimal values simultaneously. In order to mitigate these issues, we have explored using dimension-reduction techniques to reduce the size of the parameter space which must be searched and multi-objective genetic algorithms to obtain the sets of tuning parameters which provide pareto-optimal values for the objectives. These methods have enabled us to obtain improved values of the vertical emittance at the Cornell Electron Storage Ring (CESR), and to do so with greater control of orbit errors. Algorithmic tuning is a multidisciplinary endeavour, requiring expertise in beam dynamics, diagnostics, control systems and computer science, and thus a key practical problem is to formulate a common language in which experts with different specialties can communicate. We developed a solution to this problem in the form of a generic accelerator software interface that allows for rapid prototyping of optimization and control algorithms. Our interface is built on the Experimental Physics and Industrial Control Systems (EPICS) and makes possible testing control code in simulation before deployment on real accelerators, as well as deployment of third-party optimization code.

43 PARTICLE ACCELERATORS↗

SUNDIALS Multiphysics+MPIManyVector Performance Testing

In this report we document performance test results on a SUNDIALS-based multiphysics demonstration application. We aim to assess the large-scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically: 1. SUNDIALS’ new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources. 2. ARKode’s new multirate integration module, MRIStep, allowing high-order accurate calculations that subcycle "fast" processes within "slow" ones. 3. SUNDIALS’ new flexible linear solver interfaces, that allow streamlined specification of problem-specific linear solvers. 4. SUNDIALS’ new N_Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI_Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS’ generic SUNNonlinearSolver interface and accelerator-enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).

97 MATHEMATICS AND COMPUTING↗

High-Dimensional Spectral Sampling

A wide variety of applications in science and engineering depend on effectively and efficiently exploring high- dimensional parameter spaces. Common examples are uncertainty quantification, inverse problems, or various optimization approaches. In such applications, the first step is to create an initial uniform, random sampling of the space to create a baseline of knowledge. The simulations evaluated at these samples are then used for subsequent processing. In its most generic form, the goal of sampling is to produce the maximal amount of information with the minimal number of samples. However, despite being the first step and thus influencing the performance of the entire workflow, the quality of the initial sample designs is rarely analyzed or optimized. This project develops a new spectral sampling theory for analyzing and creating experiment designs in high-dimensional spaces. Further, new algorithms for estimating sample properties and sample synthesis are designed. Finally, applications of the proposed spectral sampling tools in surrogate modeling, scientific machine learning and image analysis are presented.

97 MATHEMATICS AND COMPUTING↗

Analysis of Differential Financial Impacts of LWR Load-Following Operations. Revision 1

This document reports the application of HERON (Holistic Energy Resource Optimization Network), a newly-developed RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization, to a case study conducted in collaboration with Arizona Public Service (APS). The study is a work in progress; this report presents a description of current models and corresponding results for a particular consideration. APS is anticipating several operational challenges. The first is rapid growth of Variable Renewable Energy (VRE) sources such as wind and solar power providing energy to the grid in the APS service region. This creates challenges for baseload operation of the Palo Verde Generating Station (PVGS), a nuclear power plant co-owned and operated by APS. The second challenge for APS is a change in the cooling water acquisition contract with the Sub Regional Operating Group (SROG). The existing contract is expiring soon, and a renewal is only available for a significantly higher price of water. These two challenges represent a difficulty in maintaining profitability for PVGS as local area demand is highly fluctuant and water prices may increase. Furthermore, excess production at PVGS is currently traded on the Palo Verde energy Hub (PVH), which occasionally offers very low or negative locational marginal prices (LMP) for electricity. This analysis considers the potential benefits of occasional curtailment of PVGS as a method to mitigate unprofitable production hours. Similar to load following, optional daily curtailment would allow APS to elect curtailment of PVGS as low as 70% of its nominal operating level between the hours of 08:00 and 16:00 provided local area demand is met and it is profitable to do so. This document reports a stochastic differential analysis of the Net Present Value (NPV) of PVGS as operated by NPS as a result of optional curtailment operation. To model this differential case, HERON has been designed from previous efforts in this study as a generic tool to accelerate the design and execution of RAVEN workflows for technoeconomic analysis of grid energy systems. HERON leverages both RAVEN and the CashFlow plugin for RAVEN to stochastically explore or optimize a response surface with energy grid system production capacities as inputs and mean NPV as a goal function. HERON performs these evaluations by stochastically sampling market, weather, and demand events using RAVEN’s ARMA reduced-order model capability to generate boundary conditions for dispatch optimization. The algorithm performs the dispatch optimization for each of the synthetic scenarios. For this case, HERON designs workflows to simulate behavior of the PVH, wind, solar, and APS local area demand, and use these synthetic scenarios to stochastically evaluate the differential benefit of optional curtailment for each scenario. While this study finds optional curtailment to be a potential source of revenue, the magnitude of the differential NPV is sufficiently small that it is unlikely a priority for increasing the profitability of PVGS. Further, the differential NPV shrinks throughout the simulated years, suggesting the long-term benefit of such operation given current assumptions is negligible.

14 SOLAR ENERGY↗