Shear Shock Waves Mediate Haptic Holography via Focused Ultrasound - Elastic Wave Simulations (Open Scanning Paths)

Elastic Wave Simulations - Open Scanning Paths

This dataset is part of a larger repository (DOI: 10.5281/zenodo.5248082) which houses links to the data used in the publication "Shear Shock Waves Mediate Haptic Holography via Focused Ultrasound" (Reardon et al., 2023). If you use these simulated data please cite our publication (http://www.science.org/doi/10.1126/sciadv.adf2037) and the software package k-Wave (DOI: 10.1109/ULTSYM.2014.0037).

This dataset contains the normal shear surface velocity in a cylindrical slab of tissue-like material excited by an acoustic source with a Gaussian spatial profile simulated via a pseudo-spectral numerical method. The data is provided as .mat files. The files are separated by the type of scanning path, the scanning speed of the acoustic source, and the parameters of the scanning path. Details of the simulation parameters can be found in our publication.

Line Paths - The acoustic source scanned along a linear trajectory at speeds ranging from 2 m/s to 12 m/s (scanning speed is indicated in the filename).

Zigzag Paths - The acoustic source scanned along a zigzag path on the surface of the simulated medium with x-axis scanning speed vx = 3, 4, 5, 6 m/s. At all speeds, the ultrasound focus was modulated transverse to its primary motion direction at a speed, vy, of +-2.5 m/s yielding a zigzag path (2 cm path width). The x-axis scanning speed is designated in the filename.

Letter Paths - The acoustic source scanned the a trajectory in the shape of the letter "Z." Scanning speeds ranged from 2 m/s to 12 m/s (scanning speed is designated in the filename).

Focus Control Rate - The acoustic source scanned along a linear trajectory at 7 m/s but at different focus control sample rates fc. These paths amount to a courser sampling of the linear trajectory. In lieu of updating the location of the acoustic source at each timepoint in the simulation, we specified a rate at which the location of the acoustic source would be updated. We set fc to approximately 0.7, 1.4, and 4.2 kHz (designated at the end of the filename as VeryCoarse, Coarse, and Fine, respectively). (Compare with Line_07, which has the finest path sampling and an *fc* of approximately 200 kHz.)

Data Fields

surfaceData (NxNxM) - 3D array containing the normal shear velocity of the simulated medium (in m/s) on a NxN Cartesian grid of locations at M timepoints. The simulated tissue medium was cylindrical, so locations outside the circular top surface are NaN.

sourceSignals (NxNxQ) - 3D array containing the acoustic source distribution on the NxN Cartesian grid of locations used to excite the surface of the simulated tissue medium for Q timepoints.

sourceEnvelope (Qx1) - Vector containing the amplitude envelope that was applied to sourceSignals at each timestep Q

dt - The time between adjacent timepoints in seconds (i.e. fs = 1/dt)

dx/dy - The distance between adjacent grid locations in the x/y direction of the Cartesian grid (in m)