Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix equations”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Creep rupture of polymer-matrix composites

The time-dependent creep-rupture process in graphite-epoxy laminates is examined as a function of temperature and stress level. Moisture effects are not considered. An accelerated characterization method of composite-laminate viscoelastic modulus and strength properties is reviewed. It is shown that lamina-modulus master curves can be obtained using a minimum of normally performed quality-control-type testing. Lamina-strength master curves, obtained by assuming a constant-strain-failure criterion, are presented along with experimental data, and reasonably good agreement is shown to exist between the two. Various phenomenological delayed failure models are reviewed and two (the modified rate equation and the Larson-Miller parameter method) are compared to creep-rupture data with poor results.

Brinson, H. F.↗
NASA NTRS1984

Constitutive modeling for isotropic materials

A state-of-the-art review of applicable constitutive models with selection of two for detailed comparison with a wide range of experimental tests was conducted. The experimental matrix contained uniaxial and biaxial tensile, creep, stress relaxation, and cyclic fatigue tests at temperatures to 1093 C and strain rates from .0000001 to .001/sec. Some nonisothermal cycles will also be run. The constitutive models will be incorporated into the MARC finite element structural analysis program with a demonstration computation made for advanced turbine blade configuration. In the code development work, particular emphasis is being placed on developing efficient integration algorithms for the highly nonlinear and stiff constitutive equations. Another area of emphasis is the appropriate and efficient methodology for determing constitutive constants from a minimum extent of experimental data.

Lindholm, U. S.↗
NASA NTRS1993

Preconditioned conjugate-gradient methods for low-speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations is integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the Lower-Upper Successive Symmetric Over-Relaxation iterative scheme is more efficient than a preconditioner based on Incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional Line Gauss-Seidel Relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS2011

Effect of Data Reduction and Fiber-Bridging on Mode I Delamination Characterization of Unidirectional Composites

Reliable delamination characterization data for laminated composites are needed for input in analytical models of structures to predict delamination onset and growth. The double-cantilevered beam (DCB) specimen is used to measure fracture toughness, GIc, and strain energy release rate, GImax, for delamination onset and growth in laminated composites under mode I loading. The current study was conducted as part of an ASTM Round Robin activity to evaluate a proposed testing standard for Mode I fatigue delamination propagation. Static and fatigue tests were conducted on specimens of IM7/977-3 and G40-800/5276-1 graphite/epoxies, and S2/5216 glass/epoxy DCB specimens to evaluate the draft standard "Standard Test Method for Mode I Fatigue Delamination Propagation of Unidirectional Fiber-Reinforced Polymer Matrix Composites." Static results were used to generate a delamination resistance curve, GIR, for each material, which was used to determine the effects of fiber-bridging on the delamination growth data. All three materials were tested in fatigue at a cyclic GImax level equal to 90% of the fracture toughness, GIc, to determine the delamination growth rate. Two different data reduction methods, a 2-point and a 7-point fit, were used and the resulting Paris Law equations were compared. Growth rate results were normalized by the delamination resistance curve for each material and compared to the nonnormalized results. Paris Law exponents were found to decrease by 5.4% to 46.2% due to normalizing the growth data. Additional specimens of the IM7/977-3 material were tested at 3 lower cyclic GImax levels to compare the effect of loading level on delamination growth rates. The IM7/977-3 tests were also used to determine the delamination threshold curve for that material. The results show that tests at a range of loading levels are necessary to describe the complete delamination behavior of this material.

Murri, Gretchen B.↗
NASA NTRS2017

Electromagnetic Scattering by Spheroidal Volumes of Discrete Random Medium

We use the superposition T-matrix method to compare the far-field scattering matrices generated by spheroidal and spherical volumes of discrete random medium having the same volume and populated by identical spherical particles. Our results fully confirm the robustness of the previously identified coherent and diffuse scattering regimes and associated optical phenomena exhibited by spherical particulate volumes and support their explanation in terms of the interference phenomenon coupled with the order-of-scattering expansion of the far-field Foldy equations. We also show that increasing non-sphericity of particulate volumes causes discernible (albeit less pronounced) optical effects in forward and backscattering directions and explain them in terms of the same interference/multiple-scattering phenomenon.

Electromagnetic scattering; Multi-particle groups;↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1996

Quasi-Optimal Schwarz Methods for the Conforming Spectral Element Discretization

Fast methods are proposed for solving the system K(sub N)x = b resulting from the discretization of self-adjoint elliptic equations in three dimensional domains by the spectral element method. The domain is decomposed into hexahedral elements, and in each of these elements the discretization space is formed by polynomials of degree N in each variable. Gauss-Lobatto-Legendre (GLL) quadrature rules replace the integrals in the Galerkin formulation. This system is solved by the preconditioned conjugate gradients method. The conforming finite element space on the GLL mesh consisting of piecewise Q(sub 1) elements produces a stiffness matrix K(sub h) that is spectrally equivalent to the spectral element stiffness matrix K(sub N). The action of the inverse of K(sub h) is expensive for large problems, and is therefore replaced by a Schwarz preconditioner B(sub h) of this finite element stiffness matrix. The preconditioned operator then becomes B(sub h)(exp -l)K(sub N). The technical difficulties stem from the nonregularity of the mesh. Tools to estimate the convergence of a large class of new iterative substructuring and overlapping Schwarz preconditioners are developed. This technique also provides a new analysis for an iterative substructuring method proposed by Pavarino and Widlund for the spectral element discretization.

Casarin, Mario↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS2012

Two Dimensional Heat Transfer around Penetrations in Multilayer Insulation

The objective of this task was to quantify thermal losses involving integrating MLI into real life situations. Testing specifically focused on the effects of penetrations (including structural attachments, electrical conduit/feedthroughs, and fluid lines) through MLI. While there have been attempts at quantifying these losses both analytically and experimentally, none have included a thorough investigation of the methods and materials that could be used in such applications. To attempt to quantify the excess heat load coming into the system due to the integration losses, a calorimeter was designed to study two dimensional heat transfer through penetrated MLI. The test matrix was designed to take as many variables into account as was possible with the limited test duration and system size. The parameters varied were the attachment mechanism, the buffer material (for buffer attachment mechanisms only), the thickness of the buffer, and the penetration material. The work done under this task is an attempt to measure the parasitic heat loads and affected insulation areas produced by system integration, to model the parasitic loads, and from the model produce engineering equations to allow for the determination of parasitic heat loads in future applications. The methods of integration investigated were no integration, using a buffer to thermally isolate the strut from the MLI, and temperature matching the MLI on the strut. Several materials were investigated as a buffer material including aerogel blankets, aerogel bead packages, cryolite, and even an evacuated vacuum space (in essence a no buffer condition).

Johnson, Wesley L.↗
NASA NTRS1991

Sensitivity analysis of a wing aeroelastic response

A variation of Sobieski's Global Sensitivity Equations (GSE) approach is implemented to obtain the sensitivity of the static aeroelastic response of a three-dimensional wing model. The formulation is quite general and accepts any aerodynamics and structural analysis capability. An interface code is written to convert one analysis's output to the other's input, and visa versa. Local sensitivity derivatives are calculated by either analytic methods or finite difference techniques. A program to combine the local sensitivities, such as the sensitivity of the stiffness matrix or the aerodynamic kernel matrix, into global sensitivity derivatives is developed. The aerodynamic analysis package FAST, using a lifting surface theory, and a structural package, ELAPS, implementing Giles' equivalent plate model are used.

Kapania, Rakesh K.↗
NASA NTRS1985

Iterative methods for mixed finite element equations

Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.

Nakazawa, S.↗
NASA NTRS1985

Concentration of carbon dioxide by a high-temperature electrochemical membrane cell

The performance of a molten carbonate carbon dioxide concentrator (MCCDC) cell, as a device for removal of CO2 from manned spacecraft cabins without fuel expenditure, is investigated. The test system consists of an electrochemical cell (with an Li2CO3-38 mol pct K2CO3 membrane contained in a LiAlO2 matrix), a furnace, and a flow IR analyzer for monitoring CO2. Operation of the MCCDC-driven cell was found to be suitable for the task of CO2 removal: the cell performed at extremely low CO2 partial pressures (at or above 0.1 mm Hg); cathode CO2 efficiencies of 97 percent were achieved with 0.25 CO2 inlet concentration at 19 mA sq cm, at temperatures near 873 K. Anode concentrations of up to 5.8 percent were obtained. Simple cathode and anode performance equations applied to correlate cell performance agreed well with those measured experimentally. A flow diagram for the process is included.

Kang, M. P.↗
NASA NTRS1993

A study of metal-ceramic wettability in SiC-Al using dynamic melt infiltration of SiC

Pressure-assisted infiltration with a 2014 Al alloy of plain and Cu-coated single crystal platelets of alpha silicon carbide was used to study particulate wettability under dynamic conditions relevant to pressure casting of metal-matrix composites. The total penetration length of infiltrant metal in porous compacts was measured at the conclusion of solidification as a function of pressure, infiltration time, and SiC size for both plain and Cu-coated SiC. The experimental data were analyzed to obtain a threshold pressure for the effect of melt intrusion through SiC compacts. The threshold pressure was taken either directly as a measure of wettability or converted to an effective wetting angle using the Young-Laplace capillary equation. Cu coating resulted in partial but beneficial improvements in wettability as a result of its dissolution in the melt, compared to uncoated SiC.

Asthana, R.↗
NASA NTRS1993

Preconditioned Conjugate Gradient methods for low speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS1995

The FEM-R-Matrix Approach: Use of Mixed Finite Element and Gaussian Basis Sets for Electron Molecule Collisions

For the calculation of electron molecule collision cross sections R-matrix methods automatically take advantage of the division of configuration space into an inner region (I) bounded by radius tau b, where the scattered electron is within the molecular charge cloud and the system is described by an correlated Configuration Interaction (CI) treatment in close analogy to bound state calculations, and an outer region (II) where the scattered electron moves in the long-range multipole potential of the target and efficient analytic methods can be used for solving the asymptotic Schroedinger equation plus boundary conditions.

Thuemmel, Helmar T.↗
NASA NTRS1988

Finite element (MARC) solution technologies for viscoplastic analyses

A need for development of realistic constitutive models for structural components operating at high temperatures, accompanied by appropriate solution technologies for stress/life analyses of these components is studied. Viscoplastic models provide a better description of inelastic behavior of materials, but their mathematical structure is very complex. The highly nonlinear and stiff nature of the constitutive equations makes analytical solutions difficult. Therefore, suitable solution, finite element or other numerical, technologies must be developed to make these models adaptable for better and rational designs of components. NASA-Lewis has developed several solution technologies and successfully applied them to the solution of a number of uniaxial and multiaxial problems. Some of these solution technologies are described along with the models and representative results. The solution technologies developed and presented encompass a wide range of models, such as, isotropic, anisotropic, metal matrix composites, and single crystal models.

Arya, V. K.↗
NASA NTRS1990

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗