Numerical simulation of wave propagation in frozen porous media

We propose a numerical algorithm for simulating wave propagation in frozen porous media. The original theory assumes that there is no direct contact between solid grains and ice (we include the grain-ice interaction), and predicts three compressional waves and two shear waves. The wavefield is obtained with a grid-method based on the Fourier differential operator and a 4th-order Runge-Kutta time-integration algorithm. Since the presence of slow diffusive waves makes the differential equations stiff, a time-splitting integration algorithm is used to solved the stiff part with an analytical technique.