Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗
DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2016

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗
NASA NTRS1980

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗
NASA NTRS1980

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one-dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗
NASA NTRS2006

History of Metal Veins in Acapulcoite-Lodranite Clan Meteorite GRA 95209

Graves Nunataks (GRA) 95209 has been hailed as the missing link of core formation processes in the acapulcoitelodranite parent asteroid because of the presence of a complex cm-scale metal vein network. Because the apparent liquid temperature of the metal vein (approximately 1500 C) is higher than inferred for the metamorphic grade of the meteorite, questions regarding the vein s original composition, temperature, and mechanism of emplacement have arisen. We have determined trace siderophile element compositions of metals in veins and surrounding matrix in an effort to clarify matters. We analyzed metals in GRA 95209 in a portion of thick metal vein and adjacent metal-rich (30-40 modal%), sulfide poor (less than 1%) matrix by EPMA and LA-ICP-MS for major and trace siderophile elements using methods described by [3]. We also examined metals from a metal-poor (approximately 15 modal%) and relatively sulfide-rich (2-5 modal%) region of the sample. Kamacite is the dominant metal phase in all portions of the sample. In comparison to matrix metal, vein metal contains more schreibersite and less tetrataenite, and is less commonly associated with Fe,Mn,Mg-bearing phosphates and graphite. Vein kamacite contains higher Co, P, and Cr and lower Cu and Ge. These minor variations aside, all metal types in GRA 95209 are fairly homogeneous in terms of their levels of enrichment of compatible siderophile elements (e.g. Pt, Ir, Os) relative to incompatible siderophile elements (e.g. As, Pd, Au), consistent with the loss of metal-sulfide partial melt that characterizes much of the clan. Whatever compositional differences between matrix and vein metal that may have originally existed, they have since largely co-equilibrated to similar restitic trace element compositions. We agree with [2] that metal veins, in their present state, do not represent a liquid composition. The original vein liquid was much more S-rich and emplaced at correspondingly lower liquid temperatures. Much of the Fe,Ni component solidified in cm scale conduits while S-rich melts were expelled and continued to migrate by percolation. The higher troilite content in metal poor regions of the sample results mostly from trapping of a small portion of these melts. The troilite is not remnant primary sulfide. Strong depletions of W, Mo, and especially Ga (greater than 50%, greater than 60%, and greater than 90% depletion, respectively) in metals of the metalpoor GRA 95209 lithology are localized at scales of 10-100 micrometers in the vicinity of graphite spherules. These depletions must have occurred below the temperatures at which cm-scale equilibration occurred, and future work will seek to determine their cause.

Herrin, J. S.↗
NASA NTRS2007

Method of making carbon fiber-carbon matrix reinforced ceramic composites

A method of making a carbon fiber-carbon matrix reinforced ceramic composite wherein the result is a carbon fiber-carbon matrix reinforcement is embedded within a ceramic matrix. The ceramic matrix does not penetrate into the carbon fiber-carbon matrix reinforcement to any significant degree. The carbide matrix is a formed in situ solid carbide of at least one metal having a melting point above about 1850 degrees centigrade. At least when the composite is intended to operate between approximately 1500 and 2000 degrees centigrade for extended periods of time the solid carbide with the embedded reinforcement is formed first by reaction infiltration. Molten silicon is then diffused into the carbide. The molten silicon diffuses preferentially into the carbide matrix but not to any significant degree into the carbon-carbon reinforcement. Where the composite is intended to operate between approximately 2000 and 2700 degrees centigrade for extended periods of time such diffusion of molten silicon into the carbide is optional and generally preferred, but not essential.

Williams, Brian↗
NASA NTRS1974

A recurrence matrix method for the analysis of longitudinal and torsional vibrations in non-uniform multibranch beams with variable boundary conditions

An approximate method for calculating the longitudinal and torsional natural frequencies and associated modal data of a beamlike, variable cross section multibranch structure is presented. The procedure described is the numerical integration of the first order differential equations that characterize the beam element in longitudinal motion and that satisfy the appropriate boundary conditions.

Davis, R. B.↗
DOE OSTI2021

Color halo scenario of charmonium-like hybrids

The internal structures of J PC =1 -- ,(o,1,2) -+ charmonium-like hybrids are investigated under lattice QCD in the quenched approximation. We define the Bethe-Salpeter wave function (Φn(r)) in the Coulomb gauge as the matrix element of a spatially extended hybrid-like operator ($\bar{c}cg$) between the vacuum and n-th state for J PC , with r being the spatial separation between a localized $\bar{c}c$ component and the chromomagnetic strength tensor. These wave functions exhibit some similarities for states with the aforementioned different quantum numbers, and their r-behaviors (no node for the ground states and one node for the first excited states) imply that r can be a meaningful dynamical variable for these states. Additionally, the mass splittings of the ground states and first excited states of charmonium-like hybrids in these channels are obtained for the first time to be approximately 1.2-1.4 GeV. These results do not support the flux-tube description of heavy-quarkonium-like hybrids in the Born-Oppenheimer approximation. In contrast, a charmonium-like hybrid can be viewed as a "color halo" charmonium for which a relatively localized color octet $\bar{c}c$ is surrounded by gluonic degrees of freedom, which can readily decay into a charmonium state along with one or more light hadrons. The color halo picture is compatible with the decay properties of Y (42600) and suggests LHCb and BelleII to search for (0,1,2) -+ charmonium-like hybrids in Xc0,1,2N and J/ψω(φ) final states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Subleading corrections to the S 3 free energy of necklace quiver theories dual to massive IIA

We investigate the S 3 free energy of N = 3 Chern-Simons-matter quiver gauge theories with gauge group U(N) r (r ≥ 2) where the sum of Chern-Simons levels does not vanish, beyond the leading order in the large-N limit. We take two different approaches to explore the sub-leading structures of the free energy. First we evaluate the matrix integral for the partition function in the ’t Hooft limit using a saddle point approximation. Second we use an ideal Fermi-gas model to compute the same partition function, but in the limit of fixed Chern-Simons levels. The resulting expressions for the free energy F = – log Z are then compared in the overlapping parameter regime. The Fermi-gas approach also hints at a universal 1/6 log N correction to the free energy. Since the quiver gauge theories we consider are dual to massive Type IIA theory, we expect the sub-leading correction of the planar free energy in the large ’t Hooft parameter limit to match higher-derivative corrections to the tree-level holographic dual free energy, which have not yet been fully investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Singular value decomposition and similarity renormalization group evolution of nuclear interactions

One of the main challenges for ab initio nuclear many-body theory is the growth of computational and storage costs as calculations are extended to heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform Singular Value Decompositions of nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization scheme and resolution scales. We develop and implement the Similarity Renormalization Group (SRG) evolution of the factorized interaction, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves, and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Lastly, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the In-Medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Randomized algorithms for accelerating linear algebraic computations

The project supported the development of new methodologies for performing matrix computations that form key building blocks in modern scientific computing, such as low rank approximation of matrices, and efficient representations of global operators that arise in simulations of physical phenomena.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Self-consistent Pauli corrections in Brueckner-Hartree-Fock calculations

A scheme is introduced which makes it feasible to make completely self-consistent Brueckner-Hartree-Fock (BHF) and renormalized BHF calculations for spherical, closed-shell and axially-symmetric deformed nuclei. The usual requirement or orbital self-consistency has been imposed, as well as self-consistency in the starting energies and occupation probabilities. Previously, only approximate forms were used for the Pauli operator. This approximation is removed and a method for making the necessary Pauli corrections to the reaction matrix during the approach to self-consistency is presented. A discussion of the symmetries which reduce the problem to one of manageable proportions is included.

Braley, R. C.↗
NASA NTRS1978

Thermal environment effects on strength and impact properties of boron-aluminum composites

A systematic study was conducted regarding the degradation of fracture strength and impact energy in commercial B-Al composites in both static and cyclic thermal environments. The composites used in the study contained approximately 50 vol % boron fibers, unidirectionally aligned in either a 6061 Al or 1100 Al matrix. The tensile strengths of the composites after 3000 thermal cycles as a function of upper cycle temperature are presented in graphs. The temperature at which the strengths of 6061 Al matrix, B-Al composites were significantly degraded after 3000 cycles was noticeably higher than that for the 1100 Al matrix composites. Static heating at 420 C resulted in no significant strength degradation for the 6061 Al matrix composites. In the case of 1100 matrix composites, some degradation was observed at 420 C but markedly less than in the composites cycled to 420 C.

Grimes, H. H.↗
NASA NTRS1982

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗
NASA NTRS1996

Forced response testing of an axi-centrifugal turboshaft engine

Dynamic data from tests of a T55-L-712 engine are presented. Engine stall/surge data were analyzed using digital signal processing techniques. In addition, forced response testing (system identification studies) was done at various engine speeds. Forced response testing was done using eight jet ejectors approximately equally circumferentially spaced about the compressor front face. This paper presents some preliminary results for the ground idle (approximately 60% of design speed) point. Brief descriptions of the jet injection system, the test matrix, and analysis techniques used are presented. Results of these analyses indicate a substantial transfer of energy across the compressor first stage at some frequencies and that the ejectors are effective in modifying the local flow conditions in front of the first compressor stage.

Owen, A. Karl↗
NASA NTRS1996

Forced Response Testing of an Axi-Centrifugal Turboshaft Engine

Dynamic data from tests of a T55-L,712 engine are presented. Engine stall/surge data were analyzed using digital signal processing techniques. In addition, forced response testing (system identification studies) was done at various engine speeds. Forced response testing was done using eight jet ejectors approximately equally circumferentially spaced about the compressor front face. This paper presents some preliminary results for the ground idle (approximately 60% of design speed) point. Brief descriptions of the jet injection system, the test matrix, and analysis techniques used are presented. Results of these analyses indicate a substantial transfer of energy across the compressor first stage at some frequencies and that the ejectors are effective in modifying the local flow conditions in front of the first compressor stage.

Owen, A. Karl↗