Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear operator theory”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Status of EM Pump Modeling Capability for VTR

SAS4A/SASSYS-1 is a simulation tool used to perform deterministic analyses of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. Most recently, SAS4A/SASSYS-1 has been selected as the safety analysis software for the Versatile Test Reactor (VTR), a new materials testing Sodium Fast Reactor (SFR) to be built by the Department of Energy. In order to support analysis of the VTR, which will utilize ElectroMagnetic (EM) pumps as the primary coolant pumps, the development of new EM pump modeling capabilities has been initiated for SAS4/SASSYS-1. The new physics-based EM pump model captures systems-level transient performance based on equivalent circuit theory and is highly flexible to model a wide range of electromagnetic pumps of the three-phase linear induction type subject to changes in power supply and plant conditions. Within the equivalent circuit framework, two sub-models are developed to provide options for analysts both with and without a detailed geometric characterization of their pump such that minimal information is required to achieve a realistic working model, but specific details of a particular pump design may be captured if sufficient information is provided. In addition, heat generated from pump operation and its transport through the plant system can be modeled to better capture the behavior of certain transients. This paper will outline the model and derive the governing equations of the new EM pump model. An initial comparison between the simple EM pump model and experimental data is made. This comparison shows that the equivalent circuit model can capture pump performance, with greater accuracy near the rated conditions. In order to demonstrate the transient behavior of the new model, a simple transient with an EM pump replacing the centrifugal pumps in the Advanced Burner Test Reactor is included.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Filtering states with total spin on a quantum computer

Starting from a general wave-function described on a set of spins or qubits, we propose several quantum algorithms to extract the components of this state on eigenstates of total spin S 2 and its z component S z . The method plays the role of total spin projection and gives access to the amplitudes of the initial state on a total spin basis. Here, the different algorithms have various degrees of sophistication depending on the requested tasks. They can either solely project onto the subspace with good total spin or completely uplift the degeneracy in this subspace. In the former case, when the projection on the total spin j is made, we show that the number of operations for the projection can be reduced from a quadratic to a linear dependence in j. After each measurement, the state collapses to one of the spin eigenstates that could be used for postprocessing. For this reason, we call the method total quantum spin filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗
DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Performance evaluation of a Terry GS-2 steam impulse turbine with air-water mixtures

Terry steam turbines are widely used in various industries because of their robust design. Within the nuclear power generation industry, they are used in the Reactor Core Isolation Cooling System to remove decay heat during reactor isolation events. During the Fukushima Daiichi nuclear power station disaster in Japan in 2011, the Reactor Core Isolation Cooling System and associated Terry turbine operated for over 70 hours in Unit 2; this runtime is well beyond the expected operating duration. Theories suggest the turbine was subjected to a two-phase inlet flow, which could degrade the turbine performance. In this work, an experimental test rig was constructed to test a full-scale Terry model GS-2 steam turbine under two-phase air/water flows. Steady-state efficiency and torque performance maps of the turbine were developed over a range of turbine inlet pressures (1.38–4.83 bar or 20–70 psia), air mass fractions (0.05–1.0) and rotational speeds up to 4000 RPM. Furthermore, turbine performance followed expected trends with torque varying linearly and efficiency varying quadratically with rotational speed. In addition, high-speed images of the two-phase flow entering the turbine were also analyzed to understand how changes in inlet pressure and air mass fraction affect the flow regime and homogenization. The present tests with air–water two-phase mixtures are an important step towards providing an understanding of the full-scale Terry turbine’s behavior and performance curves under two-phase conditions. The results of this work will be combined with air/water and steam/water data gathered using a small-scale Terry ZS-1 steam turbine in order to understand the scaling relationship between large and small size Terry turbines and fluid pairs. The combined data set will enable further development of analytical models over a wide range of conditions and may be used to provide technical justification for expanded use of the Terry turbines in nuclear power plant safety systems and other systems.

42 ENGINEERING↗
DOE OSTI2024

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces

Here, this paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Wind Turbine Main Bearing Rating Lives as Determined by IEC 61400-1 and ISO 281

This presentation studies the rating lives of wind turbine main bearings, as determined by the IEC 61400-1 and ISO 281 standards. A brief review of relevant bearing life theory and turbine design requirements is provided. This includes a discussion on possible shortcomings, including the existence (or not) of the bearing fatigue load limit and the validity of assuming linear damage accumulation. A detailed study is then undertaken to determine rating lives for two models of main bearing in a 1.5 MW wind turbine. Rating life assessment is carried out under different conditions, including various combinations of main bearing temperature, wind field characteristics, lubricant viscosity and contamination levels. Rating lives are found to be sufficiently above the desired 20 year design life for both bearing models under expected operating conditions. For the larger bearing, operational loads are shown to be below or close to the bearing fatigue load limit a vast majority of the time. Key sensitivities for rating life values are shown to be the temperature and contamination. Overall, the results of this study suggest that rating life assessment does not account for reported rates of main bearing failures in 1 to 3 MW wind turbines. In future work, it is recommended that efforts be undertaken to identify principal root causes of main bearing failures, possibly leading to a new application standard specific to this component. It is also recommended that impacts of partial wake impingement on main bearing rating lives are investigated.

ENGINEERING,WIND ENERGY↗
DOE OSTI2023

Wind turbine main bearing rating lives as determined by IEC 61400‐1 and ISO 281: A critical review and exploratory case study

Abstract This paper studies the rating lives of wind turbine main bearings, as determined by the IEC 61400‐1 and ISO 281 standards. A critical review of relevant bearing life theory and turbine design requirements is provided, including discussion on possible shortcomings such as the existence (or not) of the bearing fatigue load limit and the validity of assuming linear damage accumulation. A detailed exploratory case study is then undertaken to determine rating lives for two models of main bearing in a 1.5 MW wind turbine. Rating life assessment is carried out under different conditions, including various combinations of main bearing temperature, wind field characteristics, lubricant viscosity, and contamination levels. Rating lives are found to be sufficiently above the desired 20‐year design life for both bearing models under expected operating conditions. For the larger bearing, operational loads are shown to be below or close to the bearing fatigue load limit a vast majority of the time. Key sensitivities for rating life values are temperature and contamination. Overall, the results of this study suggest that an ISO 281 rating life assessment does not account for reported rates of main bearing failures in 1 to 3 MW wind turbines. It is recommended that a similar analysis be undertaken for ISO/TS 16281 rating lives, along with further efforts to identify principal root causes of main bearing failures in future work, possibly leading to a new application standard specific to this component. It is also recommended that the impacts of partial wake impingement on main bearing rating lives are investigated.

17 WIND ENERGY↗
DOE OSTI2025

Search for CP violation in events with top quarks and Z bosons at $\sqrt{s}$ = 13 and 13.6 TeV

A search for the violation of the charge-parity (CP) symmetry in the production of top quarks in association with Z bosons is presented, using events with at least three charged leptons and additional jets. The search is performed in a sample of proton-proton collision data collected by the CMS experiment at the CERN LHC in 2016–2018 at a center-of-mass energy of 13 TeV and in 2022 at 13.6 TeV, corresponding to a total integrated luminosity of 173 fb –1 . For the first time in this final state, observables that are odd under the CP transformation are employed. Also for the first time, physics-informed machine-learning techniques are used to construct these observables. While for standard model (SM) processes the distributions of these observables are predicted to be symmetric around zero, CP-violating modifications of the SM would introduce asymmetries. Two CP-odd operators $\mathcal{O}$$^{I}_{tW}$ and $\mathcal{O}$$^{I}_{tZ}$ in the SM effective field theory are considered that may modify the interactions between top quarks and electroweak bosons. The obtained results are consistent with the SM prediction within two standard deviations, and exclusion limits on the associated Wilson coefficients of –2.7 < $c$$^{I}_{tW}$ < 2.5 and –0.2 < $c$$^{I}_{tZ}$ < 2.0 and are set at 95 % confidence level. The largest discrepancy is observed in $c$$^{I}_{tZ}$ where data is consistent with positive values, with an observed local significance with respect to the SM hypothesis of 2.5 standard deviations, when only linear terms are considered.

CMS↗
DOE OSTI2022

Learning Optimal Multigrid Smoothers via Neural Networks

Multigrid methods are one of the most efficient techniques for solving large sparse linear systems arising from partial differential equations (PDEs) and graph Laplacians from machine learning applications. One of the key components of multigrid is smoothing, which aims at reducing high-frequency errors on each grid level. However, finding optimal smoothing algorithms is problem-dependent and can impose challenges for many problems. In this paper, we propose an efficient adaptive framework for learning optimized smoothers from operator stencils in the form of convolutional neural networks (CNNs). Here, the CNNs are trained on small-scale problems from a given type of PDEs based on a supervised loss function derived from multigrid convergence theories and can be applied to large-scale problems of the same class of PDEs. Numerical results on anisotropic rotated Laplacian problems and variable coefficient diffusion problems demonstrate improved convergence rates and solution time compared with classical hand-crafted relaxation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Endless Dirac nodal lines in kagome-metal Ni3In2S2

Abstract Topological semimetals are a frontier of quantum materials. In multiband electronic systems, topological band crossings can form closed curves, known as nodal lines. In the presence of spin–orbit coupling and/or symmetry-breaking operations, topological nodal lines can break into Dirac/Weyl nodes and give rise to interesting transport properties, such as the chiral anomaly and giant anomalous Hall effect. Recently, the time-reversal symmetry-breaking induced Weyl fermions are observed in a kagome-metal Co 3 Sn 2 S 2 , triggering interests in nodal-line excitations in multiband kagome systems. Here, using first-principles calculations and symmetry-based indicator theories, we find six endless nodal lines along the stacking direction of kagome layers and two nodal rings in the kagome plane in nonmagnetic Ni 3 In 2 S 2 . The linear dipsersive electronic structure, confirmed by angle-resolved photoemission spectroscopy, induces large magnetoresistance up to 2000% at 9 T. Our results establish a diverse topological landscape of multiband kagome metals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Regularization scheme dependence of the counterterms in the galaxy bias expansion

In this paper we explore how different regularization prescriptions affect the counterterms in the renormalization of the galaxy bias expansion. We work in the context of primordial local non-Gaussianity including non-linear gravitational evolution. We carry out the one-loop renormalization of the field $δ$$^{ρ}_{2}$ (i.e. the square of the matter overdensity field) up to third order in gravitational evolution. Three regularization schemes are considered and their impact on the values of the counterterms is studied. We explicitly verify that the coefficients of the non-boost invariant operators are regularization scheme independent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗
DOE OSTI2022

Analysis of the Feasibility of UAS-Based EMI Sensing for Underground Utilities Detection and Mapping

This paper investigates the feasibility of using a linear current sensing (LCS) technique integrated on an unmanned aerial system (UAS) for detecting and mapping underground infrastructure rapidly and cost-effectively. The LCS technique is based on data from a wide band of electromagnetic induction frequencies (50 kHz to 2 MHz) using a vector magnetic field gradiometer. This technique takes advantage of a slowly decaying secondary magnetic field in order to achieve greater standoff detection distance ($\frac{1}{R^2}$ vs. $\frac{1}{R^6}$ for compact metallic targets during EMI sensing, where R is the distance from a target to the sensor). These secondary magnetic fields are produced by the excite current on long conductors, allowing detection at a distance of 10 meters or more. The system operates between tens of kHz to a few MHz and uses either an active EMI source or existing EM fields to excite this linear current on a long metallic subsurface target. Once excited, these linear currents produce a secondary magnetic field that is detected with an above ground triaxial magnetic field gradiometer. By moving and tracking its geolocation, the system outputs rich datasets sufficient to support high-fidelity forward and inverse EMI models for estimating the depth and orientation of deep underground long linear metallic infrastructure. The system’s hardware and its integration to a UAS system are outlined, along with the formulation of LCS theory, and numerical and experimental data are presented. The results illustrate that the LCS technique offers large standoff detection, is adaptable to UAS, and could be used effectively for detecting deep underground infrastructure such as wires and pipes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2023

Thermomechanical Stress and Creep-Fatigue Analysis of a High-Temperature Prototype Receiver for Heating Particles

This work presents a three-dimensional (3D) thermomechanical model of a prototype-scale enclosed light trapping solar receiver for heating particles. Results of the thermoelastic model are used to estimate receiver lifetime under maximum flux conditions. A computational fluid dynamics (CFD) model is first developed to predict the temperature fields in a multi-panel assembly under steady operating conditions. Solar flux distributions on the receiver are obtained from the software package SolTrace and applied to the 3D thermal model. The subsequent particle heating is captured through a simplified 1D energy balance. Panel reradiation is considered through a surface-to-surface radiation model and natural convection loss to the surrounding air is captured in a representative fluid domain surrounding the receiver. The resulting temperature fields from the CFD analysis are used as inputs for a thermoelastic mechanical model with representative boundary conditions. With the resultant temperature and stress fields, a creep-fatigue damage and lifetime analysis is performed using the linear damage accumulation (LDA) theory. The Manson-Coffin formula and Larson Miller correlation are used to calculate the fatigue and creep, respectively. A maximum damage (corresponding to a 30-year service life) is defined for design assessment. The model was first developed and verified in detail by comparing with published results in the literature (temperature and stress profiles and distributions, and creep/fatigue damage fractions) for tubular solar receivers with supercritical carbon dioxide as the working fluid. It was then implemented to model a planar-cavity receiver with various design parameters. Specifically, three different design geometries are considered, and the results show that a maximum temperature of approximately 1200 K could be reached for each design with the given incident solar flux, with the main difference being the distribution of these temperatures. Preliminary resulting stresses for the small-scale prototype without design optimization vary from 20 MPa to 250 MPa for each design, with the maximum stresses occurring on the front face and concave geometry on the side of the panel. In future work, the developed methodology shown here will be applied to analyze a full-scale (50-150 MWth) receiver.

concentrated solar power↗
DOE OSTI2025

Linear-scaling quadruple excitations in local pair natural orbital coupled-cluster theory

Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.

Auxiliary functions↗