Search NASA⌕ Search

SEARCH · Search NASA

Results for “implicit method”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1980

Navier-Stokes calculations with a coupled strongly implicit method. II Spline deferred-corrector solutions

The coupled strongly implicit (CSIP) method described previously is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equations. 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 (equal to or less than 1000). In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗
NASA NTRS1988

Implicit methods in CFD

A class of implicit approximate factorization schemes is examined for stability and convergence characteristics. These schemes include Newton's method, factorization, and flux-vector splitting. Examples are used to show that all practical methods suffer from some limited stability or asymptotic convergence restriction. A three-dimensional factored scheme is shown which suffers from unconditional instability which can only be ameliorated by added artificial dissipation. An F3D + or - flux split scheme is described which avoids unconditional instability, but in the end has similar convergence characteristics.

Pulliam, Thomas H.↗
NASA NTRS1978

An alternating direction implicit method for the Control Data STAR-100 vector computer

An implementation of the alternating direction implicit (ADI) method for the Control Data STAR-100 computer is presented and analyzed. Two parallel algorithms, both of which are most efficient when used to solve many independent tridiagonal systems of equations, are discussed relative to their usefulness in an ADI implementation on the STAR-100 computer. It is shown that it may be desirable to alternate between the parallel algorithms as the direction of implicitness is alternated in order to eliminate the data rearrangement which would otherwise be required. The applicability of the two parallel tridiagonal solvers to several other numerical algorithms is also discussed.

Lambiotte, J. J., Jr.↗
NASA NTRS1980

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A(0)-stable linear two-step methods in conjunction with the method of approximate factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A(0)-stable method. The parameter space for which the resulting ADI schemes are second-order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗
NASA NTRS1983

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗
NASA NTRS1979

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A sub 0-stable linear two-step methods in conjunction with the method of approximation factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A sub 0-stable method. The parameter space for which the resulting ADI schemes are second order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗
NASA NTRS1997

Euler/Navier-Stokes Methods

Implicit methods have been the workhorse for the Euler and Navier-Stokes equations for the last 25 years. The ground breaking work of Dr. Joe Steger in implementing such techniques in practical Euler and Navier-Stokes codes provided the basis for all the success in this area. This presentation will highlight his contribution and technical excellence in the area of implicit methods for CFD.

Pulliam, Tom↗
NASA NTRS1982

An improved semi-implicit method for structural dynamics analysis

A semi-implicit algorithm is presented for direct time integration of the structural dynamics equations. The algorithm avoids the factoring of the implicit difference solution matrix and mitigates the unacceptable accuracy losses which plagued previous semi-implicit algorithms. This substantial accuracy improvement is achieved by augmenting the solution matrix with two simple diagonal matrices of the order of the integration truncation error.

Park, K. C.↗
NASA NTRS1981

An implicit method for the calculation of inlet flow fields

Inlet flow fields are calculated by an implicit, time marching procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation. Results are presented and discussed for an inlet of current design practice. Results include inviscid calculations performed for supercritical inlet operation with uniform and nonuniform inflow boundary conditions as well as for subcritical inlet operation with uniform inflow boundary conditions. Results for viscous calculations performed for supercritical inlet operation with uniform inflow boundary conditions are also discussed.

Biringen, S.↗
NASA NTRS1989

A fully-coupled implicit method for thermo-chemical nonequilibrium air at sub-orbital flight speeds

A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.

Park, Chul↗
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 NTRS1981

Calculation of two-dimensional inlet flow fields by an implicit method including viscous effects: User's manual

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time marching finite difference procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. The mathematical formulation of the problem and the solution algorithm are given. Numerical stability and accuracy as well as the initial and boundary conditions used are discussed. Instructions for program use and operation along with the overall program logic are also given.

Biringen, S.↗
NASA NTRS1986

A locally implicit method for fluid flow problems

The fluid flow inside the space shuttle main engine (SSME) traverses through a complex geometrical configuration. The flow is compressible, viscous, and turbulent with pockets of separated regions. Several computer codes are being developed to solve three dimensional Navier-Stokes equations with different turbulence models for analyzing the SSME internal flow. The locally implicit scheme is a computationally efficient scheme which converges rapidly in multi-grid modes for elliptic problems. It has the promise of providing a rapidly converging algorithm for steady-state viscous flow problems.

Reddy, K. C.↗
NASA NTRS1991

An implicit method for the nonlinear modelling and simulation of piezoceramic actuators displaying hysteresis

Experiments have shown that piezoceramic materials display a nonlinear relationship between the applied electric field and the actuation strain. This relationship also displays a substantial hysteresis upon reversal of the applied field. In this paper, piezoceramic actuator models are incorporated into the structural equations of motion to arrive at a set of nonlinear actively controlled structural equations of motion. A new implicit algorithm for determining the time history of the actively controlled structure is presented. The algorithm employs the trapezoidal rule for stepping the equations forward in time. The algorithm is compared to an explicit algorithm and is shown to provide greater numerical accuracy. In addition, the numerical stability and convergence characteristics are presented via example. The new algorithm displays convergence and stability properties that are comparable to the standard trapezoidal rule used for dynamic systems without hysteresis.

Leigh, Timothy D.↗