This repository includes scripts that were used to simulate the impact of cavitating bubbles in amplifying incident ultrasound waves on nearby cells.
This system requires baseline MATLAB, it was programmed on 2019, but may run on earlier versions.
The Keller-Miksis model is used to simulate the pressure amplification of an incident ultrasound signal as a function of initial bubble radius and distance from center of cavitating bubble. This is done through two steps: 1) Running Keller-Miksis simulation for all bubble radii, 2) Taking FFT of all bubble radius vs time datasets generated
Use this function, with modified parameters to represent the simulation desired, to generate the Keller-Miksis radius vs time datasets
KellerMiksis_batch.m
For testing purposes, a single Keller-Miksis simulation can be done using:
KellerMiksis_single.m
The results of each of these simulations is also output graphically by these functions as:
After the previous step, run this function, with modified parameters to represent the simulation desired, to generate the output figure
FarFieldPressure.m
The following figure is generated which shows decibel gain in the incident ultrasound frequency FFT of the pressure emissions from cavitating bubbles compared to the stimulating ultrasound signal. Initial bubble radius on the y-axis and distance from the bubble on the x-axis. In red regions the cavitating bubble’s emitted pressure has a larger FFT signal at the stimulation frequency than the ultrasound signal.
For a graph of the natural frequency of bubbles as a function of initial radius, run the following function:
NatFreq.m
Johansen, K., J.H. Song, and P. Prentice, Validity of the Keller-Miksis equation for "non-stable" cavitation and the acoustic emissions generated. 2017 Ieee International Ultrasonics Symposium (Ius), 2017.
Christopher E. Brennen, CAVITATION AND BUBBLE DYNAMICS, https://authors.library.caltech.edu/25017/4/chap4.htm#L2
Currently contributing is not suppported, please see future versions at https://github.com/drmittelstein/oncotripsy to determine whether this changes.
Please see available versions at https://github.com/drmittelstein/oncotripsy
- David Reza Mittelstein - "Modifying ultrasound waveform parameters to control, influence, or disrupt cells" Caltech Doctorate Thesis in Medical Engineering
- Acknowledgements to my colleagues in Gharib, Shapiro, and Colonius lab at Caltech who helped answer questions involved in the development of these scripts.