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

Implicit unsteady transonic airfoil calculations at supersonic freestreams

The computer code LTRAN2 has been extended to compute unsteady transonic flows about oscillating airfoils with supersonic freestreams. The LTRAN2 code uses an alternating direction implicit (ADI) algorithm to solve the two-dimensional, nonlinear, low-frequency, transonic small-disturbance (LF-TSD) equation. The modified code, LTRAN2-SS, includes a 'high-frequency' option. Steady solutions are checked against those computed by the steady TSD code, TSFOIL; unsteady computations of the linear LF-TSD equation are compared with known linear theory solutions; and new unsteady nonlinear solutions are presented. These cases include standard AGARD test cases for the NLR 7301, MBB-A3, and DO Al supercritical airfoils, as well as several NACA airfoils. The modified code enables aerodynamicists to quickly and efficiently compute transonic flows for both subsonic and supersonic freestreams, and thus resolve flutter boundaries through the full extent of the transonic dip phenomenon.

Chow, L. J.↗
NASA NTRS1982

An efficient full potential implicit method based on characteristics for analysis of supersonic flows

A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.

Shankar, V.↗
NASA NTRS1983

Implicit upwind methods for the compressible Navier-Stokes equations

A class of implicit upwind differencing methods for the compressible Navier-Stokes equations is described and applied. The methods are based on the use of local eigenvalues or wave speeds to control spatial differencing of inviscid terms and are aimed at increasing the level of accuracy and stability achievable in computation. Techniques for accelerating the rate of convergence to a steady state solution are also used. Applications to inviscid and viscous transonic flows are discussed and compared with other methods and experimental measurements. It is shown that accurate and efficient transonic airfoil calculations can be made on the Cray-l computer in less than 2 min.

Coakley, T. J.↗
NASA NTRS1982

An implicit formulation for the one-dimensional two-phase multi-interface Stefan Problem

In recent years there has been an intensified effort in the analysis of the Stefan Problem. Particular interest is directed towards determining the location of the interface between the solid and liquid phases of the substance being analyzed. A method has been developed using an implicit finite difference scheme which provides an accurate calculation of the phase change interface locations for a problem with simultaneous multiple phases. This is achieved through an iterative technique which requires the phase change energy to be correctly distributed within the nodes which are undergoing a phase transition. The new method predicts the interface locations without requiring time step size control or node size modification at the interfaces.

Williams, S. D.↗
NASA NTRS1983

Monotone implicit algorithms for the small-disturbance and full potential equations applied to transonic flows

Numerical calculations of transonic flows by potential equations typically use algorithms that change the method of calculation for regions of subsonic and supersonic flow. In this paper, implicit approximate-factorization algorithms are modified to use the monotonic switch in the type of finite-differencing that was developed by Godunov for the Euler equations. Calculations of flows over airfoils by these algorithms are compared with calculations by other methods that are in common usage. For the small-disturbance potential equation, comparisons are made with the Murman-Cole method and the monotone method of Engquist and Osher for both steady and unsteady flows. For the full potential equation, comparisons are made with the methods of Jameson and of Holst and Ballhaus for steady flows. The comparisons show that the monotone methods are more stable. For steady flows, solutions are obtained for cases where the Murman-Cole switch requires a time step over ten times smaller in order for the calculations to remain stable. These improvements are achieved with no increase in computer storage and only minor modifications in current codes.

Goorjian, P. M.↗
NASA NTRS1983

An implicit, transonic, full-potential code for cascade flow on H-grid topology

A transonic, full-potential code is developed for computing the flow through two-dimensional cascades using an H-type grid topology that employs an implicit approximate-factorization scheme. The body-conforming H-grid is generated numerically by solving Poisson's equation. The flow-solution algorithm at the coordinate mapping singularity associated with this grid is investigated using two different types of finite-difference schemes. The grid-geometry effect on these schemes is also studied by noting free-stream capturing properties. It is found that by implementing a consistent spatial differencing scheme, the mapping singularities can be resolved numerically, and the grid-geometry-induced error minimized. The code is verified by computing model cascade flow problems.

Kwak, D.↗
NASA NTRS1983

Unconditionally stable implicit-explicit algorithms for coupled thermal stress waves

An unconditionally stable implicit-explicit method is proposed for the analysis of transient coupled thermal stress waves, and the computer-implementation aspects of the method are discussed. In the method proposed here, the mechanical displacement and temperature are used as independent variables. The resulting coupled finite element matrix equations are symmetric.

Liu, W. K.↗
NASA NTRS1983

Implicit upwind methods for the compressible Navier-Stokes equations

A class of implicit upwind-differencing methods for the compressible Navier-Stokes equations is described and applied. The methods are based on the use of local eigenvalues or wave speeds to control spatial differencing of inviscid terms and are aimed at increasing the level of accuracy and stability achievable in computation. Techniques for accelerating the rate of convergence to a steady-state solution are also used. Applications to inviscid and viscous transonic flows are discussed and compared with other methods and experimental measurements. It is shown that accurate and efficient transonic airfoil calculations can be made on the Cray-1 coomputer in less than 2 min.

Coakley, T. J.↗
NASA NTRS1983

Three-dimensional implicit lambda methods

This paper derives the three dimensional lambda-formulation equations for a general orthogonal curvilinear coordinate system and provides various block-explicit and block-implicit methods for solving them, numerically. Three model problems, characterized by subsonic, supersonic and transonic flow conditions, are used to assess the reliability and compare the efficiency of the proposed methods.

Napolitano, M.↗
NASA NTRS1984

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗
NASA NTRS1984

Computation of viscous flow in planar and axisymmetric ducts by an implicit marching procedure

A streamwise marching procedure, approximately 200 times faster than a full Navier-Stokes procedure with comparable accuracy, is presented for solving problems of compressible viscous subsonic flow. Results are presented and compared with experimental data for the cases of developing turbulent flow in a circular pipe; turbulent flow in a two-dimensional S-duct; and turbulent flow in a typical subsonic diffuser. Prior to each main marching step, a preliminary marching step is taken in which the integral continuity equation and an uncoupled form of the streamwise momentum equation are solved simultaneously to obtain the viscous pressure correction. During the main marching step the equations for continuity, streamwise momentum, cross-flow momentum, and energy are solved simultaneously as a coupled system using an implicit finite-difference method, with the viscous pressure correction treated as a source term. The analysis may be used for flows with both favorable and adverse pressure gradients and to predict the location of flow separation.

Towne, C. E.↗
NASA NTRS1984

A strongly implicit procedure for the compressible Navier-Stokes equations

The strongly implicit procedure of Stone is used to obtain numerical solutions of the compressible Navier-Stokes equations in conservative form. In contrast to the spatially split Douglas-Gunn type methods, the method is shown to be numerically stable for the three-dimensional wave equation. The method is applied to a variety of external and internal two-dimensional flow problems involving shock wave boundary-layer interaction for both laminar and turbulent flows. The results are in good agreement with other methods and/or experiments. The storage penalty associated with the method is discussed and a simple, yet effective, means of minimizing the problem is presented.

Walters, R. W.↗
NASA NTRS1984

Using exact solutions to develop an implicit scheme for the baroclinic primitive equations

The exact solutions presently obtained by means of a novel method for nonlinear initial value problems are used in the development of numerical schemes for the computer solution of these problems. The method is applied to a new, fully implicit scheme on a vertical slice of the isentropic baroclinic equations. It was not possible to find a global scale phenomenon that could be simulated by the baroclinic primitive equations on a vertical slice.

Marchesin, D.↗
NASA NTRS1982

Transonic-flow computation using an explicit-implicit method

The explicit-implicit predictor corrector method of MacCormack (1981) is applied to the analysis of flows past airfoils. By comparing results obtained with different methods and meshes, it is shown that the above method provides, after certain modifications, reasonably good predictions of inviscid and viscous flows about an airfoil. Good results are also obtained for the transonic regime if the free-stream conditions are correct and if a suitable mesh is used.

Kordulla, W.↗
NASA NTRS1984

A full-implicit-continuous-Eulerian (FICE) scheme for multidimensional transient magnetohydrodynamic (MHD) flows

A full implicit continuous Eulerian (FICE) scheme is developed for solving multidimensional transient MHD flow problems. The physical system under consideration is a general case of a transient MHD flow in which an initial steady state is subject to a finite amplitude disturbance. The governing equations are described, their finite difference formulation is presented and the FICE algorithm is given. The boundary conditions are treated by classifying them into physical and computational ones. The usefulness of the FICE algorithm is demonstrated using a physical example concerning the dynamical response of the static solar atmosphere due to a representative photospheric disturbance.

Wu, S. T.↗
NASA NTRS1985

An implicit finite-difference code for a two-equation turbulence model for three-dimensional flows

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

Kaul, U. K.↗
NASA NTRS1985

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗
NASA NTRS1985

Implicit flux-split schemes for the Euler equations

Recent progress in the development of implicit algorithms for the Euler equations using the flux-vector splitting method is described. Comparisons of the relative efficiency of relaxation and spatially-split approximately factored methods on a vector processor for two-dimensional flows are made. For transonic flows, the higher convergence rate per iteration of the Gauss-Seidel relaxation algorithms, which are only partially vectorizable, is amply compensated for by the faster computational rate per iteration of the approximately factored algorithm. For supersonic flows, the fully-upwind line-relaxation method is more efficient since the numerical domain of dependence is more closely matched to the physical domain of dependence. A hybrid three-dimensional algorithm using relaxation in one coordinate direction and approximate factorization in the cross-flow plane is developed and applied to a forebody shape at supersonic speeds and a swept, tapered wing at transonic speeds.

Thomas, J. L.↗