Biomechanical Simulation using Supercomputer for Predictive Medicine
0 likes931 views
1) Biomechanical simulations using supercomputers can help predict the outcomes of new medical therapies like high intensity focused ultrasound (HIFU) for treating cancers. 2) Numerical models incorporate patient medical imaging data and simulate how ultrasound waves propagate and create focal points of high pressure for thermal ablation of tumors, accounting for interference at tissue interfaces. 3) Large supercomputers with over 640,000 cores and capable of 10 petaflops are needed to model complex fluid-structure interactions in living tissues at multiple scales using fully Eulerian methods.
![Comparison of the focal point
Without phase delay
Pmax=1.7MPa
0
2
Pressure [MPa]
With phase delay
Pmax=2.2MPa
0
2
Pressure [MPa]](https://image.slidesharecdn.com/takagisbiomechanicalsimulation-131111210456-phpapp02/85/Biomechanical-Simulation-using-Supercomputer-for-Predictive-Medicine-6-320.jpg)
Biomechanical Simulation using Supercomputer for Predictive Medicine
- 1. Biomechanical Simulation using Supercomputer for Predictive Medicine Shu TAKAGI Department of Mechanical Engineering & Department of Bioengineering The University of Tokyo K. Sugiyama (RIKEN), S. Ii (Osaka Univ.), K. Okita (Nihon Univ.) , T. Azuma ( Univ. Tokyo), Y. Matsumoto ( Univ. Tokyo )
- 2. “K (京）” computer in Kobe multi-level parallelism http://jp.fujitsu.com/solutions/hpc/brochures 10P (=1016) operations/sec = 8×104 (nodes) × 8 (cores) × 8 (operations) × 2×109 (Hz) http://www.fujitsu.com/global/about/tech/k/ 10 Peta Flops with 640,000 cores
- 3. For the Design and Development of High Intensity Focused Ultrasound Therapy (HIFU) Device
- 4. High Intensity Focused Ultrasound therapy •Prostate cancer •Breast cancer http://www.prostatecancercentre.co.uk/treatments/hifu.html HIFU therapy has been developed for the treatment of deeply-placed cancer. Brain cancer Liver cancer http://www.imasonic.com/ Displacement of focal point due to the reflection and refraction of ultrasound at the interfaces of bones Control of the focal point by an array transducer
- 6. Comparison of the focal point Without phase delay Pmax=1.7MPa 0 2 Pressure [MPa] With phase delay Pmax=2.2MPa 0 2 Pressure [MPa]
- 7. Numerical Method Suitable for Medical Image Data (Full Eulerian FSI simulations) Gaehtgens et al.(1980) Blood Cells., 6, 799.
- 8. Background Fluid-Structure coupling analysis of living body Diagnostic image (CT, MRI) Voxel data (Volume Fraction of Constituents) representing multi-component geometry without Mesh Generation (Finite-Difference or Finite-Element) Simulation on Eulerian frame
- 9. Full Eulerian approach for Fluid-Structure Interactions vs. Lagrangian Eulerian How is the two-phase distinguished ? s 1 Solid Solid 0.5 Fluid Fluid 0 by the boundary of mesh by solid volume fraction How is the solid deformation described? by the displacement of material points themselves by left Cauchy-Green deformation tensor ・ Sugiyama, Ii et al. (2011) J. Comput . Phys., 230, 596. ・ Ii, Sugiyama et al. (2011) Int. J. Numer. Meth. Fluids, 65, 150.
- 10. The Super-Fast FSI Solver: ZZ-EFSI ZZ-EFSI achieved the actual speed of 4.5 PETA FLOPS!!! Software availablle at http://www.islim.org/islim-dl_e.html
- 11. Development of Multiscale Thrombosis Simulator
- 12. Multiscale modeling of initial stage of thromobosis Blood flow (continuum mechanics) Protein(wWF)-Protein (GP1b) binding (stochastic process) •Substance diffusion •Metabolic reaction and activation •Morphology change Molecular interaction (molecular dynamics)
- 13. Adhesion of Platelets in a Shear Flow Ht = 5.5 % = 800 (s-1) t 200 (ms) Ht = 15.3 % Ht = 21.9 %
- 14. Acknowledgement Tokai University: S. Goto, S. Shiozaki, N.Tamura RIKEN: S. Noda, H. Yokota, R. Himeno Osaka University: S. Takeuchi Shanghai Jiaotong University: X. Gong, F. Liang Peking University: J. Wu York University: H. Huang References: ・ Sugiyama, el al, J. Comput. Phys., Vol.230 (2010), pp.596-627. ・ Takagi et al. J. Appl. Mech. Vol.79 (2011), 010911. ・ Ii et al., Int. J. Numerical Method. in Fluids, Vol.65 (2011). pp.43-66. ・ Ii et al., Comm. in Cmuput. Phys, Vol.12 (2013), pp.544-576. ・ Shiozaki et al. JBSE, Vol.7 (2012), pp.275-283. ・ Ii et al., J.Comput Phys, , Vol.231 (2012), pp.2328-2358.