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.
Several approximations to the Doppler broadening functions psi(x, theta) and chi(x, theta) are compared with respect to accuracy and speed of evaluation. A technique, due to A. M. Turning (1943), is shown to be at least as accurate as direct numerical quadrature and somewhat faster than Gaussian quadrature. FORTRAN 4 listings are included.
Three-dimensional finite element analysis was used to simulate progressive failure of a plain weave composite subjected to in-plane extension. The loading was parallel to one of the tow directions. The effects of various characteristics of the finite element model on predicted behavior were examined. The predicted behavior was found to be sensitive to quadrature order, mesh refinement, and the material degradation model. Also the sensitivity of the predictions to the tow waviness was studied. The predicted strength decreased considerably with increased waviness. More numerical studies and comparisons with experimental data are needed to establish reliable guidelines for accurate progressive failure prediction.
The usefulness of long-wavelength potential field anomalies in lithospheric interpretation is greatly increased with spherical Earth modeling techniques. Gauss-Legendre quadrature integration is used to calculate the anomalous potential of gravity and magnetic fields and their spatial derivatives on a spherical Earth for an arbitrary body represented by an equivalent point source distribution of gravity poles or magnetic dipoles. The distribution of equivalent point sources is determined directly from the coordinate limits of the source volume. Variable integration limits for an arbitrarily shaped body are derived from interpolation of points which approximate the body's surface envelope. The versatility of the method is enhanced by the ability to treat physical property variations within the source volume and to consider variable magnetic fields over the source and observation surface. A number of examples verify and illustrate the capabilities of the technique, including preliminary modeling of potential field signatures for Mississippi embayment crustal structure at satellite elevations.
A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.
Two iterative methods are considered, Richardson's method and a general second order method. For both methods, a variant of the method is derived for which only even numbered iterates are computed. The variant is called a leapfrog method. Comparisons between the conventional form of the methods and the leapfrog form are made under the assumption that the number of unknowns is large. In the case of Richardson's method, it is possible to express the final iterate in terms of only the initial approximation, a variant of the iteration called the grand-leap method. In the case of the grand-leap variant, a set of parameters is required. An algorithm is presented to compute these parameters that is related to algorithms to compute the weights and abscissas for Gaussian quadrature. General algorithms to implement the leapfrog and grand-leap methods are presented. Algorithms for the important special case of the Chebyshev method are also given.
We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.
In the preferred embodiment, an encrypted GPS signal is down-converted from RF to baseband to generate two quadrature components for each RF signal (L1 and L2). Separately and independently for each RF signal and each quadrature component, the four down-converted signals are counter-rotated with a respective model phase, correlated with a respective model P code, and then successively summed and dumped over presum intervals substantially coincident with chips of the respective encryption code. Without knowledge of the encryption-code signs, the effect of encryption-code sign flips is then substantially reduced by selected combinations of the resulting presums between associated quadrature components for each RF signal, separately and independently for the L1 and L2 signals. The resulting combined presums are then summed and dumped over longer intervals and further processed to extract amplitude, phase and delay for each RF signal. Precision of the resulting phase and delay values is approximately four times better than that obtained from straight cross-correlation of L1 and L2. This improved method provides the following options: separate and independent tracking of the L1-Y and L2-Y channels; separate and independent measurement of amplitude, phase and delay L1-Y channel; and removal of the half-cycle ambiguity in L1-Y and L2-Y carrier phase.
The instantaneous values of output voltages representing the wind velocity vector and the temperature at different elevations of the 250-foot meteorological tower located at NASA Wallops Flight Center are provided with the three dimensional split-film TSI Model 1080 anemometer system. The output voltages are sampled at a rate of one every 5 milliseconds, digitized and stored on digital magnetic tapes for a time period of approximately 40 minutes, with the use of a specially designed data acqusition system. A new calibration procedure permits the conversion of the digital voltages to the respective values of the temperature and the velocity components in a Cartesian coordinate system connected with the TSI probe with considerable accuracy. Power, cross, coincidence and quadrature spectra of the wind components and the temperature are obtained with the use of the fast Fourier transform. The cosine taper data window and ensemble and frequency smoothing techniques are used to provide smooth estimates of the spectral functions.
The three-dimensional equation of radiative transfer is formally solved using a Fourier-Riccati approach while calculations are performed on cloudy media embedded in a two-dimensional space. An extension to Stephens' work this study addresses the coupling between space and angle asserted by the equation of transfer. In particular, the accuracy of the computed radiation field as it is influenced by the angular resolution of the phase function and spatial discretization of the cloudy medium is discussed. The necessity of using a large number of quadrature points to calculate fluxes even when the phase function is isotropic for media exhibiting vertical and horizontal inhomogeneities is demonstrated. Effects of incorrect spatial sampling on both radiance and flux fields are also quantified by example. Radiance and flux comparisons obtained by the Fourier-Riccati model and the independent pixel approximation for inhomogeneous cloudy media illustrate the inadequacy of the latter even for tenuous clouds.
A calculation method is presented for sound propagation over an impedance discontinuity in flat ground with a homogeneous, still atmosphere. The method is based on an approximate solution to a two dimensional boundary integral equation formulation of the problem, which expresses the wave field as the solution for homogeneous ground plus an integral over half of the boundary. Through recognizing this integral as a generalized Fourier integral, asymptotic methods are applied to evaluate the part of the integral most expensive to compute by numerical quadrature. Single frequency excess attenuation results for propagation from a point source above rigid ground to a receiver above absorbing ground are discussed. The results are applied, with air attenuation and A-weighting, to a notional jet engine noise source. Simple trends are noted.
Some spectral type collocation method well suited for the approximation of fourth-order systems are proposed. The model problem is the biharmonic equation, in one and two dimensions when the boundary conditions are periodic in one direction. It is proved that the standard Gauss-Lobatto nodes are not the best choice for the collocation points. Then, a new set of nodes related to some generalized Gauss type quadrature formulas is proposed. Also provided is a complete analysis of these formulas including some new issues about the asymptotic behavior of the weights and we apply these results to the analysis of the collocation method.
Two potential coefficient fields that are complete to degree and order 360 have been computed. One field (OSU86E) excludes geophysically predicted anomalies while the other (OSU86F) includes such anomalies. These fields were computed using a set of 30' mean gravity anomalies derived from satellite altimetry in the ocean areas and from land measurements in North America, Europe, Australia, Japan and a few other areas. Where no 30' data existed, 1 deg x 1 deg mean anomaly estimates were used if available. No rigorous combination of satellite and terrestrial data was carried out. Instead advantage was taken of the adjusted anomalies and potential coefficients from a rigorous combination of the GEML2' potential coefficient set and 1 deg x 1 deg mean gravity anomalies. The two new fields were computed using a quadrature procedure with de-smoothing factors. The spectra of the new fields agree well with the spectra of the fields with 1 deg x 1 deg data out to degree 180. Above degree 180 the new fields have more power. The fields have been tested through comparison of Doppler station geoid undulations with undulations from various geopotential models. The agreement between the two types of undulations is approximately + or - 1.6 m. The use of a 360 field over a 180 field does not significantly improve the comparison. Instead it allows the comparison to be done at some stations where high frequency effects are important. In addition maps made in areas of high frequency information (such as trench areas) clearly reveal the signal in the new fields from degree 181 to 360.
Westinghouse Electric Corporation Division under NASA contract NAS3-25865 developed a proof-of-concept (POC) multichannel demultiplexer implemented as an acousto-optic radiofrequency (RF) with a spectrum analyzer. A detailed analysis of the experimental results indicate that the expected degradation caused by the acousto-optical channelizer is approximately 2.0 dB degradation at 10(exp -5) bit-error rate (BER) and 3.0 dB degradation at 10(exp -8) BER. This degradation may be quite acceptable when considering the excellent volume, mass, and power characteristics of acousto-optical channelizing relative to other technologies. In addition, system performance can be greatly improved by using digital pulse shaping in the modem and increasing the channel spacing from 40 to 45 kHz for 64 kbps quadrature phase-shift keying (QPSK) modulation.
Gauss-Legendre quadrature integration is used to calculate the anomalous potential of gravity and magnetic fields and their spatial derivatives on a spherical earth. The procedure involves representation of the anomalous source as a distribution of equivalent point gravity poles or point magnetic dipoles. The distribution of equivalent point sources is determined directly from the volume limits of the anomalous body. The variable limits of integration for an arbitrarily shaped body are obtained from interpolations performed on a set of body points which approximate the body's surface envelope. The versatility of the method is shown by its ability to treat physical property variations within the source volume as well as variable magnetic fields over the source and observation surface. Examples are provided which illustrate the capabilities of the technique, including a preliminary modeling of potential field signatures for the Mississippi embayment crustal structure at 450 km.
Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.
The effect of the adhesive layer on crack propagation in composite materials is investigated. The composite medium consists of parallel load carrying laminates and buffer strips arranged periodically and bonded with thin adhesive layers. The strips, assumed to be isotropic and linearly elastic, contain symmetric cracks of arbitrary lengths located normal to the interfaces. Two problems are considered: (1) thin adhesive layers are approximated by uncoupled tension and shear springs distributed along the interfaces of the strips for which only the case of internal cracks can be treated rigorously; (2) broken laminates and the true singular behavior in the presence of the adhesive layer are studied. The adhesive is then treated as an isotropic, linearly elastic continuum. General expressions for field quantities are obtained in terms of infinite Fourier integrals. These expressions give a system of singular integral equations in terms of the crack surface displacement derivatives. By using appropriate quadrature formulas, the integral equations reduce to a system of linear algebraic equations which are solved numerically.
The linearized global primitive equation model of Salby and Garcia (1987) was modified to study tropospheric circulations and their interaction with tropical convection. The vertical resolution was increased to approximately 1.5 km, with the effective vertical domain extending the approximately 40 km (higher altitudes in the calculation are used to ensure the radiation upper boundary condition is satisfied). A formal wave conditional instability of the second kind (CISK) parameterization was introduced. That parameterization allows low-frequency disturbances to interact with convection by organizing surface moisture. Boundary layer friction has two important effects in these calculations: it exaggerates surface convergence near the equator by driving the flow out of geostrophic balance; and it introduces a phase shift between the convergence pattern and the circulation, which drives the heating out of quadrature with the temperature field and allows a positive feedback to take place. The observational complement of this theoretical work involves Global Cloud Imagery (GCI) composited from six satellites in the International Satellite Cloud Climatology Project (ISCCP) which were simultaneously observing the Earth's cloud field. Having global coverage, horizontal resolution of half a degree, and produced synoptically every 3 hours, the GCI affords an unprecendented view of global convection. A full year (1983 to 1984) of GCI has now been successfully created. Among other phenomena, that imagery reveals a regular progression of easterly waves in the convective pattern of the Western Pacific. Those disturbances have odd symmetry about the equator and organize convection into a pattern with similar symmetry. In time-mean maps of cloud cover, these propagating disturbances produce a split Inter-Tropical Convergence Zone (ITCZ) in the western Pacific.
A theoretical investigation of the velocity profiles for laminar mixing of a high-velocity stream with a region of fluid at rest has been made assuming that the Prandtl number is unity. A method which involves only quadratures is presented for calculating the velocity profile in the mixing layer for an arbitrary value of the free-stream Mach number. Detailed velocity profiles have been calculated for free-stream Mach numbers of 0, 1, 2, 3, and 5. For each Mach number, velocity profiles are presented for both a linear and a 0.76-power variation of viscosity with absolute temperature. The calculations for a linear variation are much simpler than those for a 0.76-power variation. It is shown that by selecting the constant of proportionality in the liner approximation such that it gives the correct value for the viscosity in the high-temperature part of the mixing layer, the resulting velocity profiles are in excellent agreement with those calculated by a 0.76-power variation.