In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method with discontinuous Galerkin (DG) discretization for the physical and phase spaces for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum, and energy at the discrete level. The recently developed LoMaC low rank tensor algorithm (arXiv: 2207.00518) simultaneously evolves the macroscopic conservation laws of mass, momentum, and energy using the kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables. This paper is a generalization of our previous work, but with DG discretization to take advantage of its compactness and flexibility in handling boundary conditions and its superior accuracy in the long term. The algorithm is developed in a similar fashion as that for a finite difference scheme, by observing that the DG method can be viewed equivalently in a nodal fashion. With the nodal DG method, assuming a tensorized computational grid, one will be able to (i) derive differentiation matrices for different nodal points based on a DG upwind discretization of transport terms, and (ii) define a weighted inner product space based on the nodal DG grid points. The algorithm can be extended to the high dimensional problems by hierarchical Tucker (HT) decomposition of solution tensors and a corresponding conservative projection algorithm. In a similar spirit, the algorithm can be extended to DG methods on nodal points of an unstructured mesh, or to other types of discretization, e.g., the spectral method in velocity direction. Extensive numerical results are performed to showcase the efficacy of the method.
Wei Guo, Jannatul Ferdous Ema, Jing-Mei Qiu
. A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method with the Discontinuous Galerkin Method for the Vlasov Dynamics[J]. Communications on Applied Mathematics and Computation, 2024
, 6(1)
: 550
-575
.
DOI: 10.1007/s42967-023-00277-7
[1] Arnold, D.N., Brezzi, F., Cockburn, B., Marini, L.D.:Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Num. Anal. 39(5), 1749-1779 (2002)
[2] Cockburn, B., Luskin, M., Shu, C.-W., Süli, E.:Enhanced accuracy by post-processing for finite element methods for hyperbolic equations. Math. Comput. 72(242), 577-606 (2003)
[3] De Dios, B.A., Hajian, S.:High order and energy preserving discontinuous Galerkin methods for the Vlasov-Poisson system. arXiv:1209.4025 (2012)
[4] De Lathauwer, L., De Moor, B., Vandewalle, J.:A multilinear singular value decomposition. SIAM J. Matrix Anal. Appl. 21(4), 1253-1278 (2000)
[5] Gottlieb, S., Ketcheson, D.I., Shu, C.-W.:Strong Stability Preserving Runge-Kutta and Multistep Time Discretizations. World Scientific, NJ, United States (2011)
[6] Grasedyck, L.:Hierarchical singular value decomposition of tensors. SIAM J. Matrix Anal. Appl. 31(4), 2029-2054 (2010)
[7] Guo, W., Qiu, J.-M.:A low rank tensor representation of linear transport and nonlinear vlasov solutions and their associated flow maps. arXiv:2106.08834 (2021)
[8] Guo, W., Qiu, J.-M.:A conservative low rank tensor method for the Vlasov dynamics. arXiv:2106.08834 (2022)
[9] Guo, W., Qiu, J.-M.:A Local Macroscopic Conservative (LoMaC) low rank tensor method for the Vlasov dynamics. arXiv:2207.00518 (2022)
[10] Guo, W., Qiu, J.-M.:A low rank tensor representation of linear transport and nonlinear Vlasov solutions and their associated flow maps. J. Comput. Phys. 458, 111089 (2022)
[11] Hackbusch, W., Kühn, S.:A new scheme for the tensor representation. J. Fourier Anal. Appl. 15(5), 706-722 (2009)
[12] Mandal, J., Deshpande, S.:Kinetic flux vector splitting for Euler equations. Comput. Fluids 23(2), 447-478 (1994)
[13] Tucker, L.R.:Some mathematical notes on three-mode factor analysis. Psychometrika 31(3), 279-311 (1966)
[14] Xu, K., Martinelli, L., Jameson, A.:Gas-kinetic finite volume methods, flux-vector splitting, and artificial diffusion. J. Comput. Phys. 120(1), 48-65 (1995)
