Acoustic simulation with high frequency (ex. 1 MHz)
-
This is correct but if there is any air in that area you should remove it. The maximum grid step is calculated by looking at all the materials in the simulation and finding the one with the lowest speed of sound (air, around 300 m/s, everything else should be around 1500 m/s except for bone which is higher), and then dividing that by the simulation frequency (and some factor indicating minimum number of grid points per wavelength).
-
But when I remove the head from the grid settings, the material considered for automatic calculation of the max cell step is water (in my case, f=500 kHz, so max step = 0.3 mm). So, I suppose that in voxeling the module takes into account this cell size as the maximum. Correct?
I also tried to include the head in my simulation in order to see the difference in the grid. In this case, the automatic setting for the whole head becomes 2 mm, but in both cases I have the same total number of cells. It seems that the module considers a fixed voxelization of Duke head. That could be possible? -
I'm not sure, but the elements in the grid settings define the size of the grid and also where some grid lines should definitely be (the red lines in the grid). The actual resolution should be calculated based on the materials in the whole simulation domain within the total grid.
Did you mean to say that your default grid settings went from 0.3 mm to 2 mm?
-
I usually change the automatic settings in manual (by including the transducer and the water background and removing Duke): max step 0.3, resolution 0.45 and priority 0. In that way, all the objects in my manual settings (transducer and water background) are set at 0.3 mm. Then in voxeling I include also Duke, but in this way I don't know what is the step used for voxelizing the head.
If I include the head object also in gridding automatically S4L set its cells at 2 mm as maximum size, and the rest of the domain at 0.3 mm.
Both cases result in the same total number of cells in the whole computational domain.
For answering your question: when I include also the head as a grid object, the automatic grid settings of the ONLY head are automatically set at 2 mm (with no action by myself). -
Oh I see, but if you use manual grid settings then you should really be careful that the resolution is fine enough. I would suggest making some tests: make your resolution finer and see how much your results change. Hopefully they don't change much (maybe within a few percent).
-
Moving to manual settings, I kept the same step and resolution as the default automatic values. I only changed the priority from 50 to 0 in order to remove the red lines and reduce the number of cells I had. (that was a suggestion from a colleague of mine). Anyway, I will make some tests. Thank you for your suggestions.
-
@montanaro It is possible to do the modeling HIFU brain tumor ablation in Sim4Life.
If setting the operating frequency over 1 MHz, the computational size will be much bigger. -
-
Yes, unfortunately there are no real work arounds for that.. You need to simulate at least 10 points per wavelength (speed of sound of medium / frequency), - this is what your minimum spatial resolution should be to have results that make sense.
@montanaro Thanks for your clarification.