Hi Kevin,
Doing the long pulse is straightforward using a time-varying source in kspaceFirstOrder3D
, but k-Wave doesn't have a coupled heat-acoustic solver, and the function describing the temperature-dependence of the Grueneisen parameter would have to be somehow incorporated too.
Using k-Wave you would have to run the heat solver, kWaveDiffusion
, to get the temperature changes (incorporating heat diffusion), use a known function for the temperature-dependence of the Grueneisen parameter to get the time-dependent photoacoustic amplitude, and then use a time-varying source in kspaceFirstOrder3D
, scaled by this time-varying Grueneisen parameter, for the acoustic simulation.
This would still assume that the acoustic properties (sound speed, density) do not change with temperature, which isn't quite right.
Best wishes
Ben