Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “symmetric matrices”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Experimental investigation of the hydrodynamic forces on the shroud of a centrifugal pump impeller

Fluid-induced forces acting on a rotating impeller are known to cause rotor-dynamic problems in turbomachines. The forces generated by leakage flow along the front shroud surface of a centrifugal turbomachine impeller play an important role among these fluid-induced forces. The present research was aimed to gain a better understanding of these shroud forces. An experimental apparatus was designed and constructed to simulate the impeller shroud leakage flow. Hydrodynamic forces and steady and unsteady pressure distributions on the rotating shroud were measured as functions of eccentricity, width of shroud clearance, face seal clearance and shaft rotating speed. The forces measured from the dynamometer and manometers agreed well. The hydrodynamic force matrices were found skew-symmetric and statically unstable. This is qualitatively similar to the result of previous hydrodynamic volute force measurements. Nondimensionalized normal and tangential forces decrease slightly as Reynolds number increases. As the width of the shroud clearance decreases and/or the eccentricity increases, the hydrodynamic forces increase nonlinearly. There was some evidence found that increased front seal clearance could reduce the radial shroud forces and the relative magnitude of the destabilizing tangential force. Subharmonic pressure fluctuations were also observed which may adversely affect the behavior of the rotor system.

Zhuang, Fei↗
DOE OSTI2024

Robust Iterative Method for Symmetric Quantum Signal Processing in All Parameter Regimes

Here, this paper addresses the problem of solving nonlinear systems in the context of symmetric quantum signal processing (QSP), a powerful technique for implementing matrix functions on quantum computers. Symmetric QSP focuses on representing target polynomials as products of matrices in SU(2) that possess symmetry properties. We present a novel Newton’s method tailored for efficiently solving the nonlinear system involved in determining the phase factors within the symmetric QSP framework. Our method demonstrates rapid and robust convergence in all parameter regimes, including the challenging scenario with ill-conditioned Jacobian matrices, using standard double precision arithmetic operations. For instance, solving symmetric QSP for a highly oscillatory target function α cos(1000x) (polynomial degree ≈ 1433) takes 6 iterations to converge to machine precision when α = 0.9, and the number of iterations only increases to 18 iterations when α = 1 – 10 -9 with a highly ill-conditioned Jacobian matrix. Leveraging the matrix product state structure of symmetric QSP, the computation of the Jacobian matrix incurs a computational cost comparable to a single function evaluation. Moreover, we introduce a reformulation of symmetric QSP using real-number arithmetics, further enhancing the method’s efficiency. Extensive numerical tests validate the effectiveness and robustness of our approach, which has been implemented in the QSPPACK software package.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Discrete Kalman filtering equations of second-order form for control-structure interaction simulations

A second-order form of discrete Kalman filtering equations is proposed as a candidate state estimator for efficient simulations of control-structure interactions in coupled physical coordinate configurations as opposed to decoupled modal coordinates. The resulting matrix equation of the present state estimator consists of the same symmetric, sparse N x N coupled matrices of the governing structural dynamics equations as opposed to unsymmetric 2N x 2N state space-based estimators. Thus, in addition to substantial computational efficiency improvement, the present estimator can be applied to control-structure design optimization for which the physical coordinates associated with the mass, damping and stiffness matrices of the structure are needed instead of modal coordinates.

Park, K. C.↗
NASA NTRS1997

Graph Embedding Techniques for Bounding Condition Numbers of Incomplete Factor Preconditioning

We extend graph embedding techniques for bounding the spectral condition number of preconditioned systems involving symmetric, irreducibly diagonally dominant M-matrices to systems where the preconditioner is not diagonally dominant. In particular, this allows us to bound the spectral condition number when the preconditioner is based on an incomplete factorization. We provide a review of previous techniques, describe our extension, and give examples both of a bound for a model problem, and of ways in which our techniques give intuitive way of looking at incomplete factor preconditioners.

Guattery, Stephen↗
DOE OSTI2024

Shallow Bound States and Hints for Broad Resonances with Quark Content b ¯ c ¯ u d in B − D ¯ and B * − D ¯ Scattering from Lattice QCD

We present the first determination of the energy dependence of the B − D ¯ and B * − D ¯ isospin-0, S -wave scattering amplitudes both below and above the thresholds using lattice QCD, which allows us to investigate rigorously whether mixed bottom-charm b ¯ c ¯ u d tetraquarks exist as bound states or resonances. The scattering phase shifts are obtained using Lüscher’s method from the energy spectra in two different volumes. To ensure that no relevant energy level is missed, we use large, symmetric 7 × 7 and 8 × 8 correlation matrices that include, at both source and sink, B ( * ) − D ¯ scattering operators with the lowest three or four possible back-to-back momenta in addition to local b ¯ c ¯ u d operators. We fit the energy dependence of the extracted scattering phase shifts using effective-range expansions. We observe sharp peaks in the B ( * ) − D ¯ scattering rates close to the thresholds, which are associated with shallow bound states, either genuine or virtual, a few MeV or less below the B ( * ) − D ¯ thresholds. In addition, we find hints for resonances with masses of order 100 MeV above the thresholds and decay widths of order 200 MeV. Published by the American Physical Society 2024

Alexandrou, Constantia↗
NASA NTRS1984

System identification of analytical models of damped structures

A procedure is presented for identifying linear nonproportionally damped system. The system damping is assumed to be representable by a real symmetric matrix. Analytical mass, stiffness and damping matrices which constitute an approximate representation of the system are assumed to be available. Given also are an incomplete set of measured natural frequencies, damping ratios and complex mode shapes of the structure, normally obtained from test data. A method is developed to find the smallest changes in the analytical model so that the improved model can exactly predict the measured modal parameters. The present method uses the orthogonality relationship to improve mass and damping matrices and the dynamic equation to find the improved stiffness matrix.

Fuh, J.-S.↗
NASA NTRS1986

A generalized theory for eccentric and misalignment effects in high-pressure annular seals

High-pressure annular seal leakage and dynamic coefficients vary with eccentricity and misalignment. Recent seal leakage data with both concentric and fully eccentric alignments support the seal leakage model with surface roughness and eccentricity effects included. In this paper, the seal dynamic coefficient calculation has been generalized and allows direct calculation of the seal dynamic coefficients at any circumferential location. The generalized solution agrees with the results obtained by using the calculated values of an earlier paper and performing a coordinate transformation. The analysis results coincide with the measured data in showing that the stiffness and damping matrices of seal coefficients are not skew symmetric, and the main diagonal seal coefficients are not equal. The measured direct stiffnesses were found higher than predicted by the concentric seal theory, but this may be explained by the presence of eccentricity in the test operating mode.

Chen, W. C.↗
NASA NTRS1989

Factorization Of Positive Definite, Banded Hermitian Matrices

Report discusses application of Cholesky factorization algorithm to positive definite, banded Hermitian matrices. Begins by extending Cholesky factorization algorithm to cover uniformly-partitioned, banded, positive definite matrices of rank n that is real symmetric or Hermitian. Then two stratagems given for use of algorithm in concurrent-processing system in which N less than it has to be to enable factorization of matrix in as few serial steps as possible and where uniformly high efficiency expected from all processing elements. One of major purposes of this and related studies to maximize speedup and efficiency in system of concurrent-data-processing elements.

Salama, Moktar A.↗
NASA NTRS1989

Ritz method for transient response in systems having unsymmetric stiffness

The DMAP coding was automated to such an extent by using the device of bubble vectors, that it is useable for analyses in its present form. This feasibility study demonstrates that the Ritz Method is so compelling as to warrant coding its modules in FORTRAN and organizing the resulting coding into a new Rigid Format. Even though this Ritz technique was developed for unsymmetric stiffness matrices, it offers advantages to problems with symmetric stiffnesses. If used for the symmetric case the solution would be simplified to one set of modes, because the adjoint would be the same as the primary. Its advantage in either type of symmetry over a classical eigenvalue modal expansion is that information density per Ritz mode is far richer than per eigenvalue mode; thus far fewer modes would be needed for the same accuracy and every mode would actively participate in the response. Considerable economy can be realized in adapting Ritz vectors for modal solutions. This new Ritz capability now makes NASTRAN even more powerful than before.

Butler, Thomas G.↗
NASA NTRS1986

Concurrent Cholesky factorization of positive definite banded Hermitian matrices

First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.

Utku, S.↗
NASA NTRS1978

Extension of the tridiagonal reduction (FEER) method for complex eigenvalue problems in NASTRAN

As in the case of real eigenvalue analysis, the eigensolutions closest to a selected point in the eigenspectrum were extracted from a reduced, symmetric, tridiagonal eigenmatrix whose order was much lower than that of the full size problem. The reduction process was effected automatically, and thus avoided the arbitrary lumping of masses and other physical quantities at selected grid points. The statement of the algebraic eigenvalue problem admitted mass, damping, and stiffness matrices which were unrestricted in character, i.e., they might be real, symmetric or nonsymmetric, singular or nonsingular.

Newman, M.↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Hydrodynamic impeller stiffness, damping, and inertia in the rotordynamics of centrifugal flow pumps

The lateral hydrodynamic forces experienced by a centrifugal pump impeller performing circular whirl motions within several volute geometries were measured. The lateral forces were decomposed into: (1) time averaged lateral forces and (2) hydrodynamic force matrices representing the variation of the lateral forces with position of the impeller center. It is found that these force matrices essentially consist of equal diagonal terms and skew symmetric off diagonal terms. One consequence of this is that during its whirl motion the impeller experiences forces acting normal and tangential to the locus of whirl. Data on these normal and tangential forces are presented; it is shown that there exists a region of positive reduced whirl frequencies, within which the hydrodynamic forces can be destablizing with respect to whirl.

Jery, S.↗
DOE OSTI2026

Data Driven Correlated Noise Simulation for the ICEBERG LArTPC

Accurate electronic-noise simulation is essential for low-energy physics in liquid-argon TPCs. More realistic noise modeling allows us to better tune reconstruction algorithms and more reliably assess and optimize signal-detection thresholds. We present a data-driven noise simulation framework developed for the ICEBERG test stand for DUNE that generates synthetic noise waveforms that reproduce both (i) the measured per-channel magnitude of the Fast Fourier Transform (FFT) and (ii) frequency-dependent channel-to-channel correlations observed in ICEBERG noise data. Using a dedicated noise-only dataset, we build a compact noise model containing per-channel FFT-magnitude targets together with a small set of band-wise cross-wire color matrices. White noise is generated in the frequency domain by drawing circular-symmetric complex Gaussian coefficients with random phases and scaling them to match the measured FFT-magnitude targets, and cross-wire correlations are subsequently imposed using the stored color matrices. The model and algorithm were integrated into the LArSoft + Wire-Cell Toolkit simulation chain and validated by comparing waveform structure, frequency-domain spectra, and band-limited correlation matrices from simulated noise and ICEBERG data. This approach can be extended to other LArTPC operating conditions.

Ghosh, Avik [Iowa State U.]↗
NASA NTRS1983

Vibration control of large linear quadratic symmetric systems

Some unique properties on a class of the second order lambda matrices were found and applied to determine a damping matrix of the decoupled subsystem in such a way that the damped system would have preassigned eigenvalues without disturbing the stiffness matrix. The resulting system was realized as a time invariant velocity only feedback control system with desired poles. Another approach using optimal control theory was also applied to the decoupled system in such a way that the mode spillover problem could be eliminated. The procedures were tested successfully by numerical examples.

Jeon, G. J.↗
NASA NTRS1977

Complex eigenvalue extraction in NASTRAN by the tridiagonal reduction (FEER) method

An extension of the Tridiagonal Reduction (FEER) method to complex eigenvalue analysis in NASTRAN is described. As in the case of real eigenvalue analysis, the eigensolutions closest to a selected point in the eigenspectrum are extracted from a reduced, symmetric, tridiagonal eigenmatrix whose order is much lower than that of the full size problem. The reduction process is effected automatically, and thus avoids the arbitrary lumping of masses and other physical quantities at selected grid points. The statement of the algebraic eigenvalue problem admits mass, damping and stiffness matrices which are unrestricted in character, i.e., they may be real, complex, symmetric or unsymmetric, singular or non-singular.

Newman, M.↗
DOE OSTI2022

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1975

An interactive graphics package for the automatic node renumbering of finite element matrices

An interactive graphics software package which allows users to display the non-zero structure of large sparse symmetric materials was described and methods used to implement it as a portable FORTRAN callable subroutine were summarized. In particular, the system permits the display of the resulting matrix after reordering the rows and columns, with the reordering scheme either defined by the user or automatically generated by the program with the aim of reducing matrix bandwidth and profile. Although the primary application of the package has been to the finite element analysis of structures, it is equally well suited to the many other areas of engineering and science which use sparse matrices.

Boisvert, R. F.↗