Simulation of laser propagation in a plasma with a frequency wave equation

Authors:
S. Desroziers;F. Nataf;R. Sentis
Affiliations:
CEA, France;Laboratoire J-L-Lions, CNRS UMR 7598, Université Paris VI, Paris, France;CEA/Bruyeres, 91680 Bruyeres, France
Venue:
Journal of Computational Physics
Year:
2008

Citing 11
Cited 1

Matrix analysis

Matrix analysis
A review on the inverse of symmetric tridiagonal and block tridiagonal matrices

SIAM Journal on Matrix Analysis and Applications
A perfectly matched layer for the absorption of electromagnetic waves

Journal of Computational Physics
A domain decomposition method for the Helmholtz equation and related optimal control problems

Journal of Computational Physics
A Parallel Fast Direct Solver for Block Tridiagonal Systems with Separable Matrices of Arbitrary Dimension

SIAM Journal on Scientific Computing
Separation-of-variables as a preconditioner for an iterative Helmholtz solver

Applied Numerical Mathematics
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems
Orthogonal Eigenvectors and Relative Gaps

SIAM Journal on Matrix Analysis and Applications
A geometric optics method for high-frequency electromagnetic fields computations near fold caustics: part II. the energy

Journal of Computational and Applied Mathematics
Comparison of multigrid and incomplete LU shifted-Laplace preconditioners for the inhomogeneous Helmholtz equation

Applied Numerical Mathematics
Simulation of the Paraxial Laser Propagation Coupled with Hydrodynamics in 3D Geometry

Journal of Scientific Computing

Simulation of laser beam propagation with a paraxial model in a tilted frame

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

The aim of this work is to perform numerical simulations of the propagation of a laser in a plasma. At each time step, one has to solve a Helmholtz equation in a domain which consists in some hundreds of millions of cells. To solve this huge linear system, we use an iterative Krylov method preconditioned by a separable matrix. The corresponding linear system is solved with a block cyclic reduction method. Some enlightenments on the parallel implementation are also given. Lastly, numerical results are presented including some features concerning the scalability of the numerical method on a parallel architecture.