Engineering Papers⌕ Search

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.

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

At least 163 records · Page 9

Almost periodic solutions to difference equations

The theory of Massera and Schaeffer relating the existence of unique almost periodic solutions of an inhomogeneous linear equation to an exponential dichotomy for the homogeneous equation was completely extended to discretizations by a strongly stable difference scheme. In addition it is shown that the almost periodic sequence solution will converge to the differential equation solution. The preceding theory was applied to a class of exponentially stable partial differential equations to which one can apply the Hille-Yoshida theorem. It is possible to prove the existence of unique almost periodic solutions of the inhomogeneous equation (which can be approximated by almost periodic sequences) which are the solutions to appropriate discretizations. Two methods of discretizations are discussed: the strongly stable scheme and the Lax-Wendroff scheme.

Bayliss, A.↗

An improved solution of the simplified two-phase plume

An exact solution of the two-phase, gas-particle, solid propellant rocket exhaust plume flowfield is presented. The good agreement between the exact solution and the experimental data proves that the accuracy of the present technique is acceptable for, at least, preliminary design calculations. A good comparison between the plume impingement effects predicted by the exact solution and those by the best computer code available has also been found for various propellant formulations and rocket sizes. This shows that the present exact solution does describe the basic characteristics of the two-phase plume. In case of a parametric analysis for obtaining an optimum design, the present exact solution will be a more convenient tool than any one of the numerical techniques especially when many parameters are involved in the optimization study. Several typical results are shown herein and their physical significance is discussed.

Wang, S. Y.↗

The influence of the computational mesh on accuracy for initial value problems with discontinuous or nonunique solutions

Discontinuous, or weak, solutions of the wave equation, the inviscid form of Burgers equation, and the time-dependent, two-dimensional Euler equations are studied. A numerical method of second-order accuracy in two forms, differential and integral, is used to calculate the weak solutions of these equations for several initial value problems, including supersonic flow past a wedge, a double symmetric wedge, and a sphere. The effect of the computational mesh on the accuracy of computed weak solutions including shock waves and expansion phenomena is studied. Modifications to the finite-difference method are presented which aid in obtaining desired solutions for initial value problems in which the solutions are nonunique.

Maccormack, R. W.↗

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Thompson, J. F.↗

Inverse solutions for laminar boundary-layer flows with separation and reattachment

Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.

Carter, J. E.↗

Supersonic flow of chemically reacting gas-particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.

Penny, M. M.↗

Solution of linear systems by a singular perturbation technique

An approximate solution is obtained for a singularly perturbed system of initial valued, time invariant, linear differential equations with multiple boundary layers. Conditions are stated under which the approximate solution converges uniformly to the exact solution as the perturbation parameter tends to zero. The solution is obtained by the method of matched asymptotic expansions. Use of the results for obtaining approximate solutions of general linear systems is discussed. An example is considered to illustrate the method and it is shown that the formulas derived give a readily computed uniform approximation.

Ardema, M. D.↗

Turbulent viscous-shock-layer solutions with strong vorticity interaction

Numerical solutions of the viscous-shock-layer equations governing laminar and turbulent flows of a perfect gas and radiating and nonradiating mixtures of perfect gases in chemical equilibrium are presented for hypersonic flow over spherically blunted cones and hyperboloids. The results are compared with boundary-layer and inviscid flow-field solutions. The agreement with the inviscid flow field data is satisfactory. The agreement with boundary-layer solutions is good except in regions of strong vorticity interaction. In these flow regions, the viscous-shock-layer solutions appear to be more satisfactory than the boundary-layer solutions.

Anderson, E. C.↗

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗

Analytical Solutions for Non-Linear Differential Equations with the Help of a Digital Computer

A technique was developed with the help of a digital computer for analytic (algebraic) solutions of autonomous and nonautonomous equations. Two operational transform techniques have been programmed for the solution of these equations. Only relatively simple nonlinear differential equations have been considered. In the cases considered it has been possible to assimilate the secular terms into the solutions. For cases where f(t) is not a bounded function, a direct series solution is developed which can be shown to be an analytic function. All solutions have been checked against results obtained by numerical integration for given initial conditions and constants. It is evident that certain nonlinear differential equations can be solved with the help of a digital computer.

Cromwell, P. C.↗

Computer simulation of concentrated solid solution strengthening

The interaction forces between a straight edge dislocation moving through a three-dimensional block containing a random array of solute atoms were determined. The yield stress at 0 K was obtained by determining the average maximum solute-dislocation interaction force that is encountered by edge dislocation, and an expression relating the yield stress to the length of the dislocation and the solute concentration is provided. The magnitude of the solid solution strengthening due to solute atoms can be determined directly from the numerical results, provided the dislocation line length that moves as a unit is specified.

Kuo, C. T. K.↗

Computer simulation of the motion of a straight dislocation line in concentrated solid solutions. II

An investigation was undertaken to determine if the size and modulus interaction of a solute atom with a screw dislocation and the modulus interaction with an edge dislocation contributed to strengthening, in addition to the size interaction with an edge dislocation. The results indicate that the size interaction between solute atom and an edge dislocation accounts for most of the solid solution strengthening in f.c.c. alloys. The contribution to the yield stress from the modulus interaction with an edge dislocation is less than 15%. The interaction between a solute atom and a screw dislocation is much less than that between a solute atom and an edge dislocation.

Kuo, C. T. K.↗

Direct analytical solutions to non-uniform beam problems

The direct analytical solution to the vibration of non-uniform beams with and without discontinuities and with various boundary conditions is presented. Results are compared to results from the exact solution for certain cases where the exact solution has been obtained. It is shown that the direct solution converges to the exact solution, in fact, with 'indefinite accuracy' just as Hamilton stated that it would.

Bailey, C. D.↗

A high order accurate numerical solution algorithm for turbulent boundary layer flow

A fourth-order accurate numerical solution algorithm is derived using finite element interpolation theory for the non-linear parabolic equations governing turbulent boundary layer flow including a two-equation turbulence closure model. The results of carefully controlled numerical experiments firmly quantize for the first time performance differences between finite element and finite difference solution methodology for this type of equation. The developed algorithm takes advantage of the apparent semi-analytical formulational procedure, in establishment of a single, retarded-evaluation Jacobian matrix iterative solution algorithm. Numerical results document performance of solution economy features in terms of computer requirements and solution accuracy. The developed algorithm should find wide application in aerodynamics analysis.

Soliman, M. O.↗

Combination solution for gravity field including altimetry

A previous global solution for the earth's gravitational field and station positions employed a combination of 21 weeks of Doppler tracking data from seven satellites, a complete set of 5-deg mean gravity anomalies, and geometric constraints from the worldwide satellite triangulation network on 17 of the stations. To this is now added geoid heights over the major part of the oceans obtained from GEOS 3 altimeter observations. The solution has a form based on the simple density-layer method and is equivalent to a solution in spherical harmonics to the tenth degree and order. The rms difference of 5-deg mean geoid heights of the new solution with the GEM 10 solution is 5.6 m.

Chovitz, B. H.↗

Laminar and turbulent flow solutions with radiation and ablation injection for Jovian entry

Laminar and turbulent flow-field solutions with coupled carbon-phenolic mass injection are presented for the forebody of a probe entering a nominal Jupiter atmosphere. Solutions are obtained for a 35-degree hyperboloid and for a 45-degree spherically blunted cone using a time-dependent, finite-difference method. The radiative heating rates for the coupled laminar flow are significantly reduced as compared to the corresponding no-blowing case; however, for the coupled turbulent flow, it is found that the surface radiative heating rates are substantially increased and often exceed the corresponding no-blowing values. Turbulence is found to have no effect on the surface radiative heating rates for the no-blowing solutions. The present results are compared with the other available solutions, and some additional solutions are presented.

Kumar, A.↗

Laminar and turbulent flow solutions with radiation and ablation injection for Jovian entry

Laminar and turbulent flowfield solutions with coupled carbon-phenolic mass injection are presented for the forebody of a probe entering a nominal Jupiter atmosphere. Solutions are obtained for a 35-deg hyperboloid and for a 45-deg spherically blunted cone using a time-dependent finite-difference method. The radiative heating rates for the coupled laminar flow are significantly reduced as compared to the corresponding no-blowing case; however, for the coupled turbulent flow, it is found that the surface radiative heating rates are substantially increased and often exceed the corresponding no-blowing values. Turbulence is found to have no effect on the surface radiative heating rates for the no-blowing solutions. The present results are compared with the other available solutions and some additional solutions are presented.

Kumar, A.↗

Closed-form solutions of performability

Methods which yield closed form performability solutions for continuous valued variables are developed. The models are similar to those employed in performance modeling (i.e., Markovian queueing models) but are extended so as to account for variations in structure due to faults. In particular, the modeling of a degradable buffer/multiprocessor system is considered whose performance Y is the (normalized) average throughput rate realized during a bounded interval of time. To avoid known difficulties associated with exact transient solutions, an approximate decomposition of the model is employed permitting certain submodels to be solved in equilibrium. These solutions are then incorporated in a model with fewer transient states and by solving the latter, a closed form solution of the system's performability is obtained. In conclusion, some applications of this solution are discussed and illustrated, including an example of design optimization.

Meyer, J. F.↗