41 research outputs found
Experiments and simulations demonstrating the rapid ultrasonic rewarming of frozen beef cryovials
The development of methods to safely rewarm large volume cryopreserved
biological samples remains a barrier to the widespread adoption of
cryopreservation. Here, experiments and simulations were performed to
demonstrate that ultrasound can increase rewarming rates relative to thermal
conduction alone. An ultrasonic rewarming setup based on a custom 444 kHz
tubular piezoelectric transducer was designed, characterized, and tested with 2
mL cryovials filled with frozen ground beef. Rewarming rates were characterized
in the -20C to 5C range. Thermal conduction-based rewarming
was compared to thermal conduction plus ultrasonic rewarming, demonstrating a
ten-fold increase in rewarming rate when ultrasound was applied. The maximum
recorded rewarming rate with ultrasound was 57C per minute,
approximately 2.5 times faster than with thermal conduction alone. Coupled
acoustic and thermal simulations were developed and showed good agreement with
the heating rates demonstrated experimentally and were also used to demonstrate
spatial heating distributions with small (3C) temperature
differentials throughout the sample when the sample was below 0C. The
experiments and simulations performed in this work demonstrate the potential
for ultrasound as a rewarming method for cryopreserved tissues, as faster
rewarming rates may improve the viability of cryopreserved tissues and reduce
the time needed for cells to regain normal function.Comment: 13 pages, 11 figure
On-the-Fly Calculation of Time-Averaged Acoustic Intensity in Time-Domain Ultrasound Simulations Using a k-Space Pseudospectral Method
OBJECTIVE: This paper presents a method to calculate the average acoustic intensity during ultrasound simulation using a new approach that exploits compression of intermediate results. METHODS: One of the applications of high-intensity focused ultrasound (HIFU) simulations is the calculation of the thermal dose, which indicates the amount of tissue destroyed using a state-of-the-art k-space pseudospectral method. The thermal simulation is preceded by the calculation of the average intensity within the acoustic simulation. Due to the time staggering between the particle velocity and the acoustic pressure used in such simulations, the average intensity calculation is typically executed offline after the acoustic simulation consuming both disk space and time (the data can spread over terabytes). Our new approach calculates the average intensity during the acoustic simulation using the output coefficients of a new compression method which enables resolving the time staggering on-the-fly with huge disk space savings. To reduce RAM requirements, the article also presents a new 40-bit method for encoding compression complex coefficients. RESULTS: Experimental numerical simulations with the proposed method have shown that disk space requirements are up to 99 % lower. The simulation speed was not significantly affected by the approach and the compression error did not affect the prediction accuracy of the thermal dose. CONCLUSION: From the standpoint of supercomputers, the new approach is significantly more economical. SIGNIFICANCE: Saving computing resources increases the chances of real use of acoustic simulations in practice. The method can be applied to signals of a similar character, e.g., for electromagnetic radio waves
A Learned Born Series for Highly-Scattering Media
A new method for solving the wave equation is presented, called the learned
Born series (LBS), which is derived from a convergent Born Series but its
components are found through training. The LBS is shown to be significantly
more accurate than the convergent Born series for the same number of
iterations, in the presence of high contrast scatterers, while maintaining a
comparable computational complexity. The LBS is able to generate a reasonable
prediction of the global pressure field with a small number of iterations, and
the errors decrease with the number of learned iterations.Comment: 6 pages, 1 figur
A Helmholtz equation solver using unsupervised learning: Application to transcranial ultrasound
Transcranial ultrasound therapy is increasingly used for the non-invasive
treatment of brain disorders. However, conventional numerical wave solvers are
currently too computationally expensive to be used online during treatments to
predict the acoustic field passing through the skull (e.g., to account for
subject-specific dose and targeting variations). As a step towards real-time
predictions, in the current work, a fast iterative solver for the heterogeneous
Helmholtz equation in 2D is developed using a fully-learned optimizer. The
lightweight network architecture is based on a modified UNet that includes a
learned hidden state. The network is trained using a physics-based loss
function and a set of idealized sound speed distributions with fully
unsupervised training (no knowledge of the true solution is required). The
learned optimizer shows excellent performance on the test set, and is capable
of generalization well outside the training examples, including to much larger
computational domains, and more complex source and sound speed distributions,
for example, those derived from x-ray computed tomography images of the skull.Comment: 23 pages, 13 figure
High Resolution 3D Ultrasonic Breast Imaging by Time-Domain Full Waveform Inversion
Ultrasound tomography (UST) scanners allow quantitative images of the human
breast's acoustic properties to be derived with potential applications in
screening, diagnosis and therapy planning. Time domain full waveform inversion
(TD-FWI) is a promising UST image formation technique that fits the parameter
fields of a wave physics model by gradient-based optimization. For high
resolution 3D UST, it holds three key challenges: Firstly, its central building
block, the computation of the gradient for a single US measurement, has a
restrictively large memory footprint. Secondly, this building block needs to be
computed for each of the measurements, resulting in a massive
parallel computation usually performed on large computational clusters for
days. Lastly, the structure of the underlying optimization problem may result
in slow progression of the solver and convergence to a local minimum. In this
work, we design and evaluate a comprehensive computational strategy to overcome
these challenges: Firstly, we introduce a novel gradient computation based on
time reversal that dramatically reduces the memory footprint at the expense of
one additional wave simulation per source. Secondly, we break the dependence on
the number of measurements by using source encoding (SE) to compute stochastic
gradient estimates. Also we describe a more accurate, TD-specific SE technique
with a finer variance control and use a state-of-the-art stochastic LBFGS
method. Lastly, we design an efficient TD multi-grid scheme together with
preconditioning to speed up the convergence while avoiding local minima. All
components are evaluated in extensive numerical proof-of-concept studies
simulating a bowl-shaped 3D UST breast scanner prototype. Finally, we
demonstrate that their combination allows us to obtain an accurate 442x442x222
voxel image with a resolution of 0.5mm using Matlab on a single GPU within 24h