Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearization”

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 865 records · Page 48

CALiPER Report 21.3: Cost Effectiveness of Linear (T8) LED Lamps

Meeting performance expectations is important for driving adoption of linear LED lamps, but cost-effectiveness may be an overriding factor in many cases. Linear LED lamps cost more initially than fluorescent lamps, but energy and maintenance savings may mean that the life-cycle cost is lower. This report details a series of life-cycle cost simulations that compared a two-lamp troffer using LED lamps (38 W total power draw) or fluorescent lamps (51 W total power draw) over a 10-year study period. Variables included LED system cost ($\$$40, $\$$80, or $\$$120), annual operating hours (2,000 hours or 4,000 hours), LED installation time (15 minutes or 30 minutes), and melded electricity rate ($\$$0.06/kWh, $\$$0.12/kWh, $\$$0.18/kWh, or $\$$0.24/kWh). A full factorial of simulations allows users to interpolate between these values to aid in making rough estimates of economic feasibility for their own projects. In general, while their initial cost premium remains high, linear LED lamps are more likely to be cost-effective when electric utility rates are higher than average and hours of operation are long, and if their installation time is shorter.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Addressing Research and Development Gaps for Plasma-Material Interactions with Linear Plasma Devices

Plasma-material interactions in future fusion reactors have been identified as a knowledge gap to be dealt with before any next step device past ITER can be built. The challenges are manifold. They are related to power dissipation so that the heat fluxes to the plasma-facing components can be kept at technologically feasible levels; maximization of the lifetime of divertor plasma-facing components that allow for steadystate operation in a reactor to reach the neutron fluence required; the tritium inventory (storage) in the plasma-facing components, which can lead to potential safety concerns and reduction in the fuel efficiency; and it is related to the technology of the plasma-facing components itself, which should demonstrate structural integrity under the high temperatures and high neutron fluence. While the dissipation of power exhaust can and should be addressed in high power toroidal devices, the interaction of the plasma with the materials can be best addressed in dedicated linear devices due to their cost effectiveness and ability to address urgent research and development gaps more timely. However, new linear plasma devices are needed to investigate the PMI under fusion reactor conditions and test novel plasma-facing components. Existing linear devices are limited either in their flux, their reactor-relevant plasma transport regimes in front of the target, their fluence, or their ability to test material samples a priori exposed to high neutron fluence. The proposed Material Plasma Exposure eXperiment (MPEX) is meant to address those deficiencies and will be designed to fulfill the fusion reactor-relevant plasma parameters as well as the ability to expose a priori neutron activated materials to plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear and Nonlinear Constitutive Model for Piezoelectricity in ALEGRA-FE

This report develops and documents linear and nonlinear constitutive relations implemented in ALEGRA-FE (ferroelectric). A thermodynamic framework is created to describe the electromechanical system in the form of a free energy functional. Constitutive relations are derived by taking series expansions of the free energy functional with respect to the independent fields. First order expansion terms yield linear constitutive relations and higher order expansion terms yield non-linear constitutive relations. This document serves as supplement to Section 4 of Sandia Report SAND2013-7363, Rev 3. Methods for implementation of kinematic relations of piezoelectric models and rotation of material principal axes are covered in the supplemented report. Additional discussion on phase velocity calculation is also presented.

30 DIRECT ENERGY CONVERSION↗

Modular Linear Thermoviscoelastic Model

Time-dependent, viscoelastic responses of materials like polymers and glasses have long been studied. As such, a variety of models have been put forth to describe the behavior including simple rheological models (e.g. Maxwell, Kelvin), linear "fading memory" theories, and hereditary integral based linear thermal viscoelastic approaches as well as more recent nonlinear theories that are either integral, fictive temperature, or differential internal state variable based. The current work details a new LINEAR_THERMOVISCOELASTIC model that has been added to LAME. This formulation represents a viscoelastic theory that neglects some of the phenomenological details of the PEC/SPEC models in favor of efficiency and simplicity. Furthermore, this new model is a first step towards developing "modular" viscoelastic capabilities akin to those available with hardening descriptions for plasticity models in LAME. Specifically, multiple different (including user-defined) shift-factor forms are implemented with each being easily selected via parameter specification rather than requiring distinct material models.

36 MATERIALS SCIENCE↗

The Relative Importance of Assumed Infrasound Mechanisms on the Linear Inversion of Infrasound Time Series at the Source Physics Experiment

We invert far field infrasound data for the equivalent seismo-acoustic time domain moment tensor to assess the relative importance of two assumed seismoacoustic source mechanisms. The infrasound data were produced by a four of the underground chemical explosions that were conducted during the Source Physics Experiment (SPE). For each SPE event that we invert, we produce three set of atmospheric Green's functions: an average model based on ten years of atmospheric data, as well as two extrema models designed to maximize the variability of atmospheric conditions for the given time-of-day and day-of-year for each SPE event. To parameterize the inversion, we assume that the source of infrasonic energy results from the linear combination of explosion-induced surface spall and linear seismic-to-elastic mode conversion at the Earth's free surface. We find that the inversion yields relatively repeatable results for the estimated spall source whereas the estimated isotropic explosion source is highly variable. This suggests that the majority of the observed acoustic energy is produced by the spall source and/or our modeling of the elastic energy propagation, and data are subsequent conversion to acoustic energy via linear elastic-to-acoustic coupling at the surface, is too simplistic.

58 GEOSCIENCES↗

Linear Vlasov Solver For Microbunching Gain Estimation With Inclusion of CSR, LSC And Geometric Impedances

As is known, microbunching instability (MBI) has been one of the most challenging issues in designs of magnetic chicanes for short-wavelength free-electron lasers or linear colliders, as well as those of transport lines for recirculating or energy recovery linac machines. To more accurately quantify MBI in a single-pass system and for more complete analyses, we further extend and continue to increase the capabilities of our previously developed linear Vlasov solver [1] to incorporate more relevant impedance models into the code, including transient and steady-state free-space and/or shielding coherent synchrotron radiation (CSR) impedances, the longitudinal space charge (LSC) impedances, and the linac geometric impedances with extension of the existing formulation to include beam acceleration [2]. Then, we directly solve the linearized Vlasov equation numerically for microbunching gain amplification factor. In this study we apply this code to a beamline lattice of transport arc [3] following an upstream linac section. The resultant gain functions and spectra are presented here, and some results are compared with particle tracking simulation by ELEGANT [4]. We also discuss some underlying physics with inclusion of these collective effects and the limitation of the existing formulation. It is anticipated that this more thorough analysis can further improve the understanding of MBI mechanisms and shed light on how to suppress or compensate MBI effects in lattice designs.

Tsai, C. -Y.↗

Batched Sparse Linear Algebra (Final Report for Subcontract B648960)

This report finalizes design specifications for developing batched kernels for small tensor operations for unassembled matrix-free iterative solvers, batched solvers for partially assembled operators, and batched solvers with support for various sparse formats. The outcome of the project milestones is a set of interfaces to Batched Sparse LA solvers running on hardware accelerators for use in ECP Libraries and Applications. It is part of the development of sparse batched kernels, solvers/preconditioners as well as creating interoperability in xSDK libraries with sparse and dense batched functions to benefit ECP applications. The participants included representatives from ECP libraries (not limited to the xSDK project), applications, and vendors (AMD, Intel, and NVIDIA). Batched sparse linear algebra solvers form the new frontier for algorithmic development and performance engineering. Many applications (ECP and non-ECP alike) require simultaneous solutions of small linear systems of equations that are structurally sparse. To move towards high hardware utilization, it is important to provide these applications with appropriate interfaces to efficient batched sparse solvers running on modern hardware accelerators. We present interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the software portable between the major hardware accelerators from AMD, Intel, and NVIDIA. The presented interface specifications includes batched band, sparse iterative, and sparse direct solvers. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, SUNDIALS, and SuperLU_dist.

97 MATHEMATICS AND COMPUTING↗

Randomized Algorithms for Linear Solvers

Recently, randomized algorithms in numerical linear algebra, specifically those centered around random sketching, have gained traction in primarily theoretical research due to their potential to significantly reduce problem dimensionality at the cost of an O(1) multiplicative distortion factor. It has been assumed that this sketching can be done efficiently, but thorough investigation into how precisely to do it has been neglected. Moreover, the theory-based community has argued for sketching’s ability to reduce computational cost via complexity analysis, but has not researched how it affects the stability of the algorithms. At Sandia, efficient linear solvers that scale well on modern HPC architectures while maintaining stability are imperative for practical applications. In this LDRD, we developed a random sketching strategy that is substantially faster than existing ones, and demonstrate its superior performance in practice on a NVIDIA H100 GPU. Moreover, we show how this can be used to significantly outperform existing linear least squares solvers while improving the solver’s stability as well. Additionally, we demonstrate how this sketching strategy can be used to make a fast, stable QR factorization that can subsequently be used in s-step and block Krylov solvers. Finally, we incorporate a sketching-based block orthogonalization scheme into s-step GMRES, which is stable and faster than existing approaches on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗

Compensating for Amplifier Non-Linearity in a SEL Controller

The SEL architecture used commonly in CW superconducting linacs is a positive feedback system that requires determining the correct limits for the operating conditions to ensure stable operation. These limits are determined from the results of the amplifier calibration which characterizes the dac drive(forward power) against the cavity field. When the amplifier non-linearity precludes a simple linear fit for this characteristic, its modification to account for the amplifiers non-linearity is essential to operate the system in the SEL architecture. This process is described with test data from a 32kW SSA for the PIP-II linac.

Sankar Raman, Shrividhyaa [Fermilab]↗

Variational Quantum Linear Solver

Previously proposed quantum algorithms for solving linear systems of equations cannot be implemented in the near term due to the re quired circuit depth. Here, we propose a hybrid quantum-classical algorithm, called Variational Quantum Linear Solver (VQLS), for solving linear systems on near-term quantum computers. VQLS seeks to variationally prepare |x$\rangle$ such that A|x$\rangle$ ∝ |b$\rangle$. We derive an operationally meaningful termination condition for VQLS that allows one to guarantee that a desired solution precision ϵ is achieved. Specifically, we prove that C $⩾$ ϵ 2 /κ 2 , where C is the VQLS cost function and κ is the condition number of A. We present efficient quantum circuits to estimate C, while providing evidence for the classical hardness of its estimation. Using Rigetti’s quantum computer, we success fully implement VQLS up to a problem size of 1024 × 1024. Finally, we numerically solve nontrivial problems of size up to 2 50 × 2 50 . For the specific examples that we consider, we heuristically find that the time complexity of VQLS scales efficiently in ϵ, κ, and the system size N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Timestamp Error Detection and Estimation for PMU Data based on Linear Correlation between Relative Phase Angle and Frequency

Time synchronization is essential to synchro-phasor-based applications. However, Timestamp Error (TE) in synchrophasor data can result in application failures. This paper proposes a method for TE detection based on the linear correlation between frequency and relative phase angle. The TE converts the short-term relative phase angle from noise-like signal to one that linear with the frequency. Pearson Correlation Coefficient (PCC) is applied to measure the linear correlation and then detect the timestamp error. The time error is estimated based on the variation of frequency and relative phase angle. Case studies with actual synchrophasor data demonstrate the effectiveness of TE detection and excellent accuracy of TE estimation.

Yu, Wenpeng↗

Geometric Interpretation of a Non-Linear Extension of Quantum Mechanics

We recently introduced a particular non-linear generalization of quantum mechanics that has the property that it is exactly solvable in terms of the eigenvalues and eigenfunctions of the Hamiltonian of the usual linear quantum mechanics problem. In this paper, we suggest that the two components of the wave function represent the system described by the Hamiltonian H in two different asymptotic regions of spacetime and we show that the non-linear terms can be viewed as giving rise to gravitational effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

LINEAR POLARIZATION OF CLASS I METHANOL MASERS IN MASSIVE STAR-FORMING REGIONS

Class I methanol masers are found to be good tracers of the interaction between outflows from massive young stellar objects with their surrounding media. Although polarization observations of Class II methanol masers have been able to provide information about magnetic fields close to the central (proto)stars, polarization observations of Class I methanol masers are rare, especially at 44 and 95 GHz. We present the results of linear polarization observations of 39 Class I methanol maser sources at 44 and 95 GHz. These two lines are observed simultaneously with one of the 21 m Korean VLBI Network telescopes in single-dish mode. Approximately 60% of the observed sources have fractional polarizations of a few percent in at least one transition. This is the first reported detection of linear polarization of the 44 GHz methanol maser. The two maser transitions show similar polarization properties, indicating that they trace similar magnetic environments, although the fraction of the linear polarization is slightly higher at 95 GHz. We discuss the association between the directions of polarization angles and outflows. We also discuss some targets having different polarization properties at both lines, including DR21(OH) and G82.58+0.20, which show the 90° polarization angle flip at 44 GHz.

79 ASTRONOMY AND ASTROPHYSICS↗

Linear Polarization, Circular Polarization, and Depolarization of Gamma-ray Bursts: A Simple Case of Jitter Radiation

Linear and circular polarizations of gamma-ray bursts (GRBs) have been detected recently. We adopt a simplified model to investigate GRB polarization characteristics in this paper. A compressed two-dimensional turbulent slab containing stochastic magnetic fields is considered, and jitter radiation can produce the linear polarization under this special magnetic field topology. Turbulent Faraday rotation measure (RM) of this slab makes strong wavelength-dependent depolarization. The jitter photons can also scatter with those magnetic clumps inside the turbulent slab, and a nonzero variance of the Stokes parameter V can be generated. Furthermore, the linearly and circularly polarized photons in the optical and radio bands may suffer heavy absorptions from the slab. Thus we consider the polarized jitter radiation transfer processes. Finally, we compare our model results with the optical detections of GRB 091018, GRB 121024A, and GRB 131030A. We suggest simultaneous observations of GRB multi-wavelength polarization in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dense-gas Mass versus Star Formation Rate Relation: A Misleading Linearity?

We model the star formation relation of molecular clumps in dependence of their dense-gas mass when their volume density profile is that of an isothermal sphere (i.e., ρ{sub clump}(r)∝r{sup −2}). Dense gas is defined as gas whose volume density is higher than a threshold ρ{sub th}=700 M{sub ⊙} pc{sup −3} (i.e., HCN(1-0)-mapped gas). We divide the clump into two regions: a dense inner region (where ρ{sub clump}(r)⩾ρ{sub th}), and low-density outskirts (where ρ{sub clump}(r)<ρ{sub th}). We find that the total star formation rate of clumps scales linearly with the mass of their dense inner region, even when more than half of the clump star formation activity takes place in the low-density outskirts. We therefore emphasize that a linear star formation relation does not necessarily imply that star formation takes place exclusively in the gas whose mass is given by the star formation relation. The linearity of the star formation relation is strengthened when we account for the mass of dense fragments (e.g., cores, fibers) seeding star formation in the low-density outskirts, and which our adopted clump density profile ρ{sub clump}(r) does not resolve. We also find that the star formation relation is significantly tighter when considering the dense gas than when considering all the clump gas, as observed for molecular clouds of the Galactic plane. When the clumps have no low-density outskirts (i.e., they consist of dense gas only), the star formation relation becomes superlinear and progressively wider.

79 ASTRONOMY AND ASTROPHYSICS↗

Underlying one-step methods and nonautonomous stability of general linear methods

We generalize the theory of underlying one-step methods to strictly stable general linear methods (GLMs) solving nonautonomous ordinary differential equations (ODEs) that satisfy a global Lipschitz condition. Here, we combine this theory with the Lyapunov and Sacker-Sell spectral stability theory for one-step methods developed in [34,35,36] to analyze the stability of a strictly stable GLM solving a nonautonomous linear ODE. These results are applied to develop a stability diagnostic for the solution of nonautonomous linear ODEs by strictly stable GLMs.

97 MATHEMATICS AND COMPUTING↗

Frontal Slice Approaches for Tensor Linear Systems

Inspired by the row and column action methods for solving large-scale linear systems, in this work, we explore the use of frontal slices for solving tensor linear systems. In particular, this paper presents a novel approach for using frontal slices of a tensor $\mathcal{A}$ to solve tensor linear systems $\mathcal{A} ∗\mathcal{X} = \mathcal{B}$ where ∗ denotes the $t$-product. In addition, we consider variations of this method, including cyclic, block, and randomized approaches, each designed to optimize performance in different operational contexts. Our primary contribution lies in the development and convergence analysis of these methods. Experimental results on synthetically generated and real-world data, including applications such as image and video deblurring, demonstrate the efficacy of our proposed approaches and validate our theoretical findings.

Luo, Hengrui↗

ESPPU INPUT: C$^3$ within the "Linear Collider Vision"

The Linear Collider Vision calls for a Linear Collider Facility with a physics reach from a Higgs Factory to the TeV-scale with $e^+e^{-}$ collisions. One of the technologies under consideration for the accelerator is a cold-copper distributed-coupling linac capable of achieving high gradient. This technology is being pursued by the C$^3$ collaboration to understand its applicability to future colliders and broader scientific applications. In this input we share the baseline parameters for a C$^3$ Higgs-factory and the energy reach of up to 3 TeV in the 33 km tunnel foreseen under the Linear Collider Vision. Recent results, near-term plans and future R&D needs are highlighted.

43 PARTICLE ACCELERATORS↗