Three discontinuous Galerkin schemes for the anisotropic heat conduction equation on non-aligned grids

DOI10.1016/J.CPC.2015.10.009zbMATH Open1344.76054arXiv1505.05383OpenAlexW1866828981MaRDI QIDQ310583FDOQ310583

Authors: M. Held, M. Wiesenberger, Andreas Stegmeir 

Publication date: 8 September 2016

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: We present and discuss three discontinuous Galerkin (dG) discretizations for the anisotropic heat conduction equation on non-aligned cylindrical grids. Our most favourable scheme relies on a self-adjoint local dG (LDG) discretization of the elliptic operator. It conserves the energy exactly and converges with arbitrary order. The pollution by numerical perpendicular heat fluxes degrades with superconvergence rates. We compare this scheme with aligned schemes that are based on the flux-coordinate independent approach for the discretization of parallel derivatives. Here, the dG method provides the necessary interpolation. The first aligned discretization can be used in an explicit time-integrator. However, the scheme violates conservation of energy and shows up stagnating convergence rates for very high resolutions. We overcome this partly by using the adjoint of the parallel derivative operator to construct a second self-adjoint aligned scheme. This scheme preserves energy, but reveals unphysical oscillations in the numerical tests, which result in a decreased order of convergence. Both aligned schemes exhibit low numerical heat fluxes into the perpendicular direction. We build our argumentation on various numerical experiments on all three schemes for a general axisymmetric magnetic field, which is closed by a comparison to the aligned finite difference (FD) schemes of References [1,2]

Full work available at URL: https://arxiv.org/abs/1505.05383

Recommendations

• Finite-difference schemes for anisotropic diffusion
• Solution of 3D heat conduction equations using the discontinuous Galerkin method on unstructured grids
• A local discontinuous Galerkin method for heat-conduction equations with discontinuous coefficients
• Modelling of heat transport in magnetised plasmas using non-aligned coordinates
• scientific article; zbMATH DE number 5492345

zbMATH Keywords

discontinuous Galerkin methodsanisotropic heat conduction in magnetized plasmasfield line mapflux-coordinate independent approachnon-aligned coordinatessupport operator method

Mathematics Subject Classification ID

Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element methods applied to problems in fluid mechanics (76M10) Statistical mechanics of plasmas (82D10) Ionized gas flow in electromagnetic fields; plasmic flow (76X05)

Cites Work

• Title not available (Why is that?)
• The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems
• Superconvergence of the local discontinuous Galerkin method for elliptic problems on Cartesian grids
• Title not available (Why is that?)
• Finite-difference schemes for anisotropic diffusion
• On Locking and Robustness in the Finite Element Method
• An a priori error analysis of the local discontinuous Galerkin method for elliptic problems
• Nonlinear magnetohydrodynamics simulation using high-order finite elements.
• Support-operator finite-difference algorithms for general elliptic problems
• Finite element and higher order difference formulations for modelling heat transport in magnetised plasmas
• Modelling of heat transport in magnetised plasmas using non-aligned coordinates
• A conservative discontinuous Galerkin scheme for the 2D incompressible Navier-Stokes equations
• A flux-coordinate independent field-aligned approach to plasma turbulence simulations
• Title not available (Why is that?)
• Title not available (Why is that?)

Cited In (14)

• A finite volume flux coordinate independent approach
• New 25-point stencils with optimal accuracy for 2-D heat transfer problems. Comparison with the quadratic isogeometric elements
• Fourier-spectral element approximation of the ion-electron Braginskii system with application to tokamak edge plasma in divertor configuration
• Advances in the flux-coordinate independent approach
• Analysis of locally-aligned and non-aligned discretisation schemes for reactor-scale tokamak edge turbulence simulations
• A high-order non field-aligned approach for the discretization of strongly anisotropic diffusion operators in magnetic fusion
• An efficient high-order solver for diffusion equations with strong anisotropy on non-anisotropy-aligned meshes
• Numerical evaluation of line, surface and toroidal integrals on level sets of toroidally symmetric functions
• A partially mesh-free scheme for representing anisotropic spatial variations along field lines
• Streamline integration as a method for two-dimensional elliptic grid generation
• An efficient high-order numerical solver for diffusion equations with strong anisotropy
• GENE-X: a full-\(f\) gyrokinetic turbulence code based on the flux-coordinate independent approach
• Streamline integration as a method for structured grid generation in X-point geometry
• Directional enrichment functions for finite element solutions of transient anisotropic diffusion

Uses Software

• GRILLIX
• NIMROD

This page was built for publication: Three discontinuous Galerkin schemes for the anisotropic heat conduction equation on non-aligned grids

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q310583)