A numerical tour of wave propagation |

Another approach to improve the behavior of the discrete PML at grazing incidence consists in modifying the complex coordinate transform used classically in the PML to introduce a frequency-dependent term that implements a Butterworth-type filter in the layer. This approach has been developed for Maxwellâ€™s equations named convolutional PML (CPML) (Roden and Gedney, 2000) or complex frequency shifted-PML (CFS-PML). The key idea is that for waves whose incidence is close to normal, the presence of such a filter changes almost nothing because absorption is already almost perfect. But for waves with grazing incidence, which for geometrical reasons do not penetrate very deep in the PML, but travel there a longer way in the direction parallel to the layer, adding such a filter will strongly attenuate them and will prevent them from leaving the PML with significant energy .

Define , . Then the acoustic wave equation reads

To combine the absorbing effects into the acoustic equation, we merely need to combine two convolution terms into the above equations:

(25) |

where , are the convolution terms of and ; , are the convolution terms of and . These convolution terms can be computed via the following relation:

(26) |

where and . More details about the derivation of C-PML, the interested readers are referred to Collino and Tsogka (2001) and Komatitsch and Martin (2007).

A numerical tour of wave propagation |

2021-08-31