SEARCH · Engineering Papers
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.
Diatomic molecules inelastic collision cross sections for specific rotational transitions, discussing S matrix energy requirements for statistical analysis
A new approximate theory which links the inherent flaw concept with the theory of crack tip stress singularities at a bi-material interface was developed. Three assumptions were made: (1) the existence of inherent flaw (i.e., damage zone) at the tip of the crack, (2) a fracture of the filamentary composites initiates at a crack lying in the matrix material at the interface of the matrix/filament, and (3) the laminate fails whenever the principal load-carrying laminae fails. This third assumption implies that for a laminate consisting of 0 degree plies, cracks into matrix perpendicular to the 0 degree filaments are the triggering mechanism for the final failure. Based on this theory, a parameter bar K sub Q which is similar to the stress intensity factor for isotropic materials but with a different dimension was defined. Utilizing existing test data, it was found that bar K sub Q can be treated as a material constant. Based on this finding a fracture mechanics analysis methodology was developed. The analytical results are correlated well with test results. This new approximate theory can apply to both brittle and metal matrix composite laminates with crack or hole.
An economical algorithm is presented for predicting the position of a satellite perturbed by drag and zonal harmonics J2 through J4. Simplicity being of the essence, drag is modeled as a secular decay rate in the semimajor axis (retarded motion) with the zonal perturbations modeled from a modified version of Brouwers formulas. The algorithm is developed as an alternative on-board orbit predictor; a back up propagator requiring low energy consumption; or a ground based propagator for microcomputer applications (e.g., at the foot of an antenna). An O(J2) secular retarded state partial matrix (matrizant) is also given to employ with state estimation. The theory has been implemented in BASIC on an inexpensive microcomputer, the program occupying under 8K bytes of memory. Simulated trajectory data and real tracking data are employed to illustrate the theory's ability to accurately accommodate oblateness and drag effects.
Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.
An economical algorithm is presented for predicting the position of a satellite perturbed by drag and zonal harmonics J sub 2 through J sub 4. Simplicity being of the essence, drag is modeled as a secular decay rate in the semi-axis (retarded motion); with the zonal perturbations modeled from a modified version of the Brouwers formulas. The algorithm is developed as: an alternative on-board orbit predictor; a back up propagator requiring low energy consumption; or a ground based propagator for microcomputer applications (e.g., at the foot of an antenna). An O(J sub 2) secular retarded state partial matrix (matrizant) is also given to employ with state estimation. The theory was implemented in BASIC on an inexpensive microcomputer, the program occupying under 8K bytes of memory. Simulated trajectory data and real tracking data are employed to illustrate the theory's ability to accurately accommodate oblateness and drag effects.
The available micromechanics theories for predicting the uniaxial strengths of unidirectional composites are described. The theories presuppose that the required properties of the constituents of the composites are known and that the uniaxial loads are coincident with the principal axes of the unidirectional composite. The uniaxial strengths considered are those required to resist either in-plane loads, flexural (bending), or horizontal-beam shear. The discussion deals primarily with fiber-nonmetallic matrix composites where all the fibers are parallel and lie in the same plane. However, the micromechanics theories presented here are transferable to fiber-metallic matrix composites so long as the underlying assumptions are not violated. Some of the concepts described herein may also be applicable to particulate composites.
The theory of multiphoton ionization for an atomic system of arbitrary complexity is developed using a density matrix formalism. An expression is obtained which determines the differential N-photon ionization cross section as a function of the polarization states of the target atom and the incident radiation. The parameters which characterize the photoelectron angular distribution are related to the general reduced matrix elements for the N-photon transition. Two-photon ionization of unpolarized atoms is treated as an illustration of the use of the theory. The dependence of the multiphoton ionization cross section on the polarization state of the incident radiation, which has been observed in two- and three-photon ionization of Cs, is accounted for by the theory. Finally, the photoelectron spin polarization produced by the multiphoton ionization of unpolarized atoms, like the analogous polarization resulting from single-photon ionization, is found to depend on the circular polarization of the incident radiation.
The theory of multiphoton ionization for an atomic system of arbitrary complexity is developed using a density matrix formalism. An expression is obtained which determines the differential N-photon ionization cross section as a function of the polarization states of the target atom and the incident radiation. The parameters which characterize the photo-electron angular distribution are related to the general reduced matrix elements for the N-photon transition. Two-photon ionization of unpolarized atoms is treated as an illustration of the use of the theory. The dependence of the multiphoton ionization cross section on the polarization state of the incident radiation, which has been observed in two- and three-photon ionization of Cs, is accounted for by the theory. Finally, the photoelectron spin polarization produced by the multiphoton ionization of unpolarized atoms, like the analogous polarization resulting from single-photon ionization, is found to depend on the circular polarization of the incident radiation.
Covariant generalization of strong interaction physics using S-matrix - theoretical physics
A two-dimensional mechanical model is presented to predict the compressive strength of unidirectional fiber composites using technical beam theory and classical elasticity. First, a single fiber resting on a matrix half-plane is considered. Next, a more elaborate analysis of a uniformly laminated, unidirectional fiber composite half-plane is presented. The model configuration incorporates a free edge which introduces a buckling mode that originates at the free edge and decays into the inferior of the half-plane. It is demonstrated that for composites of low volume fraction (less than 0.3), this decay mode furnishes values of buckling strain that are below the values predicted by the Rosen (1965) model. At a higher volume fraction the buckling mode corresponds to a half wavelength that is in violation of the usual assumptions of beam theory. Causes for deviations of the model prediction froi?i existing experimental results are discussed.
A two-dimensional mechanical model is presented to predict the compressive strength of unidirectional fiber composites using technical beam theory and classical elasticity. First, a single fiber resting on a matrix half-plane is considered. Next, a more elaborate analysis of a uniformly laminated, unidirectional fiber composite half-plane is presented. The model configuration incorporates a free edge which introduces a buckling mode that originates at the free edge and decays into the interior of the half-plane. It is demonstrated that for composites of low volume fraction (less than 0.3), this decay mode furnishes values of buckling strain that are below the values predicted by the Rosen (1965) model. At a higher volume fraction the buckling mode corresponds to a half-wavelength that is in violation of the usual assumptions of beam theory. Causes for deviations of the model prediction from existing experimental results are discussed.
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.
Two efficient numerical methods for dealing with the stability of linear periodic systems are presented. Both methods combine the use of multivariable Floquet-Liapunov theory with an efficient numerical scheme for computing the transition matrix at the end of one period. The numerical properties of these methods are illustrated by applying them to the simple parametric excitation problem of a fixed end column. The practical value of these methods is shown by applying them to some helicopter rotor blade aeroelastic and structural dynamics problems. It is concluded that these methods are numerically efficient, general and practical for dealing with the stability of large periodic systems.
A strategy was developed to permit trade-offs between the number of floating point operations required and the storage requirements for the solution of certain difference problems, such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement roughly doubles the number of computations. Reducing the storage more than this tends to make the number of computations prohibitively large. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5m(2) + 1 temporary storage cells. In many cases m is a constant and quite small. For example, in solving a factored form of the three-dimensional Navier-Stokes equations, the size m of the block tridiagonals is 5. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.
Theoretical and experimental studies were conducted to determine the noise transmission into acoustic enclosures ranging from simple rectangular box models to full scale light aircraft in flight. The structural models include simple, stiffened, curved stiffened, and orthotropic panels and double wall windows. The theoretical solutions were obtained by model analysis. Transfer matrix and finite element procedures were utilized. Good agreement between theory and experiment has been achieved. An efficient acoustic add-on treatment was developed for interior noise control in a twin engine light aircraft.
In many programs solving difference equations, problem size is restricted by the number of available memory cells. A strategy has been developed to permit trade-offs between the number of floating point operations required and storage requirements for the solution of certain problems such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement will roughly double the number of computations. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5 sq m + 1 temporary storage cells. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.
The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. ScramJet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel massflow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.
The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation Indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. Scramjet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel mass flow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.