Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Finite element analysis of instability related delamination growth

Instability-related growth of through-width delaminations in laminated coupons was investigated analytically and experimentally. In the analytical effort, a geometrically nonlinear finite-element stress analysis was used, and the numerical results were compared with exact solutions. In addition, measured lateral deflections of postbuckled through-width delaminations were compared with predicted deflections. Lateral deflections, stress distributions, and strain energy release rates were calculated for various delamination lengths, delamination depths, and applied loads. Delamination growth rates observed for unidirectional graphite/epoxy bonded to aluminum were qualitatively correlated with calculated strain-energy release rates and it was found that the Mode I component of energy release dominates the growth process.

Whitcomb, J. D.↗
NASA NTRS1981

On the formation of protostellar disks

An analysis is presented of the hydrodynamic aspects of the growth of protostellar disks from the accretion of a rotating gas cloud. Streamline equations for the flow onto the disk are derived, including the special case where the disk is very thin. For the radial structure of the disk, detailed equations of conservation of mass, momentum, angular momentum, and energy are given. Exact solutions are found for the case of no internal viscosity. For the more general case where there is viscous stress, hypotheses are formulated with the objective of discovering the main qualitative characteristics of a turbulent system and to estimate the magnitude of the turbulent stress. Solutions are also found for boundary conditions, and the large and small viscosity limits. The relation of energy balance to disk thermal structure is detailed. Finally, trends in the collapse of the dust cloud are discussed in the light of the solutions given.

Cassen, P.↗
NASA NTRS1982

Criteria for the use of regression analysis for remote sensing of sediment and pollutants

An examination of limitations, requirements, and precision of the linear multiple-regression technique for quantification of marine environmental parameters is conducted. Both environmental and optical physics conditions have been defined for which an exact solution to the signal response equations is of the same form as the multiple regression equation. Various statistical parameters are examined to define a criteria for selection of an unbiased fit when upwelled radiance values contain error and are correlated with each other. Field experimental data are examined to define data smoothing requirements in order to satisfy the criteria of Daniel and Wood (1971). Recommendations are made concerning improved selection of ground-truth locations to maximize variance and to minimize physical errors associated with the remote sensing experiment.

Whitlock, C. H.↗
NASA NTRS1983

The transmission or scattering of elastic waves by an inhomogeneity of simple geometry: A comparison of theories

The extended method of equivalent inclusions is applied to study the specific wave problems: (1) the transmission of elastic waves in an infinite medium containing a layer of inhomogeneity, and (2) the scattering of elastic waves in an infinite medium containing a perfect spherical inhomogeneity. Eigenstrains are expanded as a geometric series and a method of integration based on the inhomogeneous Helmholtz operator is adopted. This study compares results, obtained by using limited number of terms in the eigenstrain expansion, with exact solutions for the layer problem and that for a perfect sphere.

Sheu, Y. C.↗
NASA NTRS1983

Calculation of potential flow past non-lifting bodies at angle of attack using axial and surface singularity methods

Two different singularity methods have been utilized to calculate the potential flow past a three dimensional non-lifting body. Two separate FORTRAN computer programs have been developed to implement these theoretical models, which will in the future allow inclusion of the fuselage effect in a pair of existing subcritical wing design computer programs. The first method uses higher order axial singularity distributions to model axisymmetric bodies of revolution in an either axial or inclined uniform potential flow. Use of inset of the singularity line away from the body for blunt noses, and cosine-type element distributions have been applied to obtain the optimal results. Excellent agreement to five significant figures with the exact solution pressure coefficient value has been found for a series of ellipsoids at different angles of attack. Solutions obtained for other axisymmetric bodies compare well with available experimental data. The second method utilizes distributions of singularities on the body surface, in the form of a discrete vortex lattice. This program is capable of modeling arbitrary three dimensional non-lifting bodies. Much effort has been devoted to finding the optimal method of calculating the tangential velocity on the body surface, extending techniques previously developed by other workers.

Shu, J. Y.↗
NASA NTRS1982

Numerical evaluation of the Rayleigh integral for planar radiators using the FFT

Rayleigh's integral formula is evaluated numerically for planar radiators of any shape, with any specified velocity in the source plane using the fast Fourier transfrom algorithm. The major advantage of this technique is its speed of computation - over 400 times faster than a straightforward two-dimensional numerical integration. The technique is developed for computation of the radiated pressure in the nearfield of the source and can be easily extended to provide, with little computation time, the vector intensity in the nearfield. Computations with the FFT of the nearfield pressure of baffled rectangular plates with clamped and free boundaries are compared with the 'exact' solution to illuminate any errors. The bias errors, introduced by the FFT, are investigated and a technique is developed to significantly reduce them.

Williams, E. G.↗
NASA NTRS1982

Potentials on large spacecraft in LEO

A system-oriented computer code is used to predict surface charging due to voltages generated within a satellite operating in the typical dense plasma environment of LEO. The use of this code is demonstrated by predicting the expansion of electric fields onto a kapton surface from a pinhole over a biased conductor in a LEO environment. The results are compared to a more-exact solution and experimental data.

Mandell, M. J.↗
NASA NTRS1983

Reduction of rms-error in shallow faceted large space antennas

This paper examines the potential for reducing root-mean-square surface error in shallow faceted reflectors by replacing flat facets with membrane facets. Exact solutions are obtained for the small lateral deflections of equilateral triangular and rectangular membranes subject to isotropic tension and parabolic edge deflections. These solutions are used to minimize the rms-error between a facet of a shallow paraboloidal surface and its approximating membrane facet. The resulting optimum placements and edge curvatures yield membrane facets which have significantly lower rms-errors than the corresponding best-fit flat facets. The rms-error reductions are about 55 percent for equilateral triangles and 25 percent to 93 percent for rectangles, depending on aspect ratio. The results suggest that the use of membrane facets conforming to curved structural members could yield reflectors with lower rms-error, or comparable error with larger facets and, hence, fewer structural members.

Fichter, W. B.↗
NASA NTRS1983

Numerical evaluation of the radiation from unbaffled, finite plates using the FFT

An iteration technique is described which numerically evaluates the acoustic pressure and velocity on and near unbaffled, finite, thin plates vibrating in air. The technique is based on Rayleigh's integral formula and its inverse. These formulas are written in their angular spectrum form so that the fast Fourier transform (FFT) algorithm may be used to evaluate them. As an example of the technique the pressure on the surface of a vibrating, unbaffled disk is computed and shown to be in excellent agreement with the exact solution using oblate spheroidal functions. Furthermore, the computed velocity field outside the disk shows the well-known singularity at the rim of the disk. The radiated fields from unbaffled flat sources of any geometry with prescribed surface velocity may be evaluated using this technique. The use of the FFT to perform the integrations in Rayleigh's formulas provides a great savings in computation time compared with standard integration algorithms, especially when an array processor can be used to implement the FFT.

Williams, E. G.↗
NASA NTRS1983

Analysis of coolant entrance boundary shape of porous region to control cooling along exit boundary

A cooled porous region has a plane surface exposed to a specified spatially varying heat flux. The coolant leaves the region through this surface, and it is desired to control the flow distribution to maintain a specified uniform surface temperature. This is accomplished by having the coolant entrance surface shaped to provide in the region the necessary variation of path length and, hence, flow resistance. The surface shape at the coolant entrance is found by solving a Cauchy boundary value problem. An exact solution is obtained that will deal with a wide variety of heating distributions for both two- and three-dimensional shapes.

Siegel, R.↗
NASA NTRS1983

Kinematics of velocity and vorticity correlations in turbulent flow

The kinematic problem of calculating second-order velocity moments from given values of the vorticity covariance is examined. Integral representation formulas for second-order velocity moments in terms of the two-point vorticity correlation tensor are derived. The special relationships existing between velocity moments in isotropic turbulence are expressed in terms of the integral formulas yielding several kinematic constraints on the two-point vorticity correlation tensor in isotropic turbulence. Numerical evaluation of these constraints suggests that a Gaussian curve may be the only form of the longitudinal velocity correlation coefficient which is consistent with the requirement of isotropy. It is shown that if this is the case, then a family of exact solutions to the decay of isotropic turbulence may be obtained which contains Batchelor's final period solution as a special case. In addition, the computed results suggest a method of approximating the integral representation formulas in general turbulent shear flows.

Bernard, P. S.↗
NASA NTRS1982

Turbulent flow simulation - A large eddy simulator's viewpoint

The state of the art of the data base on turbulent structures and numerical turbulence modelling is assessed. The failure of statistical methods to produce accurate results revealed that turbulence is more than simply random noise. Both randomness and coherence may be characteristics of turbulence, which may also feature noise that is a remnant of previous coherence. Five methods are presently used for numerically predicting turbulent flow behavior: correlations, integral methods, the k-epsilon model, and phenomenological and exact solutions to the Navier-Stokes equations. The higher level models include large eddy simulation and full simulation with numerical solutions to the Navier-Stokes equations; the two methods are applicable to experimental, rather than design work.

Ferziger, J. H.↗
NASA NTRS1983

A theory of rotating stall of multistage axial compressors. II - Finite disturbances

A small-disturbance theory of rotating stall in axial compressors is extended to finite amplitude, assuming the compressor characteristic is a parabola over the range of the disturbance. An exact solution is found which requires the operating point to be at the minimum or maximum of the parabola. If the characteristic is flat in a deep-stall regime, the previous harmonic solution applies with neither reverse flow or 'unstalling'. If the characteristic is concave upward in deep stall, the disturbance has a skewed shape, steeper at the stall-zone trailing edge as experiment shows. Propagation speed is only slightly affected by this nonlinearity. Near stall inception, negative curvature in combination with multiple stall zones can limit the nonlinear oscillation, in the manner of 'progressive' stall. If, as seems likely, lag at stall inception is negative (opposite to inertia), propagation speed exceeds 1/2 wheel speed, as experiments suggest.

Moore, F. K.↗
NASA NTRS1984

A three-dimensional nonlinear analysis of cross-ply rectangular composite plates

The results of a three-dimensional, geometrically nonlinear, finite-element analysis of the bending of cross-ply laminated anisotropic composite plates are presented. Individual laminae are assumed to be homogeneous, orthotropic and linearly elastic. A fully three-dimensional isoparametric finite element with eight nodes (i.e. linear element) and 24 degrees of freedom (three displacement components per node) is used to model the laminated plate. The finite element results of the linear analysis are found to agree very well with the exact solutions of cross-ply laminated rectangular plates under sinusoidal loading. The finite element results of the three-dimensional, geometrically nonlinear analysis are compared with those obtained by using a shear deformable, geometrically nonlinear, plate theory. It is found that the deflections predicted by the shear deformable plate theory are in fair agreement with those predicted by three-dimensional elasticity theory; however stresses were found to be not in good agreement.

Kuppusamy, T.↗
NASA NTRS1983

Stability, vibration and passive damping of partially restrained imperfect columns

A theoretical and experimental study of slender tubular columns for possible use in space structures is conducted in the presence of partial rotational end restraints. Explicit formulas are derived for computing the buckling load and the lowest natural frequency of perfectly straight uniform elastic members with rotational end restraints possessing linear moment-rotation characteristics. An exact solution in the form of a transcendental equation, and a numerical solution using second-order finite-differences are also presented. The presence of an initial imperfection is also incorporated into the numerical procedure. Vibration tests are conducted on an imperfect tubular steel member in the absence of an axial load. A damping concept consisting of a string-mass assembly is explored. Three passive damping configurations involving combinations of three lead shots were investigated. The three lead shot configurations provided considerably greater damping than the single lead shot.

Razzaq, Z.↗
NASA NTRS1983

Self-shadowing of orbiting trusses

The approach used to assess shadowing reductions on the heating of orbiting trusses involves determining the heating rates with slender member shadowing effects included and then obtaining the thermal response of the shadowed member. Steps taken to identify shadowers, find locations where shadowing occurs and calculate shadow intensity are listed. The finite element thermal structural analysis of cable stiffened space structure is delineated and the exact solution of the caternary problem is given. Typical cable surface heating rates are plotted. The structural analysis includes large deformation (nonlinear), thermal effect, and the pretension effect. Displacements and stresses are computed at different orbital positions for an orbit.

Mahaney, J.↗
NASA NTRS1984

Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

Pratt, D. T.↗
NASA NTRS1983

An eigenvalue method for solving transient heat conduction problems

The eigenvalue method, which has been used by researchers in structure mechanics, is applied to problems in heat conduction. Its formulation is decribed in terms of an examination of transient heat conduction in a square slab. Taking advantage of the availability of the exact solution, we compare the accuracy and other numerical properties of the eigenvalue method with those of existing numerical schemes. The comparsion shows that, overall, the eigenvalue method appears to be fairly attractive. Furthermore, only a few dominant eigenvalues and their corresponding eigenvectors need to be computed and retained to yield reasonably high accuracy. Greater savings are attained in the computation time for a transient problem with long time duration and a large computational domain.

Shih, T. M.↗