Some new discretization and adaptation and multigrid methods for 2-D 3-T diffusion equations

In the simulation of laser-driven implosion of a fuel capsule in inertial confinement fusion experiments, a system of two-dimensional diffusion equations coupled with electron, iron and photon temperature are widely used to approximately describe the process of energy across multiple materials and the exchange of energy among electrons, irons and photons. The numerical solution of such equations is always challenging because of its strong nonlinear phenomena and strong discontinuous interfaces. In this article, we design a symmetric finite volume method, develop the corresponding preconditioning technique, and propose a mesh adaptation algorithm based on Hessian matrix and a two-grid method. Using these new methods, we demonstrate that the energy conservation error and computation efficiency of the integrated algorithm are much better than classical method.

[1]  Pingwen Zhang,et al.  A Moving Mesh Finite Element Algorithm for Singular Problems in Two and Three Space Dimensions , 2002 .

[2]  Shi Shu,et al.  Symmetric finite volume discretizations for parabolic problems , 2003 .

[3]  J. Pasciak,et al.  Parallel multilevel preconditioners , 1990 .

[4]  Gabriel Wittum,et al.  Parallel adaptive multigrid algorithm for 2-d 3-t diffusion equations , 2004, Int. J. Comput. Math..

[5]  Enrique Zuazua,et al.  A two-grid approximation scheme for nonlinear Schrödinger equations: dispersive properties and convergence , 2005 .

[6]  Carsten Carstensen,et al.  Each averaging technique yields reliable a posteriori error control in FEM on unstructured grids. Part I: Low order conforming, nonconforming, and mixed FEM , 2002, Math. Comput..

[7]  Weizhang Huang,et al.  Variational mesh adaptation II: error estimates and monitor functions , 2003 .

[8]  Simon N. Chandler-Wilde,et al.  A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane , 2002, Numerische Mathematik.

[9]  Long Chen,et al.  Optimal anisotropic meshes for minimizing interpolation errors in Lp-norm , 2007, Math. Comput..

[10]  O. Landen,et al.  Comparisons of line emission from 2- and 3-dimensionsal simulations of ICF capsules to experiments , 2003 .

[11]  Weizhang Huang.  Mathematical Principles of Anisotropic Mesh Adaptation , 2006 .

[12]  Yunqing Huang,et al.  A SYMMETRIC FINITE VOLUME ELEMENT SCHEME ON QUADRILATERAL GRIDS AND SUPERCONVERGENCE , 2006 .

[13]  Pingwen Zhang,et al.  Moving mesh methods in multiple dimensions based on harmonic maps , 2001 .

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

[15]  Jinchao Xu,et al.  A preconditioned GMRES method for nonsymmetric or indefinite problems , 1992 .

[16]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique , 1992 .

[17]  Jinchao Xu.  Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[18]  Jinchao Xu.  A new class of iterative methods for nonselfadjoint or indefinite problems , 1992 .

[19]  Jinchao Xu,et al.  Asymptotically Exact A Posteriori Error Estimators, Part II: General Unstructured Grids , 2003, SIAM J. Numer. Anal..

[20]  Tao Tang,et al.  A multilevel successive iteration method for nonlinear elliptic problems , 2003, Math. Comput..

[21]  Mikaelian.  LASNEX simulations of the classical and laser-driven Rayleigh-Taylor instability. , 1990, Physical review. A, Atomic, molecular, and optical physics.