Hemodynamics and wall shear metrics in a pulmonary autograft: Comparing a fluid-structure interaction and computational fluid dynamics approach

OBJECTIVE

In young patients, aortic valve disease is often treated by placement of a pulmonary autograft (PA) which adapts to its new environment through growth and remodeling. To better understand the hemodynamic forces acting on the highly distensible PA in the acute phase after surgery, we developed a fluid-structure interaction (FSI) framework and comprehensively compared hemodynamics and wall shear-stress (WSS) metrics with a computational fluid dynamic (CFD) simulation.

METHODS

The FSI framework couples a prestressed non-linear hyperelastic arterial tissue model with a fluid model using the in-house coupling code CoCoNuT. Geometry, material parameters and boundary conditions are based on in-vivo measurements. Hemodynamics, time-averaged WSS (TAWSS), oscillatory shear index (OSI) and topological shear variation index (TSVI) are evaluated qualitatively and quantitatively for 3 different sheeps.

RESULTS

Despite systolic-to-diastolic volumetric changes of the PA in the order of 20 %, the point-by-point correlation of TAWSS and OSI obtained through CFD and FSI remains high (r > 0.9, p < 0.01) for TAWSS and (r > 0.8, p < 0.01) for OSI). Instantaneous WSS divergence patterns qualitatively preserve similarities, but large deformations of the PA leads to a decrease of the correlation between FSI and CFD resolved TSVI (r < 0.7, p < 0.01). Moderate co-localization between FSI and CFD is observed for low thresholds of TAWSS and high thresholds of OSI and TSVI.

CONCLUSION

FSI might be warranted if we were to use the TSVI as a mechano-biological driver for growth and remodeling of PA due to varying intra-vascular flow structures and near wall hemodynamics because of the large expansion of the PA.

N-acetyl-L-cysteine improves the performance of chronic cyclic heat-stressed finisher broilers but has no effect on tissue glutathione levels

1. It was hypothesised that dietary N-acetyl-L-cysteine (NAC) in feed, as a source of cysteine, could improve the performance of heat-stressed finisher broilers by fostering glutathione (GSH) synthesis. GSH is the most abundant intracellular antioxidant for which the sulphur amino acid cysteine is rate limiting for its synthesis.2. In the first experiment, four levels of NAC: 0, 500, 1000 and 2000 mg/kg were added to a diet with a suboptimal level of sulphur amino acids in the finisher phase. In the second experiment, NAC was compared to other sulphur amino acid sources at equal molar amounts of digestible sulphur amino acids. Birds were allocated to four groups: control, 2000 mg/kg NAC, 1479 mg/kg L-cystine, and 2168 mg/kg Ca-salt of 2-hydroxy-4-(methylthio)butanoic acid. A chronic cyclic heat stress model (temperature was increased to 34°C for 7 h daily) was initiated at 28 d of age.3. In the first experiment, growth performance and feed efficiency in the finisher phase were significantly improved by graded NAC. ADG was 88.9, 92.2, 93.7 and 97.7 g/d, and the feed-to-gain ratio was 2.18, 1.91, 1.85 and 1.81 for the 0, 500, 1000 and 2000 mg/kg NAC treatments, respectively. However, liver and heart GSH levels were not affected by NAC. On d 29, liver gene transcript of cystathionine-beta-synthase like was reduced by NAC, which suggested reduced trans-sulphuration activity. The second experiment showed that L-cystine and Ca-salt of 2-hydroxy-4-(methylthio) butanoic acid were more effective in improving performance than NAC.4. In conclusion, N-acetyl-L-cysteine improved dose-dependently growth and feed efficiency in heat-stressed finishing broilers. However, this was not associated with changes in tissue GSH levels, but more likely worked by sparing methionine and/or NAC's and cysteine's direct antioxidant properties.

A CFD study on the interplay of torsion and vortex guidance by the mitral valve on the left ventricular wash-out making use of overset meshes (Chimera technique)

Cardiovascular disease often occurs with silent and gradual alterations of cardiac blood flow that can lead to the onset of chronic pathological conditions. Image-based patient-specific Computational Fluid Dynamics (CFD) models allow for an extensive quantification of the flow field beyond the direct capabilities of medical imaging techniques that could support the clinicians in the early diagnosis, follow-up, and treatment planning of patients. Nonetheless, the large and impulsive kinematics of the left ventricle (LV) and the mitral valve (MV) pose relevant modeling challenges. Arbitrary Lagrangian-Eulerian (ALE) based computational fluid dynamics (CFD) methods struggle with the complex 3D mesh handling of rapidly moving valve leaflets within the left ventricle (LV). We, therefore, developed a Chimera-based (overset meshing) method to build a patient-specific 3D CFD model of the beating LV which includes a patient-inspired kinematic model of the mitral valve (LVMV). Simulations were performed with and without torsion. In addition, to evaluate how the intracardiac LV flow is impacted by the MV leaflet kinematics, a third version of the model without the MV was generated (LV with torsion). For all model versions, six cardiac cycles were simulated. All simulations demonstrated cycle-to-cycle variations that persisted after six cycles but were albeit marginal in terms of the magnitude of standard deviation of velocity and vorticity which may be related to the dissipative nature of the numerical scheme used. The MV was found to have a crucial role in the development of the intraventricular flow by enhancing the direct flow, the apical washout, and the propagation of the inlet jet towards the apical region. Consequently, the MV is an essential feature in the patient-specific CFD modeling of the LV. The impact of torsion was marginal on velocity, vorticity, wall shear stress, and energy loss, whereas it resulted to be significant in the evaluation of particle residence times. Therefore, including torsion could be considered in patient-specific CFD models of the LV, particularly when aiming to study stasis and residence time. We conclude that, despite some technical limitations encountered, the Chimera technique is a promising alternative for ALE methods for 3D CFD models of the heart that include the motion of valve leaflets.

Computational fluid dynamics model to predict the dynamical behavior of the cerebrospinal fluid through implementation of physiological boundary conditions

Cerebrospinal fluid (CSF) dynamics play an important role in maintaining a stable central nervous system environment and are influenced by different physiological processes. Multiple studies have investigated these processes but the impact of each of them on CSF flow is not well understood. A deeper insight into the CSF dynamics and the processes impacting them is crucial to better understand neurological disorders such as hydrocephalus, Chiari malformation, and intracranial hypertension. This study presents a 3D computational fluid dynamics (CFD) model which incorporates physiological processes as boundary conditions. CSF production and pulsatile arterial and venous volume changes are implemented as inlet boundary conditions. At the outlets, 2-element windkessel models are imposed to simulate CSF compliance and absorption. The total compliance is first tuned using a 0D model to obtain physiological pressure pulsations. Then, simulation results are compared with in vivo flow measurements in the spinal subarachnoid space (SAS) and cerebral aqueduct, and intracranial pressure values reported in the literature. Finally, the impact of the distribution of and total compliance on CSF pressures and velocities is evaluated. Without respiration effects, compliance of 0.17 ml/mmHg yielded pressure pulsations with an amplitude of 5 mmHg and an average value within the physiological range of 7-15 mmHg. Also, model flow rates were found to be in good agreement with reported values. However, when adding respiration effects, similar pressure amplitudes required an increase of compliance value to 0.51 ml/mmHg, which is within the range of 0.4-1.2 ml/mmHg measured in vivo. Moreover, altering the distribution of compliance over the four different outlets impacted the local flow, including the flow through the foramen magnum. The contribution of compliance to each outlet was directly proportional to the outflow at that outlet. Meanwhile, the value of total compliance impacted intracranial pressure. In conclusion, a computational model of the CSF has been developed that can simulate CSF pressures and velocities by incorporating boundary conditions based on physiological processes. By tuning these boundary conditions, we were able to obtain CSF pressures and flows within the physiological range.

In vitro stability and ex vivo absorption of thymol monoglucosides in the porcine gut

Thymol α-D-glucopyranoside (TαG) and thymol β-D-glucopyranoside (TβG) are believed to have different kinetic behaviours in the porcine gut than its parent aglycon thymol. However, recently, it was shown that concentrations of both glucosides decreased rapidly in the stomach and proximal small intestine following oral supplementation to piglets as did thymol. Yet, the stability of thymol glucosides in gut contents and their absorption route remains obscure. Therefore, a series of in vitro incubations were performed, simulating the impact of pH, digestive enzymes, bacterial activity and mucosal extracts on stability of these glucosides. Their absorption mechanisms were investigated using the Ussing chamber model in the presence or the absence of inhibitors of sodium-dependent glucose linked transporter 1 and lactase phlorizin hydrolase. Both glucosides remained intact at physiological pH levels in the presence of digestive enzymes. Recoveries from TαG and TβG were below 90% when incubated with small intestinal homogenates from the distal jejunum or from all sampled sites, respectively. However, no aglycon could be detected in these samples. Bacterial inoculum of the small intestine, on the other hand, hydrolysed TβG quickly with up to 44% of free aglycon appearing. TαG proved more resistant to porcine gastro-intestinal bacterial glucosidases with only trace amounts (<1%) of free thymol at the end of the incubations. Electrophysiological measurements in Ussing chambers did not suggest active transport of the glucosides. Mucosal TαG and TβG concentrations were unchanged between start and end of the absorption measurements. Additionally, no TαG and only a very limited amount of TβG were retrieved from the serosal side. Tissue associated concentrations, although marginal (<1% of luminal concentration), were mainly as intact glucoside or as aglycon for TαG and TβG, respectively. Addition of both inhibitors significantly increased the amount of intact glucosides retrieved from the mucosal tissues as compared to controls. In conclusion, bacterial hydrolysis was identified as the most important source of TβG loss, whereas TαG seemed less prone to degradation or absorption in these in vitro and ex vivo models.

Mixed impact of torsion on LV hemodynamics: A CFD study based on the Chimera technique

Image-based patient-specific Computational Fluid Dynamics (CFD) models of the Left Ventricle (LV) can be used to quantify hemodynamics-based biomarkers that can support the clinicians in the early diagnosis, follow-up and treatment planning of patients, beyond the capabilities of the current imaging modalities. We propose a workflow to build patient-specific CFD models of the LV with moving boundaries based on the Chimera technique to overcome the convergence issues previously encountered by means of the Arbitrarian Lagrangian Eulerian approach. The workflow was tested while investigating whether the torsional motion has an impact on LV fluid dynamics. Starting from 3D cine MRI scans of a healthy volunteer, six cardiac cycles were simulated in three CFD LV models: with no, physiological, and exaggerated torsion. The Chimera technique was robust in handling the impulsive motion of the LV endocardium, allowing to notice cycle-to-cycle variations in every simulated case. Torsion affected slightly velocity, vorticity, WSS. It did not affect energy loss and induced a double-sided effect in terms of residence time: the particles ejected in one beat decreased, whereas the motility of the particles remaining in the LV was affected only in the exaggerated torsion case, indicating that implementation of torsion can be discarded in case of physiological levels. Nonetheless, caution is warranted when interpreting these results given the absence of the mitral valve, the papillary muscles, and the trabeculae. The effects of the mitral valve will be evaluated within an Fluid Structure Interaction simulation framework as further development of the current model.

Patient-specific CFD models for intraventricular flow analysis from 3D ultrasound imaging: Comparison of three clinical cases

BACKGROUND

As the intracardiac flow field is affected by changes in shape and motility of the heart, intraventricular flow features can provide diagnostic indications. Ventricular flow patterns differ depending on the cardiac condition and the exploration of different clinical cases can provide insights into how flow fields alter in different pathologies.

METHODS

In this study, we applied a patient-specific computational fluid dynamics model of the left ventricle and mitral valve, with prescribed moving boundaries based on transesophageal ultrasound images for three cardiac pathologies, to verify the abnormal flow patterns in impaired hearts. One case (P1) had normal ejection fraction but low stroke volume and cardiac output, P2 showed low stroke volume and reduced ejection fraction, P3 had a dilated ventricle and reduced ejection fraction.

RESULTS

The shape of the ventricle and mitral valve, together with the pathology influence the flow field in the left ventricle, leading to distinct flow features. Of particular interest is the pattern of the vortex formation and evolution, influenced by the valvular orifice and the ventricular shape. The base-to-apex pressure difference of maximum 2mmHg is consistent with reported data.

CONCLUSION

We used a CFD model with prescribed boundary motion to describe the intraventricular flow field in three patients with impaired diastolic function. The calculated intraventricular flow dynamics are consistent with the diagnostic patient records and highlight the differences between the different cases. The integration of clinical images and computational techniques, therefore, allows for a deeper investigation intraventricular hemodynamics in patho-physiology.

Patient-specific CFD simulation of intraventricular haemodynamics based on 3D ultrasound imaging

BACKGROUND

The goal of this paper is to present a computational fluid dynamic (CFD) model with moving boundaries to study the intraventricular flows in a patient-specific framework. Starting from the segmentation of real-time transesophageal echocardiographic images, a CFD model including the complete left ventricle and the moving 3D mitral valve was realized. Their motion, known as a function of time from the segmented ultrasound images, was imposed as a boundary condition in an Arbitrary Lagrangian-Eulerian framework.

RESULTS

The model allowed for a realistic description of the displacement of the structures of interest and for an effective analysis of the intraventricular flows throughout the cardiac cycle. The model provides detailed intraventricular flow features, and highlights the importance of the 3D valve apparatus for the vortex dynamics and apical flow.

CONCLUSIONS

The proposed method could describe the haemodynamics of the left ventricle during the cardiac cycle. The methodology might therefore be of particular importance in patient treatment planning to assess the impact of mitral valve treatment on intraventricular flow dynamics.

Fluid-Structure Interaction Simulation of Prosthetic Aortic Valves: Comparison between Immersed Boundary and Arbitrary Lagrangian-Eulerian Techniques for the Mesh Representation

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations’ outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results.

Differential impact of local stiffening and narrowing on hemodynamics in repaired aortic coarctation: an FSI study

Even after successful treatment of aortic coarctation, a high risk of cardiovascular morbidity and mortality remains. Uncertainty exists on the factors contributing to this increased risk among which are the presence of (1) a residual narrowing leading to an additional resistance and (2) a less distensible zone disturbing the buffer function of the aorta. As the many interfering factors and adaptive physiological mechanisms present in vivo prohibit the study of the isolated impact of these individual factors, a numerical fluid-structure interaction model is developed to predict central hemodynamics in coarctation treatment. The overall impact of a stiffening on the hemodynamics is limited, with a small increase in systolic pressure (up to 8 mmHg) proximal to the stiffening which is amplified with increasing stiffening and length. A residual narrowing, on the other hand, affects the hemodynamics significantly. For a short segment (10 mm), the combination of a stiffening and narrowing (coarctation index 0.5) causes an increase in systolic pressure of 58 mmHg, with 31 mmHg due to narrowing and an additional 27 mmHg due to stiffening. For a longer segment (25 mm), an increase in systolic pressure of 50 mmHg is found, of which only 9 mmHg is due to stiffening.

An animal-specific FSI model of the abdominal aorta in anesthetized mice

Recent research has revealed that angiotensin II-induced abdominal aortic aneurysm in mice can be related to medial ruptures occurring in the vicinity of abdominal side branches. Nevertheless a thorough understanding of the biomechanics near abdominal side branches in mice is lacking. In the current work we present a mouse-specific fluid-structure interaction (FSI) model of the abdominal aorta in ApoE(-/-) mice that incorporates in vivo stresses. The aortic geometry was based on contrast-enhanced in vivo micro-CT images, while aortic flow boundary conditions and material model parameters were based on in vivo high-frequency ultrasound. Flow waveforms predicted by FSI simulations corresponded better to in vivo measurements than those from CFD simulations. Peak-systolic principal stresses at the inner and outer aortic wall were locally increased caudal to the celiac and left lateral to the celiac and mesenteric arteries. Interestingly, these were also the locations at which a tear in the tunica media had been observed in previous work on angiotensin II-infused mice. Our preliminary results therefore suggest that local biomechanics play an important role in the pathophysiology of branch-related ruptures in angiotensin-II infused mice. More elaborate follow-up research is needed to demonstrate the role of biomechanics and mechanobiology in a longitudinal setting.

Influence of valve size, orientation and downstream geometry of an aortic BMHV on leaflet motion and clinically used valve performance parameters

The aim of this study was to reconcile some of our own previous work and the work of others to generate a physiologically realistic numerical simulation environment that allows to virtually assess the performance of BMHVs. The model incorporates: (i) a left ventricular deformable model to generate a physiological inflow to the aortic valve; (ii) a patient-specific aortic geometry (root, arch and descending aorta); (iii) physiological pressure and flow boundary conditions. We particularly studied the influence of downstream geometry, valve size and orientation on leaflet kinematics and functional indices used in clinical routine. Compared to the straight tube geometry, the patient-specific aorta leads to a significant asynchronous movement of the valve, especially during the closing of the valve. The anterior leaflet starts to close first, impacts the casing at the closed position and remains in this position. At the same time, the posterior leaflet impacts the pivoting mechanisms at the fully open position. At the end of systole, this leaflet subsequently accelerates to the closed position, impacting the casing with an angular velocity of approximately -477 rad/s. The valve size greatly influences the transvalvular pressure gradient (TPG), but does not change the overall leaflet kinematics. This is in contrast to changes in valve orientation, where changing valve orientation induces large differences in leaflet kinematics, but the TPG remains approximately the same.

Validation of a numerical FSI simulation of an aortic BMHV by in vitro PIV experiments

In this paper, a validation of a recently developed fluid-structure interaction (FSI) coupling algorithm to simulate numerically the dynamics of an aortic bileaflet mechanical heart valve (BMHV) is performed. This validation is done by comparing the numerical simulation results with in vitro experiments. For the in vitro experiments, the leaflet kinematics and flow fields are obtained via the particle image velocimetry (PIV) technique. Subsequently, the same case is numerically simulated by the coupling algorithm and the resulting leaflet kinematics and flow fields are obtained. Finally, the results are compared, revealing great similarity in leaflet motion and flow fields between the numerical simulation and the experimental test. Therefore, it is concluded that the developed algorithm is able to capture very accurately all the major leaflet kinematics and dynamics and can be used to study and optimize the design of BMHVs.

Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Changes in the pig small intestinal mucosal glutathione kinetics after weaning

Glutathione (GSH) serves as a major endogenous antioxidant and its kinetics have been poorly described in the weaned pig. This study was to assess the effect of birth weight, sex, and days postweaning on the small intestine (SI) mucosal GSH kinetics. At weaning (18.8 ± 0.44 d) 34 pairs of intrauterine growth restricted (IUGR) and normal birth weight sex-matched littermates were selected and fed a starter diet ad libitum until 1 h before sampling at 0, 2, 5, 12 and 28 d postweaning. Mucosa was collected from 2 SI sites, at 5% and at 75% of total length, to determine GSH and glutathione disulfide (GSSG), glutathione peroxidase (GSH-Px) and malondialdehyde (MDA), and plasm GSH-Px and MDA. At both 5 and 75% of total length, the GSH-Px activity and GSH concentrations increased gradually with increasing days postweaning to peak at day 12 (P < 0.05). The GSH-Px activity and GSH concentrations at 5% of SI length were consistently higher as compared to 75% of SI length (e.g., at day 12, 43.2 and 28.9 units/mg protein and 21.5 and 15.4 μmol/g protein, respectively). The GSSG:GSH ratio at 5% of total length was 2-fold higher at day 5 compared to all other days (P < 0.05), possibly indicating that the mucosal redox balance was disturbed in that time window. The higher GSH-Px activity, GSH content, and GSSG:GSH ratio in the proximal SI might illustrate the higher need for antioxidant action at that site. Plasma MDA and GSH-Px activity followed a comparable pattern as in the small intestine.

Velocity profiles in the human ductus venosus: a numerical fluid structure interaction study

The veins distributing oxygenated blood from the placenta to the fetal body have been given much attention in clinical Doppler velocimetry studies, in particular the ductus venosus. The ductus venosus is embedded in the left liver lobe and connects the intra-abdominal portion of the umbilical vein (IUV) directly to the inferior vena cava, such that oxygenated blood can bypass the liver and flow directly to the fetal heart. In the current work, we have developed a mathematical model to assist the clinical assessment of volumetric flow rate at the inlet of the ductus venosus. With a robust estimate of the velocity profile shape coefficient (VC), the volumetric flow rate may be estimated as the product of the time-averaged cross-sectional area, the time-averaged cross-sectional maximum velocity and the VC. The time average quantities may be obtained from Doppler ultrasound measurements, whereas the VC may be estimated from numerical simulations. The mathematical model employs a 3D fluid structure interaction model of the bifurcation formed by the IUV, the ductus venosus and the left portal vein. Furthermore, the amniotic portion of the umbilical vein, the right liver lobe and the inferior vena cava were incorporated as lumped model boundary conditions for the fluid structure interaction model. A hyperelastic material is used to model the structural response of the vessel walls, based on recently available experimental data for the human IUV and ductus venous. A parametric study was constructed to investigate the VC at the ductus venosus inlet, based on a reference case for a human fetus at 36 weeks of gestation. The VC was found to be [Formula: see text] (Mean [Formula: see text] SD of parametric case study), which confirms previous studies in the literature on the VC at the ductus venosus inlet. Additionally, CFD simulations with rigid walls were performed on a subsection of the parametric case study, and only minor changes in the predicted VCs were observed compared to the FSI cases. In conclusion, the presented mathematical model is a promising tool for the assessment of ductus venosus Doppler velocimetry.

Modeling hemodynamics in vascular networks using a geometrical multiscale approach: numerical aspects

On the one hand the heterogeneity of the circulatory system requires the use of different models in its different compartments, featuring different assumptions on the spatial degrees of freedom. On the other hand, the mutual interactions between its compartments imply that these models should preferably not be considered separately. These requirements have led to the concept of geometrical multiscale modeling, where the main idea is to couple 3D models with reduced 1D and/or 0D models. As such detailed information on the flow field in a specific region of interest can be obtained while accounting for the global circulation. However, the combination of models with different mathematical features gives rise to many difficulties such as the assignment of boundary conditions at the interface between two models and the development of robust coupling algorithms, as the subproblems are usually solved in a partitioned way. This review aims to give an overview of the most important aspects concerning 3D-1D-0D coupled models. In addition, some applications are presented in order to illustrate the potentialities of these coupled models.

Magnetic resonance imaging of the sternal extremity of the clavicle in forensic age estimation: towards more sound age estimates

OBJECTIVES

The development of the sternal extremity of the clavicle was examined prospectively with magnetic resonance imaging (MRI) for forensic bone age assessment. The objectives were to examine whether these results were comparable with the literature on plain radiography and CT, to provide age prediction intervals that correctly reflect the high variability in clavicular growth among individuals and to investigate whether MRI of the clavicle can be used to differentiate between individuals younger or older than 18 years.

METHODS

The clavicles of 220 volunteers (16-26 years) were examined with 3-T MRI and evaluated according to the Schmeling and Kreitner classification (Schmeling, Int J Legal Med 118:5-8, 2004; Kreitner, Eur Radiol 8:1116-1122, 1998). An additional hand/wrist radiograph was taken and evaluated according to Greulich and Pyle (1959). After a descriptive analysis, a multivariate ordinal regression model was fitted and embedded in a Bayesian framework based on Thevissen et al. (Int J Legal Med 124:35-42, 2009).

RESULTS

The descriptive results were concordant with the literature, although the Kreitner classification is recommended and simultaneous evaluation of the hand is considered a basic requirement (Schmeling, Int J Legal Med 118:5-8, 2004). The 95 % credibility interval for both genders with bilateral stage IV is 20-26 years. The corresponding estimated probability of being younger than 18 years is 0.8 % in females and 0.2 % in males.

CONCLUSION

MRI of the sternal extremity of the clavicle can be used to differentiate between being younger or older than 18, but a simultaneous evaluation of the hand/wrist is essential. Future evaluation of the predictive performance of the model, using comparable but larger reference samples, is necessary to validate these results.

The upstream boundary condition influences the leaflet opening dynamics in the numerical FSI simulation of an aortic BMHV

In this paper, the influence of the upstream boundary condition in the numerical simulation of an aortic bileaflet mechanical heart valve (BMHV) is studied. Three three-dimensional cases with different upstream boundary conditions are compared. The first case consists of a rigid straight tube with a velocity profile at its inlet. In the second case, the upstream geometry is a contracting left ventricle (LV), positioned symmetrically with respect to the valve. In the last case, the LV is positioned asymmetrical with respect to the valve. The cases are used to simulate the same three-dimensional BMHV. The change in time of the LV volume is calculated such that the flow rate through the valve is identical in each case. The opening dynamics of the BMHV are modelled using fluid-structure interaction. The simulations show that differences occur in the leaflet movement of the three cases. In particular, with the asymmetric LV, one of the leaflets impacts the blocking mechanism at its open position with a 34% higher velocity than when using the velocity profile, and with an 88% higher velocity than in the symmetric LV case. Therefore, when simulating such an impact, the upstream boundary condition needs to be chosen carefully.

Accuracy of carotid strain estimates from ultrasonic wall tracking: a study based on multiphysics simulations and in vivo data

We used a multiphysics model to assess the accuracy of carotid strain estimates derived from a 1-D ultrasonic wall tracking algorithm. The presented tool integrates fluid-structure interaction (FSI) simulations with an ultrasound simulator (Field II), which allows comparison of the ultrasound (US) images with a ground truth. Field II represents tissue as random points on which US waves reflect and whose position can be updated based on the flow field and vessel wall deformation from FSI. We simulated the RF-signal of a patient-specific carotid bifurcation, including the blood pool as well as the vessel wall and surrounding tissue. Distension estimates were obtained from a wall tracking algorithm using tracking points at various depths within the wall, and further processed to assess radial and circumferential strain. The simulated data demonstrated that circumferential strain can be estimated with reasonable accuracy (especially for the common carotid artery and at the lumen-intima and media-adventitia interface), but the technique does not allow to reliably assess intra-arterial radial strain. These findings were supported by in vivo data of 10 healthy adults, showing similar circumferential and radial strain profiles throughout the arterial wall. We concluded that these deviations are present due to the complex 3-D vessel wall deformation, the presence of specular reflections and, to a lesser extent, the spatially varying beam profile, with the error depending on the phase in the cardiac cycle and the scanning location.

FSI simulation of asymmetric mitral valve dynamics during diastolic filling

In this article, we present a fluid-structure interaction algorithm accounting for the mutual interaction between two rigid bodies. The algorithm was used to perform a numerical simulation of mitral valve (MV) dynamics during diastolic filling. In numerical simulations of intraventricular flow and MV motion, the asymmetry of the leaflets is often neglected. In this study the MV was rendered as two rigid, asymmetric leaflets. The 2D simulations incorporated the dynamic interaction of blood flow and leaflet motion and an imposed subject-specific, transient left ventricular wall movement obtained from ultrasound recordings. By including the full Jacobian matrix in the algorithm, the speed of the simulation was enhanced by more than 20% compared to using a diagonal Jacobian matrix. Furthermore, our results indicate that important features of the flow field may not be predicted by the use of symmetric leaflets or in the absence of an adequate model for the left atrium.

The prolonged P(300) latency in recently detoxified alcohol-dependent patients is related to activation of the inflammatory response system

The aims of this study were to examine the late components of the auditory event-related potentials (AERPs), i.e. N(100), P(200) and P(300), in recently detoxified alcohol-dependent patients compared to normal controls and to investigate whether there is a relationship between alterations in these AERPs and signs of activation of the inflammatory response system (IRS). The study subjects consisted of 14 healthy volunteers and 14 recently detoxified alcohol-dependent patients. All subjects performed a two-tone auditory discrimination task, using a standard "oddball" paradigm. The alcohol-dependent patients had their blood sampled to examine IRS markers, such as erythrocyte sedimentation rate (ESR), serum copper concentrations and the number of leukocytes. The P(300) latency was significantly greater in recently detoxified alcohol-dependent patients than in normal controls. There were significant correlations between the P(300) latency and the ESR (r = 0.84, p = 0.009), serum copper concentrations (r = 0.73, p = 0.01) and number of monocytes (r = 0.71, p = 0.006). It is concluded that the P(300) latency is prolonged in detoxified, chronic alcohol-dependent patients and is positively related to indicators of IRS activation. It is hypothesized that activation of the IRS may play a role in the delayed P(300) latency in recently detoxified, alcohol-dependent patients.

Linear CT-scan measurements in alcohol-dependent patients with and without delirium tremens

The aims of the present study were to examine whether chronic alcohol dependence and the development of delirium tremens are characterized by changes in linear CT measurements of brain liquor spaces and intracranial distances indicative of prefrontal atrophy, and frontal (sub)cortical or temporal (sub)cortical atrophy. Toward this end linear measurements were performed in 47 alcohol dependent patients with and without a history of delirium tremens and in 10 healthy volunteers using CT-scanning. The following linear measurements were calculated: (1) the Evans ratio; (2) the cella media index (CMI); (3) the maximum width of the third ventricle; (4) the maximum width of the fourth ventricle; (5) the maximum frontal subarachnoid space (MFSS); (6) the maximum width of the anterior interhemispheric fissure (MIF), and (7) the maximum width of the Sylvian fissure (MSF). The alcoholics were divided into subgroups according to the Münchner Alkoholismus Test (MALT) and the presence of delirium tremens. The MFSS of the alcohol-dependent patients was significantly larger than that of the controls. The MIF and MSF of high MALT scorers were significantly larger than those of low scorers and controls. Alcohol-dependent patients with a known history of delirium tremens had significantly larger MIF and MSF than did patients without delirium tremens and controls. The results suggest that alcohol dependence is characterized by prefrontal atrophy, and that frontal cortical and temporal (sub)cortical atrophy may be related to the development of delirium tremens.

Detection of human papilloma virus (HPV) genomes by the primed in situ (PRINS) labelling technique

Primed in situ Labelling, a technique based on primer mediated DNA synthesis, has become a useful tool in cytogenetics, especially for chromosome mapping, banding and the investigation of sequence organization in fresh metaphase preparations. Its application in the routine surgical pathology laboratory has been hampered by the fact that the technique did not work on paraffin-embedded, formalin-fixed tissue. We investigated cervical biopsies (n = 20) with morphological signs of HPV infection and found that the PRINS method is at least as sensitive as a classical in situ hybridization assay for detecting HPV DNA in paraffin-embedded, formalin-fixed tissue. In all investigated cases (n = 20), HPV DNA was found by both methods. The PRINS method was able to demonstrate HPV DNA not only in superficial koilocytotic squamous cells but also in non-koilocytotic cells in the deeper spinous cell layers, and even in some basal cells. We describe an economical protocol using conventional consensus HPV oligonucleotide DNA primers. The described method is rapid (approximately 3 hours) and easy to perform for screening and subtyping HPV infection in the routine surgical pathology laboratory.

Lower activity of serum peptidases in abstinent alcohol-dependent patients

This study examines i) the activity of serum prolyl endopeptidase (PEP) and dipeptidlyl peptidase IV (DPP IV) in detoxified alcohol-dependent patients without liver disease versus normal controls, and ii) the relationships between serum DPP IV and PEP activity and the production of cytokines or cytokine receptors, such as interleukin-6 (IL-6), tumor necrosis factor-alpha (TNF-alpha), interferon-y (IFN-y), IL-1 receptor antagonist (IL-1RA), and IL-10, and granulocyte-macrophage colony stimulatory factor (GM-CSF). Alcohol-dependent patients had significantly lower serum PEP and DPP IV activity than normal controls. We found that 58.3% and 50.0% of the alcohol-dependent patients, respectively, had PEP and DPP IV activities, which were lower than the mean control values minus 2 SD. There were significant inverse correlations between lowered serum DPP IV and PEP activity and the increased production of IL-6, INF-gamma, IL-IRA, IL-10, and GM-CSF. These results show that lower serum DPP IV and PEP activity may be related to the pathophysiology of alcohol dependence.