Dragomir‐Daescu, Dan

AbstractComputational fluid dynamics (CFD) has been studied as a tool for the stratification of aneurysm rupture risk. We performed CFD analysis in a patient operated on for a ruptured anterior communicating artery aneurysm. The point of rupture was identified during surgery. The aneurysm and blood vessels were segmented from computed tomography angiography to prepare a model for simulations. We found that the streamlines showed a concentrated inflow jet directed straight at the rupture point, and high wall shear stress was found at the point of rupture in the aneurysm sac. Thus specific local hemodynamics may be indicative of the aneurysm rupture site.

Quantitative computed tomography-based finite-element analysis (QCT/FEA) has become increasingly popular in an attempt to understand and possibly reduce vertebral fracture risk. It is known that scanning acquisition settings affect Hounsfield units (HU) of the CT voxels. Material properties assignments in QCT/FEA, relating HU to Young's modulus, are performed by applying empirical equations. The purpose of this study was to evaluate the effect of QCT scanning protocols on predicted stiffness values from finite-element models. One fresh frozen cadaveric torso and a QCT calibration phantom were scanned six times varying voltage and current and reconstructed to obtain a total of 12 sets of images. Five vertebrae from the torso were experimentally tested to obtain stiffness values. QCT/FEA models of the five vertebrae were developed for the 12 image data resulting in a total of 60 models. Predicted stiffness was compared to the experimental values. The highest percent difference in stiffness was approximately 480% (80 kVp, 110 mAs, U70), while the lowest outcome was ∼1% (80 kVp, 110 mAs, U30). There was a clear distinction between reconstruction kernels in predicted outcomes, whereas voltage did not present a clear influence on results. The potential of QCT/FEA as an improvement to conventional fracture risk prediction tools is well established. However, it is important to establish research protocols that can lead to results that can be translated to the clinical setting.

With the increased availability of computational resources, the past decade has seen a rise in the use of computational fluid dynamics (CFD) for medical applications. There has been an increase in the application of CFD to attempt to predict the rupture of intracranial aneurysms, however, while many hemodynamic parameters can be obtained from these computations, to date, no consistent methodology for the prediction of the rupture has been identified. One particular challenge to CFD is that many factors contribute to its accuracy; the mesh resolution and spatial/temporal discretization can alone contribute to a variation in accuracy. This failure to identify the importance of these factors and identify a methodology for the prediction of ruptures has limited the acceptance of CFD among physicians for rupture prediction. The International CFD Rupture Challenge 2013 seeks to comment on the sensitivity of these various CFD assumptions to predict the rupture by undertaking a comparison of the rupture and blood-flow predictions from a wide range of independent participants utilizing a range of CFD approaches. Twenty-six groups from 15 countries took part in the challenge. Participants were provided with surface models of two intracranial aneurysms and asked to carry out the corresponding hemodynamics simulations, free to choose their own mesh, solver, and temporal discretization. They were requested to submit velocity and pressure predictions along the centerline and on specified planes. The first phase of the challenge, described in a separate paper, was aimed at predicting which of the two aneurysms had previously ruptured and where the rupture site was located. The second phase, described in this paper, aims to assess the variability of the solutions and the sensitivity to the modeling assumptions. Participants were free to choose boundary conditions in the first phase, whereas they were prescribed in the second phase but all other CFD modeling parameters were not prescribed. In order to compare the computational results of one representative group with experimental results, steady-flow measurements using particle image velocimetry (PIV) were carried out in a silicone model of one of the provided aneurysms. Approximately 80% of the participating groups generated similar results. Both velocity and pressure computations were in good agreement with each other for cycle-averaged and peak-systolic predictions. Most apparent “outliers” (results that stand out of the collective) were observed to have underestimated velocity levels compared to the majority of solutions, but nevertheless identified comparable flow structures. In only two cases, the results deviate by over 35% from the mean solution of all the participants. Results of steady CFD simulations of the representative group and PIV experiments were in good agreement. The study demonstrated that while a range of numerical schemes, mesh resolution, and solvers was used, similar flow predictions were observed in the majority of cases. To further validate the computational results, it is suggested that time-dependent measurements should be conducted in the future. However, it is recognized that this study does not include the biological aspects of the aneurysm, which needs to be considered to be able to more precisely identify the specific rupture risk of an intracranial aneurysm.

Osteoporosis is characterized by bony material loss and decreased bone strength leading to a significant increase in fracture risk. Patient-specific quantitative computed tomography (QCT) finite element (FE) models may be used to predict fracture under physiological loading. Material properties for the FE models used to predict fracture are obtained by converting grayscale values from the CT into volumetric bone mineral density (vBMD) using calibration phantoms. If there are any variations arising from the CT acquisition protocol, vBMD estimation and material property assignment could be affected, thus, affecting fracture risk prediction. We hypothesized that material property assignments may be dependent on scanning and postprocessing settings including voltage, current, and reconstruction kernel, thus potentially having an effect in fracture risk prediction. A rabbit femur and a standard calibration phantom were imaged by QCT using different protocols. Cortical and cancellous regions were segmented, their average Hounsfield unit (HU) values obtained and converted to vBMD. Estimated vBMD for the cortical and cancellous regions were affected by voltage and kernel but not by current. Our study demonstrated that there exists a significant variation in the estimated vBMD values obtained with different scanning acquisitions. In addition, the large noise differences observed utilizing different scanning parameters could have an important negative effect on small subregions containing fewer voxels.

ABSTRACT Although a large number of studies have addressed the age-related changes in bone mineral density (BMD), there is a paucity of data for the assessment of femoral strength loss with age in both genders. We determined the variation of strength with age in femurs of women and men by mechanical tests on a cohort of 100 cadaveric femurs. In addition, the age-related neck BMD loss in our cadaveric cohort was found to be similar with BMD loss of four published population-based studies. Given the strong correlation found in our cadaveric study between BMD and femoral strength, we also estimated the femoral strength of the four populations based on their reported neck BMDs. Our study showed that men's femurs in our cadaveric cohort were stronger than women's femurs by about 800 N at the same BMD level, and by 1750 N at the same age. The strength differences were not explained satisfactorily by the size difference between men's and women's bones. Similar to the findings of clinical studies, the BMD values of men at all ages were larger than that of women. The age-related loss rates in BMD and strength were not statistically different between the two genders of our cadaveric cohort. After normalization, strength decreased more than 40% faster than BMD. On average, men reached a certain BMD value about 16 years later than women, and for strength about 23 years later, which may explain the higher rate of hip fracture in postmenopausal women. In patient population cohorts men reached a similar BMD value about 16 to 25 years later than women, whereas for estimated strength, sometimes more than 40 years later. © 2015 American Society for Bone and Mineral Research.

We present a quantified description of the fluid flow and a novel flowline-based meshing technique to create adaptive grids for Computational Fluid Dynamics (CFD) simulations in patient-specific intracranial aneurysms. The adaptive grid density is obtained such that it captures the fine geometrical details of the flow with high grid density, while smoother flow characteristics are calculated with a coarser grid density. The correlation between the topological characteristics of the flow and the element size of the adaptive grid results in a practical mathematical formula for calculating the element size using only one uniform base mesh and a user defined implementation in CFD post processors.

The Modeling and Analysis group in the Division of Engineering at Mayo Clinic, Rochester works in the area of numerical modeling and experimental validation of cerebral aneurysms. In our modeling we use patient-specific geometries from imaging studies provided by Dr. Kallmes. The CFD analysis is performed with commercial packages Mimics (Materialise, Leuven, Belgium), for segmentation and ANSYS (ANSYS Inc. Canonsburg, PA) for meshing, simulation and post-processing. The experiments are conducted using a tomographic Particle Image Velocimetry system from LaVision Inc. The group is led by Dan Dragomir-Daescu, PhD and David F. Kallmes, MD and consists of two postdoctoral fellows and three engineers. For this project, CFD simulations were performed by postdoctoral fellow Simona Hodis, using the finite volume based solver Fluent.

We present a computational fluid dynamics (CFD) analysis of the hemodynamic environment of an anterior communicating artery that spontaneously ruptured immediately following three-dimensional rotational angiography. Subsequent digital subtraction angiography allowed for the localization of the point of rupture within the aneurysm dome. CFD analysis demonstrated a concentrated jet that impinged directly at the site of rupture. Peak systolic pressure and wall shear stress were both maximal near the rupture location.

Convection induced by repetitive compression of porous tissue scaffolds enhances solute transport inside the scaffold. Our previous experiments have shown that pore size, shape, and orientation with respect to strain direction greatly influence loading-induced solute transport. The objective of this study was to develop a computational model of deformation-induced solute transport in porous tissue scaffolds, which included the pore geometry of the scaffold. This geometry consisted of a cubic scaffold with single channel in the middle of the scaffold, immersed in a fluid reservoir. Cylindrical pores with circular or elliptic cross section, and spheroid pores were modeled. The scaffold was cyclically compressed from one side, causing fluid motion and dispersion of solute inside the scaffold pore. Scaffold deformation was solved using the finite element method, and fluid flow and solute transport were solved using the finite volume method. The distortion of the scaffold–fluid interface was transferred as a boundary condition to the fluid flow solver. Both convection and diffusion were included in the computations. The solute transport rates in the different scaffold pore geometries agreed well with our previous experimental results obtained with X-ray microimaging. This model will be used to explore transport properties of a spectrum of novel scaffold designs.

The objective of this study was to investigate the influence of pore geometry on the transport rate and depth after repetitive mechanical deformation of porous scaffolds for tissue engineering applications. Flexible cubic imaging phantoms with pores in the shape of a circular cylinder, elliptic cylinder, and spheroid were fabricated from a biodegradable polymer blend using a combined 3D printing and injection molding technique. The specimens were immersed in fluid and loaded with a solution of a radiopaque solute. The solute distribution was quantified by recording 20 μm pixel-resolution images in an X-ray microimaging scanner at selected time points after intervals of dynamic straining with a mean strain of 8.6 ± 1.6% at 1.0 Hz. The results show that application of cyclic strain significantly increases the rate and depth of solute transport, as compared to diffusive transport alone, for all pore shapes. In addition, pore shape, pore size, and the orientation of the pore cross-sectional asymmetry with respect to the direction of strain greatly influence solute transport. Thus, pore geometry can be tailored to increase transport rates and depths in cyclically deformed scaffolds, which is of utmost importance when thick, metabolically functional tissues are to be engineered.