Higher-order variational finite difference schemes for solving 3-D paraxial wave equations using splitting techniques

Eliane Bécache (1) , Francis Collino (1) , Patrick Joly (1)

1 INRIA Rocquencourt

Numerical schemes for solving 3-D paraxial equations are constructed using splitting techniques. The solution can be reduced to a series of 2-D paraxial equations in each direction of splitting. The discretization along the depth is based on higher-order conservative schemes. The discretization along the transverse variables is based on higher-order finite difference variational schemes. Numerical experiments illustrate the advantages of higher-order schemes, which are much less dispersive, even for a small number of discretization points per wavelength.

Domaines: Analyse numérique [math.NA]

Eliane Bécache, Francis Collino, Patrick Joly. Higher-order variational finite difference schemes for solving 3-D paraxial wave equations using splitting techniques. Wave Motion, 2000, ⟨10.1016/S0165-2125(99)00038-4⟩. ⟨hal-01856321⟩