arrow
Volume 29, Issue 3
A Simple Solver for the Two-Fluid Plasma Model Based on PseudoSpectral Time-Domain Algorithm

Benoit Morel, Remo Giust, Kazem Ardaneh & Francois Courvoisier

Commun. Comput. Phys., 29 (2021), pp. 955-978.

Published online: 2021-01

Export citation
  • Abstract

We present a solver of 3D two-fluid plasma model for the simulation of short-pulse laser interactions with plasma. This solver resolves the equations of the two-fluid plasma model with ideal gas closure. We also include the Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral Time-Domain (PSTD) method to solve Maxwell's curl equations. We use a Strang splitting to integrate Euler equations with source term: while Euler equations are solved with a composite scheme mixing Lax-Friedrichs and Lax-Wendroff schemes, the source term is integrated with a fourth-order Runge-Kutta scheme. This two-fluid plasma model solver is simple to implement because it only relies on finite difference schemes and Fast Fourier Transforms. It does not require spatially staggered grids. The solver was tested against several well-known problems of plasma physics. Numerical simulations gave results in excellent agreement with analytical solutions or with previous results from the literature.

  • AMS Subject Headings

35L02, 35L03

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-29-955, author = {Morel , BenoitGiust , RemoArdaneh , Kazem and Courvoisier , Francois}, title = {A Simple Solver for the Two-Fluid Plasma Model Based on PseudoSpectral Time-Domain Algorithm}, journal = {Communications in Computational Physics}, year = {2021}, volume = {29}, number = {3}, pages = {955--978}, abstract = {

We present a solver of 3D two-fluid plasma model for the simulation of short-pulse laser interactions with plasma. This solver resolves the equations of the two-fluid plasma model with ideal gas closure. We also include the Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral Time-Domain (PSTD) method to solve Maxwell's curl equations. We use a Strang splitting to integrate Euler equations with source term: while Euler equations are solved with a composite scheme mixing Lax-Friedrichs and Lax-Wendroff schemes, the source term is integrated with a fourth-order Runge-Kutta scheme. This two-fluid plasma model solver is simple to implement because it only relies on finite difference schemes and Fast Fourier Transforms. It does not require spatially staggered grids. The solver was tested against several well-known problems of plasma physics. Numerical simulations gave results in excellent agreement with analytical solutions or with previous results from the literature.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2020-0117}, url = {http://global-sci.org/intro/article_detail/cicp/18572.html} }
TY - JOUR T1 - A Simple Solver for the Two-Fluid Plasma Model Based on PseudoSpectral Time-Domain Algorithm AU - Morel , Benoit AU - Giust , Remo AU - Ardaneh , Kazem AU - Courvoisier , Francois JO - Communications in Computational Physics VL - 3 SP - 955 EP - 978 PY - 2021 DA - 2021/01 SN - 29 DO - http://doi.org/10.4208/cicp.OA-2020-0117 UR - https://global-sci.org/intro/article_detail/cicp/18572.html KW - Two-fluid plasma model, 3D Hydrodynamic code, Lax-Wendroff, composite scheme, PSTD, laser-plasma interaction. AB -

We present a solver of 3D two-fluid plasma model for the simulation of short-pulse laser interactions with plasma. This solver resolves the equations of the two-fluid plasma model with ideal gas closure. We also include the Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral Time-Domain (PSTD) method to solve Maxwell's curl equations. We use a Strang splitting to integrate Euler equations with source term: while Euler equations are solved with a composite scheme mixing Lax-Friedrichs and Lax-Wendroff schemes, the source term is integrated with a fourth-order Runge-Kutta scheme. This two-fluid plasma model solver is simple to implement because it only relies on finite difference schemes and Fast Fourier Transforms. It does not require spatially staggered grids. The solver was tested against several well-known problems of plasma physics. Numerical simulations gave results in excellent agreement with analytical solutions or with previous results from the literature.

Morel , BenoitGiust , RemoArdaneh , Kazem and Courvoisier , Francois. (2021). A Simple Solver for the Two-Fluid Plasma Model Based on PseudoSpectral Time-Domain Algorithm. Communications in Computational Physics. 29 (3). 955-978. doi:10.4208/cicp.OA-2020-0117
Copy to clipboard
The citation has been copied to your clipboard