Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Implicit finite-difference computations of unsteady transonic flows about airfoils, including the treatment of irregular shock-wave motions

A computer code, LTRAN2, has been constructed that efficiently computes low-frequency unsteady transonic flows about airfoils in motion. The code solves the two-dimensional, nonlinear, low-frequency, small-disturbance equation by an alternating-direction implicit (ADI) algorithm. First, as a check on the code, LTRAN2 is used to solve the linear, low-frequency, small-disturbance equation. The results are compared with known answers from linear theory. Then nonlinear results are presented. Tijdeman's experimental observations of shock wave motions resulting from airfoil flap oscillations are qualitatively reproduced by LTRAN2. These computations compare favorably with those of Magnus and Yoshihara and are obtained in substantially less computer time. Also the underlying theory of the governing equation and the construction of the algorithm are discussed.

Ballhaus, W. F.↗
NASA NTRS1977

An implicit shock-fitting scheme for unsteady transonic flow computations

The alternating-direction implicit scheme developed by NASA Ames for unsteady transonic flows has been modified to include a shock-fitting algorithm as well as an analytically stretched coordinate system. The shock-fitting procedure treats shock waves as discontinuities normal to the free stream. Improvements in shock position and the unsteady pressure distributions are obtained by this modification. The various types of shock motion observed experimentally by Tijdeman are well simulated in calculations using the modified computational scheme. The method of detecting shock wave formation and the procedure for fitting a moving shock wave are illustrated. Results for a pulsating parabolic arc airfoil and for an NACA 64A006 airfoil with oscillating quarter-chord flap are presented and discussed.

Yu, N. J.↗
NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗
NASA NTRS1977

Implicit finite difference simulation of flow about arbitrary geometries with application to airfoils

Finite difference procedures are used to solve either the Euler equations or the 'thin layer' Navier-Stokes equations subject to arbitrary boundary conditions. An automatic grid generation program is employed, and because an implicit finite difference algorithm for the flow equations is used, time steps are not severely limited when grid points are finely distributed. Computational efficiency and compatibility to vectorized computer processors is maintained by use of approximate factorization techniques. Computed results for both inviscid and viscous flow about airfoils are described and compared to various known solutions.

Steger, J. L.↗
NASA NTRS1977

Application of an axisymmetric implicit blunt body procedure - Computation of solar wind flows past terrestrial planets

A recently developed axisymmetric implicit unsteady Euler equation solver has been applied to the calculation of solar wind flows past terrestrial planets. The primary emphasis of this paper is on the application of that computational technique to the calculation of the flow field properties of steady, axisymmetric, supersonic/hypersonic flows past blunt-nosed shapes characteristic of interplanetary solar wind flows past both magnetic and nonmagnetic planets. The purpose is to demonstrate both the accuracy and flexibility of this new technique as well as to present a novel application of such supersonic blunt-body procedures.

Chaussee, D. S.↗
NASA NTRS1977

Study of the application of an implicit model-following flight controller to lift-fan VTOL aircraft

An implicit model-following flight controller is proposed. This controller is relatively simple in concept: it provides an input/output relationship that is approximately that of any selected second order system; it provides good gust alleviation; and it is self-trimming. The flight controller was applied to all axes of a comprehensive mathematical model of a lift-fan V/STOL transport. Power management controls and displays were designed to match the various modes of control provided by the flight controller. A piloted simulation was performed using a six degree of freedom simulator. The fixed-operating-point handling qualities throughout the powered lift flight envelope received pilot ratings of 3-1/2 or better. Approaches and vertical landings in IFR zero-zero conditions received pilot ratings varying from 2-1/2 to 4 depending on the type of approach and weather conditions.

Merrick, V. K.↗
NASA NTRS1978

On improving the iterative convergence properties of an implicit approximate-factorization finite difference algorithm

The iterative convergence properties of an approximate-factorization implicit finite-difference algorithm are analyzed both theoretically and numerically. Modifications to the base algorithm were made to remove the inconsistency in the original implementation of artificial dissipation. In this way, the steady-state solution became independent of the time-step, and much larger time-steps can be used stably. To accelerate the iterative convergence, large time-steps and a cyclic sequence of time-steps were used. For a model transonic flow problem governed by the Euler equations, convergence was achieved with 10 times fewer time-steps using the modified differencing scheme. A particular form of instability due to variable coefficients is also analyzed.

Desideri, J. A.↗
NASA NTRS1978

An implicit algorithm for the conservative transonic full potential equation using an arbitrary mesh

A new, implicit approximate factorization (AF) algorithm designed to solve the conservative full-potential equation for the transonic flow past arbitrary airfoils has been developed. The new algorithm uses an upwind bias of the density coefficient to provide stability in supersonic regions. This allows the simple two- and three-banded matrix form of the AF scheme to be retained over the entire flow field, even in regions of supersonic flow. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Airfoil pressure distributions have been computed and are in good agreement with independent results.

Holst, T. L.↗
NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗
NASA NTRS1978

Navier-Stokes calculations with a coupled strongly implicit method. Part 2: Spline solutions

A coupled strongly implicit method is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equation. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered. In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗
NASA NTRS1979

An implicit algorithm for the conservative, transonic full-potential equation with effective rotated differencing

A new differencing scheme for the conservative full potential equation which effectively simulates rotated differencing is presented. The scheme was implemented by an appropriate upwind bias of the density coefficient along coordinate directions. A fast, fully implicit, approximate factorization iteration scheme was then used to solve the resulting difference equations. Solutions for a number of traditionally difficult transonic airfoil test cases are presented.

Holst, T. L.↗
NASA NTRS1978

On the accuracy and convergence of implicit numerical integration of finite element generated ordinary differential equations

A study of accuracy and convergence of linear functional finite element solution to linear parabolic and hyperbolic partial differential equations is presented. A variable-implicit integration procedure is employed for the resultant system of ordinary differential equations. Accuracy and convergence is compared for the consistent and two lumped assembly procedures for the identified initial-value matrix structure. Truncation error estimation is accomplished using Richardson extrapolation.

Baker, A. J.↗
NASA NTRS1980

Calculation of two-dimensional inlet flow fields in a supersonic free stream by an implicit marching code with nonorthogonal mesh generation: User's manual

An implicit, shock-capturing finite-difference code which is used to calculate two-dimensional inlet flow fields in a supersonic free stream is explained. The Euler equations are subjected to general nonorthogonal transformation and a body-fitted coordinate system is employed. The mathematical formulation of the problem is given along with the numerical algorithm. Initial and boundary conditions, numerical stability, program limitations, and accuracy is discussed. An overall program logic as well as instructions for program use and operation are also furnished.

Biringen, S. H.↗
NASA NTRS1980

A diagonal form of an implicit approximate-factorization algorithm with application to a two dimensional inlet

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic free stream flow about and through inlets. The modification transforms the coupled system of equations into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D. S.↗
NASA NTRS1980

Calculation of inlet flow fields by an implicit technique

Two-dimensional inlet flow fields in a supersonic free stream are calculated by an implicit, shock-capturing, finite-difference method. The Euler equations are subjected to a general curvilinear transformation and a body-fitted coordinate system is employed. The method is used to solve supercritical, critical, and subcritical flow fields which are simulated by prescribing appropriate conditions at the inlet outflow boundary. Results are presented for a drooped-cowl inlet.

Biringen, S.↗
NASA NTRS1980

Implicit computations of unsteady transonic flow governed by the full-potential equation in conservation form

An alternating-direction implicit algorithm is presented for solving the conservative, full-potential equation for unsteady, transonic flow. A new development is the time-linearization of the density function. This linearization reduces the solution process from one of solving a system of two equations at each mesh point to one of solving a single equation. Two sample cases are computed. First, a one-dimensional traveling shock wave is computed and compared with the analytic solution. Second, a two-dimensional case is computed of a flow field that results from a thickening and subsequently thinning airfoil. The resulting flow field, which includes a traveling shock wave, is compared to the flow field obtained from the low-frequency, small-disturbance, transonic equation.

Goorjian, P. M.↗
NASA NTRS1978

An efficient explicit-implicit-characteristic method for solving the compressible Navier-Stokes equations

Explicit, implicit, and characteristic finite-difference methods are applied to solve model equations representative of the compressible Navier-Stokes equations. An approach is then formulated for solving the Navier-Stokes equation at high Reynolds numbers. The approach has drastically reduced the computation time required to obtain viscous flow solutions. Computational results for shock wave separated flows are presented.

Maccormack, R. W.↗
NASA NTRS1980

Implicit model following and parameter identification of unstable aircraft

A transformation in the s-plane is described which has utility in implicit model-following optimal control design application and in estimation or parameter identification problems. The objective of the transformation is, for the control problem, to achieve an unstable closed-loop system, and, for the estimation problem, to alleviate algorithm convergence problems that may arise in identifying unstable systems. For the control problem, the transformation is a shift along the real (sigma) axis of the plant and model poles and zeros. This transformation is shown to be equivalent to a modified performance index but offers the advantage of compatibility with existing optimal control solution algorithms. For the estimation problem, the data are multiplied by an exponential function and the assumed measurement and process noise covariances are appropriately modified. Examples of both control and estimation applications are presented.

Lebacqz, J. V.↗