SEARCH · Engineering Papers

Results for “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 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Numerical solutions for viscous and potential flow about arbitrary two-dimensional bodies using body-fitted coordinate systems

A procedure for numerical solution of the time-dependent, incompressible Navier-Stokes equations for the flow about arbitrarily shaped two-dimensional bodies is given. This solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having a coordinate line coincident with the body contour regardless of its shape. The implicit solution utilizes the vorticity-stream function formulation with a false-position iterative adjustment of the surface vorticity in satisfaction of the no-slip boundary condition. Excellent agreement with the Blasius boundary layer solution is obtained for a semi-infinite flat plate. Results are presented for Reynolds numbers up to 2000 for several airfoils and a cambered rock.

Thames, F. C.↗

NASA NTRS1978

Numerical solution of the Navier-Stokes equations for super-sonic flows with strong shocks

The numerical solution of the full Navier-Stokes Equations for viscous flows with high Mach numbers and a strong detached bow shock was obtained. Two dimensional flows around a circular cylinder, and a circular cylinder with an aft-body in the form of a fairing, were considered. The solution of the compressible N.S. equations was accomplished by the method of finite differences. An implicit scheme of solution, the S.O.R., was used with the optimum acceleration parameters determined by trial and error. The tensor notation was used in writing the N-S Equations transformed into general curvilinear coordinates. The equations for the generation of the coordinate system were solved, followed by the solution of the N.S. equations, at the end of a set of given number of time steps. "Wiggles", constituted the one major problem that needed to be overcome. These oscillations give rise to quantities such as negative temperatures, which ultimately caused the computational program to break down. Certain dissipative finite-difference schemes damped these oscillations.

Devarayalu, K.↗

NASA NTRS1977

Integration-free interval doubling for Riccati equation solutions

Starting with certain identities obtained by Reid (1972) and Redheffer (1962) for general matrix Riccati equations (RE's), we give various algorithms for the case of constant coefficients. The algorithms are based on two ideas - first, relate the RE solution with general initial conditions to anchored RE solutions; and second, when the coefficients are constant, the anchored solutions have a basic shift-invariance property. These ideas are used to construct an integration-free, superlinearly convergent iterative solution to the algebraic RE. Preliminary numerical experiments show that our algorithms, arranged in square-root form, provide a method that is numerically stable and appears to be competitive with other methods of solving the algebraic RE.

Sidhu, G. S.↗

NASA NTRS1977

Solution of transient optimization problems by using an algorithm based on nonlinear programming

A new algorithm is presented for solution of dynamic optimization problems which are nonlinear in the state variables and linear in the control variables. It is shown that the optimal control is bang-bang. A nominal bang-bang solution is found which satisfies the system equations and constraints, and influence functions are generated which check the optimality of the solution. Nonlinear optimization (gradient search) techniques are used to find the optimal solution. The algorithm is used to find a minimum time acceleration for a turbofan engine.

Teren, F.↗

NASA NTRS1977

Recursive analytical solution describing artificial satellite motion perturbed by an arbitrary number of zonal terms

An analytical first order solution has been developed which describes the motion of an artificial satellite perturbed by an arbitrary number of zonal harmonics of the geopotential. A set of recursive relations for the solution, which was deduced from recursive relations of the geopotential, was derived. The method of solution is based on Von-Zeipel's technique applied to a canonical set of two-body elements in the extended phase space which incorporates the true anomaly as a canonical element. The elements are of Poincare type, that is, they are regular for vanishing eccentricities and inclinations. Numerical results show that this solution is accurate to within a few meters after 500 revolutions.

Mueller, A. C.↗

NASA NTRS1978

A discussion of the solution for the motion of Pluto

The semianalytical solution for Pluto given by Nacozy and Diehl (1978) is compared with the numerical solution obtained by Williams and Benson (1971). The effect of the Pluto-Uranus near resonance is discussed along with how it can be incorporated into the solution. Details are given on the calculation of certain quantities in the long-term solution. A calculation of librational periods of the perihelion, eccentricity, and inclination is performed. The calculation is based on the method of Giacaglia and Hori (1968).

Nacozy, P. E.↗

NASA NTRS1980

Newton's method: A link between continuous and discrete solutions of nonlinear problems

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

Thurston, G. A.↗

NASA NTRS1980

Iterative solutions for one-dimensional diffusion with time varying surface composition and composition-dependent diffusion coefficient

Solutions are given for one-dimensional diffusion problems with a time varying surface composition and also a composition dependent diffusion coefficient. The most general solution does not require special mathematical functions to fit the variation in surface composition or D(C). In another solution, a series expansion may be used to fit the time dependent surface concentration. These solutions make use of iterative calculations that converge rapidly and are highly stable. Computer times are much shorter than that required for finite difference calculations and can efficiently make use of interactive graphics terminals. Existing gas carburization data were used to provide an illustration of an iterative approach with a time varying carbon composition at the free surface.

Chow, M.↗

NASA NTRS1980

A rapid method for the approximate determination of nonlinear solutions Application to aerodynamic flows

A method for determining highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameters, is developed and evaluated. The procedure consists of defining a unit perturbation by employing two or more nonlinear solutions which differ from one another by a nominal change in some geometric or flow parameter, and then using that unit perturbation to predict a family of related nonlinear solutions over a range of parameter variation. Coordinate straining is used in determining the unit perturbation to account for the movement of discontinuities and maxima of high-gradient regions due to the perturbation. Although the procedure is generally applicable, results are presented here for nonlinear aerodynamic applications. Attention is focused in particular on transonic flows which are strongly supercritical and exhibit large surface shock movement over the parametric range studied; and on subsonic flows which display large pressure variations in the stagnation and peak suction pressure regions. Flows past both isolated airfoils and compressor cascades involving a variety of flow and geometry parameter changes are considered. Comparisons with the corresponding 'exact' nonlinear solutions indicate a remarkable accuracy and range of validity of such a procedure. Computational time is trivial.

Stahara, S. S.↗

NASA NTRS1980

The solution of variable-step implicit difference equations for dynamic systems analysis

A semi-implicit iterative technique is given for the solution of implicit matrix difference equations which arise in variable-step direct time integration of dynamic systems equations. For all practical purposes, the technique avoids repeated refactoring of the stepsize-dependent solution matrix and achieves the solution convergence for engineering accuracy requirements within a permissible number of iterations. Extension of the present method to nonlinear problems treated by the pseudo-force procedure is straightforward. For nonlinear solution procedures based on the tangent modulus formulation, the technique requires modification to make use of the tangent stiffness matrix. It is noted that the question of whether to include the nonlinear terms within the iteration loop has not yet been resolved.

Park, K. C.↗

NASA NTRS1981

Solutions of the Helmholtz equation with boundary conditions for force-free magnetic fields

It is shown that the solution, with one ignorable coordinate, for the Taylor minimum energy state (resulting in a force-free magnetic field) in either a straight cylindrical or a toroidal geometry with arbitrary cross section can be reduced to the solution of either an inhomogeneous Helmholtz equation or a Grad-Shafranov equation with simple boundary conditions. Standard Green's function theory is, therefore, applicable. Detailed solutions are presented for the Taylor state in toroidal and cylindrical domains having a rectangular cross section. The focus is on solutions corresponding to the continuous eigenvalue spectra. Singular behavior at 90 deg corners is explored in detail.

Rasband, S. N.↗

NASA NTRS1981

Effect of zone size on the convergence of exact solutions for diffusion in single phase systems with planar, cylindrical or spherical geometry

Exact solutions for diffusion in single phase binary alloy systems with constant diffusion coefficient and zero-flux boundary condition have been evaluated to establish the optimum zone size of applicability. Planar, cylindrical and spherical interface geometry, and finite, singly infinite, and doubly infinite systems are treated. Two solutions are presented for each geometry, one well suited to short diffusion times, and one to long times. The effect of zone-size on the convergence of these solutions is discussed. A generalized form of the diffusion solution for doubly infinite systems is proposed.

Unnam, J.↗

NASA NTRS1982

Numerical solution of potential flow about arbitrary 2-dimensional multiple bodies

A procedure for the finite-difference numerical solution of the lifting potential flow about any number of arbitrarily shaped bodies is given. The solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having coordinate lines coincident with the contours of all bodies in the field, regardless of their shapes and number. The effects of all numerical parameters involved are analyzed and appropriate values are recommended. Comparisons with analytic solutions for single Karman-Trefftz airfoils and a circular cylinder pair show excellent agreement. The technique of application of the boundary-fitted coordinate systems to the numerical solution of partial differential equations is illustrated.

Thompson, J. F.↗

NASA NTRS1982

Relaxation solution of the full Euler equations

A numerical procedure for the relaxation solution of the full steady Euler equations is described. By embedding the Euler system in a second order surrogate system, central differencing may be used in subsonic regions while retaining matrix forms well suited to iterative solution procedures and convergence acceleration techniques. Hence, this method allows the development of stable, fully conservative differencing schemes for the solution of quite general inviscid flow problems. Results are presented for both subcritical and shocked supercritical internal flows. Comparisons are made with a standard time dependent solution algorithm.

Johnson, G. M.↗

NASA NTRS1982

Influence of corrosive solutions on microhardness and chemistry of magnesium oxide /001/ surfaces

X-ray photoelectron spectroscopy analyses and hardness experiments were conducted on cleaved magnesium oxide /001/ surfaces. The magnesium oxide bulk crystals were cleaved to specimen size along the /001/ surface, and indentations were made on the cleaved surface in corrosive solutions containing HCl, NaOH, or HNO3 and in water without exposing the specimen to any other environment. The results indicated that chloride (such as MgCl2) and sodium films are formed on the magnesium oxide surface as a result of interactions between an HCl-containing solution and a cleaved magnesium oxide surface. The chloride films soften the magnesium oxide surface. In this case microhardness is strongly influenced by the pH value of the solution. The lower the pH, the lower the microhardness. Sodium films, which are formed on the magnesium oxide surface exposed to an NaOH containing solution, do not soften the magnesium oxide surface.

Ishigaki, H.↗

NASA NTRS1982

Integrating matrix solution of the hybrid state vector equations for beam vibration

A simple, versatile, and efficient computational technique has been developed for dynamic analysis of linear elastic beam and rod type of structures. Moreover, the method provides a rather general solution approach for two-point boundary value problems that are described by a single independent spatial variable. For structural problems, the method is implemented by a mixed state vector formulation of the differential equations, combined with an integrating matrix solution procedure. Highly accurate solutions are easily achieved with this approach. Example solutions are given for beam vibration problems including discontinuous stiffness and mass parameters, elastic restraint boundary conditions, concentrated inertia loading, and rigid body modes

Lehman, L. L.↗

NASA NTRS1982

Solution to the backward-Kolmogorov equation for a nonstationary oscillation problem

The transition probability density function of a Markovian approximation of the response amplitude of an oscillator under nonstationary excitation is determined in an analytical form. A solution is presented for the associated backward-Kolmogorov equation by transforming the equation into a form amenable to solution by the method of the separation of variables. This procedure results in a boundary value problem which is then solved by using an infinite series of Laguerre polynomials. It is found that the infinite series solution is equivalent to a closed-form solution involving a Bessel function.

Solomos, G. P.↗

NASA NTRS1982

Time-dependent solutions of multimode convection equations

Truncated modal equations are used to study the time evolution of thermal convection. In the Boussinesq approximation these nonlinear equations are obtained by expanding the fluctuating velocity and temperature fields in a finite set of planforms of the horizontal coordinates. Numerical studies dealing with two or three modes with triad interactions are discussed. Rich time dependence was found in these cases: periodic and aperiodic solutions can be obtained, along with various steady solutions. Three-mode solutions reproduce the qualitative appearance of spoke-pattern convection as observed in experiments at high Prandtl numbers. Though the values of the periods of the time-dependent solutions do not agree with those of the experiments, their variation with Rayleigh number compares favorably. Except at the highest Rayleigh number considered (10,000,000), the theoretical Nusselt numbers agree well with experiment.

Toomre, J.↗