Search for: All records

We consider a model convection-diffusion problem and present useful connections between the finite differences and finite element discretization methods. We introduce a general upwinding Petrov-Galerkin discretization based on bubble modification of the test space and connect the method with the general upwinding approach used in finite difference discretization. We write the finite difference and the finite element systems such that the two corresponding linear systems have the same stiffness matrices, and compare the right hand side load vectors for the two methods. This new approach allows for improving well known upwinding finite difference methods and for obtaining new error estimates. We prove that the exponential bubble Petrov-Galerkin discretization can recover the interpolant of the exact solution. As a consequence, we estimate the closeness of the related finite difference solutions to the interpolant. The ideas we present in this work, can lead to building efficient new discretization methods for multidimensional convection dominated problems. 

For a model convection-diffusion problem, we obtain new error estimates for a general upwinding finite element discretization based on bubble modification of the test space. The key analysis tool is finding representations of the optimal norms on the trial spaces at the continuous and discrete levels. We analyze and compare three methods: the standard linear discretization, the saddle point least square and the upwinding Petrov-Galerkin methods. We conclude that the bubble upwinding Petrov-Galerkin method is the most performant discretization for the one dimensional model. Our results for the model convection-diffusion problem can be extended for creating new and efficient discretizations for the multidimensional cases. 

We consider a model convection-diffusion problem and present our recent analysis and numerical results regarding mixed finite element formulation and discretization in the singular perturbed case when the convection term dominates the problem. Using the concepts of optimal norm and saddle point reformulation, we found new error estimates for the case of uniform meshes. We compare the standard linear Galerkin discretization to a saddle point least square discretization that uses quadratic test functions, and explain the non-physical oscillations of the discrete solutions. We also relate a known upwinding Petrov–Galerkin method and the stream-line diffusion discretization method, by emphasizing the resulting linear systems and by comparing appropriate error norms. The results can be extended to the multidimensional case in order to find efficient approximations for more general singular perturbed problems including convection dominated models. 

We consider a preconditioning technique for mixed methods with a conforming test space and a nonconforming trial space. Our method is based on the classical saddle point disccretization theory for mixed methods and the theory of preconditioning symmetric positive definite operators. Efficient iterative processes for solving the discrete mixed formulations are proposed and choices for discrete compatible spaces are provided. For discretization, a basis is needed only for the test spaces and assembly of a global saddle point system is avoided. We provide approximation properties for the discretization and iteration errors and also provide a sharp estimate for the convergence rate of the proposed algorithm in terms of the condition number of the elliptic preconditioner and the discrete inf− sup and sup− sup constants of the pair of discrete spaces. We focus on applications to elliptic PDEs with discontinuous coefficients. Numerical results for two and three dimensional domains are included to support the proposed method.