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,225 records · Page 68

A combined finite element-boundary element formulation for solution of axially symmetric bodies

A new method is presented for the computation of electromagnetic scattering from axially symmetric bodies. To allow the simulation of inhomogeneous cross sections, the method combines the finite element and boundary element techniques. Interior to a fictitious surface enclosing the scattering body, the finite element method is used which results in a sparce submatrix, whereas along the enclosure the Stratton-Chu integral equation is enforced. By choosing the fictitious enclosure to be a right circular cylinder, most of the resulting boundary integrals are convolutional and may therefore be evaluated via the FFT with which the system is iteratively solved. In view of the sparce matrix associated with the interior fields, this reduces the storage requirement of the entire system to O(N) making the method attractive for large scale computations. The details of the corresponding formulation and its numerical implementation are described.

Collins, Jeffrey D.↗

The prediction of pseudo-resonance positions in the Schwinger variational principle

The Schwinger Variational Principle is appled to s-wave electron-hydrogen atom scattering. It is shown computationally, that, consistent with a paper by B. Apagyi, P. Levay, and K. Ladanyi, there are pseudo-resonances at the static exchange level of approximation, but not at the static level. The T-matrix was employed, as well as, the K-matrix version of the Schwinger Principle, with a real Slater basis, and obtained the same results in both. The origin of the pseudo-resonances as resulting from singularities in the separable potential that is effectively employed in the Lippmann-Schwinger equation from which the Schwinger Variational Principle can be derived. The determination of the pseudo-resonance parameters from the separable potential is computationally inexpensive and may be used to predict the pseudo-resonance parameters for the scattering calculations so that they may be avoided.

Weatherford, Charles A.↗

Nonlinear programming extensions to rational function approximations of unsteady aerodynamics

This paper deals with approximating unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft. Two methods of formulating these approximations are extended to include both the same flexibility in constraining them and the same methodology in optimizing nonlinear parameters as another currently used 'extended least-squares' method. Optimal selection of 'nonlinear' parameters is made in each of the three methods by use of the same nonlinear (nongradient) optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is of lower order than that required when no optimization of the nonlinear terms is performed. The free 'linear' parameters are determined using least-squares matrix techniques on a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from the different approaches are described, and results are presented which show comparative evaluations from application of each of the extended methods to a numerical example. The results obtained for the example problem show a significant (up to 63 percent) reduction in the number of differential equations used to represent the unsteady aerodynamic forces in linear time-invariant equations of motion as compared to a conventional method in which nonlinear terms are not optimized.

Tiffany, Sherwood H.↗

A triangular thin shell finite element: Linear analysis

The formulation of the linear stiffness matrix for a doubly-curved triangular thin shell element, using a modified potential energy principle, is described. The strain energy component of the potential energy is expressed in terms of displacements and displacement gradients by use of consistent Koiter strain-displacement equations. The element inplane and normal displacement fields are approximated by complete cubic polynomials. The interelement displacement admissibility conditions are met in the global representation by imposition of constraint conditions on the interelement boundaries; the constraints represent the modification of the potential energy. Errors due to the nonzero strains under rigid body motion are shown to be of small importance for practical grid refinements through performance of extensive comparison analyses.

Thomas, G. R.↗

Light scattering by tenuous particles - A generalization of the Rayleigh-Gans-Rocard approach

We consider scattering by arbitrarily shaped particles that satisfy two conditions: (1) that the polarizability of the particle relative to the ambient medium be small compared to 1 and (2) that the phase shift introduced by the particle be less than 2. We solve the integro-differential equation proposed by Shifrin by using the method of successive iterations and then applying a Fourier transform. For the second iteration, results are presented that accurately describe scattering by a broad class of particles. The phase function and other elements of the scattering matrix are shown to be in excellent agreement with Mie theory for spherical scatterers.

Acquista, C.↗

Eigenvalues and stability problems of rotors

The essential theoretical results of the application of a developed transfer matrix method to the free transverse vibration of a rotor are shown. Gyroscopic and shear effects, rotary inertia, and external and internal damping as well as the influence of sleeve bearings and rotor supports are taken into consideration. The eigenvalues of the motion equations of the rotor are searched by using a modified determinant method.

Walczyk, Z.↗

Quasinormal modes of Reissner-Nordstrom black holes

A matrix-eigenvalue algorithm is presented for accurately computing the quasi-normal frequencies and modes of charged static blackholes. The method is then refined through the introduction of a continued-fraction step. The approach should generalize to a variety of nonseparable wave equations, including the Kerr-Newman case of charged rotating blackholes.

Leaver, Edward W.↗

Fully Deuterated Aliphatic Hydrocarbons Obtained From Iron Carbide Treated with DCl and D2O

According to Oparin, Mendeleev thought that the origin of petroleum was the result of the hydrolysis of iron carbides by superheated steam under pressure from the deep interior of the Earth through geological formations where the metal carbides exist. As early as 1877, Mendeleev described the reaction leading to the synthesis of hydrocarbons according to the general equation 3Fe(sub m)C(sub n) + mH2O yields mFe3O4 + C(sub 3n)H(sub 8m). Other experimental studies on the production of hydrocarbons from cast iron have been reported. Because of the possibility that hydrocarbons may have been trapped within the carbon matrix of the cast iron, which usually has a high content of carbon, we have studied the reaction of pure iron carbide with deuterium chloride and deuterated water. This was done in order to distinguish any newly formed deuterated hydrocarbons from any possible impurities of trapped hydrocarbons. The experiments were carried out by simply allowing iron carbide to react with concentrated deuterium chloride in D2O. The volatile hydrocarbon fraction examined by gas chromatography-mass spectrometry (GC/MS), using a Finnigan 1020/OWA instrument. contained low molecular weight hydrocarbons in a range C3 to C7. Lower molecular weight hydrocarbons were not detected by GC/MS because the MS scanning mode was preset above mass 40 to exclude components of air. The identified hydrocarbons are similar to those obtained under prebiotic conditions using high frequency discharge. The hydrocarbons found in common were propane, butane, pentane, 3-methylpentane, hexane, and heptane. The percent yields decline with increasing carbon number (propane 11%, n-heptane 1%). Similar results have been obtained by the direct treatment of metal carbides by pulse laser vaporization mass spectrometry. These results show that the hydrolysis of iron carbides may have been a significant source of hydrocarbons on the primitive Earth. There appears to be a predominance of straight chain isomers, which may have been important in the prebiotic synthesis of membrane-forming compounds.

Marquez, C.↗

An Upper Bound on Orbital Debris Collision Probability When Only One Object has Position Uncertainty Information

Upper bounds on high speed satellite collision probability, P (sub c), have been investigated. Previous methods assume an individual position error covariance matrix is available for each object. The two matrices being combined into a single, relative position error covariance matrix. Components of the combined error covariance are then varied to obtain a maximum P (sub c). If error covariance information for only one of the two objects was available, either some default shape has been used or nothing could be done. An alternative is presented that uses the known covariance information along with a critical value of the missing covariance to obtain an approximate but useful P (sub c) upper bound. There are various avenues along which an upper bound on the high speed satellite collision probability has been pursued. Typically, for the collision plane representation of the high speed collision probability problem, the predicted miss position in the collision plane is assumed fixed. Then the shape (aspect ratio of ellipse), the size (scaling of standard deviations) or the orientation (rotation of ellipse principal axes) of the combined position error ellipse is varied to obtain a maximum P (sub c). Regardless as to the exact details of the approach, previously presented methods all assume that an individual position error covariance matrix is available for each object and the two are combined into a single, relative position error covariance matrix. This combined position error covariance matrix is then modified according to the chosen scheme to arrive at a maximum P (sub c). But what if error covariance information for one of the two objects is not available? When error covariance information for one of the objects is not available the analyst has commonly defaulted to the situation in which only the relative miss position and velocity are known without any corresponding state error covariance information. The various usual methods of finding a maximum P (sub c) do no good because the analyst defaults to no knowledge of the combined, relative position error covariance matrix. It is reasonable to think, given an assumption of no covariance information, an analyst might still attempt to determine the error covariance matrix that results in an upper bound on the P (sub c). Without some guidance on limits to the shape, size and orientation of the unknown covariance matrix, the limiting case is a degenerate ellipse lying along the relative miss vector in the collision plane. Unless the miss position is exceptionally large or the at-risk object is exceptionally small, this method results in a maximum P (sub c) too large to be of practical use. For example, assuming that the miss distance is equal to the current ISS alert volume along-track (+ or -) distance of 25 kilometers and that the at-risk area has a 70 meter radius. The maximum (degenerate ellipse) P (sub c) is about 0.00136. At 40 kilometers, the maximum P (sub c) would be 0.00085 which is still almost an order of magnitude larger than the ISS maneuver threshold of 0.0001. In fact, a miss distance of almost 340 kilometers is necessary to reduce the maximum P (sub c) associated with this degenerate ellipse to the ISS maneuver threshold value. Such a result is frequently of no practical value to the analyst. Some improvement may be made with respect to this problem by realizing that while the position error covariance matrix of one of the objects (usually the debris object) may not be known the position error covariance matrix of the other object (usually the asset) is almost always available. Making use of the position error covariance information for the one object provides an improvement in finding a maximum P (sub c) which, in some cases, may offer real utility. The equations to be used are presented and their use discussed.

Frisbee, Joseph H., Jr.↗

Two-dimensional computer simulation of EMVJ and grating solar cells under AMO illumination

A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.

Gray, J. L.↗

A three-dimensional pressure flux-split RNS application to sub/supersonic flow in inlets and ducts

The reduced Navier-Stokes (RNS) formulation is combined with a pressure based flux split procedure for the computational analysis of three dimensional flow. Only lowest order diffusion, that is required to satisfy the no-slip boundary condition is retained in this approximation. The governing equations are appropriately differenced such that the physical boundary conditions are closely coupled and combined to form a closed discrete system. No numerical or characteristic type boundary conditions are required. A global pressure relaxation procedure, i.e. multi-sweep PNS, is considered. A sparse matrix direct solver is applied to the crossplane solution. At the outflow boundary only the pressure or pressure gradient is prescribed. The resulting formulation has been applied to the solution of three-dimensional subsonic/supersonic internal flows in inlets and ducts.

Khosla, P. K.↗

Tensor-GMRES method for large sparse systems of nonlinear equations

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Feng, Dan↗

Composite Nanomechanics: A Mechanistic Properties Prediction

A unique mechanistic theory is described to predict the properties of nanocomposites. The theory is based on composite micromechanics with progressive substructuring down to a nanoscale slice of a nanofiber where all the governing equations are formulated. These equations have been programmed in a computer code. That computer code is used to predict 25 properties of a mononanofiber laminate. The results are presented graphically and discussed with respect to their practical significance. Most of the results show smooth distributions. Results for matrix-dependent properties show bimodal through-the-thickness distribution with discontinuous changes from mode to mode.

Chamis, Christos C.↗

Hybrid-dual-fourier tomographic algorithm for a fast three-dimensionial optical image reconstruction in turbid media

A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model.

Alfano, Robert R.↗

Vibration of prestressed periodic lattice structures

Equations are developed for vibration of general lattice structures that have repetitive geometry. The method of solution is an extension of a previous paper for buckling of similar structures. The theory is based on representing each member of the structure with the exact dynamic stiffness matrix and taking advantage of the repetitive geometry to obtain an eigenvalue problem involving the degrees-of-freedom at a single node in the lattice. Results are given for shell-and beam-like lattice structures and for rings stiffened with tension cables and a central mast. The variation of frequency with external loading and the effect of local member vibration on overall modes is shown.

Anderson, M. S.↗

Variability simulations with a steady, linearized primitive equations model

Solutions of the steady, primitive equations on a sphere, linearized about a zonally symmetric basic state are computed for the purpose of simulating monthly mean variability in the troposphere. The basic states are observed, winter monthly mean, zonal means of zontal and meridional velocities, temperatures and surface pressures computed from the 15 year NMC time series. A least squares fit to a series of Legendre polynomials is used to compute the basic states between 20 H and the equator, and the hemispheres are assumed symmetric. The model is spectral in the zonal direction, and centered differences are employed in the meridional and vertical directions. Since the model is steady and linear, the solution is obtained by inversion of a block, pente-diagonal matrix. The model simulates the climatology of the GFDL nine level, spectral general circulation model quite closely, particularly in middle latitudes above the boundary layer. This experiment is an extension of that simulation to examine variability of the steady, linear solution.

Kinter, J. L., III↗

Upwind relaxation methods for the Navier-Stokes equations using inner iterations

A supersonic and subsonic test problem is solved using an upwind line relaxation algorithm for the Navier-Stokes equations which employs inner iterations. In the supersonic test case, the inner iteration method was found to require only about 65 percent of the CPU time which was required by the most efficient line relaxation method without inner iterations. In the subsonic test case, poor matrix conditioning resulted in less efficiency than line relaxation methods employing a more conventional CPU-saving strategy.

Taylor, Arthur C., III↗

Spectral simulation of unsteady compressible flow past a circular cylinder

An unsteady compressible viscous wake flow past a circular cylinder was successfully simulated using spectral methods. A new approach in using the Chebyshev collocation method for periodic problems is introduced. It was further proved that the eigenvalues associated with the differentiation matrix are purely imaginary, reflecting the periodicity of the problem. It was been shown that the solution of a model problem has exponential growth in time if improper boundary conditions are used. A characteristic boundary condition, which is based on the characteristics of the Euler equations of gas dynamics, was derived for the spectral code. The primary vortex shedding frequency computed agrees well with the results in the literature for Mach = 0.4, Re = 80. No secondary frequency is observed in the power spectrum analysis of the pressure data.

Don, Wai-Sun↗