Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 559 records · Page 31

Fracture and contact problems for an elastic wedge

The plane elastostatic contact problem for an infinite elastic wedge of arbitrary angle is discussed. The medium is loaded through a frictionless rigid wedge of a given symmetric profile. Using the Mellin transform formulation the mixed boundary value problem is reduced to a singular integral equation with the contact stress as the unknown function. With the application of the results to the fracture of the medium in mind, the main emphasis in the study has been on the investigation of the singular nature of the stress state around the apex of the wedge and on the determination of the contact pressure.

Erdogan, F.↗

The ion optics of a two grid electron-bombardment thruster

A detailed experimental investigation has been performed to determine the ion beam divergence of an electron-bombardment ion thruster as a function of grid geometry changes. The results show that, to a good approximation, each geometrical grid parameter independently affects one aspect of grid set performance. These observations are used to develop a graphical technique for predicting the ion beam divergence of an arbitrary ion source and grid geometry combination. The usefulness of this technique is demonstrated by comparing predicted ion beam divergence of the 30-cm diameter Engineering Model ion thruster with independent experimental determinations. Good agreement is shown between predicted and experimental results.

Aston, G.↗

Reflections, diffractions, and surface waves for an interior impedance wedge of arbitrary angle

The asymptotic-impedance wedge solution for plane-wave illumination at normal incidence is examined for interior wedge diffraction. An efficient method for calculating the diffraction coefficient for arbitrary wedge angle is presented. The asymptotic solution isolates the incident, singly reflected, multiply reflected, diffracted, surface-wave, and associated-surface-wave transition fields. Multiply reflected fields (of any order) from the exact solution arise as ratios of auxiliary Maliuzhinets functions; however, by using properties of these functions, the representation can be reduced to products of reflection coefficients, much more efficient for calculation. A surface-wave transition field is added to the surface wave to retain continuity of the total field at the surface wave boundaries. This formulation is equally valid for both exterior and interior wedges with uniform but different impedances on each face, for both soft and hard polarizations.

Griesser, Timothy↗

Coupled elastic response of open section laminated composite beams subject to generalized beam loading

A generalized beam theory for arbitrary open section laminated composite beams subjected to generalized beam loading is developed using a strength of materials approach. At present, the theory is limited to statistically determinate beams. Solution of the general system of equations yields strain and curvature distributions as functions of the axial coordinate x and the contour coordinate s. Average stresses on a given cross section or ply stresses in a local coordinate system are determined using these distributions. Data obtained for thin-walled open section composite beams reveal that elastic coupling occurs among curvatures, shearing strains, and axial strain, depending on the laminate stiffness parameters. Laminate ply layup does not affect the location of the shear center in the thin-wall approximation, but may have a significant impact for thicker walled open sections. The average axial stress distribution is found to be distinctly different from the isotropic distribution, but the average shear stress distribution is the same as that in the isotropic case.

Zvarick, Albert G.↗

QPOs from Random X-ray Bursts around Rotating Black Holes

We continue our earlier studies of quasi-periodic oscillations (QPOs) in the power spectra of accreting, rapidly-rotating black holes that originate from the geometric 'light echoes' of X-ray flares occurring within the black hole ergosphere. Our present work extends our previous treatment to three-dimensional photon emission and orbits to allow for arbitrary latitudes in the positions of the distant observers and the X-ray sources in place of the mainly equatorial positions and photon orbits of the earlier consideration. Following the trajectories of a large number of photons we calculate the response functions of a given geometry and use them to produce model light curves which we subsequently analyze to compute their power spectra and autocorrelation functions. In the case of an optically-thin environment, relevant to advection-dominated accretion flows, we consistently find QPOs at frequencies of order of approximately kHz for stellar-mass black hole candidates while order of approximately mHz for typical active galactic nuclei (approximately equal to 10(exp 7) solar mass) for a wide range of viewing angles (30 degrees to 80 degrees) from X-ray sources predominantly concentrated toward the equator within the ergosphere. As in out previous treatment, here too, the QPO signal is produced by the frame-dragging of the photons by the rapidly-rotating black hole, which results in photon 'bunches' separated by constant time-lags, the result of multiple photon orbits around the hole. Our model predicts for various source/observer configurations the robust presence of a new class of QPOs, which is inevitably generic to curved spacetime structure in rotating black hole systems.

Kukumura, Keigo↗

Unsteady flow past wings having sharp-edge separation

A vortex-lattice technique is developed to model unsteady, incompressible flow past thin wings. This technique predicts the shape of the wake as a function of time; thus, it is not restricted by planform, aspect ratio, or angle of attack as long as vortex bursting does not occur and the flow does not separate from the wing surface. Moreover, the technique can be applied to wings of arbitrary curvature undergoing general motion; thus, it can treat rigid-body motion, arbitrary wing deformation, gusts in the freestream, and periodic motions. Numerical results are presented for low-aspect rectangular wings undergoing a constant-rate, rigid-body rotation about the trailing edge. The results for the unsteady motion are compared with those predicted by assuming quasi-steady motion. The present results exhibit hysteretic behavior.

Atta, E. H.↗

Shear flow aerodynamics - Lifting surface theory

A lifting surface theory based on a parallel shear flow model is presented for steady, incompressible flows. The theory is intended to account approximately for the presence of a boundary layer. The method of Fourier transforms is used to calculate the pressure on a surface of infinite extent and arbitrary contour. Immediately above the surface is a region of sheared flow (the boundary layer), outside of which the flow velocity is constant. The Fourier transform of the pressure on this surface is used to derive the shear flow equivalent to the kernel function of classical potential flow lifting surface theory. The kernel function provides an integral relation between the upwash at a given point on the surface and the pressure everywhere on the surface. This relation is treated as an integral equation for the pressure, and is solved numerically. Computations are presented for the lift and pitching moment on a flat plate in two-dimensional flow, and for flat, rectangular wings of aspect ratio 1, 2, and 5. As expected, the shear layer decreases the lift curve slope; however, the shear layer (whose thickness is constant along the wing chord) has little effect on the center of pressure.

Ventres, C. S.↗

Navier-Stokes calculations with a coupled strongly implicit method. Part 2: Spline solutions

A coupled strongly implicit method is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equation. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered. In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗

Continuous analysis of stresses from arbitrary surface loads on a half space

A new form of elemental surface load on a half space is introduced, presuming a quasi-pyramidal variation of load which is doubly linear in each of four rectangular parts of a surface rectangle. Approximations of arbitrary load distributions by sums of such elements are continuous, piecewise linear in two directions and well adaptable. The loads may be normal or tangential. The explicit solutions obtained for all stress and displacement components due to each elemental load involve only elementary functions, are free of the discontinuities which arise with stepwise elements, and are suitable for computing. Some illustrative stress distributions are presented for elemental loads and for multiple pyramidal loads involving both normal and tangential loads. The value of the load continuity in the more complicated analyses of surface cracks is also illustrated.

Bell, J. C.↗

Navier-Stokes calculations with a coupled strongly implicit method. II Spline deferred-corrector solutions

The coupled strongly implicit (CSIP) method described previously is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equations. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered (equal to or less than 1000). In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗

A method for computing the kernel of the downwash integral equation for arbitrary complex frequencies

For the design of active controls to stabilize flight vehicles, which requires the use of unsteady aerodynamics that are valid for arbitrary complex frequencies, algorithms are derived for evaluating the nonelementary part of the kernel of the integral equation that relates unsteady pressure to downwash. This part of the kernel is separated into an infinite limit integral that is evaluated using Bessel and Struve functions and into a finite limit integral that is expanded in series and integrated termwise in closed form. The developed series expansions gave reliable answers for all complex reduced frequencies and executed faster than exponential approximations for many pressure stations.

Desmarais, R. N.↗

Likelihood Analysis for Mega Pixel Maps

The derivation of cosmological parameters from astrophysical data sets routinely involves operations counts which scale as O(N(exp 3) where N is the number of data points. Currently planned missions, including MAP and Planck, will generate sky maps with N(sub d) = 10(exp 6) or more pixels. Simple "brute force" analysis, applied to such mega-pixel data, would require years of computing even on the fastest computers. We describe an algorithm which allows estimation of the likelihood function in the direct pixel basis. The algorithm uses a conjugate gradient approach to evaluate X2 and a geometric approximation to evaluate the determinant. Monte Carlo simulations provide a correction to the determinant, yielding an unbiased estimate of the likelihood surface in an arbitrary region surrounding the likelihood peak. The algorithm requires O(N(sub d)(exp 3/2) operations and O(Nd) storage for each likelihood evaluation, and allows for significant parallel computation.

Kogut, Alan J.↗

Notes on SAW Tag Interrogation Techniques

We consider the problem of interrogating a single SAW RFID tag with a known ID and known range in the presence of multiple interfering tags under the following assumptions: (1) The RF propagation environment is well approximated as a simple delay channel with geometric power-decay constant alpha >/= 2. (2) The interfering tag IDs are unknown but well approximated as independent, identically distributed random samples from a probability distribution of tag ID waveforms with known second-order properties, and the tag of interest is drawn independently from the same distribution. (3) The ranges of the interfering tags are unknown but well approximated as independent, identically distributed realizations of a random variable rho with a known probability distribution f(sub rho) , and the tag ranges are independent of the tag ID waveforms. In particular, we model the tag waveforms as random impulse responses from a wide-sense-stationary, uncorrelated-scattering (WSSUS) fading channel with known bandwidth and scattering function. A brief discussion of the properties of such channels and the notation used to describe them in this document is given in the Appendix. Under these assumptions, we derive the expression for the output signal-to-noise ratio (SNR) for an arbitrary combination of transmitted interrogation signal and linear receiver filter. Based on this expression, we derive the optimal interrogator configuration (i.e., transmitted signal/receiver filter combination) in the two extreme noise/interference regimes, i.e., noise-limited and interference-limited, under the additional assumption that the coherence bandwidth of the tags is much smaller than the total tag bandwidth. Finally, we evaluate the performance of both optimal interrogators over a broad range of operating scenarios using both numerical simulation based on the assumed model and Monte Carlo simulation based on a small sample of measured tag waveforms. The performance evaluation results not only provide guidelines for proper interrogator design, but also provide some insight on the validity of the assumed signal model. It should be noted that the assumption that the impulse response of the tag of interest is known precisely implies that the temperature and range of the tag are also known precisely, which is generally not the case in practice. However, analyzing interrogator performance under this simplifying assumption is much more straightforward and still provides a great deal of insight into the nature of the problem.

Barton, Richard J.↗

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

The NASA Multiscale Analysis Tool (NASMAT) was developed recently to allow a wide variety of multiscale analysis problems to be effectively and efficiently solved. The architecture of NASMAT was established specifically to enable parallelized, “plug-and-play” functionality to reduce the complexity associated with adding new features to the code in the future and to allow end users to rapidly implement and evaluate user-defined capabilities. Additionally, the tool utilizes recursive data structures and subroutines to allow for an arbitrary number of length scales when performing multiscale analyses of heterogeneous materials. These features permit the rapid integration of user-defined capabilities (e.g., a material model, micromechanics approach, or failure theory) at all stages within a NASMAT calculation while leveraging built-in techniques where needed. Additionally, these features allow NASMAT to both be called from an external program as well as call an external program. This paper specifically focuses on the multiscale integration and interoperability of NASMAT with other analysis techniques through an illustrative, multiscale analysis of a 3D woven polymer matrix composite (PMC).

NASMAT↗

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

Multiple relaxation and inhomogeneous broadening in resonance enhanced Raman scattering - Application to tunable infrared generation

The solutions for the imaginary susceptibility of the Raman field transition with arbitrary relaxation rates and field strengths are examined for differing sets of relaxation rates with emphasis on alkali metal vapors which have spontaneous emission dominated relaxation. The model is further expanded to include Doppler broadening and used to predict the peak gain as a function of detuning for a frequency doubled alexandrite laser-pumped cesium vapor gain cell.

Ryan, J. C.↗