An Efficient High-order Numerical Solver for Diffusion Equations with Strong Anisotropy

In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. We demonstrate the accuracy produced by the high-order scheme and develop an efficient preconditioning technique to solve the corresponding linear system, which is robust to the mesh size and anisotropy of the problem. Several numerical tests are provided to validate the accuracy and efficiency of the proposed algorithm.

[1]  Markus Held,et al.  Three discontinuous Galerkin schemes for the anisotropic heat conduction equation on non-aligned grids , 2015, Comput. Phys. Commun..

[2]  Jinchao Xu,et al.  The auxiliary space method and optimal multigrid preconditioning techniques for unstructured grids , 1996, Computing.

[3]  Jacob B. Schroder,et al.  Smoothed aggregation solvers for anisotropic diffusion , 2012, Numer. Linear Algebra Appl..

[4]  E. LESTER SMITH,et al.  AND OTHERS , 2005 .

[5]  S. V. Nepomnyaschikh,et al.  Decomposition and Fictitious Domains Methods for Elliptic Boundary Value Problems , 2017 .

[6]  Min Tang,et al.  Uniformly Convergent Scheme for Strongly Anisotropic Diffusion Equations with Closed Field Lines , 2018, SIAM J. Sci. Comput..

[7]  Luca F. Pavarino Courant.  Domain Decomposition Algorithms for the p-version Finite Element Method for Elliptic Problems , 1992 .

[8]  Weizhang Huang,et al.  Maximum principle for the finite element solution of time‐dependent anisotropic diffusion problems , 2013 .

[9]  Benedict Dingfelder,et al.  A locally field-aligned discontinuous Galerkin method for anisotropic wave equations , 2020, J. Comput. Phys..

[10]  R. Eymard,et al.  3D Benchmark on Discretization Schemes for Anisotropic Diffusion Problems on General Grids , 2008 .

[11]  Bobby Philip,et al.  Adaptive algebraic smoothers , 2012, J. Comput. Appl. Math..

[12]  Luca F. Pavarino.  Additive Schwarz methods for thep-version finite element method , 1993 .

[13]  Béatrice Rivière,et al.  Discontinuous Galerkin Methods for Second-Order Elliptic PDE with Low-Regularity Solutions , 2011, J. Sci. Comput..

[14]  Barry Koren,et al.  Finite-difference schemes for anisotropic diffusion , 2014, J. Comput. Phys..

[15]  A. Ern,et al.  A discontinuous Galerkin method with weighted averages for advection–diffusion equations with locally small and anisotropic diffusivity , 2008 .

[16]  M. V. Umansky,et al.  On Numerical Solution of Strongly Anisotropic Diffusion Equation on Misaligned Grids , 2005 .

[17]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[18]  Michael Griebel,et al.  Tensor product type subspace splittings and multilevel iterative methods for anisotropic problems , 1995, Adv. Comput. Math..

[19]  Weizhang Huang,et al.  Discrete maximum principle and a Delaunay-type mesh condition for linear finite element approximations of two-dimensional anisotropic diffusion problems , 2010, 1008.0562.

[20]  Jinchao Xu,et al.  Analysis of a two-level method for anisotropic diffusion equations on aligned and nonaligned grids , 2013, Numer. Linear Algebra Appl..

[21]  Jonathan J. Hu,et al.  A new smoothed aggregation multigrid method for anisotropic problems , 2007, Numer. Linear Algebra Appl..

[22]  HEN,et al.  Convergence Analysis of V-Cycle Multigrid Methods for Anisotropic Elliptic Equations , 2011 .

[23]  Paola F. Antonietti,et al.  Schwarz domain decomposition preconditioners for discontinuous Galerkin approximations of elliptic problems: non-overlapping case , 2007 .

[24]  Weizhang Huang,et al.  An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems , 2010, J. Comput. Phys..

[25]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[26]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[27]  Patrick Tamain,et al.  A high-order non field-aligned approach for the discretization of strongly anisotropic diffusion operators in magnetic fusion , 2020, Comput. Phys. Commun..

[28]  Blanca Ayuso de Dios,et al.  Uniformly Convergent Iterative Methods for Discontinuous Galerkin Discretizations , 2009, J. Sci. Comput..

[29]  Ludmil T. Zikatanov,et al.  An algebraic multilevel method for anisotropic elliptic equations based on subgraph matching , 2012, Numer. Linear Algebra Appl..

[30]  P. B. Snyder,et al.  BOUT++: A framework for parallel plasma fluid simulations , 2008, Comput. Phys. Commun..

[31]  Sibylle Günter,et al.  Finite element and higher order difference formulations for modelling heat transport in magnetised plasmas , 2007, J. Comput. Phys..

[32]  Uri Shumlak,et al.  Spectral element spatial discretization error in solving highly anisotropic heat conduction equation , 2010, Comput. Phys. Commun..

[33]  Rusty C. Holleman,et al.  Numerical diffusion for flow‐aligned unstructured grids with application to estuarine modeling , 2013 .

[34]  M. Ottaviani,et al.  An alternative approach to field-aligned coordinates for plasma turbulence simulations , 2010, 1002.0748.

[35]  Ludmil T. Zikatanov,et al.  A Uniform Additive Schwarz Preconditioner for High-Order Discontinuous Galerkin Approximations of Elliptic Problems , 2017, J. Sci. Comput..

[36]  S. C. Brenner,et al.  Convergence of Multigrid Algorithms for Interior Penalty Methods , 2005 .

[37]  Guillaume Latu,et al.  Comparison of Numerical Solvers for Anisotropic Diffusion Equations Arising in Plasma Physics , 2015, J. Sci. Comput..

[38]  Stephen C. Jardin,et al.  A triangular finite element with first-derivative continuity applied to fusion MHD applications , 2004 .

[39]  Weizhang Huang,et al.  Conditioning of finite element equations with arbitrary anisotropic meshes , 2012, Math. Comput..

[40]  Panayot S. Vassilevski,et al.  Algebraic Multilevel Preconditioning of Anisotropic Elliptic Problems , 1994, SIAM J. Sci. Comput..

[41]  STEVE SCHAFFER,et al.  A Semicoarsening Multigrid Method for Elliptic Partial Differential Equations with Highly Discontinuous and Anisotropic Coefficients , 1998, SIAM J. Sci. Comput..

[42]  Jacques Liandrat,et al.  A new conservative finite-difference scheme for anisotropic elliptic problems in bounded domain , 2020, J. Comput. Phys..

[43]  Sibylle Günter,et al.  Modelling of heat transport in magnetised plasmas using non-aligned coordinates , 2005 .