Ultrafast Ultrasound Imaging

Among medical imaging modalities, such as computed tomography (CT) and magnetic resonance imaging (MRI), ultrasound imaging stands out due to its temporal resolution. Owing to the nature of medical ultrasound imaging, it has been used for not only observation of the morphology of living organs but also functional imaging, such as blood flow imaging and evaluation of the cardiac function. Ultrafast ultrasound imaging, which has recently become widely available, significantly increases the opportunities for medical functional imaging. Ultrafast ultrasound imaging typically enables imaging frame-rates of up to ten thousand frames per second (fps). Due to the extremely high temporal resolution, this enables visualization of rapid dynamic responses of biological tissues, which cannot be observed and analyzed by conventional ultrasound imaging. This Special Issue includes various studies of improvements to the performance of ultrafast ultrasoun

Level Set Methods for MRE Image Processing and Analysis

Ph.DDOCTOR OF PHILOSOPH

The LifeV library: engineering mathematics beyond the proof of concept

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar

Overcoming conventional modeling limitations using image- driven lattice-boltzmann method simulations for biophysical applications

The challenges involved in modeling biological systems are significant and push the boundaries of conventional modeling. This is because biological systems are distinctly complex, and their emergent properties are results of the interplay of numerous components/processes. Unfortunately, conventional modeling approaches are often limited by their inability to capture all these complexities. By using in vivo data derived from biomedical imaging, image-based modeling is able to overcome this limitation. In this work, a combination of imaging data with the Lattice-Boltzmann Method for computational fluid dynamics (CFD) is applied to tissue engineering and thrombogenesis. Using this approach, some of the unanswered questions in both application areas are resolved. In the first application, numerical differences between two types of boundary conditions: “wall boundary condition” (WBC) and “periodic boundary condition” (PBC), which are commonly utilized for approximating shear stresses in tissue engineering scaffold simulations is investigated. Surface stresses in 3D scaffold reconstructions, obtained from high resolution microcomputed tomography images are calculated for both boundary condition types and compared with the actual whole scaffold values via image-based CFD simulations. It is found that, both boundary conditions follow the same spatial surface stress patterns as the whole scaffold simulations. However, they under-predict the absolute stress values approximately by a factor of two. Moreover, it is found that the error grows with higher scaffold porosity. Additionally, it is found that the PBC always resulted in a lower error than the WBC. In a second tissue engineering study, the dependence of culture time on the distribution and magnitude of fluid shear in tissue scaffolds cultured under flow perfusion is investigated. In the study, constructs are destructively evaluated with assays for cellularity and calcium deposition, imaged using µCT and reconstructed for CFD simulations. It is found that both the shear stress distributions within scaffolds consistently increase with culture time and correlate with increasing levels of mineralized tissues within the scaffold constructs as seen in calcium deposition data and µCT reconstructions. In the thrombogenesis application, detailed analysis of time lapse microscopy images showing yielding of thrombi in live mouse microvasculature is performed. Using these images, image-based CFD modeling is performed to calculate the fluid-induced shear stresses imposed on the thrombi’s surfaces by the surrounding blood flow. From the results, estimates of the yield stress (A critical parameter for quantifying the extent to which thrombi material can resist deformation and breakage) are obtained for different blood vessels. Further, it is shown that the yielding observed in thrombi occurs mostly in the outer shell region while the inner core remains intact. This suggests that the core material is different from the shell. To that end, we propose an alternative mechanism of thrombogenesis which could help explain this difference. Overall, the findings from this work reveal that image-based modeling is a versatile approach which can be applied to different biomedical application areas while overcoming the difficulties associated with conventional modeling

Data-driven simulation for augmented surgery

International audienceTo build an augmented view of an organ during surgery, it is essential to have a biomechanical model with appropriate material parameters and boundary conditions , able to match patient specific properties. Adaptation to the patient's anatomy is obtained by exploiting the image-rich context specific to our application domain. While information about the organ shape, for instance, can be obtained preoper-atively, other patient-specific parameters can only be determined intraoperatively. To this end, we are developing data-driven simulations, which exploit information extracted from a stream of medical images. Such simulations need to run in real-time. To this end we have developed dedicated numerical methods, which allow for real-time computation of finite element simulations. The general principle consists in combining finite element approaches with Bayesian methods or deep learning techniques, that allow to keep control over the underlying computational model while allowing for inputs from the real world. Based on a priori knowledge of the mechanical behavior of the considered organ, we select a constitutive law to model its deformations. The predictive power of such constitutive law highly depends on the knowledge of the material parameters and A. Mendizaba

FluTAS: A GPU-accelerated finite difference code for multiphase flows

We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-B\'enard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows