Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical solution”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Numerical grid generation; Proceedings of the Symposium on Numerical Generation of Curvilinear Coordinate Systems and Their Use in the Numerical Solution of Partial Differential Equations, Nashville, TN, April 13-16, 1982

General curvilinear coordinate systems are considered along with the error induced by coordinate systems, basic differential models for coordinate generation, elliptic grid generation, conformal grid generation, algebraic grid generation, orthogonal grid generation, patched coordinate systems, and solid mechanics applications of boundary fitted coordinate systems. Attention is given to coordinate system control and adaptive meshes, the application of body conforming curvilinear grids for finite difference solution of external flow, the use of solution adaptive grids in solving partial differential equations, adaptive gridding for finite difference solutions to heat and mass transfer problems, and the application of curvilinear coordinate generation techniques to the computation of internal flows. Other topics explored are related to the solution of nonlinear water wave problems using boundary-fitted coordinate systems, the numerical modeling of estuarine hydrodynamics on a boundary-fitted coordinate system, and conformal grid generation for multielement airfoils.

Thompson, J. F.↗
NASA NTRS1972

Flow field analysis of aircraft configurations using a numerical solution to the three-dimensional unified supersonic/hypersonic small disturbance equations, part 1

The unified small disturbance equations are numerically solved using the well-known Lax-Wendroff finite difference technique. The method allows complete determination of the inviscid flow field and surface properties as long as the flow remains supersonic. Shock waves and other discontinuities are accounted for implicity in the numerical method. This technique was programed for general application to the three-dimensional case. The validity of the method is demonstrated by calculations on cones, axisymmetric bodies, lifting bodies, delta wings, and a conical wing/body combination. Part 1 contains the discussion of problem development and results of the study. Part 2 contains flow charts, subroutine descriptions, and a listing of the computer program.

Gunness, R. C., Jr.↗
NASA NTRS1990

Evolution of midplate hotspot swells: Numerical solutions

The evolution of midplate hotspot swells on an oceanic plate moving over a hot, upwelling mantle plume is numerically simulated. The plume supplies a Gaussian-shaped thermal perturbation and thermally-induced dynamic support. The lithosphere is treated as a thermal boundary layer with a strongly temperature-dependent viscosity. The two fundamental mechanisms of transferring heat, conduction and convection, during the interaction of the lithosphere with the mantle plume are considered. The transient heat transfer equations, with boundary conditions varying in both time and space, are solved in cylindrical coordinates using the finite difference ADI (alternating direction implicit) method on a 100 x 100 grid. The topography, geoid anomaly, and heat flow anomaly of the Hawaiian swell and the Bermuda rise are used to constrain the models. Results confirm the conclusion of previous works that the Hawaiian swell can not be explained by conductive heating alone, even if extremely high thermal perturbation is allowed. On the other hand, the model of convective thinning predicts successfully the topography, geoid anomaly, and the heat flow anomaly around the Hawaiian islands, as well as the changes in the topography and anomalous heat flow along the Hawaiian volcanic chain.

Liu, Mian↗
NASA NTRS1979

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗
NASA NTRS1986

Numerical solutions of three-dimensional Navier-Stokes equations for closed-bluff bodies

Methods of grid adaption are reviewed and a method is developed with the capability of adaption to several flow variables. This method is based on a variational approach and is an algebraic method which does not require the solution of partial differential equations. The method has also been formulated in such a way that there is no need for any matrix inversion. The method is used in conjunction with the calculation of hypersonic flow over a blunt nose. The equations of motion are the compressible Navier-Stokes equations where all viscous terms are retained. They are solved by the MacCormack time-splitting method and a movie has been produced which shows simultaneously the transient behavior of the solution and the grid adaption. The results are compared with the experimental and other numerical results.

Abolhassani, Jamshid S.↗
NASA NTRS1992

Effect of airfoil (trailing-edge) thickness on the numerical solution of panel methods based on the Dirichlet boundary condition

The practical limit of airfoil thickness ratio for which acceptable engineering results are obtainable with the Dirichlet boundary-condition-based numerical methods is investigated. This is done by studying the effect of thickness on the calculated pressure distribution near the trailing edge and by comparing the aerodynamic coefficients with available exact solutions. The first objective of this study, owing to the wide use of such computational methods, is to demonstrate the numerical symptoms that occur when the body or wing thickness approaches zero and to increase the awareness of potential users of these methods. Additionally, an effort is made to obtain the practical limits of the trailing-edge thickness where such problems will appear in the flow solution, and to propose some possible cures for very thin airfoils or those with cusped trailing edges.

Yon, Steven↗
NASA NTRS1977

Numerical solution for the temperature distribution in a cooled guide vane blade of a radial gas turbine

A two dimensional finite difference numerical technique is presented to determine the temperature distribution of an internal cooled blade of radial turbine guide vanes. A simple convection cooling is assumed inside the guide vane blade. Such cooling has relatively small cooling effectiveness at the leading edge and at the trailing edge. Heat transfer augmentation in these critical areas may be achieved by using impingement jets and film cooling. A computer program is written in FORTRAN IV for IBM 370/165 computer.

Hosny, W. M.↗
NASA NTRS1982

Numerical solutions for the spin-up of a stratified fluid

The model of Warn-Varnas et al. (1978) is used to numerically examine the spin-up flow of a thermally stratified fluid in a cylinder with an insulating side wall, and comparison of the results with the laser-Doppler measurements of Lee (1975) shows excellent agreement. It is shown that flow gradients are created in the interior of the fluid during the meridional circulation spin-up phase, and that the azimuthal flow decayed faster than has been predicted by Wallin (1969). It is established that viscous diffusion in the interior, arising from the interior-flow gradients, is the cause of the discrepancy with Wallin's theory.

Hyun, J. M.↗
NASA NTRS1975

Numerical solutions of supersonic and hypersonic laminar flows over a two-dimensional compression corner

An efficient time-splitting, second-order accurate, numerical scheme is used to solve the complete Navier-Stokes equations for supersonic and hypersonic laminar flow over a two-dimensional compression corner. A fine, exponentially stretched mesh spacing is used in the region near the wall for resolving the viscous layer. Good agreement is obtained between the present computed results and experimental measurement for a Mach number of 14.1, a Reynolds number of 104,000, and wedge angles of 15, 18, and 24 deg. The details of the pressure variation across the boundary layer are given, and a correlation between the leading edge shock and the peaks in surface pressure and heat transfer is observed.

Hung, C. M.↗
NASA NTRS1978

Numerical solution of flow problems using body-fitted coordinate systems

The technique of boundary-fitted coordinate systems is based on a method of automatic numerical generation of a general curvilinear coordinate system having a coordinate line coincident with each boundary of a general multi-connected region containing any number of arbitrarily shaped bodies. Once the curvilinear coordinate system is generated, any partial differential system of interest can be solved on this coordinate system by transforming the equations and solving the resulting system in finite difference approximation on the rectangular transformed plane. This method of automatic body-fitted curvilinear coordinate generation is used to construct finite-difference solutions of the full, time dependent Navier-Stokes equations for the unsteady viscous flow about arbitrary two-dimensional airfoils, or any other two-dimensional bodies. Finally, initial results for three-dimensional applications are also presented.

Thompson, J. F.↗
NASA NTRS1959

Numerical Solution of the Flow of a Perfect Gas Over A Circular Cylinder at Infinite Mach Number

A solution for the two-dimensional flow of an inviscid perfect gas over a circular cylinder at infinite Mach number is obtained by numerical methods of analysis. Nonisentropic conditions of curved shock waves and vorticity are included in the solution. The analysis is divided into two distinct regions, the subsonic region which is analyzed by the relaxation method of Southwell and the supersonic region which was treated by the method of characteristics. Both these methods of analysis are inapplicable on the sonic line which is therefore considered separately. The shapes of the sonic line and the shock wave are obtained by iteration techniques. The striking result of the solution is the strong curvature of the sonic line and of the other lines of constant Mach number. Because of this the influence of the supersonic flow on the sonic line is negligible. On comparison with Newtonian flow methods, it is found that the approximate methods show a larger variation of surface pressure than is given by the present solution.

Hamaker, Frank M.↗
NASA NTRS1973

Direct numerical solution of three-dimensional equations containing elliptic operators.

A direct three-dimensional elliptic solver is presented for application in a wide class of numerical methods for solving partial differential equations in physics and engineering. The derived algorithm and FORTRAN code implement Buzbee, Golub and Nielson's proposed extension of Buneman's Cyclic-Reduction Poisson solver to three dimensions. Both a 'most direct' cyclic reduction and a revised method (to eliminate roundoff error difficulties) are derived. Tests on an IBM 360/67 computer, using various optional combinations of subroutines, showed significant differences in accuracy and computing time, with the optimum subroutine combination depending on mesh size.

Martin, E. D.↗
NASA NTRS1978

Numerical solutions and laser-Doppler measurements of spin-up

The spin-up flow in a cylinder of homogeneous fluid has been examined both experimentally and numerically. A series of laser-Doppler measurements was made of the zonal flow over a range of Ekman numbers and Rossby numbers at various locations in the interior of the flow. These measurements exceed previous ones in accuracy. The weak inertial modes excited by the impulsive start are detectable. The numerical simulations used the primitive equations in axisymmetric form and employed finite-difference techniques on both constant and variable grids. The number of grid points necessary to resolve the Ekman layers was determined. A thorough comparison of the simulations and the experimental measurements is made which includes the details of the amplitude and frequency of the inertial modes. Agreement to within the experimental tolerance is achieved. Analytical results for conditions identical to those in the experiments are not available but some similar linear and nonlinear theories are also compared with the experiments.

Warn-Varnas, A.↗
NASA NTRS1979

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗
NASA NTRS1961

Numerical Solutions for Supersonic Flow of an Ideal Gas Around Blunt Two-Dimensional Bodies

The method described is an inverse one; the shock shape is chosen and the solution proceeds downstream to a body. Bodies blunter than circular cylinders are readily accessible, and any adiabatic index can be chosen. The lower limit to the free-stream Mach number available in any case is determined by the extent of the subsonic field, which in turn depends upon the body shape. Some discussion of the stability of the numerical processes is given. A set of solutions for flows about circular cylinders at several Mach numbers and several values of the adiabatic index is included.

Fuller, Franklyn B.↗
NASA NTRS1972

On the numerical solution of two-dimensional, laminar compressible flows with imbedded shock waves.

The complete, time-dependent Navier-Stokes equations are expressed in conservation form and solved by employing an explicit finite difference numerical technique which incorporates artificial viscosity terms of the form first suggested by Rusanov for numerical stability in the vicinity of shock waves. Surface boundary conditions are developed in a consistent and unique manner through the use of a physically oriented extrapolation procedure. From numerical experimentation an extended range for the explicit stability parameter is established. Also employed is an additional convergence parameter which relates incremental spatial steps. Convergence of the transient solution to a steady state flow was obtained after 400 to 500 time steps.

Goodrich, W. D.↗