Numerical Conservation Properties of H(div)-Conforming Least-Squares Finite Element Methods for the Burgers Equation
Numerical Conservation Properties of H(div)-Conforming Least-Squares Finite Element Methods for the Burgers Equation
复制标题
Burgers方程的H(div)相容最小二乘有限元方法的数值守恒性质
DOI:
--
复制
发表时间:
2005
影响因子:
3.1
通讯作者:
Luke N. Olson
中科院分区:
文献类型:
--
作者:
H. Sterck;T. Manteuffel;S. McCormick;Luke N. Olson
Least-squares finite element methods (LSFEMs) for the inviscid Burgers equation are studied. The scalar nonlinear hyperbolic conservation law is reformulated by introducing the flux vector, or the associated flux potential, explicitly as additional dependent variables. This reformulation highlights the smoothness of the flux vector for weak solutions, namely, $ff(u) in H({
m div},Omega)$. The standard least-squares (LS) finite element (FE) procedure is applied to the reformulated equations using H(div)-conforming FE spaces and a Gauss--Newton nonlinear solution technique. Numerical results are presented for the one-dimensional Burgers equation on adaptively refined space-time domains, indicating that the H(div)-conforming FE methods converge to the entropy weak solution of the conservation law. The H(div)-conforming LSFEMs do not satisfy a discrete exact conservation property in the sense of Lax and Wendroff. However, weak conservation theorems that are analogous to the Lax--Wendroff theorem for conservative finite difference methods are proved for the H(div)-conforming LSFEMs. These results illustrate that discrete exact conservation in the sense of Lax and Wendroff is not a necessary condition for numerical conservation but can be replaced by minimization in a suitable continuous norm.