Multiphase fluid‐solid coupled analysis of shock‐bubble‐stone interaction in shockwave lithotripsy. (13th January 2017)

• Record Type: Journal Article
• Authors: Wang, Kevin G.
• Is Part Of: International journal for numerical methods in biomedical engineering. Volume 33:Number 10(2017:Oct.)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2017-01-13
• Subjects: cavitation -- embedded boundary method -- FIVER -- fluid‐solid interaction -- multiphase flow -- shock‐induced fracture -- shockwave lithotripsy; Biomedical engineering -- Periodicals; Imaging systems in medicine -- Periodicals; Numerical analysis -- Periodicals; Engineering mathematics -- Periodicals; 610.28
• DOI: 10.1002/cnm.2855
• Languages: English
• ISSNs: 2040-7939

Abstract

Abstract: A novel multiphase fluid‐solid–coupled computational framework is applied to investigate the interaction of a kidney stone immersed in liquid with a lithotripsy shock wave (LSW) and a gas bubble near the stone. The main objective is to elucidate the effects of a bubble in the shock path to the elastic and fracture behaviors of the stone. The computational framework couples a finite volume 2‐phase computational fluid dynamics solver with a finite element computational solid dynamics solver. The surface of the stone is represented as a dynamic embedded boundary in the computational fluid dynamics solver. The evolution of the bubble surface is captured by solving the level set equation. The interface conditions at the surfaces of the stone and the bubble are enforced through the construction and solution of local fluid‐solid and 2‐fluid Riemann problems. This computational framework is first verified for 3 example problems including a 1D multimaterial Riemann problem, a 3D shock‐stone interaction problem, and a 3D shock‐bubble interaction problem. Next, a series of shock‐bubble‐stone–coupled simulations are presented. This study suggests that the dynamic response of a bubble to LSW varies dramatically depending on its initial size. Bubbles with an initial radius smaller than a threshold collapse within 1 μ s after the passage of LSW, whereas larger bubbles do not. For a typical LSW generated by an electrohydraulic lithotripter ( p m a x = 35.0MPa, p m i n =− 10.1MPa), this threshold is approximately 0.12mm. Moreover, this study suggests that a noncollapsing bubble imposes a negative effect on stone fracture as it shields part of the LSW from the stone. On the other hand, a collapsing bubble may promote fracture on the proximal surface of the stone, yet hinder fracture from stone interior. Abstract : A 3D computational fluid dynamics (CFD)‐computational solid dynamics (CSD) coupled computational framework is applied to investigate the interaction of model kidney stones immersed in liquid with a lithotripsy shock wave (LSW) and a gas bubble near the stone. The simulation results suggest that bubbles smaller than a certain threshold may collapse violently during the process, thereby promoting fracture on stone surface, yet hindering fracture in the interior.