Review of recent results using computational fluid dynamics simulations in patients receiving mechanical assist devices for end-stage heart failure

Many end-stage heart failure patients are not eligible to undergo heart transplantation due to organ shortage, and even those under consideration for transplantation might suffer long waiting periods. A better understanding of the hemodynamic impact of left ventricular assist devices (LVAD) on the cardiovascular system is therefore of great interest. Computational fluid dynamics (CFD) simulations give the opportunity to study the hemodynamics in this patient population using clinical imaging data such as computed tomographic angiography. This article reviews a recent study series involving patients with pulsatile and constant-flow LVAD devices in which CFD simulations were used to qualitatively and quantitatively assess blood flow dynamics in the thoracic aorta, demonstrating its potential to enhance the information available from medical imaging.

Comparison of morphological and rheological conditions between conventional and eversion carotid endarterectomy using computational fluid dynamics - a pilot study

PURPOSE To compare postoperative morphological and rheological conditions after eversion carotid endarterectomy versus conventional carotid endarterectomy using computational fluid dynamics. BASIC METHODS Hemodynamic metrics (velocity, wall shear stress, time-averaged wall shear stress and temporal gradient wall shear stress) in the carotid arteries were simulated in one patient after conventional carotid endarterectomy and one patient after eversion carotid endarterectomy by computational fluid dynamics analysis based on patient specific data. PRINCIPAL FINDINGS Systolic peak of the eversion carotid endarterectomy model showed a gradually decreased pressure along the stream path, the conventional carotid endarterectomy model revealed high pressure (about 180 Pa) at the carotid bulb. Regions of low wall shear stress in the conventional carotid endarterectomy model were much larger than that in the eversion carotid endarterectomy model and with lower time-averaged wall shear stress values (conventional carotid endarterectomy: 0.03-5.46 Pa vs. eversion carotid endarterectomy: 0.12-5.22 Pa). CONCLUSIONS Computational fluid dynamics after conventional carotid endarterectomy and eversion carotid endarterectomy disclosed differences in hemodynamic patterns. Larger studies are necessary to assess whether these differences are consistent and might explain different rates of restenosis in both techniques.

Computational fluid dynamics modelling in cardiovascular medicine.

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards 'digital patient' or 'virtual physiological human' representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

Rapid diagnosis of experimental meningitis by bacterial heat production in cerebrospinal fluid

BACKGROUND: Calorimetry is a nonspecific technique which allows direct measurement of heat generated by biological processes in the living cell. We evaluated the potential of calorimetry for rapid detection of bacterial growth in cerebrospinal fluid (CSF) in a rat model of bacterial meningitis. METHODS: Infant rats were infected on postnatal day 11 by direct intracisternal injection with either Streptococcus pneumoniae, Neisseria meningitidis or Listeria monocytogenes. Control animals were injected with sterile saline or heat-inactivated S. pneumoniae. CSF was obtained at 18 hours after infection for quantitative cultures and heat flow measurement. For calorimetry, 10 microl and 1 microl CSF were inoculated in calorimetry ampoules containing 3 ml trypticase soy broth (TSB). RESULTS: The mean bacterial titer (+/- SD) in CSF was 1.5 +/- 0.6 x 108 for S. pneumoniae, 1.3 +/- 0.3 x 106 for N. meningitidis and 3.5 +/- 2.2 x 104 for L. monocytogenes. Calorimetric detection time was defined as the time until heat flow signal exceeded 10 microW. Heat signal was detected in 10-microl CSF samples from all infected animals with a mean (+/- SD) detection time of 1.5 +/- 0.2 hours for S. pneumoniae, 3.9 +/- 0.7 hours for N. meningitidis and 9.1 +/- 0.5 hours for L. monocytogenes. CSF samples from non-infected animals generated no increasing heat flow (<10 microW). The total heat was the highest in S. pneumoniae ranging from 6.7 to 7.5 Joules, followed by L. monocytogenes (5.6 to 6.1 Joules) and N. meningitidis (3.5 to 4.4 Joules). The lowest detectable bacterial titer by calorimetry was 2 cfu for S. pneumoniae, 4 cfu for N. meningitidis and 7 cfu for L. monocytogenes. CONCLUSION: By means of calorimetry, detection times of <4 hours for S. pneumoniae and N. meningitidis and <10 hours for Listeria monocytogenes using as little as 10 microl CSF were achieved. Calorimetry is a new diagnostic method allowing rapid and accurate diagnosis of bacterial meningitis from a small volume of CSF.

Modifying peripheral IV catheters with side holes and side slits results in favorable changes in fluid dynamic properties during the injection of iodinated contrast material

OBJECTIVE: The purpose of this study was to compare a standard peripheral end-hole angiocatheter with those modified with side holes or side slits using experimental optical techniques to qualitatively compare the contrast material exit jets and using numeric techniques to provide flow visualization and quantitative comparisons. MATERIALS AND METHODS: A Schlieren imaging system was used to visualize the angiocatheter exit jet fluid dynamics at two different flow rates. Catheters were modified by drilling through-and-through side holes or by cutting slits into the catheters. A commercial computational fluid dynamics package was used to calculate numeric results for various vessel diameters and catheter orientations. RESULTS: Experimental images showed that modifying standard peripheral IV angiocatheters with side holes or side slits qualitatively changed the overall flow field and caused the exiting jet to become less well defined. Numeric calculations showed that the addition of side holes or slits resulted in a 9-30% reduction of the velocity of contrast material exiting the end hole of the angiocatheter. With the catheter tip directed obliquely to the wall, the maximum wall shear stress was always highest for the unmodified catheter and was always lowest for the four-side-slit catheter. CONCLUSION: Modified angiocatheters may have the potential to reduce extravasation events in patients by reducing vessel wall shear stress.

A Nonstationary Space-Time Gaussian Process Model for Partially Converged Simulations

In the context of expensive numerical experiments, a promising solution for alleviating the computational costs consists of using partially converged simulations instead of exact solutions. The gain in computational time is at the price of precision in the response. This work addresses the issue of fitting a Gaussian process model to partially converged simulation data for further use in prediction. The main challenge consists of the adequate approximation of the error due to partial convergence, which is correlated in both design variables and time directions. Here, we propose fitting a Gaussian process in the joint space of design parameters and computational time. The model is constructed by building a nonstationary covariance kernel that reflects accurately the actual structure of the error. Practical solutions are proposed for solving parameter estimation issues associated with the proposed model. The method is applied to a computational fluid dynamics test case and shows significant improvement in prediction compared to a classical kriging model.

Comparing the demonstration of freedom from Trichinella infection of domestic pigs by traditional and risk-based surveillance

Traditionally, the routine artificial digestion test is applied to assess the presence of Trichinella larvae in pigs. However, this diagnostic method has a low sensitivity compared to serological tests. The results from artificial digestion tests in Switzerland were evaluated over a time period of 15 years to determine by when freedom from infection based on these data could be confirmed. Freedom was defined as a 95% probability that the prevalence of infection was below 0.0001%. Freedom was demonstrated after 12 years at the latest. A new risk-based surveillance approach was then developed based on serology. Risk-based surveillance was also assessed over 15 years, starting in 2010. It was shown that by using this design, the sample size could be reduced by at least a factor of 4 when compared with the traditional testing regimen, without lowering the level of confidence in the Trichinella-free status of the pig population.

[Acute "difficult" gastrointestinal bleeding]

Gastrointestinal bleeding with its point of origin outside the reach of conventional gastro- and colonoscopy represents an extraordinary diagnostic and therapeutic challenge. Bleeding may originate from the small bowel distal to the duodenojejunal junction (middle gastrointestinal bleeding) or from the biliary tree (haemobilia) or from the pancreatic ductal system (haemosuccus pancreaticus). This particular type of gastrointestinal bleeding is often intermittend and caused by a variety of different pathologies. Angiography is the diagnostic method of choice for further investigation. It allows precise localization of the bleeding site and simultaneous interventional therapy (embolization/coiling). The importance of further diagnostic modalities such as scintigraphy, capsule endoscopy, push-enteroscopy and double-balloon-enteroscopy is discussed.

Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography

The alveolated structure of the pulmonary acinus plays a vital role in gas exchange function. Three-dimensional (3D) analysis of the parenchymal region is fundamental to understanding this structure-function relationship, but only a limited number of attempts have been conducted in the past because of technical limitations. In this study, we developed a new image processing methodology based on finite element (FE) analysis for accurate 3D structural reconstruction of the gas exchange regions of the lung. Stereologically well characterized rat lung samples (Pediatr Res 53: 72-80, 2003) were imaged using high-resolution synchrotron radiation-based X-ray tomographic microscopy. A stack of 1,024 images (each slice: 1024 x 1024 pixels) with resolution of 1.4 mum(3) per voxel were generated. For the development of FE algorithm, regions of interest (ROI), containing approximately 7.5 million voxels, were further extracted as a working subunit. 3D FEs were created overlaying the voxel map using a grid-based hexahedral algorithm. A proper threshold value for appropriate segmentation was iteratively determined to match the calculated volume density of tissue to the stereologically determined value (Pediatr Res 53: 72-80, 2003). The resulting 3D FEs are ready to be used for 3D structural analysis as well as for subsequent FE computational analyses like fluid dynamics and skeletonization.

Hemodynamics and Flow Characteristics of a New Dialysis Port

Renal replacement therapy by hemodialysis requires a permanent vascular access. Implantable ports offer a potential alternative to standard vascular access strategies although their development is limited both in number and extent. We explored the fluid dynamics within two new percutaneous bone-anchored dialysis port prototypes, both by in vitro experiments and computer simulation. The new port is to be fixed to bone and allows the connection of a dialysis machine to a central venous catheter via a built-in valve. We found that the pressure drop induced by the two ports was between 20 and 50 mmHg at 500 ml/min, which is comparable with commercial catheter connectors (15–80 mmHg). We observed the formation of vortices in both geometries, and a shear rate in the physiological range (<10,000s-1), which is lower than maximal shear rates reported in commercial catheters (up to 13,000s-1). A difference in surface shear rate of 15% between the two ports was obtained.

A longitudinal study of Type-B aortic dissection and endovascular repair scenarios: computational analyses

Conservative medical treatment is commonly first recommended for patients with uncomplicated Type-B aortic dissection (AD). However, if dissection-related complications occur, endovascular repair or open surgery is performed. Here we establish computational models of AD based on radiological three-dimensional images of a patient at initial presentation and after 4-years of best medical treatment (BMT). Computational fluid dynamics analyses are performed to quantitatively investigate the hemodynamic features of AD. Entry and re-entries (functioning as entries and outlets) are identified in the initial and follow-up models, and obvious variations of the inter-luminal flow exchange are revealed. Computational studies indicate that the reduction of blood pressure in BMT patients lowers pressure and wall shear stress in the thoracic aorta in general, and flattens the pressure distribution on the outer wall of the dissection, potentially reducing the progressive enlargement of the false lumen. Finally, scenario studies of endovascular aortic repair are conducted. The results indicate that, for patients with multiple tears, stent-grafts occluding all re-entries would be required to effectively reduce inter-luminal blood communication and thus induce thrombosis in the false lumen. This implicates that computational flow analyses may identify entries and relevant re-entries between true and false lumen and potentially assist in stent-graft planning.

State-of-the-art aortic imaging: Part II - applications in transcatheter aortic valve replacement and endovascular aortic aneurysm repair.

Transcatheter aortic valve replacement (TAVR) as well as thoracic and abdominal endovascular aortic repair (TEVAR and EVAR) rely on accurate pre- and postprocedural imaging. This review article discusses the application of imaging, including preprocedural assessment and measurements as well as postprocedural imaging of complications. Furthermore, the exciting perspective of computational fluid dynamics (CFD) based on cross-sectional imaging is presented. TAVR is a minimally invasive alternative for treatment of aortic valve stenosis in patients with high age and multiple comorbidities who cannot undergo traditional open surgical repair. Given the lack of direct visualization during the procedure, pre- and peri-procedural imaging forms an essential part of the intervention. Computed tomography angiography (CTA) is the imaging modality of choice for preprocedural planning. Routine postprocedural follow-up is performed by echocardiography to confirm treatment success and detect complications. EVAR and TEVAR are minimally invasive alternatives to open surgical repair of aortic pathologies. CTA constitutes the preferred imaging modality for both preoperative planning and postoperative follow-up including detection of endoleaks. Magnetic resonance imaging is an excellent alternative to CT for postoperative follow-up, and is especially beneficial for younger patients given the lack of radiation. Ultrasound is applied in screening and postoperative follow-up of abdominal aortic aneurysms, but cross-sectional imaging is required once abnormalities are detected. Contrast-enhanced ultrasound may be as sensitive as CTA in detecting endoleaks.

The Nature of Scientific Proof in the Age of Simulations. Is numerical mimicry a third way of establishing truth?

Is numerical mimicry a third way of establishing truth? Kevin Heng received his M.S. and Ph.D. in astrophysics from the Joint Institute for Laboratory Astrophysics (JILA) and the University of Colorado at Boulder. He joined the Institute for Advanced Study in Princeton from 2007 to 2010, first as a Member and later as the Frank & Peggy Taplin Member. From 2010 to 2012 he was a Zwicky Prize Fellow at ETH Z¨urich (the Swiss Federal Institute of Technology). In 2013, he joined the Center for Space and Habitability (CSH) at the University of Bern, Switzerland, as a tenure-track assistant professor, where he leads the Exoplanets and Exoclimes Group. He has worked on, and maintains, a broad range of interests in astrophysics: shocks, extrasolar asteroid belts, planet formation, fluid dynamics, brown dwarfs and exoplanets. He coordinates the Exoclimes Simulation Platform (ESP), an open-source set of theoretical tools designed for studying the basic physics and chemistry of exoplanetary atmospheres and climates (www.exoclime.org). He is involved in the CHEOPS (Characterizing Exoplanet Satellite) space telescope, a mission approved by the European Space Agency (ESA) and led by Switzerland. He spends a fair amount of time humbly learning the lessons gleaned from studying the Earth and Solar System planets, as related to him by atmospheric, climate and planetary scientists. He received a Sigma Xi Grant-in-Aid of Research in 2006