Application of High Dimensional B-Spline Interpolation in Solving the Gyro-Kinetic Vlasov Equation Based on Semi-Lagrangian Method. (6th July 2017)

• Record Type: Journal Article
• Authors: Xiao, Xiaotao; Ye, Lei; Xu, Yingfeng; Wang, Shaojie
• Is Part Of: Communications in computational physics. Volume 22:Number 3(2017:Sep.)
• Page Start: 789
• Page End: 802
• Publication Date: 2017-07-06
• Subjects: 65D05, -- 65Z05, -- 68U20, -- 82D10; High dimensional interpolation, -- semi-Lagrangian, -- gyro-kinetic simulation, -- numerical Lie transform; Mathematical physics -- Data processing -- Periodicals; Physics -- Data processing -- Periodicals; 530.150285
• DOI: 10.4208/cicp.OA-2016-0092
• Languages: English
• ISSNs: 1815-2406

Abstract

Abstract: The computation efficiency of high dimensional (3D and 4D) B-spline interpolation, constructed by classical tensor product method, is improved greatly by precomputing the B-spline function. This is due to the character of NLT code, i.e. only the linearised characteristics are needed so that the unperturbed orbit as well as values of the B-spline function at interpolation points can be precomputed at the beginning of the simulation. By integrating this fixed point interpolation algorithm into NLT code, the high dimensional gyro-kinetic Vlasov equation can be solved directly without operator splitting method which is applied in conventional semi-Lagrangian codes. In the Rosenbluth-Hinton test, NLT runs a few times faster for Vlasov solver part and converges at about one order larger time step than conventional splitting code.