scholarly journals Numerical Analysis of Least-Squares Group Finite Element Method for Coupled Burgers' Problem

2021 ◽  
Vol 18 (4(Suppl.)) ◽  
pp. 1521
Author(s):  
Najat Jalil Noon

In this paper, a least squares group finite element method for solving coupled Burgers' problem in   2-D is presented. A fully discrete formulation of least squares finite element method is analyzed, the backward-Euler scheme for the time variable is considered, the discretization with respect to space variable is applied as biquadratic quadrangular elements with nine nodes for each element. The continuity, ellipticity, stability condition and error estimate of least squares group finite element method are proved.  The theoretical results  show that the error estimate of this method is . The numerical results are compared with the exact solution and other available literature when the convection-dominated case to illustrate the efficiency of the proposed method that are solved through implementation in MATLAB R2018a.

Author(s):  
Rajeev Kumar ◽  
Brian H. Dennis

The least-squares finite element method (LSFEM), which is based on minimizing the l2-norm of the residual, has many attractive advantages over Galerkin finite element method (GFEM). It is now well established as a proper approach to deal with the convection dominated fluid dynamic equations. The least-squares finite element method has a number of attractive characteristics such as the lack of an inf-sup condition and the resulting symmetric positive system of algebraic equations unlike GFEM. However, the higher continuity requirements for second-order terms in the governing equations force the introduction of additional unknowns through the use of an equivalent first-order system of equations or the use of C1 continuous basis functions. These additional unknowns lead to increased memory and computing time requirements that have prevented the application of LSFEM to large-scale practical problems, such as three-dimensional compressible viscous flows. A simple finite element method is proposed that employs a least-squares method for first-order derivatives and a Galerkin method for second order derivatives, thereby avoiding the need for additional unknowns required by pure a LSFEM approach. When the unsteady form of the governing equations is used, a streamline upwinding term is introduced naturally by the leastsquares method. Resulting system matrix is always symmetric and positive definite and can be solved by iterative solvers like pre-conditioned conjugate gradient method. The method is stable for convection-dominated flows and allows for equalorder basis functions for both pressure and velocity. The stability and accuracy of the method are demonstrated with preliminary results of several benchmark problems solved using low-order C0 continuous elements.


Author(s):  
Rajeev Kumar ◽  
Brian H. Dennis

The least-squares finite element method (LSFEM), based on minimizing the l2-norm of the residual is now well established as a proper approach to deal with the convection dominated fluid dynamic equations. The least-squares finite element method has a number of attractive characteristics such as the lack of an inf-sup condition and the resulting symmetric positive system of algebraic equations unlike Galerkin finite element method (GFEM). However, the higher continuity requirements for second-order terms in the governing equations force the introduction of additional unknowns through the use of an equivalent first-order system of equations or the use of C1 continuous basis functions. These additional unknowns lead to increased memory and computational requirements that have limited the application of LSFEM to large-scale practical problems. A novel finite element method is proposed that employs a least-squares method for first-order derivatives and a Galerkin method for second order derivatives, thereby avoiding the need for additional unknowns required by a pure LSFEM approach. When the unsteady form of the governing equations is used, a streamline upwinding term is introduced naturally by the least-squares method. Resulting system matrix is always symmetric and positive definite and can be solved by iterative solvers like pre-conditioned conjugate gradient method. The method is stable for convection-dominated flows and allows for equal-order basis functions for both pressure and velocity. The method has been successfully applied here to solve complex buoyancy-driven flow with Boussinesq approximation in a square cavity with differentially heated vertical walls using low-order C0 continuous elements.


2012 ◽  
Vol 2012 ◽  
pp. 1-24 ◽  
Author(s):  
Yang Liu ◽  
Hong Li ◽  
Jinfeng Wang ◽  
Wei Gao

A new positive definite expanded mixed finite element method is proposed for parabolic partial integrodifferential equations. Compared to expanded mixed scheme, the new expanded mixed element system is symmetric positive definite and both the gradient equation and the flux equation are separated from its scalar unknown equation. The existence and uniqueness for semidiscrete scheme are proved and error estimates are derived for both semidiscrete and fully discrete schemes. Finally, some numerical results are provided to confirm our theoretical analysis.


Sign in / Sign up

Export Citation Format

Share Document