Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix 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 1,621 records · Page 90

A coupled-mode theory for multiwaveguide systems satisfying the reciprocity theorem and power conservation

Two sets of coupled-mode equations for multiwaveguide systems are derived using a generalized reciprocity relation; one set for a lossless system, and the other for a general lossy or lossless system. The second set of equations also reduces to those of the first set in the lossless case under the condition that the transverse field components are chosen to be real. Analytical relations between the coupling coefficients are shown and applied to the coupling of mode equations. It is shown analytically that these results satisfy exactly both the reciprocity theorem and power conservation. New orthogonal relations between the supermodes are derived in matrix form, with the overlap integrals taken into account.

Chuang, Shun-Lien↗

Unveiling the sea: universality of the transverse momentum dependent quark distributions at small x

Within the Colour Glass Condensate effective theory, we demonstrate that back-to-back dijet correlations in dilute-dense collisions involving a small-quark from the nuclear target can be factorised in terms of universal transverse momentum dependent distributions (TMDs) for the sea quarks. Two building blocks are needed to construct all these TMDs at the operator level: the sea quark TMD operator which appears in semi-inclusive Deep-Inelastic Scattering (SIDIS) or in the Drell-Yan process and the elastic-matrix for a quark-antiquark dipole. Compared to SIDIS, the saturation effects are stronger for dijet production in forward proton-nucleus collisions, due to additional scattering in the initial and final state, effectively resulting in a larger value for the nuclear saturation momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simultaneous stabilization and simultaneous pole placement by nonswitching dynamic compensation

The 'simultaneous stabilization problem' is defined and theorems are proposed for its solution. The problem consists in answering the question: given an r-tuple G sub 1(s), G sub r(s) of p x m proper transfer functions, does there exist a compensator K(s) such that the closed loop systems G sub 1(s) (I+K(s)G sub 1(s)) (-1), G sub r(s) (I+K(s) G sub r(s)) (-1) are (internally) stable. This question arises in reliability theory, where G sub 2(s), G sub r(s) represents a plant G sub 1(s) operating in various modes of failure and K(s) is a nonswitching stabilizing compensator. It is important in the stability analysis and design of a plant which can be switched into various operating modes. The simultaneous stabilization problem can also apply to the stabilization of a nonlinear system which is linearized at several equilibria. Conditions are defined for pole placement and the generalized Sylvestor matrix is discussed.

Ghosh, B. K.↗

Simultaneous stabilization and simultaneous pole-placement by nonswitching dynamic compensation

The 'simultaneous stabilization problem' is defined and theorems are proposed for its solution. The problem consists in answering the question: given an r-tuple G sub 1(s), G sub r(s) of p x m proper transfer functions, does there exist a compensator K(s) such that the closed loop systems G sub 1(s) (I+K(s)G sub 1(s)) (-1), G sub r(s) (I+K(s) G sub r(s)) (-1) are (internally) stable. This question arises in reliability theory, where G sub 2(s), G sub r(s) represents a plant G sub 1(s) operating in various modes of failure and K(s) is a nonswitching stabilizing compensator. It is important in the stability analysis and design of a plant which can be switched into various operating modes. The simultaneous stabilization problem can also apply to the stabilization of a nonlinear system which is linearized at several equilibria. Conditions are defined for pole placement and the generalized Sylvestor matrix is discussed. Previously announced in STAR as N82-31031

Ghosh, B. K.↗

An information theory of image gathering

Shannon's mathematical theory of communication is extended to image gathering. Expressions are obtained for the total information that is received with a single image-gathering channel and with parallel channels. It is concluded that the aliased signal components carry information even though these components interfere with the within-passband components in conventional image gathering and restoration, thereby degrading the fidelity and visual quality of the restored image. An examination of the expression for minimum mean-square-error, or Wiener-matrix, restoration from parallel image-gathering channels reveals a method for unscrambling the within-passband and aliased signal components to restore spatial frequencies beyond the sampling passband out to the spatial frequency response cutoff of the optical aperture.

Fales, Carl L.↗

The 𝑚=2 amplituhedron and the hypersimplex: Signs, clusters, tilings, Eulerian numbers

The hypersimplex Δ <#comment/> k + 1 , n \Delta _{k+1,n} is the image of the positive Grassmannian G r k + 1 , n ≥ <#comment/> 0 Gr^{\geq 0}_{k+1,n} under the moment map. It is a polytope of dimension n − <#comment/> 1 n-1 in R n \mathbb {R}^n . Meanwhile, the amplituhedron A n , k , 2 ( Z ) \mathcal {A}_{n,k,2}(Z) is the projection of the positive Grassmannian G r k , n ≥ <#comment/> 0 Gr^{\geq 0}_{k,n} into the Grassmannian G r k , k + 2 Gr_{k,k+2} under a map Z ~ <#comment/> \tilde {Z} induced by a positive matrix Z ∈ <#comment/> M a t n , k + 2 > 0 Z\in Mat_{n,k+2}^{>0} . Introduced in the context of scattering amplitudes , it is not a polytope, and has full dimension 2 k 2k inside G r k , k + 2 Gr_{k,k+2} . Nevertheless, there seem to be remarkable connections between these two objects via T-duality , as conjectured by Łukowski, Parisi, and Williams [Int. Math. Res. Not. (2023)]. In this paper we use ideas from oriented matroid theory, total positivity, and the geometry of the hypersimplex and positroid polytopes to obtain a deeper understanding of the amplituhedron. We show that the inequalities cutting out positroid polytopes —images of positroid cells of G r k + 1 , n ≥ <#comment/> 0 Gr^{\geq 0}_{k+1,n} under the moment map—translate into sign conditions characterizing the T-dual Grasstopes —images of positroid cells of G r k , n ≥ <#comment/> 0 Gr^{\geq 0}_{k,n} under Z ~ <#comment/> \tilde {Z} . Moreover, we subdivide the amplituhedron into chambers , just as the hypersimplex can be subdivided into simplices, with both chambers and simplices enumerated by the Eulerian numbers. We use these properties to prove the main conjecture of Łukowski, Parisi, and Williams [Int. Math. Res. Not. (2023)]: a collection of positroid polytopes is a tiling of the hypersimplex if and only if the collection of T-dual Grasstopes is a tiling of the amplituhedron A n , k , 2 ( Z ) \mathcal {A}_{n,k,2}(Z) for all Z Z . Moreover, we prove Arkani-Hamed–Thomas–Trnka’s conjectural sign-flip characterization of A n , k , 2 \mathcal {A}_{n,k,2} , and Łukowski–Parisi–Spradlin–Volovich’s conjectures on m = 2 m=2 cluster adjacency and on positroid tiles for A n , k , 2 \mathcal {A}_{n,k,2} (images of 2 k 2k -dimensional positroid cells which map injectively into A n , k , 2 \mathcal {A}_{n,k,2} ). Finally, we introduce new cluster structures in the amplituhedron.

Parisi, Matteo↗

Considering nonlocality in the optical potentials within eikonal models

Background: For its simplicity, the eikonal method is the tool of choice to analyze nuclear reactions at high energies (E > 100 MeV/nucleon), including knockout reactions. However, so far, the effective interactions used in this method are assumed to be fully local. Purpose: Given the recent studies on nonlocal optical potentials, in this work we assess whether nonlocality in the optical potentials is expected to impact reactions at high energies and then explore different avenues for extending the eikonal method to include nonlocal interactions. Method: We compare angular distributions obtained for nonlocal interactions (using the exact R-matrix approach for elastic scattering and the adiabatic distorted wave approximation for transfer) with those obtained using their local-equivalent interactions. Results: Our results show that transfer observables are significantly impacted by nonlocality in the high-energy regime. Because knockout reactions are dominated by stripping (transfer to inelastic channels), nonlocality is expected to have a large effect on knockout observables too. Three approaches are explored for extending the eikonal method to nonlocal interactions, including an iterative method and a perturbation theory. Conclusions: None of the derived extensions of the eikonal model provide a good description of elastic scattering. Here, this paper suggests that nonlocality removes the formal simplicity associated with the eikonal model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling fission dynamics at the barrier in a discrete-basis formalism

A configuration-interaction model is presented for the barrier region of induced fission. The configuration space is composed of seniority-zero configurations constructed from self-consistent mean-field wave functions. The Hamiltonian matrix elements between configurations include diabatic and pairing interactions between particles. Other aspects of the Hamiltonian are treated statistically, guided by phenomenological input of compound-nucleus transmission coefficients. In this exploratory study the configuration space is restricted to neutron excitations only. A key observable calculated in the model is the fission-to-capture branching ratio. We find that both pairing and diabatic interactions are important for achieving large branching to the fission channels. In accordance with the transition-state theory of fission, the calculated branching ratio is found to be quite insensitive to the fission decay widths of the pre-scission configurations. Furthermore, the barrier-top dynamics appear to be quite different from transition-state theory in that the transport is distributed over many excited configurations at the barrier top.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simple robust control laws for robot manipulators. Part 2: Adaptive case

A new class of asymptotically stable adaptive control laws is introduced for application to the robotic manipulator. Unlike most applications of adaptive control theory to robotic manipulators, this analysis addresses the nonlinear dynamics directly without approximation, linearization, or ad hoc assumptions, and utilizes a parameterization based on physical (time-invariant) quantities. This approach is made possible by using energy-like Lyapunov functions which retain the nonlinear character and structure of the dynamics, rather than simple quadratic forms which are ubiquitous to the adaptive control literature, and which have bound the theory tightly to linear systems with unknown parameters. It is a unique feature of these results that the adaptive forms arise by straightforward certainty equivalence adaptation of their nonadaptive counterparts found in the companion to this paper (i.e., by replacing unknown quantities by their estimates) and that this simple approach leads to asymptotically stable closed-loop adaptive systems. Furthermore, it is emphasized that this approach does not require convergence of the parameter estimates (i.e., via persistent excitation), invertibility of the mass matrix estimate, or measurement of the joint accelerations.

Bayard, D. S.↗

Cross Sections for Electron Impact Excitation of Ions Relevant to Planetary Atmospheres Observation

The goal of this research grant was to calculate accurate oscillator strengths and electron collisional excitation strengths for inelastic transitions in atomic species of relevance to Planetary Atmospheres. Large scale configuration-interaction atomic structure calculations have been performed to obtain oscillator strengths and transition probabilities for transitions among the fine-structure levels and R-matrix method has been used in the calculations of electron-ion collision cross sections of C II, S I, S II, S III, and Ar II. A number of strong features due to ions of sulfur have been detected in the spectra of Jupiter satellite Io. The electron excitation cross sections for the C II and S II transitions are studied in collaboration with the experimental atomic physics group at the Jet Propulsion Laboratory. There is excellent agreement between experiment and theory which provide an accurate and broad-base test of the ability of theoretical methods used in the calculation of atomic processes. Specifically, research problems have been investigated for: electron impact excitation cross sections of C II: electron impact excitation cross sections of S III; energy levels and oscillator strengths for transitions in S III; collision strengths for electron collisional excitation of S II; electron impact excitation of inelastic transitions in Ar II; oscillator strengths of fine-structure transitions in neutral sulfur; cross sections for inelastic scattering of electrons from atomic nitrogen; and excitation of atomic ions by electron impact.

Tayal, Swaraj S.↗

Baryon acoustic oscillations in thin redshift shells from BOSS DR12 and eBOSS DR16 galaxies

ABSTRACT In an age of large astronomical data sets and severe cosmological tensions, the case for model independent analyses is compelling. We present a set of 14 baryon acoustic oscillations measurements in thin redshift shells with $3\,\mathrm{ per} \,\mathrm{ cent}$ precision that were obtained by analysing BOSS DR12 and eBOSS DR16 galaxies in the redshift range 0.32 < z < 0.66. Thanks to the use of thin shells, the analysis is carried out using just redshifts and angles so that the fiducial model is only introduced when considering the mock catalogues, necessary for the covariance matrix estimation and the pipeline validation. We compare our measurements, with and without supernova data, to the corresponding constraints from Planck 2018, finding good compatibility. A Monte Python module for this likelihood is available at github.com/ranier137/angularBAO.

79 ASTRONOMY AND ASTROPHYSICS↗

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM’s Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally invariant system on a quantum computer is presented, which we call “scalable circuits ADAPT-VQE” (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined with use of classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and is shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wave functions are found to become independent of the number of spatial sites, L . This allows a controlled extrapolation of the circuits, determined with use of small or modest-sized systems, to arbitrarily large L . The circuits for the Schwinger model are determined on lattices up to L = 14 (28 qubits) with the Qiskit classical simulator, and are subsequently scaled up to prepare the L = 50 (100 qubits) vacuum on IBM’s 127-superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introduction of an improved error-mitigation technique, which we call “operator decoherence renormalization”, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical matrix product state simulations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Applications of velocity potential function to acoustic duct propagation and radiation from inlets using finite element theory

A finite element velocity potential program was developed to study acoustic wave propagation in complex geometries. For irrotational flows, relatively low sound frequencies, and plane wave input, the finite element solutions showed significant effects of inlet curvature and flow gradients on the attenuation of a given acoustic liner in a realistic variable area turbofan inlet. The velocity potential approach can not be used to estimate the effects of rotational flow on acoustic propagation, since the potential acoustic disturbances propagate at the speed of the media in sheared flow. Approaches are discussed that are being considered for extending the finite element solution to include the far field, as well as the internal portion of the duct. A new matrix partitioning approach is presented that can be incorporated in previously developed programs to allow the finite element calculation to be marched into the far field. The partitioning approach provided a large reduction in computer storage and running times.

Baumeister, K. J.↗

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high k modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental verification and practical application of torquewhirl theory of rotordynamic instability

A theory developed by Vance in 1978 to explain the destabilizing effect of torque on a whirling rotor was experimentally verified. The measurements made on a specially designed test apparatus are described. New computer models were also developed to investigate the effect of torque on rotordynamic stability of multidisk flexible rotor bearing systems. The effect of torque was found to be most pronounced when the system is already marginally stable from other influences. The modifications required to include torque in a typical shaft transfer matrix are described, and results are shown which identify the type of rotor design most sensitive to load torque.

Vance, J. M.↗

Off axis loading of silicon carbide fibers in a titanium aluminide matrix

Silicon carbide SCS-6 fibers in a Ti24A1 + 11Nb matrix were subjected to off axis loading in a 'thin-slice' pushout test, resulting in various combinations of shear, radial compression, and tension along the fibers as a function of orientation angle. The load necessary for debonding decreased as the orientation angle increased, whereas the average frictional sliding stress after 60 s of sliding remained relatively constant for orientation angles less than 30 deg. Analyses of the specimen bending stresses and of the contact stresses by finite element modeling and thin plate theory are presented.

Kennefick, C. M.↗

Compensator improvement for multivariable control systems

A theory and the associated numerical technique are developed for an iterative design improvement of the compensation for linear, time-invariant control systems with multiple inputs and multiple outputs. A strict constraint algorithm is used in obtaining a solution of the specified constraints of the control design. The result of the research effort is the multiple input, multiple output Compensator Improvement Program (CIP). The objective of the Compensator Improvement Program is to modify in an iterative manner the free parameters of the dynamic compensation matrix so that the system satisfies frequency domain specifications. In this exposition, the underlying principles of the multivariable CIP algorithm are presented and the practical utility of the program is illustrated with space vehicle related examples.

Mitchell, J. R.↗