{"title":"An Implicit, Asymptotic-Preserving and Energy-Charge-Conserving Method for the Vlasov-Maxwell System Near Quasi-Neutrality","authors":"Chuwen Ma, Shi Jin","doi":"10.4208/cicp.oa-2023-0133","DOIUrl":null,"url":null,"abstract":"An implicit, asymptotic-preserving and energy-charge-conserving (APECC)\nParticle-In-Cell (PIC) method is proposed to solve the Vlasov-Maxwell (VM) equations in the quasi-neutral regime. Charge conservation is enforced by particle orbital\naveraging and fixed sub-time steps. The truncation error depending on the number of sub-time steps is further analyzed. The temporal discretization is chosen by\nthe Crank-Nicolson method to conserve the discrete energy exactly. The key step in\nthe asymptotic-preserving iteration for the nonlinear system is based on a decomposition of the current density deduced from the Vlasov equation in the source of the\nMaxwell model. Moreover, we show that the convergence is independent of the quasineutral parameter. Extensive numerical experiments show that the proposed method\ncan achieve asymptotic preservation and energy-charge conservation.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"93 1","pages":""},"PeriodicalIF":2.6000,"publicationDate":"2024-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Communications in Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.4208/cicp.oa-2023-0133","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"PHYSICS, MATHEMATICAL","Score":null,"Total":0}
引用次数: 0
Abstract
An implicit, asymptotic-preserving and energy-charge-conserving (APECC)
Particle-In-Cell (PIC) method is proposed to solve the Vlasov-Maxwell (VM) equations in the quasi-neutral regime. Charge conservation is enforced by particle orbital
averaging and fixed sub-time steps. The truncation error depending on the number of sub-time steps is further analyzed. The temporal discretization is chosen by
the Crank-Nicolson method to conserve the discrete energy exactly. The key step in
the asymptotic-preserving iteration for the nonlinear system is based on a decomposition of the current density deduced from the Vlasov equation in the source of the
Maxwell model. Moreover, we show that the convergence is independent of the quasineutral parameter. Extensive numerical experiments show that the proposed method
can achieve asymptotic preservation and energy-charge conservation.
期刊介绍:
Communications in Computational Physics (CiCP) publishes original research and survey papers of high scientific value in computational modeling of physical problems. Results in multi-physics and multi-scale innovative computational methods and modeling in all physical sciences will be featured.