Cerebral aneurysm flow diverter modeled as a thin inhomogeneous porous medium in hemodynamic simulations.

Rapid and accurate simulation of cerebral aneurysm flow modifications by flow diverters (FDs) can help improving patient-specific intervention and predicting treatment outcome. However, when FD devices are explicitly represented in computational fluid dynamics (CFD) simulations, flow around the stent wires must be resolved, leading to high computational cost. Classic porous medium (PM) methods can reduce computational expense but cannot capture the inhomogeneous FD wire distribution once implanted on a cerebral artery and thus cannot accurately model the post-stenting aneurysmal flow. We report a novel approach that models the FD flow modification as a thin inhomogeneous porous medium (iPM). It improves over the classic PM approaches in two ways. First, the FD is more appropriately treated as a thin screen, which is more accurate than the classic 3D-PM-based Darcy-Forchheimer relation. Second, pressure drop is calculated cell-by-cell using the local FD geometric parameters across an inhomogeneous PM. We applied the iPM technique to simulating the post-stenting hemodynamics of three patient-specific aneurysms. To test its accuracy and speed, we compared the results from the iPM technique against CFD simulations with explicit FD devices. The iPM CFD ran 500% faster than the explicit CFD while achieving 94%-99% accuracy; thus, iPM is a promising clinical bedside modeling tool to assist endovascular interventions with FD and stents.

Improving accuracy for finite element modeling of endovascular coiling of intracranial aneurysm.

BACKGROUND: Computer modeling of endovascular coiling intervention for intracranial aneurysm could enable a priori patient-specific treatment evaluation. To that end, we previously developed a finite element method (FEM) coiling technique, which incorporated simplified assumptions. To improve accuracy in capturing real-life coiling, we aimed to enhance the modeling strategies and experimentally test whether improvements lead to more accurate coiling simulations. METHODS: We previously modeled coils using a pre-shape based on mathematical curves and mechanical properties based on those of platinum wires. In the improved version, to better represent the physical properties of coils, we model coil pre-shapes based on how they are manufactured, and their mechanical properties based on their spring-like geometric structures. To enhance the deployment mechanics, we include coil advancement to the aneurysm in FEM simulations. To test if these new strategies produce more accurate coil deployments, we fabricated silicone phantoms of 2 patient-specific aneurysms in duplicate, deployed coils in each, and quantified coil distributions from intra-aneurysmal cross-sections using coil density (CD) and lacunarity (L). These deployments were simulated 9 times each using the original and improved techniques, and CD and L were calculated for cross-sections matching those in the experiments. To compare the 2 simulation techniques, Euclidean distances (dMin, dMax, and dAvg) between experimental and simulation points in standardized CD-L space were evaluated. Univariate tests were performed to determine if these distances were significantly different between the 2 simulations. RESULTS: Coil deployments using the improved technique agreed better with experiments than the original technique. All dMin, dMax, and dAvg values were smaller for the improved technique, and the average values across all simulations for the improved technique were significantly smaller than those from the original technique (dMin: p = 0.014, dMax: p = 0.013, dAvg: p = 0.045). CONCLUSION: Incorporating coil-specific physical properties and mechanics improves accuracy of FEM simulations of endovascular intracranial aneurysm coiling.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements.

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally-derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound [Formula: see text] over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Size-dependent strength of dental adhesive systems.

OBJECTIVE: The aim of this study is to explain the influence of peripheral interface stress singularities on the testing of tensile bond strength. The relationships between these theoretically predicted singularities and the effect of specimen size on the measured bond strength are evaluated. METHODS: Finite element method (FEM) and boundary element method (BEM) analyses of microtensile bond strength test specimens were performed and the presence of localized high stress concentrations and singularities was analyzed. The specimen size effect predicted by the models was compared to previously published experimental data. RESULTS: FEM analysis of single-material trimmed hour-glass versus cast cylindrical specimens showed different theoretical stress distributions, with the dumbbell or cylindrical specimens showing a more homogeneous distribution of the stress on the critical symmetry plane. For multi-material specimens, mathematical singularities at the free edge of the bonded interface posed a computational challenge that resulted in mesh-dependence in the standard FEM analysis. A specialized weighted-traction BEM analysis, designed to eliminate mesh-dependence by capturing the effect of the singularity, predicted a specimen size effect that corresponds to that published previously in the literature. SIGNIFICANCE: The results presented here further support the attention to specimen dimensions that has already broadened the empirical use of the microtensile test methods. FEM and BEM analyses that identify stress concentrations and especially marginal stress singularities must be accounted for in reliable bonding strength assessments. Size-dependent strength variations generally attributed to the effects of flaw distributions throughout the interfacial region are not as relevant as the presence of singularities at bonded joint boundaries - as revealed by both FEM and BEM analyses, when interpreted from a generalized fracture mechanics perspective. Furthermore, this size-dependence must be considered when evaluating or designing dental adhesive systems.

Computer modeling of deployment and mechanical expansion of neurovascular flow diverter in patient-specific intracranial aneurysms.

Flow diverter (FD) is an emerging neurovascular device based on self-expandable braided stent for treating intracranial aneurysms. Variability in FD outcome has underscored a need for investigating the hemodynamic effect of fully deployed FD in patient-specific aneurysms. Image-based computational fluid dynamics, which can provide important hemodynamic insight, requires accurate representation of FD in deployed states. We developed a finite element analysis (FEA) based workflow for simulating mechanical deployment of FD in patient-specific aneurysms. We constructed FD models of interlaced wires emulating the Pipeline Embolization Device, using 3D finite beam elements to account for interactions between stent strands, and between the stent and other components. The FEA analysis encompasses all steps that affect the final deployed configuration including stent crimping, delivery and expansion. Besides the stent, modeling also includes key components of the FD delivery system such as microcatheter, pusher, and distal coil. Coordinated maneuver of these components allowed the workflow to mimic clinical operation of FD deployment and to explore clinical strategies. The workflow was applied to two patient-specific aneurysms. Parametric study indicated consistency of the deployment result against different friction conditions, but excessive intra-stent friction should be avoided. This study demonstrates for the first time mechanical modeling of braided FD stent deployment in cerebral vasculature to produce realistic deployed configuration, thus paving the way for accurate CFD analysis of flow diverters for reliable prediction and optimization of treatment outcome.