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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

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.↗

NASA NTRS1983

Finite element formulations for acoustical radiation

Finite and infinite element techniques are applied to linear acoustical problems involving infinite anechoic boundaries. Theory is presented for a simple one dimensional model based on Webster's horn equation. Results are then presented both for the one dimensional model and for two axisymmetric test cases. Comparisons with exact solutions indicate that both the infinite element and wave envelope schemes are effective in correctly predicting the near field. The wave envelope scheme is also shown to be capable of resolving the far field radiation pattern.

Astley, R. J.↗

NASA NTRS1984

A theoretical study of the effects of interacting grain boundaries on electron-beam-induced currents

The influence of several grain boundaries on the electron-beam-induced current in polycrystalline silicon is investigated. The simple geometry contemplated here, consisting of two plane parallel grain boundaries possessing different effective, albeit constant, surface recombination velocities, intersecting a planar junction at right angles allows for an exact solution of the corresponding boundary value problem. It is shown that if the distance between the grain boundaries w becomes less than the diffusion length for minority carriers L within the grain, the grain boundary with the smaller effective recombination velocity becomes less and less visible the smaller the ratio w/L.

Luke, K. L.↗

NASA NTRS1984

Propagation and radiation of sound from flanged circular ducts with circumferentially varying wall admittances. I Semi-infinite ducts. II - Finite ducts with sources

Sound propagation in infinite, semiinfinite, and finite circular ducts with circumferentially varying wall admittances is investigated analytically. The infinite case is considered, and an example demonstrates the effects of wall-admittance distribution on dispersion characteristics and mode shapes. An exact solution is obtained for the semiinfinite case, a circular duct with a flanged opening: sidelobe suppression and circumferential-mode energy scattering leading to radiated-field asymmetry are found. A finite duct system with specified hard-walled pressure sources is examined in detail, evaluating reflection coefficients, transmission losses, and radiated-field directivity. Graphs and diagrams are provided, and the implications of the results obtained for the design of aircraft-turbofan inlet liners are discussed.

Fuller, C. R.↗

NASA NTRS1984

Method for calculating multidimensional electric fields in photovoltaic modules

A finite element method for evaluating the electrical isolation characteristics of photovoltaic modules was developed; its accuracy was verified by comparison with an exact solution for a geometry similar to that of solar cells. Tests on a square test coupon, employed in electrical isolation tests, and a group of disc-shaped solar cells illustrated the finite element method's usefulness in evaluating module encapsulation designs. Finite element models had to avoid adjacent large and small elements and elements with large aspect ratios, and the NASTRAN output had to be curve fitted to calculate the maximum field. Geometric limits were indicated: cells with very sharp edges, and cells much thinner or thicker than the dielectric pottant layer.

Kallis, J. M.↗

NASA NTRS1983

A multi-item maintenance center inventory model for low-demand reparable items

In many military and commercial contexts, complex equipment undergoes scheduled maintenance overhauls at regular intervals during which all failed components are replaced. Failure to have replacements on hand for all failed parts requires emergency measures at premium cost. When reparable parts are highly reliable and expensive, both holding and shortage costs are high. This model determines the reparable parts inventory for a maintenance center under three alternative criteria: (1) maximizing job-completion rate subject to constraint on total holding costs, (2) minimizing total holding costs plus expected job noncompletion costs, and (3) minimizing total holding costs subject to a required minimum job-completion rate. Exact solutions may be obtained using dynamic programming. Approximate solutions, found easily by marginal analysis, have readily computed bounds on possible error. The solution methods for the three formulations are illustrated in a simple example.

Schaefer, M. K.↗

NASA NTRS1983

Planar dynamics of a uniform beam with rigid bodies affixed to the ends

The planar dynamics of a uniform elastic beam subject to a variety of geometric and natural boundary conditions and external excitations were analyzed. The beams are inextensible and capable of small transverse bending deformations only. Classical beam vibration eigenvalue problems for a cantilever with tip mass, a cantilever with tip body and an unconstrained beam with rigid bodies at each are examined. The characteristic equations, eigenfunctions and orthogonality relations for each are derived. The forced vibration of a cantilever with tip body subject to base acceleration is analyzed. The exact solution of the governing nonhomogeneous partial differential equation with time dependent boundary conditions is presented and compared with a Rayleigh-Ritz approximate solution. The arbitrary planar motion of an elastic beam with rigid bodies at the ends is addressed. Equations of motion are derived for two modal expansions of the beam deflection. The motion equations are cast in a first order form suitable for numerical integration. Selected FORTRAN programs are provided.

Storch, J.↗

NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗

NASA NTRS1985

Deceleration of a supersonic flow behind a curved shock wave with isentropic precompression

Three-dimensional supersonic flows of an ideal fluid in the neighborhood of bodies formed by being cut out along the streamlines of an axisymmetric flow are investigated. The flow consists of a region of isoentropic compression and a region of vortex flow. An exact solution with variable entropy is used to describe the flow in the vortex region. In the continuous flow region an approximate solution is constructed by expanding the solution in a series in a small parameter. The effect of the shape of the excision and the vorticity of the flow on compression of the jet and and the total pressure loss coefficient is studied.

Dulov, V. G.↗

NASA NTRS1985

A higher-order shear deformation theory of laminated elastic shells

A higher-order shear deformation theory of elastic shells is developed for shells laminated of orthotropic layers. The theory is a modification of the Sanders' theory and accounts for parabolic distribution of the transverse shear strains through thickness of the shell and tangential stress-free boundary conditions on the boundary surfaces of the shell. The Navier-type exact solutions for bending and natural vibration are presented for cylindrical and spherical shells under simply supported boundary conditions.

Reddy, J. N.↗