Approximate factorization implicit methods

A totally different approach to solving the Navier-Stokes equations is made in alternating direction implicit (ADI) methods (Briley and McDonald, 1975) and approximate factorization implicit (AFI) methods (Beam and Warming, 1977, 1978). These methods apply an approximate spatial factorization technique to avoid the inversion of huge banded matrices and are computationally very efficient. 

The fractional step, or time splitting, concept is more a generic operator splitting approach than a particular solution method [30, 211, 124, 92, 49]. It is essentially an approximate factorization of the methods applied to the different operators in an equation or a set of equations. The overall set of operators can be solved explicitly, implicitly or by a combination of both implicit and explicit discretization schemes. 

A variety of explicit (Dufort-Frankel, Lax-Wendroff, Runge-Kutta) and implicit (approximate factorization, LU-SGS) or hybrid schemes have been employed for integration in time. Because of the complexity of the incompressible Navier-Stokes equations, stability analyses to determine critical time steps are difficult. As a general rule, the allowable time step for an explicit method is proportional to the ratio of the smallest grid size to the largest convective velocity (or the wave propagation speed for an artificial compressibility method). 

Since Sg satisfies a system of equations identical with that satisfied by X, except that Rg replaces H, a standard method for improving the result of a direct method is to do the replacement and solve as before for the correction. This obtained correction will not, of course, be the true Sg) but only an approximation Sg, and it will have been obtained as a result of a set of operations that is equivalent to the formation of CaRg, in accordance with Eq. (2-3). The Ca = C is not known explicitly, but is defined implicitly (see the methods of factorization below). 

In our approach [1, 2] termed the dynamic method the complex susceptibility x = x — ix" is determined by a law of undamped motion of a dipole in a given potential well and by dissipation mechanism often described as stosszahlansatz in the underlying kinetic or Boltzmann equation. In this review we shall refer to this (dynamic) method as the ACF method, since it is actually based on calculation of the spectrum of the dipolar autocorrelation function (ACF). Actually we use a one-particle approximation, in which the form of an employed potential well (being in many cases rectangular or close to it) is taken a priori. Correlation of the particles coordinates is characterized implicitly by the Kirkwood correlation factor g, its value being taken from the experimental data. The ACF method is simple and effective, because we do not employ the stochastic equations of motions. This feature distinguishes our method from other well-known approaches—for example, from those described in books [13, 14]. 

An analysis of chemical desorption has recently been published (Chem.Eng.Sci., 21 0980)), which is based on a number of simplifying assumptions the film theory model is assumed, the diffusivities of all species are taken to be equal to each other, and in the solution of the differential equations an approximation which is second order with respect to distance from the gas-liquid interface is used this approximation was introduced as early as 1948 by Van Krevelen and Hoftizer. However, the assumptions listed above are not at all drastic, and two crucial elements are kept in the analysis reversibility of the chemical reactions and arbitrary chemical mechanisms and stoichiometry.The result is a methodology for developing, for any given chemical mechanism, a highly nonlinear, implicit, but algebraic equation for the calculation of the rate enhancement factor as a function of temperature, bulk-liquid composition, interface gas partial pressure and physical mass transfer coefficient The method of solution is easily gene ralized to the case of unequal diffusivities and corrections for differences between the film theory and the penetration theory models can be calculated. 

---