Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Thrust chamber performance using Navier-Stokes solution

The viscous, axisymmetric flow in the thrust chamber of the space shuttle main engine (SSME) was computed on the CRAY 205 computer using the general interpolants method (GIM) code. Results show that the Navier-Stokes codes can be used for these flows to study trends and viscous effects as well as determine flow patterns; but further research and development is needed before they can be used as production tools for nozzle performance calculations. The GIM formulation, numerical scheme, and computer code are described. The actual SSME nozzle computation showing grid points, flow contours, and flow parameter plots is discussed. The computer system and run times/costs are detailed.

Chan, J. S.↗
NASA NTRS1985

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from a source surface, given the proper boundary conditions. The results of this work; describing the mathematical formalism that was developed, the numerical techniques employed, and the stability criteria developed for these numerical schemes are presented. An analytical solution is used for a test case; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent ( 5%).

Wu, S. T.↗
NASA NTRS1985

Absorbing boundary conditions for exterior problems

Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary condition at infinity and be accurate with respect to the interior numerical scheme. Such boundary conditions are commonly referred to as absorbing boundary conditions. A treatment is given of these boundary conditions for wave-like equations.

Hariharan, S. I.↗
NASA NTRS1985

A procedure for the calculation of supersonic flows with strong viscous-inviscid interaction

The present paper is concerned with the calculation of strong viscous-inviscid interactions in two-dimensional laminar supersonic flows with and without separation. The equations solved are the so-called parabolized Navier-Stokes equations. The streamwise pressure gradient term is written as a combination of a forward and a backward difference to provide a path for upstream propagation of information. Global iteration is utilized to repeatedly update the pressure field from an initial guess until convergence is achieved. The numerical scheme employed is a new alternating direction explicit (ADE) procedure which is used as an alternative to the more difficult to program multigrid strategy to accelerate convergence. Results are presented for flows past two flat plate related bodies.

Barnett, M.↗
NASA NTRS1985

Contributions to the finite element solution of the fan noise radiation problem

The radiation of fan generated noise to the far field from a nacelle of realistic geometry is investigated using the finite element method. Several innovations have been introduced to minimize the computational requirements and create a highly efficient numerical scheme. The innovations include: (1) formulation of the problem in terms of velocity potential and density in such a way that no inlet mean flow velocity derivatives are required in the field equations, (2) the use of 'wave envelope' elements in an outer region permitting a grid much coarser than would be used for conventional finite elements, (3) the use of a mesh which deforms with an increase of forward flight speed so that mesh lines are always lines of constant phase and rays for a point source, permitting the use of wave envelope elements and simple boundary conditions for any case of forward velocity, (4) an efficient scheme for introducing the noise source via modal amplitude coefficients, and (5) the use of a frontal solution technique which for physically realistic problems drastically reduces the active storage requirements. The finite element scheme is outlined, as are the specific details of the innovations. Results are given for cases where comparable experimental data are available.

Eversman, W.↗
NASA NTRS1985

A three-dimensional boundary-layer analysis including heat-transfer and blade-rotation effects

An analysis has been developed for calculating three-dimensional viscous flows over arbitrary bodies through the solution of complete three-dimensional boundary-layer equations. The governing equations include the effect of blade rotation and are written in nonorthogonal surface-fitted coordinates to allow maximum flexibility. A coupled finite-difference numerical scheme is used for solving the governing equations. The accuracy of this scheme has been assessed through comparisons with known analytical/experimental results for several test cases. The present analysis has then been applied to compute viscous flows over a turbine endwall and airfoil suction surface, and fairly good correlation has been obtained with available experimental data for surface heat transfer and viscous streamline skewing.

Vatsa, V. N.↗
NASA NTRS1985

Contact stresses in pin-loaded orthotropic plates

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Hyer, M. W.↗
NASA NTRS1985

Simulation of turbulent mixing layers

A hybrid numerical scheme, combined with a scalar convection scheme, is used to compute turbulent mixing layers. The hybrid scheme is used to solve the vorticity equation by dividing the domain into two regions, applying a Lagrangian method in one region and an Eulerian method in the other. The scalar convection scheme is used to track an interface which separates a dyed fluid from a dye-free fluid, thereby distinguishing the fluid that starts on the high speed side from the fluid that starts on the low speed side. A two-dimensional mixing layer has been simulated using this new method; the roll up of a vortex sheet into structures that interact and pair is observed. The growth of the momentum thickness and the development of the mean turbulence statistics are found to be consistent with experimental observations. Forcing of the layer at moderate and low frequencies shows a dramatic effect on the growth of the momentum thickness and the entrainment ratio. Low frequency forcing enhances the growth of the layer; moderate frequency forcing inhibits the growth.

Mansour, N. N.↗
NASA NTRS1986

Transonic aeroelasticity of wings with tip stores

The presence of tip stores influences both aerodynamic and aeroelastic performances of wings. Such effects are more pronounced in the transonic regime. In this study, transonic aeroelasticity of wings with tip stores is studied for the first time by a theoretical method using the unsteady-small disturbance transonic aerodynamic equations coupled with modal structural equations of motion. The aerodynamic and structural equations of motion are simultaneously integrated by a time-accurate numerical scheme. To validate the tip store simulation, aeroelastic computations are made for a typical rectangular wing with a tip store and results are correlated with available wind tunnel data for the corresponding wing without a tip store at various flight conditions. Aeroelastic computations are also made for a typical fighter wing with a tip store. Present computations show that it is important to account for the aerodynamics of the tip store, particularly in the transonic regime where the tip store can make the wing aeroelastically unstable.

Guruswamy, G. P.↗
NASA NTRS1986

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗
NASA NTRS1987

Numerical simulation of transonic propeller flow using a 3-dimensional small disturbance code employing novel helical coordinates

The numerical simulation of three-dimensional transonic flow about propeller blades is discussed. The equations for the unsteady potential flow about propellers is given for an arbitrary coordinate system. From this the small disturbance form of the equation is derived for a new helical coordinate system. The new coordinate system is suited to propeller flow and allows cascade boundary conditions to be applied straightforward. A numerical scheme is employed which solves the steady flow as an asymptotic limit of unsteady flow. Solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an eight bladed cascade. Both high and low advance ratio cases are given which include a lifting case as well as nonlifting cases. The nonlifting cases are compared to solutions from a Euler code.

Snyder, Aaron↗
NASA NTRS1986

Computer code for gas-liquid two-phase vortex motions: GLVM

A computer program aimed at the phase separation between gas and liquid at zero gravity, induced by vortex motion, is developed. It utilizes an explicit solution method for a set of equations describing rotating gas-liquid flows. The vortex motion is established by a tangential fluid injection. A Lax-Wendroff two-step (McCormack's) numerical scheme is used. The program can be used to study the fluid dynamical behavior of the rotational two-phase fluids in a cylindrical tank. It provides a quick/easy sensitivity test on various parameters and thus provides the guidance for the design and use of actual physical systems for handling two-phase fluids.

Yeh, T. T.↗
NASA NTRS1986

Some aspects of vortex flows determined from the thin-layer Navier-Stokes equations

Vortex flows caused by a swept wing are described in some detail. The thin-shear-layer Navier-Stokes equations with an algebraic eddy-viscosity turbulence model are solved. An implicit, factored numerical scheme is used. On the upper surface of the wing, flow is primarily along the free-stream and spanwise direction, respectively, at small and large angles of attack.

Mehta, Unmeel↗
NASA NTRS1986

A model of inner cometary ionospheres

A hydrodynamical model to describe the movement of the thermalized charged components in the inner ionosphere of comet Halley is presented. Photoelectrons are included by applying a two-stream-type approach. The numerical scheme describes shock transitions in a natural way. Solutions are obtained for a number of different assumptions concerning electron heating rates but all show that the electron temperature increases sharply where the collisional electron neutral coupling becomes unimportant. This temperature increase is accompanied by an increase in the plasma pressure and in its associated polarization electric field, and causes the plasma flow to go subsonic. In certain cases this transition occurs as an inner shock which may explain the observed ion pile-up in Halley's comet.

Korosmezey, A.↗
NASA NTRS1987

An implementation of a barotropic quasigeostrophic model of ocean circulation on the MPP

The implementation on the Massively Parallel Processor (MPP) of a barotropic quasigeostrophic model of ocean circulation is discussed. The mathematical model, including scalings and boundary conditions is discussed. The numerical scheme, which uses compact differencing is also discussed. The implementation of this model on the MPP is then presented. Finally, some performance results are given and compared to results obtained using the VPS-32 and one processor of a CRAY-2.

Grosch, C. E.↗
NASA NTRS1986

Reynolds stress closure in jet flows using wave models

A wave model was developed for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Numerical methods are under development to solve the nonseparable Rayleigh equation. A model problem is being used to assist in the algorithm development. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique are being examined.

Morris, P. J.↗
NASA NTRS1985

Nonlinear acoustic wave propagation in atmosphere. Absorbing boundary conditions for exterior problems

Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary conditions at infinity and be accurate with respect to the interior numerical scheme. The treatment of these boundary conditions for wave-like equations is discussed.

Hariharan, S. I.↗
NASA NTRS1987

On the interaction of small and large eddies in two dimensional turbulent flows

Some results concerning the interaction of small and large eddies to two dimensional turbulent flows are presented. It is shown that the amplitude of small structures decays exponentially to a small value, and from this is inferred a simplified interaction law of small and large eddies. Beside their intrinsic interest for the understanding of the physics of turbulence, these results lead to new numerical schemes to be studied in a separate work.

Foias, C.↗