Iterative solution of elliptic problems by approximate factorization

AbstractAn iterative method for the numerical solution of singularly perturbed second-order linear elliptic problems is presented. It is a defect correction iteration in which the approximate operator is the product of two first-order operators, which is readily inverted numerically. The approximate operator is generated by formal asymptotic factorization of the original operator. Hence this is a QUasi Analytic Defect correction iteration (QUAD). Both its continuous and discrete versions are analyzed in one dimension. The scheme is extended to a variety of two dimensional operators and it is analyzed for a model advection-diffusion equation. Numerical calculations show the effectiveness of the scheme over a wide range of values of the small parameter

Proceedings for the ICASE Workshop on Heterogeneous Boundary Conditions

Domain Decomposition is a complex problem with many interesting aspects. The choice of decomposition can be made based on many different criteria, and the choice of interface of internal boundary conditions are numerous. The various regions under study may have different dynamical balances, indicating that different physical processes are dominating the flow in these regions. This conference was called in recognition of the need to more clearly define the nature of these complex problems. This proceedings is a collection of the presentations and the discussion groups

A rapidly converging domain decomposition algorithm for a time delayed parabolic problem with mixed type boundary conditions exhibiting boundary layers.

[EN]A rapidly converging domain decomposition algorithm is introduced for a time delayed parabolic problem with mixed type boundary conditions exhibiting boundary layers. Firstly, a space-time decomposition of the original problem is considered. Subsequently, an iterative process is proposed, wherein the exchange of information to neighboring subdomains is accomplished through the utilization of piecewise-linear interpolation. It is shown that the algorithm provides uniformly convergent numerical approximations to the solution. Our analysis utilizes some novel auxiliary problems, barrier functions, and a new maximum principle result. More importantly, we showed that only one iteration is needed for small values of the perturbation parameter

Numerical solution of singularly perturbed convection-diffusion-reaction problems with two small parameters

Preprint version, the final publication is available at Springer via http://dx.doi.org/10.1007/s10543-015-0559-8This paper discusses the numerical solution of 1-D convection-diffusion-reaction problems that are singularly perturbed with two small parameters using a new mesh-adaptive upwind scheme that adapts to the boundary layers. The meshes are generated by the equidistribution of a special positive monitor function. Uniform, parameter independent convergence is shown and holds even in the limit that the small parameters are zero. Numerical experiments are presented that illustrate the theoretical findings, and show that the new approach has better accuracy compared with current methods.DFG, SFB 1029, Substantial efficiency increase in gas turbines through direct use of coupled unsteady combustion and flow dynamic