SEARCH · Search NASA

Results for “Implicit methods”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1979

A two dimensional model of viscous slot flow

The effects of viscosity on flows through slots are investigated. Using Briley's formulation of the Alternating Direction Implicit method together with the Thompson transformation, numerical solution of viscous incompressible flow through a slot is simulated, for Reynolds number 1000 and openness ratio one tenth. The problem is formulated in terms of stream function and vorticity. A wall iteration of Israeli type is employed to iteratively induce the no slip condition on velocity.

Cooke, C. H.↗

NASA NTRS1980

Asymptotic features of shock-wave boundary-layer interaction

A semi-implicit method is applied to solve the Navier-Stokes equations numerically and to evaluate the features of the free-interaction phenomenon that occurs when a shock wave impinges on a Blasius boundary layer. Comparisons are made with predictions of the triple-deck theory and experiment. Results include pressure and skin-friction distribution in the free-interaction region for various values of Reynolds number.

Hussaini, M. Y.↗

NASA NTRS1981

Analysis of turbulent flow about an isolated airfoil using a time dependent Navier-Stokes procedure

The procedure solves the Navier-Stokes equations by the consistently split linearized block implicit method of Briley and McDonald in a body fitted coordinate system. The procedure is described and results are presented for flow about an airfoil whose incidence changes from 6 degrees to 19 degrees at a Reynolds number of one million and Mach number of 0.2. In addition, the unsteady flow about an airfoil held at a constant 19 degree incidence is examined and compared to data.

Shamroth, S. J.↗

NASA NTRS1981

Accelerated finite-volume calculation of transonic potential flows

A fully conservative finite-volume algorithm is used to calculate transonic potential flows past isolated airfoils and through two-dimensional channels. The difference equations are solved by a multigrid technique which uses an alternating direction implicit method as a smoothing algorithm. The finite-volume formulation provides a framework within which to treat flows past complicated geometries, while the multigrid/alternating direction scheme provides rapid convergence of the solution to very small residuals.

Jameson, A.↗

NASA NTRS1983

The NCOREL computer program for 3D nonlinear supersonic potential flow computations

An innovative computational technique (NCOREL) was established for the treatment of three dimensional supersonic flows. The method is nonlinear in that it solves the nonconservative finite difference analog of the full potential equation and can predict the formation of supercritical cross flow regions, embedded and bow shocks. The method implicitly computes a conical flow at the apex (R = 0) of a spherical coordinate system and uses a fully implicit marching technique to obtain three dimensional cross flow solutions. This implies that the radial Mach number must remain supersonic. The cross flow solutions are obtained by using type dependent transonic relaxation techniques with the type dependency linked to the character of the cross flow velocity (i.e., subsonic/supersonic). The spherical coordinate system and marching on spherical surfaces is ideally suited to the computation of wing flows at low supersonic Mach numbers due to the elimination of the subsonic axial Mach number problems that exist in other marching codes that utilize Cartesian transverse marching planes.

Siclari, M. J.↗

NASA NTRS1985

Numerical simulation of flow through scramjet inlets using a three-dimensional Navier-Stokes code

A three-dimensional Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either the fully explicit or explicit-implicit method due to MacCormack. An algebraic two-layer eddy viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects which are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flow field. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Kumar, A.↗

NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗

NASA NTRS1986

Numerical solution of 2-D thrust reversing and thrust vectoring nozzle

The flowfield within and around two dimensional thrust reversing and thrust vectoring nozzles has been calculated using a new unfactored implicit method with a multiple zone grid. Computations are done for fully deployed thrust reversing nozzles, partially deployed thrust reversing nozzles with thrust vectoring, and a nozzle transitioning from partially to fully deployed. Agreement with available experimental data is good.

Imlay, S.↗

NASA NTRS1987

Efficient solutions to the Euler equations for supersonic flow with embedded subsonic regions

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two dimensions is described. Convergence of the basic algorithm to the steady state is quadratic for fully supersonic flows and is linear for other flows. This is in contrast to the block alternating direction implicit methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented herein is easily coupled with methods to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, and yields a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, Robert W.↗

NASA NTRS1987

Numerical analysis of flow through scramjet engine inlets

A set of computer programs has been developed to analyze flow through supersonic combustion ramjet (scramjet) inlets. These programs solve either the two- or three-dimensional Euler/Navier-Stokes equations in full conservation form by MacCormack's explicit or explicit-implicit method. An algebraic two-layer eddy viscosity model is used for turbulent flow calculations. The programs are operational on Control Data CYBER-200 series vector-processing computer system and have been optimized to take maximum advantage of the vector processing capability of the system. Since their development, the programs have been extensively verified and used to analyze a number of very complex inlet configurations. In this paper, results are presented from two-dimensional, quasi-three-dimensional, and three-dimensional analyses of the inlet flow field to illustrate the use of the programs.

Kumar, Ajay↗

NASA NTRS1987

Numerical analysis of flow through scramjet engine inlets

A set of computer programs was developed to analyze flow through supersonic combustion ramjet (scramjet) inlets. These programs solve either the two or three dimensional Euler/Navier-Stokes equations in full conservation form by MacCormack's explicit or explicit-implicit method. An algebraic two-layer eddy viscosity model is used for turbulent flow calculations. The programs are operational in Control Data CYBER-200 series vector-processing computer system and were optimized to take maximum advantage of the vector processing capability of the system. Since their development, the programs were extensively verified and used to analyze a number of very complex inlet configurations. Results are presented from two-dimensional, quasi-three-dimensional, and three-dimensional analyses of the inlet flow field to illustrate the use of the program.

Kumar, Ajay↗

NASA NTRS1988

Unsteady aerodynamics of rotor-stator interaction in a turbine stage

The two-dimensional, compressible, unsteady, thin-layer Navier-Stokes equations have been used to study the turbulent flow field of a turbine rotor/stator configuration in the high pressure oxidizer turbopump of the Space Shuttle main engine. The calculations were performed on a system of patched and overlaid grids. The grid system consists of 'O-' and 'H-' grids, with the inner 'O-'grid enclosing the blade surface for an accurate resolution of the leading and trailing edges, while the outer 'H-'grid makes the treatment of boundary conditions easier. The integration scheme used is an iterative, factored, implicit method with numerical fluxes evaluated by the third-order accurate upwind-biased Osher scheme. Computed results in the form of pressure contours, Mach number contours, time-averaged surface pressure, unsteady pressure amplitude and unsteady velocity vectors are presented. Many flow phenomena, such as wake cutting by rotor airfoils and evolution of unsteady passage vortices, are discussed.

Yang, R.-J.↗

NASA NTRS1988

Accuracy of a class of concurrent algorithms for transient finite element analysis

The accuracy of a new class of concurrent procedures for transient finite element analysis is examined. A phase error analysis is carried out which shows that wave retardation leading to unacceptable loss of accuracy may occur if a Courant condition based on the dimensions of the subdomains is violated. Numerical tests suggest that this Courant condition is conservative for typical structural applications and may lead to a marked increase in accuracy as the number of subdomains is increased. Theoretical speed-up ratios are derived which suggest that the algorithms under consideration can be expected to exhibit a performance superior to that of globally implicit methods when implemented on parallel machines.

Ortiz, Michael↗

NASA NTRS1989

Computation of aircraft component flow fields at transonic Mach numbers using a three-dimensional Navier-Stokes algorithm

A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.

Shrewsbury, George D.↗

NASA NTRS1988

Element-by-element and implicit-explicit finite element formulations for computational fluid dynamics

Preconditioner algorithms to reduce the computational effort in FEM analyses of large-scale fluid-dynamics problems are presented. A general model problem is constructed on the basis of the convection-diffusion equation and the two-dimensional vorticity/stream-function formulation of the Navier-Stokes equations; this problem is then analyzed using element-by-element, implicit-explicit, and adaptive implicit-explicit approximation schemes. Numerical results for the two-dimensional advection and rigid-body rotation of a cosine hill, flow past a circular cylinder, and driven cavity flow are presented in extensive graphs and shown to be in good agreement with those obtained using implicit methods.

Tezduyar, T. E.↗

NASA NTRS1989

Numerical simulation of flow through a two-strut scramjet inlet

A three-dimensional, Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either a fully explicit or explicit-implicit method. An algebraic, two-layer eddy-viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects that are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flowfield. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Kumar, Ajay↗

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 NTRS1990

Extension of a three-dimensional viscous wing flow analysis user's manual: VISTA 3-D code

Three-dimensional unsteady viscous effects can significantly influence the performance of fixed and rotary wing aircraft. These effects are important in both flows about helicopter rotors in forward flight and flows about three-dimensional (swept and tapered) supercritical wings. A computational procedure for calculating such flow field was developed. The procedure is based upon an alternating direction technique employing the Linearized Block Implicit method for solving three-dimensional viscous flow problems. In order to demonstrate the viability of this method, two- and three-dimensional problems are computed. These include the flow over a two-dimensional NACA 0012 airfoil under steady and oscillating conditions, and the steady, skewed, three-dimensional flow on a flat plate. Although actual three-dimensional flows over wings were not obtained, the ground work was laid for considering such flows. In this report a description of the computer code is given.

Weinberg, Bernard C.↗