Hi,

I have a question about the relationship between elastic simulation time and coefficients.

For instance, when I use below constants:

-----------------------------------------

% water

cp1 = 1540; % compressional wave speed [m/s]

cs1 = 0; % shear wave speed [m/s]

rho1 = 1000; % density [kg/m^3]

alpha0_p1 = 0; % compressional wave absorption [dB/(MHz^2 cm)]

alpha0_s1 = 0; % shear wave absorption [dB/(MHz^2 cm)]

% skull

cp2 = 3000; % compressional wave speed [m/s]

cs2 = 1500; % shear wave speed [m/s]

rho2 = 1850; % density [kg/m^3]

alpha0_p2 = 3; % compressional wave absorption [dB/(MHz^2 cm)]

alpha0_s2 = 6; % shear wave absorption [dB/(MHz^2 cm)]

% create the time array

cfl = 0.13;

t_end = 40e-6;

-----------------------------------------

estimated simulation time 1hours 15min 41.4135s...

-----------------------------------------

When I change the skull's attenuation values to alpha0_p2 = 8.83, alpha0_s2 = 19.5 (belongs to your paper(Modeling power law absorption and dispersion in viscoelastic solids using a split-field

and the fractional Laplacian)). Unfortunately, simulation gives a black screen. I think I have to change the cfl number. But, when I decrease the cfl coefficient, simulation time takes much longer.

How can I minimize the simulation time under these circumstances?

Best,

Michael