The Least-Squares Finite Element Method (LSFEM) has been around for decades. Recent developments have demonstrated both theoretical completeness and numerical stability of the optimal LSFEM in incompressible fluid dynamics applications. Along with the wider applications of LSFEM, researchers and practitioners have observed, in certain ranges of CFD problems, some different levels of non-conservative behavior pending on the least-squares formulations. Specifically noticeable terms are the mass conservation in incompressible flow calculations. In this paper, numerical results of different levels of conservation inconsistency are presented. By carefully examining the optimal LSFEM equations and the solution process, it is observed that the equal-weighting nature of the LSFEM formulation is the numerical culprit of such conservation problem. As the additional variables in multi-physics simulations are included, the non-conservative behavior becomes more pronounced. A numerical weighting method is proposed to solve the least-squares first order differential equations by selectively imposing higher weighting on the critical mass conservation equation based on the actual physical properties of the problem. Improved numerical results using such weighted optimal least-squares method are presented to demonstrate the applicability to general problems. Several possible refinement approaches are also discussed.
展开▼