Analyzing and Modeling the Performance of the HemeLB Lattice-Boltzmann Simulation Environment

We investigate the performance of the HemeLB lattice-Boltzmann simulator for cerebrovascular blood flow, aimed at providing timely and clinically relevant assistance to neurosurgeons. HemeLB is optimised for sparse geometries, supports interactive use, and scales well to 32,768 cores for problems with ~81 million lattice sites. We obtain a maximum performance of 29.5 billion site updates per second, with only an 11% slowdown for highly sparse problems (5% fluid fraction). We present steering and visualisation performance measurements and provide a model which allows users to predict the performance, thereby determining how to run simulations with maximum accuracy within time constraints.Comment: Accepted by the Journal of Computational Science. 33 pages, 16 figures, 7 table

Lattice-Boltzmann simulations of cerebral blood flow

Computational haemodynamics play a central role in the understanding of blood behaviour in the cerebral vasculature, increasing our knowledge in the onset of vascular diseases and their progression, improving diagnosis and ultimately providing better patient prognosis. Computer simulations hold the potential of accurately characterising motion of blood and its interaction with the vessel wall, providing the capability to assess surgical treatments with no danger to the patient. These aspects considerably contribute to better understand of blood circulation processes as well as to augment pre-treatment planning. Existing software environments for treatment planning consist of several stages, each requiring significant user interaction and processing time, significantly limiting their use in clinical scenarios. The aim of this PhD is to provide clinicians and researchers with a tool to aid in the understanding of human cerebral haemodynamics. This tool employs a high performance fluid solver based on the lattice-Boltzmann method (coined HemeLB), high performance distributed computing and grid computing, and various advanced software applications useful to efficiently set up and run patient-specific simulations. A graphical tool is used to segment the vasculature from patient-specific CT or MR data and configure boundary conditions with ease, creating models of the vasculature in real time. Blood flow visualisation is done in real time using in situ rendering techniques implemented within the parallel fluid solver and aided by steering capabilities; these programming strategies allows the clinician to interactively display the simulation results on a local workstation. A separate software application is used to numerically compare simulation results carried out at different spatial resolutions, providing a strategy to approach numerical validation. This developed software and supporting computational infrastructure was used to study various patient-specific intracranial aneurysms with the collaborating interventionalists at the National Hospital for Neurology and Neuroscience (London), using three-dimensional rotational angiography data to define the patient-specific vasculature. Blood flow motion was depicted in detail by the visualisation capabilities, clearly showing vortex fluid ow features and stress distribution at the inner surface of the aneurysms and their surrounding vasculature. These investigations permitted the clinicians to rapidly assess the risk associated with the growth and rupture of each aneurysm. The ultimate goal of this work is to aid clinical practice with an efficient easy-to-use toolkit for real-time decision support

Lattice-Boltzmann interactive blood flow simulation pipeline.

PURPOSE:Cerebral aneurysms are one of the prevalent cerebrovascular disorders in adults worldwide and caused by a weakness in the brain artery. The most impressive treatment for a brain aneurysm is interventional radiology treatment, which is extremely dependent on the skill level of the radiologist. Hence, accurate detection and effective therapy for cerebral aneurysms still remain important clinical challenges. In this work, we have introduced a pipeline for cerebral blood flow simulation and real-time visualization incorporating all aspects from medical image acquisition to real-time visualization and steering. METHODS:We have developed and employed an improved version of HemeLB as the main computational core of the pipeline. HemeLB is a massive parallel lattice-Boltzmann fluid solver optimized for sparse and complex geometries. The visualization component of this pipeline is based on the ray marching method implemented on CUDA capable GPU cores. RESULTS:The proposed visualization engine is evaluated comprehensively and the reported results demonstrate that it achieves significantly higher scalability and sites updates per second, indicating higher update rate of geometry sites' values, in comparison with the original HemeLB. This proposed engine is more than two times faster and capable of 3D visualization of the results by processing more than 30 frames per second. CONCLUSION:A reliable modeling and visualizing environment for measuring and displaying blood flow patterns in vivo, which can provide insight into the hemodynamic characteristics of cerebral aneurysms, is presented in this work. This pipeline increases the speed of visualization and maximizes the performance of the processing units to do the tasks by breaking them into smaller tasks and working with GPU to render the images. Hence, the proposed pipeline can be applied as part of clinical routines to provide the clinicians with the real-time cerebral blood flow-related information

Weighted decomposition in high-performance lattice-Boltzmann simulations: Are some lattice sites more equal than others?

Obtaining a good load balance is a significant challenge in scaling up lattice-Boltzmann simulations of realistic sparse problems to the exascale. Here we analyze the effect of weighted decomposition on the performance of the HemeLB lattice-Boltzmann simulation environment, when applied to sparse domains. Prior to domain decomposition, we assign wall and in/outlet sites with increased weights which reflect their increased computational cost. We combine our weighted decomposition with a second optimization, which is to sort the lattice sites according to a space filling curve. We tested these strategies on a sparse bifurcation and very sparse aneurysm geometry, and find that using weights reduces calculation load imbalance by up to 85 %, although the overall communication overhead is higher than some of our runs.This work has received funding from the CRESTA and MAPPER projects within the EC-FP7 (ICT-2011.9.13) under Grant Agreements nos. 287703 and 261507, and from EPSRC Grants EP/I017909/1 (www.2020science.net) and EP/I034602/1

Planning for steerable needles in neurosurgery

The increasing adoption of robotic-assisted surgery has opened up the possibility to control innovative dexterous tools to improve patient outcomes in a minimally invasive way. Steerable needles belong to this category, and their potential has been recognised in various surgical fields, including neurosurgery. However, planning for steerable catheters' insertions might appear counterintuitive even for expert clinicians. Strategies and tools to aid the surgeon in selecting a feasible trajectory to follow and methods to assist them intra-operatively during the insertion process are currently of great interest as they could accelerate steerable needles' translation from research to practical use. However, existing computer-assisted planning (CAP) algorithms are often limited in their ability to meet both operational and kinematic constraints in the context of precise neurosurgery, due to its demanding surgical conditions and highly complex environment. The research contributions in this thesis relate to understanding the existing gap in planning curved insertions for steerable needles and implementing intelligent CAP techniques to use in the context of neurosurgery. Among this thesis contributions showcase (i) the development of a pre-operative CAP for precise neurosurgery applications able to generate optimised paths at a safe distance from brain sensitive structures while meeting steerable needles kinematic constraints; (ii) the development of an intra-operative CAP able to adjust the current insertion path with high stability while compensating for online tissue deformation; (iii) the integration of both methods into a commercial user front-end interface (NeuroInspire, Renishaw plc.) tested during a series of user-controlled needle steering animal trials, demonstrating successful targeting performances. (iv) investigating the use of steerable needles in the context of laser interstitial thermal therapy (LiTT) for maesial temporal lobe epilepsy patients and proposing the first LiTT CAP for steerable needles within this context. The thesis concludes with a discussion of these contributions and suggestions for future work.Open Acces

Modeling and Simulation in Engineering

This book provides an open platform to establish and share knowledge developed by scholars, scientists, and engineers from all over the world, about various applications of the modeling and simulation in the design process of products, in various engineering fields. The book consists of 12 chapters arranged in two sections (3D Modeling and Virtual Prototyping), reflecting the multidimensionality of applications related to modeling and simulation. Some of the most recent modeling and simulation techniques, as well as some of the most accurate and sophisticated software in treating complex systems, are applied. All the original contributions in this book are jointed by the basic principle of a successful modeling and simulation process: as complex as necessary, and as simple as possible. The idea is to manipulate the simplifying assumptions in a way that reduces the complexity of the model (in order to make a real-time simulation), but without altering the precision of the results

Developmental delays and subcellular stress as downstream effects of sonoporation

Posters: no. 2Control ID: 1672434OBJECTIVES: The biological impact of sonoporation has often been overlooked. Here we seek to obtain insight into the cytotoxic impact of sonoporation by gaining new perspectives on anti-proliferative characteristics that may emerge within sonoporated cells. We particularly focused on investigating the cell-cycle progression kinetics of sonoporated cells and identifying organelles that may be stressed in the recovery process. METHODS: In line with recommendations on exposure hardware design, an immersion-based ultrasound platform has been developed. It delivers 1 MHz ultrasound pulses (100 cycles; 1 kHz PRF; 60 s total duration) with 0.45 MPa peak negative pressure to a cell chamber that housed HL-60 leukemia cells and lipid-shelled microbubbles at a 10:1 cell-tobubble ratio (for 1e6/ml cell density). Calcein was used to facilitate tracking of sonoporated cells with enhanced uptake of exogenous molecules. The developmental trend of sonoporated cells was quantitatively analyzed using BrdU/DNA flow cytometry that monitors the cell population’s DNA synthesis kinetics. This allowed us to measure the temporal progression of DNA synthesis of sonoporated cells. To investigate whether sonoporation would upset subcellular homeostasis, post-exposure cell samples were also assayed for various proteins using Western blot analysis. Analysis focus was placed on the endoplasmic reticulum (ER): an important organelle with multi-faceted role in cellular functioning. The post-exposure observation time spanned between 0-24 h. RESULTS: Despite maintaining viability, sonoporated cells were found to exhibit delays in cell-cycle progression. Specifically, their DNA synthesis time was lengthened substantially (for HL-60 cells: 8.7 h for control vs 13.4 h for the sonoporated group). This indicates that sonoporated cells were under stress: a phenomenon that is supported by our Western blot assays showing upregulation of ER-resident enzymes (PDI, Ero1), ER stress sensors (PERK, IRE1), and ER-triggered pro-apoptotic signals (CHOP, JNK). CONCLUSIONS: Sonoporation, whilst being able to facilitate internalization of exogenous molecules, may inadvertently elicit a cellular stress response. These findings seem to echo recent calls for reconsideration of efficiency issues in sonoporation-mediated drug delivery. Further efforts would be necessary to improve the efficiency of sonoporation-based biomedical applications where cell death is not desirable.postprin

A study on the change in plasma membrane potential during sonoporation

Posters: no. 4Control ID: 1680329OBJECTIVES: There has been validated that the correlation of sonoporation with calcium transients is generated by ultrasound-mediated microbubbles activity. Besides calcium, other ionic flows are likely involved in sonoporation. Our hypothesis is the cell electrophysiological properties are related to the intracellular delivery by ultrasound and microbubbles. In this study, a real-time live cell imaging platform is used to determine whether plasma membrane potential change is related to the sonoporation process at the cellular level. METHODS: Hela cells were cultured in DMEM supplemented with 10% FBS in Opticell Chamber at 37 °C and 5% CO2, and reached 80% confluency before experiments. The Calcein Blue-AM, DiBAC4(3) loaded cells in the Opticell chamber filled with PI solution and Sonovue microbubbles were immerged in a water tank on a inverted fluorescence microscope. Pulsed ultrasound (1MHz freq., 20 cycles, 20Hz PRF, 0.2-0.5MPa PNP) was irradiated at the angle of 45° to the region of interest for 1s.The real-time fluorescence imaging for different probes was acquired by a cooled CCD camera every 20s for 10min. The time-lapse fluorescence images were quantitatively analyzed to evaluate the correlation of cell viability, intracellular delivery with plasma membrane potential change. RESULTS: Our preliminary data showed that the PI fluorescence, which indicated intracellular delivery, was immediately accumulated in cells adjacent to microbubbles after exposure, suggesting that their membranes were damaged by ultrasound-activated microbubbles. However, the fluorescence reached its highest level within 4 to 6 minutes and was unchanged thereafter, indicating the membrane was gradually repaired within this period. Furthermore, using DIBAC4(3), which detected the change in the cell membrane potential, we found that the loss of membrane potential might be associated with intracellular delivery, because the PI fluorescence accumulation was usually accompanied with the change in DIBAC4 (3) fluorescence. CONCLUSIONS: Our study suggests that there may be a linkage between the cell membrane potential change and intracellular delivery mediated by ultrasound and microbubbles. We also suggest that other ionic flows or ion channels may be involved in the cell membrane potential change in sonoporation. Further efforts to explore the cellular mechanism of this phenomenon will improve our understanding of sonoporation.postprin

How sonoporation disrupts cellular structural integrity: morphological and cytoskeletal observations

Posters: no. 1Control ID: 1672429OBJECTIVES: In considering sonoporation for drug delivery applications, it is essential to understand how living cells respond to this puncturing force. Here we seek to investigate the effects of sonoporation on cellular structural integrity. We hypothesize that the membrane morphology and cytoskeletal behavior of sonoporated cells under recovery would inherently differ from that of normal viable cells. METHODS: A customized and calibrated exposure platform was developed for this work, and the ZR-75-30 breast carcinoma cells were used as the cell model. The cells were exposed to either single or multiple pulses of 1 MHz ultrasound (pulse length: 30 or 100 cycles; PRF: 1kHz; duration: up to 60s) with 0.45 MPa spatial-averaged peak negative pressure and in the presence of lipid-shelled microbubbles. Confocal microscopy was used to examine insitu the structural integrity of sonoporated cells (identified as ones with exogenous fluorescent marker internalization). For investigations on membrane morphology, FM 4-64 was used as the membrane dye (red), and calcein was used as the sonoporation marker (green); for studies on cytoskeletal behavior, CellLight (green) and propidium iodide (red) were used to respectively label actin filaments and sonoporated cells. Observation started from before exposure to up to 2 h after exposure, and confocal images were acquired at real-time frame rates. Cellular structural features and their temporal kinetics were quantitatively analyzed to assess the consistency of trends amongst a group of cells. RESULTS: Sonoporated cells exhibited membrane shrinkage (decreased by 61% in a cell’s cross-sectional area) and intracellular lipid accumulation (381% increase compared to control) over a 2 h period. The morphological repression of sonoporated cells was also found to correspond with post-sonoporation cytoskeletal processes: actin depolymerization was observed as soon as pores were induced on the membrane. These results show that cellular structural integrity is indeed disrupted over the course of sonoporation. CONCLUSIONS: Our investigation shows that the biophysical impact of sonoporation is by no means limited to the induction of membrane pores: e.g. structural integrity is concomitantly affected in the process. This prompts the need for further fundamental studies to unravel the complex sequence of biological events involved in sonoporation.postprin