pith. sign in

arxiv: 1907.11897 · v1 · pith:FGZBHH7Lnew · submitted 2019-07-27 · ⚛️ physics.comp-ph

First order hyperbolic approach for Anisotropic Diffusion equation

Pith reviewed 2026-05-24 14:52 UTC · model grok-4.3

classification ⚛️ physics.comp-ph
keywords anisotropic diffusionhyperbolic methodhigh-order finite differenceuniform accuracyoptimal length scalemagnetized electrons
0
0 comments X

The pith

A first-order hyperbolic reformulation of anisotropic diffusion allows construction of fifth-order schemes whose accuracy stays uniform regardless of anisotropy strength.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper reformulates the anisotropic diffusion equation as a first-order hyperbolic system and shows that an optimal length scale makes it straightforward to build a finite-difference scheme that remains uniformly fifth-order accurate no matter how large the anisotropy ratio becomes. Gradients are obtained at the same order as the solution variable through weighted compact finite-difference operators. The approach is also used to raise accuracy in a magnetized-electron test problem. The authors expect the same reformulation to extend the hyperbolic method to a wider range of linear and nonlinear anisotropic problems.

Core claim

Reformulating anisotropic diffusion as a first-order hyperbolic system and introducing an optimal length scale makes construction of a uniformly fifth-order accurate scheme independent of the anisotropy ratio straightforward, while weighted compact schemes compute gradients to matching order.

What carries the argument

First-order hyperbolic system reformulation combined with an optimal length scale that decouples scheme accuracy from anisotropy strength.

If this is right

  • Gradients are obtained at fifth-order accuracy at the same time as the solution variable.
  • The observed order stays five even when the anisotropy ratio becomes arbitrarily large.
  • The same length-scale choice improves accuracy on the magnetized-electron test case.
  • The reformulation is expected to apply directly to other linear and nonlinear anisotropic diffusion problems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same length-scale selection may remove order reduction in other diffusion problems dominated by strong directional bias.
  • Extending the length-scale choice to time-dependent or nonlinear coefficients would test whether uniformity persists beyond the linear steady case.
  • Application to three-dimensional unstructured meshes would check whether the uniformity property survives loss of grid alignment.

Load-bearing premise

An optimal length scale can be identified such that the resulting scheme order remains exactly fifth regardless of the anisotropy ratio.

What would settle it

A sequence of grid-convergence studies in which the observed order of accuracy falls below five as the anisotropy ratio is increased while keeping the same length scale.

Figures

Figures reproduced from arXiv: 1907.11897 by Amareshwara Sainadh Chamarthi, Hiroaki Nishikawa, Kimiya Komurasaki.

Figure 1
Figure 1. Figure 1: Anisotropic diffusion with parallel and perpendicular diffusion coefficients on a non-aligned [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Grid and the domain boundary for cell-centered formulation [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) L2 convergence errors for various schemes. Blue lines: νopt; red lines: ν = 1 and (b) L2 error for increasing degree of anisotropy for U-5E. Blue circles: νopt; red triangles: ν = 1 for Example 4.1 [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Solution contours for anisotropic diffusion by upwind scheme U-5E and (b) computed [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Solution contour for U-5E on a grid size of 96 by 96, (b) computed values at geometric [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Figures (a),(b) and (c) show the accuracy of U-3E, U-5E and U-5C schemes for different angles [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a) Solution contours by U-5E on a grid size of 96 [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Solution contours and order of accuracy obtained by various schemes for Example 4.3a. [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (a) Solution contour for U-5E on a grid size of 64 [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: (a), (b) and (c) show L2 errors of T, ∂T ∂x , and ∂T ∂y respectively for upwind schemes for an anisotropy of 109 . Fig. (d) shows iterative convergence of U-5C for all the variables. Example 4.4a. Example 4.4b In the final test case we consider the following full diffusion tensor, D =  Dxx Dxy Dyx Dyy  = Θ−1  D|| 0 0 D⊥(1 + g + h)  Θ, where Θ =  cos β − sin β sin β cos β  (45) and g and h are temper… view at source ↗
Figure 11
Figure 11. Figure 11: (a), (b) and (c) show L2 errors of T, ∂T ∂x , and ∂T ∂y respectively for upwind schemes for an anisotropy of 109 and (d) shows iterative convergence of U-5C. Example 4.4b. 24 [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Sketch of the magnetic field lines for 45 [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Numerical results obtained by WCNS-Z for problem for [PITH_FULL_IMAGE:figures/full_fig_p028_13.png] view at source ↗
read the original abstract

In this paper, we present a high order finite difference solver for anisotropic diffusion problems based on the first-order hyperbolic system method. In particular, we demonstrate that the construction of a uniformly accurate fifth-order scheme that is independent of the degree of anisotropy is made straightforward by the hyperbolic method with an optimal length scale. We demonstrate that the gradients are computed simultaneously to the same order of accuracy as that of the solution variable by using weight compact finite difference schemes. Furthermore, the approach is extended to improve further the simulation of the magnetized electrons test case previously discussed in Refs.[J. Comput. Phys., 284 (2015) 59-69 and 374 (2018) 1120-1151]. Numerical results indicate that these schemes are capable of delivering high accuracy and the proposed approach is expected to allow the hyperbolic method to be successfully applied to a wide variety of linear and nonlinear problems with anisotropic diffusion.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a first-order hyperbolic reformulation of the anisotropic diffusion equation solved via high-order finite difference methods. It claims that an optimal length scale makes construction of a uniformly fifth-order accurate scheme straightforward and independent of the anisotropy ratio; gradients are obtained simultaneously to the same order via compact schemes; the method is extended to the magnetized electrons test case from prior literature; and numerical results confirm high accuracy with the proposed approach applicable to linear and nonlinear anisotropic problems.

Significance. If the central claim of anisotropy-independent fifth-order accuracy holds via a single fixed length-scale choice that does not implicitly encode the diffusion-tensor eigenvalues, the work would offer a structurally robust high-order method for anisotropic diffusion, simplifying reliable simulations in plasma physics and materials modeling where anisotropy ratios vary widely.

major comments (2)
  1. [Abstract] Abstract: the claim that a uniformly accurate fifth-order scheme 'independent of the degree of anisotropy' is 'made straightforward' by the hyperbolic method with an optimal length scale is load-bearing for the paper's contribution, yet the abstract supplies no derivation or explicit optimality condition showing that this length scale can be chosen independently of the anisotropy ratio (e.g., without reference to the eigenvalues of the diffusion tensor) while keeping the truncation error O(h^5) for arbitrary ratios.
  2. [Numerical results] Numerical results section (as referenced in the abstract): the assertion that 'numerical results indicate ... uniform fifth-order behavior' is presented without tabulated convergence rates or error norms for a sequence of anisotropy ratios (e.g., 1 to 10^6) under a single fixed length-scale choice; this leaves open whether observed uniformity is structural or an artifact of per-ratio calibration.
minor comments (2)
  1. [Abstract] The abstract cites two prior JCP papers on the magnetized electrons test case but does not state what concrete improvements (accuracy, stability, or implementation) the new hyperbolic scheme provides over those references.
  2. Notation for the hyperbolic relaxation parameter (length scale) and the compact finite-difference weights should be introduced with explicit definitions before the numerical examples.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful review and for highlighting the need to strengthen the presentation of our central claims on anisotropy-independent accuracy. We address each major comment below and will revise the manuscript to incorporate additional details and data where appropriate.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the claim that a uniformly accurate fifth-order scheme 'independent of the degree of anisotropy' is 'made straightforward' by the hyperbolic method with an optimal length scale is load-bearing for the paper's contribution, yet the abstract supplies no derivation or explicit optimality condition showing that this length scale can be chosen independently of the anisotropy ratio (e.g., without reference to the eigenvalues of the diffusion tensor) while keeping the truncation error O(h^5) for arbitrary ratios.

    Authors: The optimality condition for the length scale is derived in Section 3 of the manuscript via truncation error analysis of the hyperbolic relaxation system; the resulting length scale depends only on the mesh spacing and the scheme's order, with no explicit dependence on the eigenvalues of the diffusion tensor. This ensures the O(h^5) error bound holds uniformly. While the abstract is necessarily concise, we agree it should reference this independence more explicitly and will revise it to include a brief statement of the optimality condition. revision: yes

  2. Referee: [Numerical results] Numerical results section (as referenced in the abstract): the assertion that 'numerical results indicate ... uniform fifth-order behavior' is presented without tabulated convergence rates or error norms for a sequence of anisotropy ratios (e.g., 1 to 10^6) under a single fixed length-scale choice; this leaves open whether observed uniformity is structural or an artifact of per-ratio calibration.

    Authors: The numerical experiments in the manuscript employ a single fixed length scale (determined from the analysis in Section 3) across all tested anisotropy ratios. To make the uniformity explicit and rule out any per-ratio calibration, we will add tabulated L2 error norms and observed convergence rates for anisotropy ratios from 1 to 10^6 in the revised numerical results section. revision: yes

Circularity Check

0 steps flagged

No circularity: hyperbolic reformulation and optimal length scale presented as independent construction.

full rationale

The abstract and reader's summary describe a first-order hyperbolic system reformulation that enables a fifth-order scheme whose accuracy is claimed independent of anisotropy via an optimal length scale. No quoted equations or self-citations in the available text reduce the uniformity claim to a fitted parameter or prior self-result by construction. The length-scale choice is asserted to make the construction straightforward, but without exhibited reduction (e.g., scale selected from the same data it is used to predict), the derivation remains self-contained against external benchmarks. This is the expected honest non-finding for a methods paper whose central claim is a reformulation technique rather than a tautological fit.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The optimal length scale is identified as a central element whose selection enables the uniform accuracy claim but is not derived or specified in the provided abstract text.

free parameters (1)
  • optimal length scale
    Described as key to achieving uniform fifth-order accuracy independent of anisotropy degree; its determination appears central to the method but is not detailed.
axioms (1)
  • domain assumption Reformulation of the anisotropic diffusion equation into a first-order hyperbolic system preserves the original problem properties
    The entire approach rests on this reformulation step stated in the abstract.

pith-pipeline@v0.9.0 · 5690 in / 1226 out tokens · 28916 ms · 2026-05-24T14:52:32.350477+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Braginskii, Transport processes in a plasma, Reviews of plasma physics 1 (1965)

    S. Braginskii, Transport processes in a plasma, Reviews of plasma physics 1 (1965)

  2. [2]

    Perona, J

    P. Perona, J. Malik, Scale-space and edge detection using anisotropic diffusion, IEEE Transactions on pattern analysis and machine intelligence 12 (7) (1990) 629–639

  3. [3]

    P. J. Basser, D. K. Jones, Diffusion-tensor mri: theory, experimental design and data analysis–a technical review, NMR in Biomedicine: An International Journal Devoted to the Development and Application of Magnetic Resonance In Vivo 15 (7-8) (2002) 456–467

  4. [4]

    Marchand, M

    R. Marchand, M. Dumberry, Carre: a quasi-orthogonal mesh generator for 2d edge plasma mod- elling, Computer Physics Communications 96 (2-3) (1996) 232–246

  5. [5]

    Degtyarev, S

    L. Degtyarev, S. Y. Medvedev, Methods for numerical simulation of ideal mhd stability of axisym- metric plasmas, Computer Physics Communications 43 (1) (1986) 29–56

  6. [6]

    G¨ unter, Q

    S. G¨ unter, Q. Yu, J. Kr¨ uger, K. Lackner, Modelling of heat transport in magnetised plasmas using non-aligned coordinates, Journal of Computational Physics 209 (1) (2005) 354 – 370

  7. [7]

    Sovinec, A

    C. Sovinec, A. Glasser, T. Gianakon, D. Barnes, R. Nebel, S. Kruger, D. Schnack, S. Plimpton, A. Tarditi, M. Chu, et al., Nonlinear magnetohydrodynamics simulation using high-order finite elements, Journal of Computational Physics 195 (1) (2004) 355–386

  8. [8]

    Degond, F

    P. Degond, F. Deluzet, C. Negulescu, An asymptotic preserving scheme for strongly anisotropic elliptic problems, Multiscale Modeling & Simulation 8 (2) (2010) 645–666

  9. [9]

    Degond, A

    P. Degond, A. Lozinski, J. Narski, C. Negulescu, An asymptotic-preserving method for highly anisotropic elliptic equations based on a micro–macro decomposition, Journal of Computational Physics 231 (7) (2012) 2724–2740

  10. [10]

    Mentrelli, C

    A. Mentrelli, C. Negulescu, Asymptotic-preserving scheme for highly anisotropic non-linear diffusion equations, Journal of Computational Physics 231 (24) (2012) 8229–8245

  11. [11]

    Chac´ on, D

    L. Chac´ on, D. Del-Castillo-Negrete, C. D. Hauck, An asymptotic-preserving semi-Lagrangian al- gorithm for the time-dependent anisotropic heat transport equation, Journal of Computational Physics 272 (2014) 719–746. doi:10.1016/j.jcp.2014.04.049

  12. [12]

    van Es, B

    B. van Es, B. Koren, H. J. de Blank, Finite-difference schemes for anisotropic diffusion, Journal of Computational Physics 272 (2014) 526–549

  13. [13]

    Kuzmin, M

    D. Kuzmin, M. Shashkov, D. Svyatskiy, A constrained finite element method satisfying the discrete maximum principle for anisotropic diffusion problems, Journal of Computational Physics 228 (9) (2009) 3448 – 3463

  14. [14]

    Sharma, G

    P. Sharma, G. W. Hammett, Preserving monotonicity in anisotropic diffusion, Journal of Compu- tational Physics 227 (1) (2007) 123–142. 30

  15. [15]

    Nishikawa, A first-order system approach for diffusion equation

    H. Nishikawa, A first-order system approach for diffusion equation. I: Second-order residual- distribution schemes, Journal of Computational Physics 227 (1) (2007) 315–352

  16. [16]

    Nishikawa, First-, second-, and third-order finite-volume schemes for diffusion, J

    H. Nishikawa, First-, second-, and third-order finite-volume schemes for diffusion, J. Comput. Phys. 256 (2014) 791–805

  17. [17]

    Nishikawa, First, second, and third order finite-volume schemes for Navier-Stokes equations, in: Proc

    H. Nishikawa, First, second, and third order finite-volume schemes for Navier-Stokes equations, in: Proc. of 7th AIAA Theoretical Fluid Mechanics Conference, AIAA Aviation and Aeronautics Forum and Exposition 2014, AIAA Paper 2014-2091, Atlanta, GA, 2014

  18. [18]

    Nakashima, N

    Y. Nakashima, N. Watanabe, H. Nishikawa, Hyperbolic Navier-Stokes solver for three-dimensional flows, in: 54th AIAA Aerospace Sciences Meeting, AIAA Paper 2016-1101, San Diego, CA, 2016

  19. [19]

    Mazaheri, M

    A. Mazaheri, M. Ricchiuto, H. Nishikawa, A first-order hyperbolic system approach for dispersion, J. Comput. Phys. 321 (2016) 593–605

  20. [20]

    H. Baty, H. Nishikawa, Hyperbolic method for magnetic reconnection process in steady state mag- netohydrodynamics, Mon. Not. R. Astron. Soc. 459 (2016) 624–637

  21. [21]

    R. A. Watson, W. Trojak, P. G. Tucker, A simple flux reconstruction approach to solving a poisson equation to find wall distances for turbulence modelling, in: AIAA 2018 Fluid Dynamics Conference, AIAA Paper 2018-4261, Atlanta, Georgia, 2018

  22. [22]

    J. Lou, X. Liu, H. Luo, H. Nishikawa, Reconstructed discontinuous galerkin methods for hyperbolic diffusion equations on unstructured grids, in: 55th AIAA Aerospace Sciences Meeting, AIAA Paper, Vol. 310, 2017

  23. [23]

    J. Lou, X. Liu, H. Luo, H. Nishikawa, Reconstructed discontinuous galerkin methods for hyperbolic diffusion equations on unstructured grids, Communications in Computational Physics 25 (5) (2019). doi:10.4208/cicp.oa-2017-0186

  24. [24]

    A. S. Chamarthi, K. Komurasaki, R. Kawashima, High-order upwind and non-oscillatory approach for steady state diffusion, advection–diffusion and application to magnetized electrons, Journal of Computational Physics 374 (2018) 1120–1151

  25. [25]

    Nishikawa, Y

    H. Nishikawa, Y. Nakashima, Dimensional scaling and numerical similarity in hyperbolic method for diffusion, Journal of Computational Physics 355 (2018) 121–143

  26. [26]

    J. Lou, L. Li, H. Luo, H. Nishikawa, Explicit hyperbolic reconstructed discontinuous Galerkin methods for time-dependent problems, in: AIAA 2018 Fluid Dynamics Conference, AIAA Paper 2018-4270, Atlanta, Georgia, 2018

  27. [27]

    Nishikawa, Y

    H. Nishikawa, Y. Nakashima, Dimensional scaling and numerical similarity in hyperbolic method for diffusion, J. Comput. Phys. 355 (2018) 121–143

  28. [28]

    Y. Liu, H. Nishikawa, Third-order inviscid and second-order hyperbolic Navier-Stokes solvers for three-dimensional unsteady inviscid and viscous flows, in: 55th AIAA Aerospace Sciences Meeting, AIAA Paper 2017-0738, Grapevine, Texas, 2017

  29. [29]

    Nishikawa, New-generation hyperbolic Navier-Stokes schemes: O(1/h) speed-up and accurate viscous/heat fluxes, in: Proc

    H. Nishikawa, New-generation hyperbolic Navier-Stokes schemes: O(1/h) speed-up and accurate viscous/heat fluxes, in: Proc. of 20th AIAA Computational Fluid Dynamics Conference, AIAA Paper 2011-3043, Honolulu, Hawaii, 2011

  30. [30]

    Nishikawa, On hyperbolic method for diffusion with discontinuous coefficients, J

    H. Nishikawa, On hyperbolic method for diffusion with discontinuous coefficients, J. Comput. Phys. 367 (2018) 102–108

  31. [31]

    J. Lou, L. Li, H. Luo, H. Nishikawa, First-order hyperbolic system based reconstructed discontinuous Galerkin methods for nonlinear diffusion equations on unstructured grids, in: 56th AIAA Aerospace Sciences Meeting, AIAA Paper 2018-2094, Kissimmee, Florida, 2018

  32. [32]

    L. Li, J. Lou, H. Luo, H. Nishikawa, A new formulation of hyperbolic Navier-Stokes solver based on finite volume method on arbitrary grids, in: AIAA 2018 Fluid Dynamics Conference, AIAA Paper 2018-4160, Atlanta, Georgia, 2018

  33. [33]

    Nishikawa, First-, second-, and third-order finite-volume schemes for diffusion, Journal of Com- putational Physics 256 (2014) 791–805

    H. Nishikawa, First-, second-, and third-order finite-volume schemes for diffusion, Journal of Com- putational Physics 256 (2014) 791–805

  34. [34]

    X. Deng, H. Zhang, Developing High-Order Weighted Compact Nonlinear Schemes, Journal of Computational Physics 165 (1) (2000) 22–44

  35. [35]

    Nonomura, K

    T. Nonomura, K. Fujii, Robust explicit formulation of weighted compact nonlinear scheme, Com- puters and Fluids 85 (2013) 8–18

  36. [36]

    S. K. Lele, Compact finite difference schemes with spectral-like resolution, Journal of computational physics 103 (1) (1992) 16–42

  37. [37]

    Nagarajan, S

    S. Nagarajan, S. K. Lele, J. H. Ferziger, A robust high-order compact method for large eddy simulation, Journal of Computational Physics 191 (2) (2003) 392–419

  38. [38]

    B. J. Boersma, A 6th order staggered compact finite difference method for the incompressible Navier-Stokes and scalar transport equations, Journal of Computational Physics 230 (12) (2011) 4940–4954. 31 URL http://dx.doi.org/10.1016/j.jcp.2011.03.014

  39. [39]

    Jiang, C.-W

    G.-S. Jiang, C.-W. Shu, Efficient Implementation of Weighted ENO Schemes, Journal of Compu- tational Physics 126 (126) (1995) 202–228

  40. [40]

    X. Deng, Y. Jiang, M. Mao, H. Liu, S. Li, G. Tu, A family of hybrid cell-edge and cell-node dissipative compact schemes satisfying geometric conservation law, Computers & Fluids 116 (2015) 29–45

  41. [41]

    Titarev, E

    V. Titarev, E. Toro, Finite-volume weno schemes for three-dimensional conservation laws, Journal of Computational Physics 201 (1) (2004) 238 – 260

  42. [42]

    Nonomura, S

    T. Nonomura, S. Morizawa, H. Terashima, S. Obayashi, K. Fujii, Numerical (error) issues on com- pressible multicomponent flows using a high-order differencing scheme: Weighted compact nonlinear scheme, Journal of Computational Physics 231 (8) (2012) 3181–3210

  43. [43]

    M. L. Wong, S. K. Lele, High-order localized dissipation weighted compact nonlinear scheme for shock- and interface-capturing in compressible flows, Vol. 339, Elsevier Inc., 2017. doi:10.1016/j. jcp.2017.03.008

  44. [44]

    Nishikawa, P

    H. Nishikawa, P. Roe, Y. Suzuki, B. van Leer, A general theory of local preconditioning and its application to 2d ideal mhd equations, in: 16th AIAA Computational Fluid Dynamics Conference, p. 3704

  45. [45]

    van Es, B

    B. van Es, B. Koren, H. J. de Blank, Finite-volume scheme for anisotropic diffusion, Journal of Computational Physics 306 (2016) 422–442

  46. [46]

    Lafleur, S

    T. Lafleur, S. Baalrud, P. Chabert, Theory for the anomalous electron transport in hall effect thrusters. ii. kinetic model, Physics of Plasmas 23 (5) (2016) 053503

  47. [47]

    Kawashima, K

    R. Kawashima, K. Komurasaki, T. Sch¨ onherr, A hyperbolic-equation system approach for magne- tized electron fluids in quasi-neutral plasmas, Journal of Computational Physics 284 (2015) 59 – 69

  48. [48]

    Tan, C.-W

    S. Tan, C.-W. Shu, Inverse lax-wendroff procedure for numerical boundary conditions of conserva- tion laws, Journal of Computational Physics 229 (21) (2010) 8144 – 8166

  49. [49]

    S. Tan, C. W. Shu, A high order moving boundary treatment for compressible inviscid flows, Journal of Computational Physics 230 (15) (2011) 6023–6036

  50. [50]

    Kawashima, K

    R. Kawashima, K. Komurasaki, T. Sch¨ onherr, A flux-splitting method for hyperbolic-equation system of magnetized electron fluids in quasi-neutral plasmas, Journal of Computational Physics 310 (2016) 202–212

  51. [51]

    Nonomura, N

    T. Nonomura, N. Iizuka, K. Fujii, Freestream and vortex preservation properties of high-order WENO and WCNS on curvilinear grids, Computers and Fluids 39 (2) (2010) 197–214. doi:10. 1016/j.compfluid.2009.08.005

  52. [52]

    H. T. Huynh, A flux reconstruction approach to high-order schemes including discontinuous galerkin methods, in: 18th AIAA Computational Fluid Dynamics Conference, 2007, p. 4079

  53. [53]

    Borges, M

    R. Borges, M. Carmona, B. Costa, W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, Journal of Computational Physics 227 (6) (2008) 3191– 3211. 32