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

Elastic Wave Simulations - Closed 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.

Circle Paths - The acoustic source was scanned at a constant linear speed along a circular trajectories with two different diameters - 1 cm and 3 cm (indicated in the filename) and for at least 2 pattern repetitions. The linear scanning speed ranged from 2 to 20 m/s and is designated in the filename.

Square Paths - The acoustic source was scanned at a constant speed along square trajectories with two different edge lengths - 1 cm and 3 cm (indicated in the filename) and for at least 2 pattern repetitions. The scan speed ranged from 2 m/s to 10 m/s and is designated in the filename.

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

nCycles - Number of pattern repetitions

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)