scholarly journals Effective Source Term Discretizations for Higher Accuracy Finite Volume Discretization of Parabolic Equations

Author(s):  
Yaw Kyei

A finite volume method is applied to develop space-time discretizations for parabolic equations based on an equation error method.A space-time expansion of the local equation error based on flux integral formulation of the equation is first designed using a desiredframework of neighboring quadrature points for the solution and local source terms. The quadrature weights are then determined through aminimization process for the error which constitutes all local compact fluxes about each centroid within the computational domain.In utilizing a local source term distribution to account for diffusive fluxes, the right minimizing quadrature weights and collocationpoints including subgrid points for the source terms may be determined and optimized for higher accuracies as well as robust higher-ordercomputational convergence. The resulting local residuals form a more complete description of the truncation errors which are then utilizedto assess the computational performances of the resulting schemes. The effectiveness of the discretization method is demonstrated by theresults and analysis of the schemes.

Author(s):  
T Thomas ◽  
C Pfrommer ◽  
R Pakmor

Abstract We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh Arepo code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretised using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of Arepo. The interaction of CRs and gyroresonant Alfvén waves is described by short-timescale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magneto-hydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetised discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magneto-hydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.


Sign in / Sign up

Export Citation Format

Share Document