bioheat_transfer()
Description: This function uses the Bioheat Transfer Equation to calculate the temperature rise at each point in a medium.
Usage
temperature_rise = bioheat_transfer(input_pressure, medium, dx, dy, dz,iterations);
Arguments
- input_pressure A 3D input pressure field with the pressure at each point in Pa
- medium A FOCUS medium struct
- dx The distance between adjacent points in x in m
- dy The distance between adjacent points in y in m
- dz The distance between adjacent points in z in m
- iterations The number of iterations to use in the calculations
Output Parameters
- temperature_rise The temperature rise at each point in the original pressure field in degrees °C.
Notes
Note that the number of iterations required for a stable temperature field will depend on the spatial step size. The number of iterations should be increased until successive results do not significantly differ to make sure that the result is correct.