Mathematical modeling of two-phase compressible fluid filtration based on modified adaptive method of minimum amendments

The work objective is to build and investigate the modified adaptive method of minimum amendments (MAMMA) which is destined for the numerical simulation of the two-phase compressible fluid filtration in porous media. This approach allows overcoming the known use limitations of other methods of the finite-difference equations solution, such as: crucial differential pressures acting on the oil-and-water bearing formation; and the compressibility of the medium at the considerable gas content in the oil phase. An approximation method - an explicit one for defining the function of water saturation, and an implicit one for the pressure function computation - is selected as the research basis. When setting the initial boundary value problem and its sampling, the process of the two-phase compressible fluid filtration in the space-dimensional domain with the lateral area bounded below by the subface of stratum, and above - by the bed top, is considered. A two-layer iterative method of the variational type - a modified method of minimal amendments adapted for solving finite-difference equations of the two-phase compressible fluid with a non-selfadjoint operator under the most general assumptions on the properties of the grid-problem operator is built. It is shown that a MAMMA has the asymptotic convergence rate characteristic of the “classical” alternate triangular method that does not use the Chebyshev acceleration technique and can be applied to the problems with a self-adjoint operator. Numerical experiments have confirmed the high efficiency of MAMMA. It is established that to achieve the specified accuracy, the number of iterations at the MAMMA reduces to 3-20 times as compared to the method of Seidel and the overrelaxation method

Sufficient conditions for convergence of positive solutions to linearized two-dimensional sediment transport problem

Introduction. The sediment transportation is one of the major processes that define the magnitude and back surface changing rate for water bodies. The most used prognostic studies in this field are based on the mathematical models that allow reducing, and in some cases, eliminating the expensive and often environmentally burdensome experiments. Spatially one-dimensional models are mainly used to predict changes in bottom topography. For actual coastal systems with irregular coastal configuration, the flow vector is generally not orthogonal to the tangent line for the coastline at each of its points. It also may not coincide with the wind stress vector. Therefore, to solve lots of practically important problems associated with the prediction of the dynamics of the back surface of water basins, it is necessary to use spatially two-dimensional models of sediment transportation and effective numerical methods of their implementation. Materials and Methods. The authors (A.I. Sukhinov, A.E. Chistyakov, E.A. Protsenko, and V.V. Sidoryakina) have earlier proposed a spatially two-dimensional model of sediment transport that satisfies the basic conservation laws (of material balance and momentum) which is a quasilinear equation of parabolic type. The linear difference schemes are constructed and studied; model and some practically important tasks are solved. However, the theoretical study on the proximity of solutions for the original nonlinear initial-boundary value problem and the linearized continuous task - on which basis a discrete model (difference scheme) was built - remained in the shadow. The study of the linearized problem correctness and the determination of sufficient conditions for positivity of solutions are of special interest because only positive solutions to this sediment transport problem have physical value within the framework of the considered models. Research Results. The investigated nonlinear two-dimensional model of sediment transport in the coastal zone of shallow reservoirs takes account of the following physically significant factors and parameters: soil porosity; critical value of the tangent stress at which load transport starts; turbulent interaction; dynamically varying of the bottom geometry; wind currents; and bottom friction. Linearization is carried out on the time grid - nonlinear coefficients of the parabolic equation are taken for the previous time grid step. Next, a chain of tasks connected by the initial data - final solutions of the linearized mixed Cauchy problems chain on a uniform time grid is constructed, and thus, the linearization for the initial 2D nonlinear model is carried out in large. Earlier, the authors have proved the existence and uniqueness of the solution to a linear tasks chain. Prior estimate of the proximity of the linearized problem chain solution to the initial non-linear task solution has been also obtained. The conditions of its solution positivity and their convergence to the nonlinear sediment transport problem are defined in the norm of the Hilbert space L1 with the rate O(τ) where τ is a time step. Discussion and Conclusions. The obtained research results of the spatially two-dimensional nonlinear sediment transport model can be used for predicting the nonlinear hydrodynamic processes, improving their accuracy and reliability due to the availability of new accounting functionality of physically important factors, including the specification of the boundary conditions