My setup includes a set of source points consisting 16 planes, each of 10x10 voxels placed on top of an imaginary sphere. I have also a set of sensors from Z = 0 to Z = top of the Sphere passing through its center. The sensor recordings are used to compute the propagation delay and from there the speed.
When there is only Water as the medium, the propagation is matching to 1500 m/s which I set for water.
In the next experiment I add another layer of medium using below lines:
medium.sound_speed = 1500*ones(kgrid.Nx, kgrid.Ny, kgrid.Nz);
medium.sound_speed(Volume==1) = 2900;
medium.density = 1000*ones(kgrid.Nx, kgrid.Ny, kgrid.Nz);
medium.density(Volume==1) = 2150;
medium.alpha_power = 1.5;
medium.alpha_coeff = 0.75*ones(kgrid.Nx, kgrid.Ny, kgrid.Nz);
medium.alpha_coeff(Volume==1) = 20.4307;
medium.BonA = 6*ones(kgrid.Nx, kgrid.Ny, kgrid.Nz);
medium.BonA(Volume==1) = 10;
in which Volume is a logical variable of size (Nx,Ny,Nz) and indicates an object which is designated using 1s and all other entries are 0.
The problem is after adding this layer of new material, the propagation speed slows down even before wave reach the new layer. The time signal peaks reach to the sensors, even for those in the first medium is exactly doubled and even though the sound speed for the second medium is set much higher than the first medium.
I have investigated this issue and if I comment this line:
medium.sound_speed(Volume==1) = 2900;
and leave all other lines there, the issue does not appear. Though this means I have neglected the second layer.
I appreciate if K-wave team can have a look at this issue, and I wonder if anybody else has measured propagation delays in heterogeneous mediums to see if it matches to what it should be.
Thanks for your times,
Maryam