Parallel Iterative Methods with Factored Preconditioning Matrices for Solving Discrete Elliptic Equations

Parallel approximate factorization method for solving discrete elliptic equations

Applications of the Lanczos method

A variable gravity model of the otolith membrane

A new technique for solving the large system of linear algebraic equations associated with implicit differencing of multidimensio nal partial differential equations is presented. The coefficient matrix of the equations is factored, and then approximations to certain terms in the matrix are obtained from series expansions. The resulting system of equations is solved easily. The method is developed and demonstrated using a simple representative two-dimensional equation. Very good results are obtained when one direction is dominant. MPLICIT finite-difference schemes for the solution of multidimensional partial differential equations are usually stable and therefore applicable to a large class of problems. However, they are difficult to implement and may require an excessive amount of computer storage and time. The long computing time arises from the need to solve the large system of linear algebraic equations that result from the differencing. The computing time can be reduced significantly by approximating the coefficient matrix of the linear equations with a matrix that produces a system of equations that are relatively easy to solve. Among such methods are the alternating direction method (ADI)l used by Beam and Warming2 and Stone's strongly implicit method,3 which has been tested by Linetal.4 In this paper, a new technique for solving the large system of linear algebraic equations associated with implict differencing of multidimensional partial differential equations is presented. This method, called the pseudo-elimination method (PE), is shown to be faster than Stone's method for certain problems. The method is directly applicable to linear and linearized nonlinear systems of parabolic or elliptic partial differential equations. In order to discuss the method, a simple linear partial differential equation will be used; however, it should be kept in mind that the PE method is applicable to much more complicated problems. The question of whether the method will work when applied to difficult problems is not addressed. The scope of this paper is limited to presenting the method and illustrating, via a simple problem, that the method has some merit and deserves further study.

The purpose of this paper is to derive a means of scaling the dependent variables of a system of linear differential equations so that the normalized system accurately reflects the dependence of the system on initial data. The theory is developed for the special case of a system of first order hyperbolic partial differential equations in two independent variables. Several other systems of differential equations, to which the same normalization could be applied, are also indicated. The scaling is equivalent to row and column scaling commonly used in linear algebra. For the case of positive matrices (that is, for the case where the diagonalizing matrix and its inverse have no zero elements), the theory for the “1” and “$\infty $ ” norms has been thoroughly developed by Bauer [1]. The assumption of positive matrices is too stringent for some applications. Thus, the theory is generalized to show that for irreducible matrices a best scaling exists (for both the “1” and “ $\infty $” norms), while for simply reducible matrices a best scaling exists for at least one of these two norms. Several examples are given.

Introduction. Computer Program Packages. Typical Problems. Classification of Equations. Discrete Methods. Finite Differences and Computational Molecules. Finite Difference Operators. Method of Weighted Residuals. Finite Elements. Method of Lines. Errors. Stability and Convergence. Irregular Boundaries. Choice of Discrete Network. Dimensionless Forms. Parabolic Equations: Introduction. Properties of a Simple Explicit Method. Fourier Stability Method. Implicit Methods. Additional Stability Considerations. Matrix Stability Analysis. Extension of Matrix Stability Analysis. Consistency, Stability, And Convergence. Pure Initial Value Problems. Variable Coefficients. Examples of Equations with Variable Coefficients. General Concepts of Error Reduction. Methods of Lines (MOL) for Parabolic Equations. Weighted Residuals and the Method of Lines. Bubnov-Galerkin Scheme for Parabolic Equations. Finite Elements and Parabolic Equations. Hermite Basis. Finite Elements and Parabolic Equations. General Basis Fucntion. Finite Elements and Parabolic Equations. Special Basis Functions. Explicity Finite Difference Methods for Nonlinear Problems. Further Applications on One Dimentions. Asymmetric Approximations. Elliptic Equations: Introduction. Simple Finite Difference Schemes. Direct Methods. Iterative Methods. Linear Elliptic Equations. Some Poit Iterative Methods. Convergence of Point Iterative Methods. Rates of Convergence. Accelerations. Conjugate Gradient Method. Extensions of SOR. Auliative Examples of Over-Relaxation. Other Point Iterative Methods. Block Iterative Methods. Alternating Direction Methods. Summary of ADI Results. Triangular Elements. Boundary Element Method (BEM). Spectral Methods. Some Nonlinear Examples. Hyperbolic Equations: Introduction. The Quasilinear System. Introductory Examples. Method of Characteristics. Constant States and Simple Waves. Typical Application of Characteristics. Finite Differences for First-Order Equations. Lax-Wendroff Methods and Other Algorithms Dissipation and Dispersion. Explicity Finite Difference Methods. Attenuation. Implicit Methods for Second-Order Equations. Time Quasilinear Examples. Simultaneous First-Order Equations. Explicit Methods. An Implicity Method for First-Order Equations. Hybrid Methods for First-Order Equations. Finite Elements and the Wave Equation. Spectral Methods and Periodic Systems. Gas Dyunamics in One Space Variable. Eulerian Difference Equations. Lagrangian Difference Equations. Hopscotch Methods for Conservation Laws. Explicity-Implicity Schemes for Conservation Laws. Special Topics: Introduction. Singularities. Shocks. Eigenvale Problems. Parabolic Equations in Segeral Space Variables. Additional Comments on Elliptic Equations. Hyperbolic Equations in Higher Dimensions. Mixed Systems. Higher Order Equations in Elasticity and Vibrations. Computational Fluid Mechanics. Stream Function. Vorticity Method for Fluid Mechanics. Primitive Variable Methods for Fluid Mechanics. Vector Potential Methods for Fluid Mechanics. Introduction to Monte Carlo Mehtods. Fast Fourier Transform and Applications. Method of Fractional Steps. Applications of Group Theory in Computation. Computational Ocean Acoustics. Enclosure Methods. Chapter References. Author Index. Subject Index.

Based on investigating the eigenvalue distribution of the preconditioned matrices that arise in the finite-element discretization of indefinite elliptic problems with preconditioning matrices, constructed on the basis of the principal coercive part of the considered elliptic problems, it is proved that the Krylov subspace CG-type methods like MINRES converge with a rate of convergence determined only by the stability properties of the discrete problem in terms of a priori estimates and the accuracy of the preconditioning matrix with respect to the principal coercive part of the considered finite-element stiffness matrix. This is also demonstrated by numerical experiments. The method is of practical interest if there are only a few negative eigenvalues.

The work presents the analysis of the free boundary value problem related to the one-dimensional invasion model of new species in biofilm reactors. In the framework of continuum approach to mathematical modelling of biofilm growth, the problem consists of a system of non-linear hyperbolic partial differential equations governing the microbial species growth and a system of semi-linear elliptic partial differential equations describing the substrate trends. The model is completed with a system of elliptic partial differential equations governing the diffusion and reaction of planktonic cells, which are able to switch their mode of growth from planktonic to sessile when specific environmental conditions are found. Two systems of non-linear differential equations for the substrate and planktonic cells mass balance within the bulk liquid are also considered. The free boundary evolution is governed by a differential equation that accounts for detachment. The qualitative analysis is performed and a uniqueness and existence result is presented. Furthermore, two special models of biological and engineering interest are discussed numerically. The invasion of Anammox bacteria in a constituted biofilm inhabiting the deammonification units of the wastewater treatment plants is simulated. Numerical simulations are run to evaluate the influence of the colonization process on biofilm structure and activity.

Parallel iterative regularization methods for solving systems of ill-posed equations

CHAPTER 9 - SYSTEMS OF DIFFERENTIAL EQUATIONS

We considered a non-stationary optimal control system with delay. Non-stationary optimal control systems are described by systems of differential or differential-algebraic equations. Variable coefficients do not allow, in general, to construct a solution to such systems in quadratures. Numerical or asymptotic methods are used to solve such systems. In this paper, asymptotic methods are used, in particular, the Feshchenko-Shkil' method for integrating of singularly perturbed systems and the Wasow’s method for systems with an unstable spectrum. In this paper, we construct a transformation that reduces the optimal delay control system to a system that does not contain terms without delay. This transformation makes it possible to integrate the system by the method of steps without solving systems of differential equations at each step. The system of equations obtained as a result of the transformation of the original system is somewhat easier to study in terms of building a solution. However, the problem of optimizing the control of both systems requires both a separate mathematical study and clarification of the practical reality of spectrum instability in such systems. If the instability of the spectrum is caused by the degeneracy of the main matrix, this leads to the unboundedness of the system solution as the small parameter approaches zero. The aforementioned growth of the solution can create emergency situations in real systems. Future research concerns solving the problem of optimal control of systems with an unstable spectrum and studying the question of the reality and physical meaning of turning points in specific systems. In future studies, it is also planned to construct an asymptotic representation of the solution of a control system with delays and degeneracies. Keywords: optimal control system, systems of differential equations with delay, asymptotic solution, turning point.

In the article we obtained sufficient conditions of the existence of the nonlinear Noetherian boundary value problem solution for the system of differential-algebraic equations which are widely used in mechanics, economics, electrical engineering, and control theory. We studied the case of the nondegenerate system of differential algebraic equations, namely: the differential algebraic system that is solvable relatively to the derivative. In this case, the nonlinear system of differential algebraic equations is reduced to the system of ordinary differential equations with an arbitrary continuous function. The studied nonlinear differential-algebraic boundary-value problem in the article generalizes the numerous statements of the non-linear non-Gath boundary value problems considered in the monographs of А.М. Samoilenko, E.A. Grebenikov, Yu.A. Ryabov, A.A. Boichuk and S.M. Chuiko, and the obtained results can be carried over matrix boundary value problems for differential-algebraic systems. The obtained results in the article of the study of differential-algebraic boundary value problems, in contrast to the works of S. Kempbell, V.F. Boyarintsev, V.F. Chistyakov, A.M. Samoilenko and A.A. Boychuk, do not involve the use of the central canonical form, as well as perfect pairs and triples of matrices. To construct solutions of the considered boundary value problem, we proposed the iterative scheme using the method of simple iterations. The proposed solvability conditions and the scheme for finding solutions of the nonlinear Noetherian differential-algebraic boundary value problem, were illustrated with an example. To assess the accuracy of the found approximations to the solution of the nonlinear differential-algebraic boundary value problem, we found the residuals of the obtained approximations in the original equation. We also note that obtained approximations to the solution of the nonlinear differential-algebraic boundary value problem exactly satisfy the boundary condition.

Parallel diagonally implicit Runge-Kutta-Nyström methods

Boundary-fitted computational meshes are generated by solving a system of elliptic partial differential equations. The grid points at the boundary surfaces are arbitrarily distributed. The interior grid distribution is controlled by a technique which insures a grid line orthogonality at the boundaries. The grid is automatically adjusted while the grid line orthogonality at the boundaries is achieved in an iterative process. This technique is demonstrated for three different configurations. Desirable cell size and shape are obtained for the boundary shapes which normally cause poor grid distribution.

In this paper, our main goal is to present a non-classical MPFA method coupled with a high order finite volume method for the simulation of oil water displacements in heterogeneous and anisotropic petroleum reservoirs using general polygonal meshes. The governing equations are solved using the IMPES (IMplicit Pressure and Explicit Saturation) strategy, where the elliptic pressure equation is discretized by a linear Multipoint Flux Approximation method using Harmonic points (MPFA-H) capable to handle strongly heterogeneous and anisotropic media. Besides, to approximate the advective term that characterizes the hyperbolic saturation equation, we use the Multidimensional Optimal Order Detection (MOOD) method. This technique is based in an “a posteriori” limitation procedure, i.e., the limitation procedure is done after the calculation of “candidate solutions” only where necessary to guarantee certain physical properties such as monotonicity. To show the potential of our finite volume formulation, we solve some benchmark problems found in literature.