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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗

NASA NTRS1985

Finite element analysis of partly wrinkled membranes

The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.

Miller, R. K.↗

NASA NTRS1991

An exact algebraic solution of the infimum in H-infinity optimization with output feedback

This paper presents a simple and noniterative procedure for the computation of the exact value of the infimum in the standard H-infinity-optimal control with output feedback. The problem formulation is general and does not place any restrictions on the direct feedthrough terms between the control input and the controlled output variables, and between the disturbance input and the measurement output variables. The method is applicable to systems that satisfy (1) the transfer function from the control input to the controlled output is right-invertible and has no invariant zeros on the j(w) axis and, (2) the transfer function from the disturbance to the measurement output is left-invertible and has no invariant zeros on the j(w) axis. A set of necessary and sufficient conditions for the solvability of H-infinity-almost disturbance decoupling problem via measurement feedback with internal stability is also given.

Chen, Ben M.↗

NASA NTRS1984

Alfven waves in a thermally stratified fluid

The properties of Alfven waves propagating along a uniform horizontal field in a highly conducting incompressible medium in the presence of strong convective instability are examined in the Boussinesq approximation. In particular, it is sought to determine whether there are exact solutions to the dynamical equations in the presence of convective forces. It is shown that a class of exact solutions of arbitrary amplitude, but of limited form, which may be of some physical interest, does exist. For large amplitudes, any mixtures of polarization states are shown to cause scattering into new modes.

Parker, E. N.↗

NASA NTRS1984

A refined nonlinear theory of plates with transverse shear deformation

A higher-order shear deformation theory of plates accounting for the von Karman strains is presented. The theory contains the same dependent unknowns as in the Hencky-Mindlin type first-order shear deformation theory and accounts for parabolic distribution of the transverse shear strains through the thickness of the plate. Exact solutions of simply supported plates are obtained using the linear theory and the results are compared with the exact solutions of three-dimensional elasticity theory, the first order shear deformation theory, and the classical plate theory. The present theory predicts the deflections, stresses, and frequencies more accurately when compared to the first-order theory and the classical plate theory.

Reddy, J. N.↗

NASA NTRS1991

Analytical solution for boundary heat fluxes from a radiating rectangular medium

Reference is made to the work of Shah (1979) which demonstrated the possibility of partially integrating the radiative equations analytically to obtain an 'exact' solution. Shah's solution was given as a double integration of the modified Bessel function of order zero. Here, it is shown that the 'exact' solution for a rectangular region radiating to cold black walls can be conveniently derived, and expressed in simple form, by using an integral function, Sn, analogous to the exponential integral function appearing in plane-layer solutions.

Siegel, R.↗

NASA NTRS1998

Robust Multigrid Smoothers for Three Dimensional Elliptic Equations with Strong Anisotropies

We discuss the behavior of several plane relaxation methods as multigrid smoothers for the solution of a discrete anisotropic elliptic model problem on cell-centered grids. The methods compared are plane Jacobi with damping, plane Jacobi with partial damping, plane Gauss-Seidel, plane zebra Gauss-Seidel, and line Gauss-Seidel. Based on numerical experiments and local mode analysis, we compare the smoothing factor of the different methods in the presence of strong anisotropies. A four-color Gauss-Seidel method is found to have the best numerical and architectural properties of the methods considered in the present work. Although alternating direction plane relaxation schemes are simpler and more robust than other approaches, they are not currently used in industrial and production codes because they require the solution of a two-dimensional problem for each plane in each direction. We verify the theoretical predictions of Thole and Trottenberg that an exact solution of each plane is not necessary and that a single two-dimensional multigrid cycle gives the same result as an exact solution, in much less execution time. Parallelization of the two-dimensional multigrid cycles, the kernel of the three-dimensional implicit solver, is also discussed. Alternating-plane smoothers are found to be highly efficient multigrid smoothers for anisotropic elliptic problems.

Llorente, Ignacio M.↗

NASA NTRS1982

Three-dimensional finite-element analysis of layered composite plates

Results are presented for an investigation of the three-dimensional, geometrically nonlinear, finite-element analysis of the bending of laminated anisotropic composite plates. The individual laminae are treated as homogeneous, transversely isotropic, and linearly elastic. A fully three-dimensional isoparametric finite element with eight modes (i.e., linear element) and 24 degrees of freedom (three displacement components per node) is used. The numerical results obtained using this linear analysis are compared with the exact solutions given in Pagano (1969, 1970). It is found that the results of the linear analysis converge to the exact solution as the mesh is refined.

Putcha, N. S.↗

NASA NTRS1980

Two-stream approximations to radiative transfer in planetary atmospheres - A unified description of existing methods and a new improvement

Existing two-stream approximations to radiative transfer theory for particulate media are shown to be represented by identical forms of coupled differential equations if the intensity is replaced by integrals of the intensity over hemispheres. One set of solutions thus suffices for all methods and provides convenient analytical comparisons. The equations also suggest modifications of the standard techniques so as to duplicate exact solutions for thin atmospheres and thus permit accurate determinations of the effects of typical aerosol layers. Numerical results for the plane albedos of plane-parallel atmospheres are given for conventional and modified Eddington approximations, conventional and modified two-point quadrature schemes, the hemispheric-constant method and the delta-function method, all for comparison with accurate discrete-ordinate solutions. A new two-stream approximation is introduced that reduces to the modified Eddington approximation in the limit of isotropic phase functions and to the exact solution in the limit of extreme anisotropic scattering. Comparisons of plane albedos and transmittances show the new method to be generally superior over a wide range of atmospheric conditions (including cloud and aerosol layers), especially in the case of nonconservative scattering.

Meador, W. E.↗

NASA NTRS1976

Use of variational methods in the determination of wind-driven ocean circulation

Simple polynomial approximations and a variational approach were used to predict wind-induced circulation in rectangular ocean basins. Stommel's and Munk's models were solved in a unified fashion by means of the proposed method. Very good agreement with exact solutions available in the literature was shown to exist. The method was then applied to more complex situations where an exact solution seems out of the question.

Gelos, R.↗

NASA NTRS1986

An entropy maximization problem related to optical communication

In relation to a problem in optical communication, the paper considers the general problem of maximizing the entropy of a stationary radom process that is subject to an average transition cost constraint. By using a recent result of Justesen and Hoholdt, an exact solution to the problem is presented and a class of finite state encoders that give a good approximation to the exact solution is suggested.

Mceliece, R. J.↗

NASA NTRS2008

Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis

Two Reynolds-averaged Navier-Stokes computer codes - one unstructured and one structured - are applied to two workshop cases (for the 3rd Workshop on CFD Uncertainty Analysis, held at Instituto Superior Tecnico, Lisbon, in October 2008) for the purpose of uncertainty analysis. The Spalart-Allmaras turbulence model is employed. The first case uses the method of manufactured solution and is intended as a verification case. In other words, the CFD solution is expected to approach the exact solution as the grid is refined. The second case is a validation case (comparison against experiment), for which modeling errors inherent in the turbulence model and errors/uncertainty in the experiment may prevent close agreement. The results from the two computer codes are also compared. This exercise verifies that the codes are consistent both with the exact manufactured solution and with each other. In terms of order property, both codes behave as expected for the manufactured solution. For the backward facing step, CFD uncertainty on the finest grid is computed and is generally very low for both codes (whose results are nearly identical). Agreement with experiment is good at some locations for particular variables, but there are also many areas where the CFD and experimental uncertainties do not overlap.

Rumsey, C. L.↗

NASA NTRS1973

Auroral ion velocity distributions using a relaxation model.

Calculation of ion velocity distributions for a weakly-ionized plasma subjected to crossed electric and magnetic fields for application to the auroral ionosphere. By replacing the Boltzmann collision integral with a simple relaxation model, an exact solution to Boltzmann's equation could be obtained. This solution has the advantage over a series expansion in that all the higher-order velocity moments are inherent in it. The exact solution is particularly advantageous when studying large departures of the distribution from its Maxwellian form, because these departures are caused by the higher velocity moments. In general, however, a simple relaxation model can only be used to obtain qualitative information on the distribution function. Consequently, it is possible to determine when the higher-order velocity moments affect the ion velocity distribution and the nature of their effect, but it is not possible to obtain accurate quantitative results.

St-Maurice, J.-P.↗

NASA NTRS1976

Features of sound propagation through and stability of a finite shear layer

The plane wave propagation, the stability and the rectangular duct mode problems of a compressible inviscid linearly sheared parallel, but otherwise homogeneous flow, are shown to be governed by Whittaker's equation. The exact solutions for the perturbation quantities are essentially Whittaker M-functions. A number of known results are obtained as limiting cases of exact solutions. For the compressible finite thickness shear layer it is shown that no resonances and no critical angles exist for all Mach numbers, frequencies and shear layer velocity profile slopes except in the singular case of the vortex sheet.

Koutsoyannis, S. P.↗

NASA NTRS1994

Parasitic instabilities in magnetized, differentially rotating disks

Velikhov, Chandrasekhar, Balbus & Hawley have discovered a MHD instability of differentially rotating fluids that may explain the effective viscosity of accretion disks. If the unperturbed magnetic energy density is much less than the gas pressure, this `magnetorotational' instability (MRI) arises on small scales and is approximately incompressible. If the unperturbed field is suffiently, strong, the MRI is suppressed. We therefore ask whether the MRI mechanism can amplify even very weak initial fields until they reach equipartition with the gas pressure. We show that in the total incompressible limit, the MRI modes are exact solutions of the nonlinear fluid equations, even if the perturbed magnetic field is much larger than the unperturbed field. Also, we present a new exact solution in the opposite limit that the magnetic pressure is much larger than that of the gas. On the other hand, we show that the incompressible MRI modes are themselves subject to parasitic instabilities with instantaneous growth rates proportional to the MRI amplitude. Some of the parasitic instabilities are of the Kelvin-Helmholtz type, while others are less familiar. The eigenfunctions of the latter group suggest that they may promote rapid reconnection of the MRI field. Thus, parasitic modes may stop MRI growth at subequipartition amplitudes if it developes from a sufficiently weak initital field.

Goodman, Jeremy↗

NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗

NASA NTRS1986

Using Hough harmonics to validate and assess nonlinear shallow-water models

The implementation of a technique for locating programming errors in shallow-water codes, establishing the correctness of the code, and assessing the performance of the numerical model under various flow conditions is described. The right-hand side of the differential equations is modified in such a way that the exact solution of the nonlinear initial-value problem is known, so that the truncation errors of the numerical scheme can be studied in detail. The exact solution is prescribed to be any linear combination of Hough harmonics which propagate in time according to their natural frequencies.

Dee, Dick P.↗

NASA NTRS1978

Characterization of acoustic disturbances in linearly sheared flows

The equation describing the plane wave propagation, the stability, or the rectangular duct mode characteristics in a compressible inviscid linearly sheared parallel, but otherwise homogeneous flow, is shown to be reducible to Whittaker's equation. The resulting solutions, which are real, viewed as functions of two variables, depend on a parameter and an argument, the values of which have precise physical meanings depending on the problem. The exact solutions in terms of Whittaker functions are used to obtain a number of known results of plane wave propagation and stability in linearly sheared flows as limiting cases in which the speed of sound goes to infinity (incompressible limit) or the shear layer thickness, or wave number, goes to zero (vortex sheet limit). The usefulness of the exact solutions is then discussed in connection with the problems of plane wave propagation and the stability of a finite thickness layer with a linear velocity profile.

Koutsoyannis, S. P.↗