In this paper we review a number of auxiliary space based preconditioners for the second order definite and semi-definite Maxwell problems discretized with the lowest order Nedelec finite elements. We discuss the para...In this paper we review a number of auxiliary space based preconditioners for the second order definite and semi-definite Maxwell problems discretized with the lowest order Nedelec finite elements. We discuss the parallel implementation of the most promising of these methods, the ones derived from the recent Hiptmair-Xu （HX） auxiliary space decomposition [Hiptmair and Xu, SIAM J. Numer. Anal., 45 （2007）, pp. 2483-2509]. An extensive set of numerical experiments demonstrate the scalability of our implementation on large-scale H（curl） problems.展开更多
We consider H（curl, Ω）-elliptic variational problems on bounded Lipschitz polyhedra and their finite element Galerkin discretization by means of lowest order edge elements. We assume that the underlying tetrahedral...We consider H（curl, Ω）-elliptic variational problems on bounded Lipschitz polyhedra and their finite element Galerkin discretization by means of lowest order edge elements. We assume that the underlying tetrahedral mesh has been created by successive local mesh refinement, either by local uniform refinement with hanging nodes or bisection refinement. In this setting we develop a convergence theory for the the so-called local multigrid correction scheme with hybrid smoothing. We establish that its convergence rate is uniform with respect to the number of refinement steps. The proof relies on corresponding results for local multigrid in a H^1 （Ω）-context along with local discrete Helmholtz-type decompositions of the edge element space.展开更多
A partially orthonormal basis is constructed with better conditioning properties for tetrahedral H（curl）-conforming Nedelec elements. The shape functions are classified into several categories with respect to their ...A partially orthonormal basis is constructed with better conditioning properties for tetrahedral H（curl）-conforming Nedelec elements. The shape functions are classified into several categories with respect to their topological entities on the reference 3-simplex. The basis functions in each category are constructed to achieve maximum orthogonaiity. The numerical study on the matrix conditioning shows that for the mass and quasi-stiffness matrices, and in a logarithmic scale the condition number grows linearly vs. order of approximation up to order three. For each order of approximation, the condition number of the quasi-stiffness matrix is about one order less than the corresponding one for the mass matrix. Also, up to order six of approximation the conditioning of the mass and quasi- stiffness matrices with the proposed basis is better than the corresponding one with the Ainsworth-Coyle basis Internat. J. Numer. Methods. Engrg., 58：2103-2130, 2003. except for order four with the quasi-stiffness matrix. Moreover, with the new basis the composite matrix μM ＋ S has better conditioning than the Ainsworth-Coyle basis for a wide range of the parameter μ.展开更多
We construct a well-conditioned hierarchical basis for triangular H(curl)-conforming elements with selected orthogonality.The basis functions are grouped into edge and interior functions,and the later is further group...We construct a well-conditioned hierarchical basis for triangular H(curl)-conforming elements with selected orthogonality.The basis functions are grouped into edge and interior functions,and the later is further grouped into normal and bubble functions.In our construction,the trace of the edge shape functions are orthonormal on the associated edge.The interior normal functions,which are perpendicular to an edge,and the bubble functions are both orthonormal among themselves over the reference element.The construction is made possible with classic orthogonal polynomials,viz.,Legendre and Jacobi polynomials.For both the mass matrix and the quasi-stiffness matrix,better conditioning of the new basis is shown by a comparison with the basis previously proposed by Ainsworth and Coyle[Comput.Methods.Appl.Mech.Engrg.,190(2001),6709-6733].展开更多
We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin(IPDG-H)method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equ...We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin(IPDG-H)method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations.The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain.It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method.The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals.Within a unified framework for adaptive finite element methods,we prove the reliability of the estimator up to a consistency error.The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.展开更多
基金This work performed under the auspices of the U.S.Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.UCRL-JRNL-237306
文摘In this paper we review a number of auxiliary space based preconditioners for the second order definite and semi-definite Maxwell problems discretized with the lowest order Nedelec finite elements. We discuss the parallel implementation of the most promising of these methods, the ones derived from the recent Hiptmair-Xu （HX） auxiliary space decomposition [Hiptmair and Xu, SIAM J. Numer. Anal., 45 （2007）, pp. 2483-2509]. An extensive set of numerical experiments demonstrate the scalability of our implementation on large-scale H（curl） problems.
基金supported in part by China NSF under the grant 60873177by the National Basic Research Project under the grant 2005CB321702
文摘We consider H（curl, Ω）-elliptic variational problems on bounded Lipschitz polyhedra and their finite element Galerkin discretization by means of lowest order edge elements. We assume that the underlying tetrahedral mesh has been created by successive local mesh refinement, either by local uniform refinement with hanging nodes or bisection refinement. In this setting we develop a convergence theory for the the so-called local multigrid correction scheme with hybrid smoothing. We establish that its convergence rate is uniform with respect to the number of refinement steps. The proof relies on corresponding results for local multigrid in a H^1 （Ω）-context along with local discrete Helmholtz-type decompositions of the edge element space.
文摘A partially orthonormal basis is constructed with better conditioning properties for tetrahedral H（curl）-conforming Nedelec elements. The shape functions are classified into several categories with respect to their topological entities on the reference 3-simplex. The basis functions in each category are constructed to achieve maximum orthogonaiity. The numerical study on the matrix conditioning shows that for the mass and quasi-stiffness matrices, and in a logarithmic scale the condition number grows linearly vs. order of approximation up to order three. For each order of approximation, the condition number of the quasi-stiffness matrix is about one order less than the corresponding one for the mass matrix. Also, up to order six of approximation the conditioning of the mass and quasi- stiffness matrices with the proposed basis is better than the corresponding one with the Ainsworth-Coyle basis Internat. J. Numer. Methods. Engrg., 58：2103-2130, 2003. except for order four with the quasi-stiffness matrix. Moreover, with the new basis the composite matrix μM ＋ S has better conditioning than the Ainsworth-Coyle basis for a wide range of the parameter μ.
基金The research was supported in part by a DOE grant 304(DEFG0205ER25678)a NSFC grant(10828101).
文摘We construct a well-conditioned hierarchical basis for triangular H(curl)-conforming elements with selected orthogonality.The basis functions are grouped into edge and interior functions,and the later is further grouped into normal and bubble functions.In our construction,the trace of the edge shape functions are orthonormal on the associated edge.The interior normal functions,which are perpendicular to an edge,and the bubble functions are both orthonormal among themselves over the reference element.The construction is made possible with classic orthogonal polynomials,viz.,Legendre and Jacobi polynomials.For both the mass matrix and the quasi-stiffness matrix,better conditioning of the new basis is shown by a comparison with the basis previously proposed by Ainsworth and Coyle[Comput.Methods.Appl.Mech.Engrg.,190(2001),6709-6733].
基金The work of the first author has been supported by the German Na-tional Science Foundation DFG within the Research Center MATHEON and by the WCU program through KOSEF(R31-2008-000-10049-0).The other authors acknowledge sup-port by the NSF grant DMS-0810176.1
文摘We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin(IPDG-H)method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations.The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain.It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method.The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals.Within a unified framework for adaptive finite element methods,we prove the reliability of the estimator up to a consistency error.The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.
