The influence of data shape acquisition process and geometric accuracy of the mandible for numerical simulation

Computer-aided technologies have allowed new 3D modelling capabilities and engineering analyses based on experimental and numerical simulation. It has enormous potential for product development, such as biomedical instrumentation and implants. However, due to the complex shapes of anatomical structures, the accuracy of these technologies plays an important key role for adequate and accurate finite element analysis (FEA).

The objective of this study was to determine the influence of the geometry variability between two digital models of a human model of the mandible. Two different shape acquisition techniques, CT scan and 3D laser scan, were assessed. A total of 130 points were controlled and the deviations between the measured points of the physical and 3D virtual models were assessed.

The results of the FEA study showed a relative difference of 20% for the maximum displacement and 10% for the maximum strain between the two geometries.     

Optimizing cardiac material parameters with a genetic algorithm

Determining the unknown material parameters of intact ventricular myocardium can be challenging due to highly nonlinear material behavior. Previous studies combining a gradient-search optimization procedure with finite element analysis (FEA) were limited to two-dimensional (2D) models or simplified three-dimensional (3D) geometries. Here we present a novel scheme to estimate unknown material parameters for ventricular myocardium by combining a genetic algorithm (GA) with nonlinear finite element analysis. This approach systematically explores the domain of the material parameters. The objective function to minimize was the error between simulated strain data and finite element model strains. The proposed scheme was validated for a 2D problem using a realistic material law for ventricular myocardium. Optimized material parameters were generally within 0.5% of the true values. To demonstrate the robustness of the new scheme, unknown material parameters were also determined for a realistic 3D heart model with an exponential hyperelastic material law. When using strains from two material points, the algorithm converged to parameters within 5% of the true values. We conclude that the proposed scheme is robust when estimating myocardial material parameters in 2D and 3D models.     

Use of the photoelastic method and finite element analysis in the assessment of wall strain in abdominal aortic aneurysm models

Abdominal aortic aneurysm (AAA) is a significant health problem. Current clinical rupture-risk relies primarily on the maximum diameter of the AAA and also growth rate. However, AAAs are a patient-specific problem and recently, numerical tools have been employed to estimate rupture-potential. Alternatively, experimental assessment of AAA biomechanics receives less attention, yet, rigorous validation of numerical tools is required prior to clinical acceptance. This paper examines the use of the photoelastic method to assess wall strain and its validation using finite element analysis (FEA) in a small number of patient-specific AAA models. Experimental models were manufactured in-house using the injection-moulding procedure together with a commercially available photoelastic material. The material was mechanically characterised prior to testing, with models examined under three loading regimes (80, 120 and 160mmHg). Each experimental model was imaged using computed tomography (CT) and reconstructed in three dimensions (3D) for numerical analyses. Experimental wall strain was measured and numerical wall strain calculated with finite element analysis (FEA). Results were qualitatively and quantitatively compared. There was good qualitative agreement between the experimental and numerical methods, with similar trends apparent throughout all models at all pressures. Overall, acceptable percentage errors between the techniques were observed for all models. Median errors of -6.5%, -0.4% and 3.9% for the models at 80, 120 and 160mmHg pressures, respectively, were determined. The photoelastic method is a very useful experimental tool that provides instant, easy to interpret, information regarding wall strain. The technique is useful for validation of numerical AAA studies.     

A geometric approach to study the contact mechanisms in the patellofemoral joint of normal versus patellofemoral pain syndrome subjects

The biomechanics of the patellofemoral (PF) joint is complex in nature, and the aetiology of such manifestations of PF instability as patellofemoral pain syndrome (PFPS) is still unclear. At this point, the particular factors affecting PFPS have not yet been determined. This study has two objectives: (1) The first is to develop an alternative geometric method using a three-dimensional (3D) registration technique and linear mapping to investigate the PF joint contact stress using an indirect measure: the depth of virtual penetration (PD) of the patellar cartilage surface into the femoral cartilage surface. (2) The second is to develop 3D PF joint models using the finite element analysis (FEA) to quantify in vivo cartilage contact stress and to compare the peak contact stress location obtained from the FE models with the location of the maximum PD. Magnetic resonance images of healthy and PFPS subjects at knee flexion angles of 15°, 30° and 45° during isometric loading have been used to develop the geometric models. The results obtained from both approaches demonstrated that the subjects with PFPS show higher PD and contact stresses than the normal subjects. Maximum stress and PD increase with flexion angle, and occur on the lateral side in healthy and on the medial side in PFPS subjects. It has been concluded that the alternative geometric method is reliable in addition to being computationally efficient compared with FEA, and has the potential to assess the mechanics of PFPS with an accuracy similar to the FEA.     

Finite element analysis of peri-implant bone volume affected by stresses around Morse taper implants: effects of implant positioning to the bone crest

Abstract

Objectives: The purpose of the present study was to evaluate the distribution and magnitude of stresses through the bone tissue surrounding Morse taper dental implants at different positioning relative to the bone crest. Materials and Methods: A mandibular bone model was obtained from a computed tomography scan. A three-dimensional (3D) model of Morse taper implant-abutment systems placed at the bone crest (equicrestal) and 2?mm bellow the bone crest (subcrestal) were assessed by finite element analysis (FEA). FEA was carried out on axial and oblique (45°) loading at 150 N relatively to the central axis of the implant. The von Mises stresses were analysed considering magnitude and volume of affected peri-implant bone. Results: On vertical loading, maximum von Mises stresses were recorded at 6-7?MPa for trabecular bone while values ranging from 73 up to 118?MPa were recorded for cortical bone. On oblique loading at the equiquestral or subcrestal positioning, the maximum von Mises stresses ranged from 15 to 21?MPa for trabecular bone while values at 150?MPa were recorded for the cortical bone. On vertical loading, >99.9vol.% cortical bone volume was subjected to a maximum of 2?MPa while von Mises stress values at 15?MPa were recorded for trabecular bone. On oblique loading, >99.9vol.% trabecular bone volume was subjected to maximum stress values at 5?MPa, while von Mises stress values at 35?MPa were recorded for >99.4vol.% cortical bone. Conclusions: Bone volume-based stress analysis revealed that most of the bone volume (>99% by vol) was subjected to significantly lower stress values around Morse taper implants placed at equicrestal or subcrestal positioning. Such analysis is commentary to the ordinary biomechanical assessment of dental implants concerning the stress distribution through peri-implant sites.     

Spatiotemporal representation of 3D hand trajectory based on beta-elliptic models

The aim of this paper was to model the hand trajectory during grasping by an extension in 3D of the 2D written language beta-elliptic model. The interest of this model is that it takes into account both geometric and velocity information. The method relies on the decomposition of the task space trajectories in elementary bricks. The latter is characterized by a velocity profile modelled with beta functions and a geometry modelled with elliptic shapes. A data base of grasping movements has been constructed and the errors of reconstruction were assessed (distance and curvature) considering two variations of the beta-elliptic model (‘quarter ellipse’ and ‘two tangents points’ method). The results showed that the method based on two tangent points outperforms the quarter ellipse method with average and maximum relative errors of 2.73% and 8.62%, respectively, and a maximum curvature error of 9.26% for the former. This modelling approach can find interesting application to characterize the improvement due to a rehabilitation or teaching process by a quantitative measurement of hand trajectory parameters.     

Estimation of 3D shape, internal density and mechanics of proximal femur by combining bone mineral density images with shape and density templates

Measurement of bone mineral density (BMD) by dual-energy X-ray absorptiometry (DXA) alone is only a moderate predictor of fracture risk. Finite element analysis (FEA) of bone mechanics, based on DXA images, may improve the prediction of fracture risk. We developed a method to estimate the 3D shape and density distribution of the proximal femur, using a 2D BMD image and a femur shape template. Proximal femurs of eighteen human cadavers were imaged using computed tomography and divided into two sets (N = 9 + 9). The template was created from the samples in first set by using 3D generalized Procrustes analysis and thin-plate splines. Subsequently, the template and 2D BMD image were utilized to estimate the shape and internal density distribution of the femurs in the second set. Finally, FEA was conducted based on the original and the estimated bone models to evaluate the effect of geometrical and density distributional errors on the mechanical strength. The volumetric errors induced by the estimation itself were low (<1.4%). In the estimation of bones in the second set, the mean distance difference between the estimated and the original bone surfaces was 0.80 ± 0.19 mm, suggesting feasible estimation of the femoral shape. The mean absolute error in voxel-by-voxel BMD was 120±8 mg cm?3. In FEA, the stiffness of the proximal femur differed by -7±16% between the original and estimated bones. The present method, in comparison with methods used in previous studies, improved the prediction of the geometry, the BMD distribution and the mechanical characteristics of the proximal femur. Potentially, the proposed method could ultimately improve the determination of bone fracture risk.     

Additive Manufacturing of Anatomical Models from Computed Tomography Scan Data

The purpose of the study presented here was to investigate the manufacturability of human anatomical models from Computed Tomography (CT) scan data via a 3D desktop printer which uses fused deposition modelling (FDM) technology. First, Digital Imaging and Communications in Medicine (DICOM) CT scan data were converted to 3D Standard Triangle Language (STL) format by using InVaselius digital imaging program. Once this STL file is obtained, a 3D physical version of the anatomical model can be fabricated by a desktop 3D FDM printer. As a case study, a patient’s skull CT scan data was considered, and a tangible version of the skull was manufactured by a 3D FDM desktop printer. During the 3D printing process, the skull was built using acrylonitrile-butadiene-styrene (ABS) co-polymer plastic. The printed model showed that the 3D FDM printing technology is able to fabricate anatomical models with high accuracy. As a result, the skull model can be used for preoperative surgical planning, medical training activities, implant design and simulation to show the potential of the FDM technology in medical field. It will also improve communication between medical stuff and patients. Current result indicates that a 3D desktop printer which uses FDM technology can be used to obtain accurate anatomical models.     

Semi-automated mitral valve morphometry and computational stress analysis using 3D ultrasound

In vivo human mitral valves (MV) were imaged using real-time 3D transesophageal echocardiography (rt-3DTEE), and volumetric images of the MV at mid-systole were analyzed by user-initialized segmentation and 3D deformable modeling with continuous medial representation, a compact representation of shape. The resulting MV models were loaded with physiologic pressures using finite element analysis (FEA). We present the regional leaflet stress distributions predicted in normal and diseased (regurgitant) MVs. Rt-3DTEE, semi-automated leaflet segmentation, 3D deformable modeling, and FEA modeling of the in vivo human MV is tenable and useful for evaluation of MV pathology.     

Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study

Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm², which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions.     

Breast volume assessment based on 3D surface geometry: verification of the method using MR imaging

Differences in breast volume and contour are subjectively estimated by surgeons. 3D surface imaging using 3D scanners provides objective breast volume quantification, but precision and accuracy of the method requires verification. Breast volumes of five test individuals were assessed using a 3D surface scanner. Magnetic resonance imaging (MRI) reference volumes were obtained to verify and compare the 3D scan measurements. The anatomical thorax wall curvature was segmented using MRI data and compared to the interpolated curvature of the posterior breast volume delimitation of 3D scan data. MRI showed higher measurement precision, mean deviation (expressed as percentage of volume) of 1.10+/-0.34% compared to 1.63+/-0.53% for the 3D scanner. Mean MRI [right (left) breasts: 638 (629)+/-143 (138) cc] and 3D scan [right (left) breasts: 493 (497)+/-112 (116) cc] breast volumes significantly correlated [right (left) breasts: r=0.982 (0.977), p=0.003 (0.004)]. The posterior thorax wall of the 3D scan model showed high agreement with the MRI thorax wall curvature [mean positive (negative) deviation: 0.33 (-0.17)+/-0.37 cm]. High correspondence and correlation of 3D scan data with MRI-based verifications support 3D surface imaging as sufficiently precise and accurate for breast volume measurements.     

Can sulci protect the brain from traumatic injury?

The influence of sulci in dynamic finite element simulations of the human head has been investigated. First, a detailed 3D FE model was constructed based on an MRI scan of a human head. A second model with a smoothed brain surface was created based on the same MRI scan as the first FE model. These models were validated against experimental data to confirm their human-like dynamic responses during impact. The validated FE models were subjected to several acceleration impulses and the maximum principle strain and strain rate in the brain were analyzed. The results suggested that the inclusion of sulci should be considered for future FE head models as it alters the strain and strain distribution in an FE model.     

Identifying the Functional Flexion-extension Axis of the Knee: An In-Vivo Kinematics Study

Purpose

This study aimed to calculate the flexion-extension axis (FEA) of the knee through in-vivo knee kinematics data, and then compare it with two major anatomical axes of the femoral condyles: the transepicondylar axis (TEA) defined by connecting the medial sulcus and lateral prominence, and the cylinder axis (CA) defined by connecting the centers of posterior condyles.

Methods

The knee kinematics data of 20 healthy subjects were acquired under weight-bearing condition using bi-planar x-ray imaging and 3D-2D registration techniques. By tracking the vertical coordinate change of all points on the surface of femur during knee flexion, the FEA was determined as the line connecting the points with the least vertical shift in the medial and lateral condyles respectively. Angular deviation and distance among the TEA, CA and FEA were measured.

Results

The TEA-FEA angular deviation was significantly larger than that of the CA-FEA in 3D and transverse plane (3.45° vs. 1.98°, p < 0.001; 2.72° vs. 1.19°, p = 0.002), but not in the coronal plane (1.61° vs. 0.83°, p = 0.076). The TEA-FEA distance was significantly greater than that of the CA-FEA in the medial side (6.7 mm vs. 1.9 mm, p < 0.001), but not in the lateral side (3.2 mm vs. 2.0 mm, p = 0.16).

Conclusion

The CA is closer to the FEA compared with the TEA; it can better serve as an anatomical surrogate for the functional knee axis.     

A novel approach to in vivo mitral valve stress analysis

Three-dimensional (3-D) echocardiography allows the generation of anatomically correct and time-resolved geometric mitral valve (MV) models. However, as imaged in vivo, the MV assumes its systolic geometric configuration only when loaded. Customarily, finite element analysis (FEA) is used to predict material stress and strain fields rendered by applying a load on an initially unloaded model. Therefore, this study endeavors to provide a framework for the application of in vivo MV geometry and FEA to MV physiology, pathophysiology, and surgical repair. We hypothesize that in vivo MV geometry can be reasonably used as a surrogate for the unloaded valve in computational (FEA) simulations, yielding reasonable and meaningful stress and strain magnitudes and distributions. Three experiments were undertaken to demonstrate that the MV leaflets are relatively nondeformed during systolic loading: 1) leaflet strain in vivo was measured using sonomicrometry in an ovine model, 2) hybrid models of normal human MVs as constructed using transesophageal real-time 3-D echocardiography (rt-3DE) were repeatedly loaded using FEA, and 3) serial rt-3DE images of normal human MVs were used to construct models at end diastole and end isovolumic contraction to detect any deformation during isovolumic contraction. The average linear strain associated with isovolumic contraction was 0.02 ± 0.01, measured in vivo with sonomicrometry. Repeated loading of the hybrid normal human MV demonstrated little change in stress or geometry: peak von Mises stress changed by <4% at all locations on the anterior and posterior leaflets. Finally, the in vivo human MV deformed minimally during isovolumic contraction, as measured by the mean absolute difference calculated over the surfaces of both leaflets between serial MV models: 0.53 ± 0.19 mm. FEA modeling of MV models derived from in vivo high-resolution truly 3-D imaging is reasonable and useful for stress prediction in MV pathologies and repairs.     

Testing different 3D techniques using geometric morphometrics: Implications for cranial fluctuating asymmetry in humans

This study aimed to test the performance of 3D digitizer, CT scanner, and surface scanner in detecting cranial fluctuating asymmetry. Sets of 32 landmarks (6 in the midline and 13 bilateral) were acquired from 14 archeological crania using a 3D digitizer, and from 3D models generated from a CT scanner and surface scanner using Viewbox 4. Levels of shape variation were analyzed in MorphoJ using Procrustes analysis of variance and Principal component analysis. Intra-observer error accounted for 1.7%, 1.8%, and 4.5% of total shape variation for 3D digitizer, CT scanner, and surface scanner respectively. Fluctuating asymmetry accounted for 15%–16% of total shape variation. Variation between techniques accounted for 18% of total shape variation. We found a higher level of missing landmarks in our surface scan data than for both 3D digitizer and CT scanner data, and both 3D model-based techniques sometimes obscured taphonomic damage. All three 3D techniques are appropriate for measuring cranial fluctuating asymmetry. We advise against combining data collected with different techniques.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

In vivo serial MRI-based models and statistical methods to quantify sensitivity and specificity of mechanical predictors for carotid plaque rupture: location and beyond

It has been hypothesized that mechanical risk factors may be used to predict future atherosclerotic plaque rupture. Truly predictive methods for plaque rupture and methods to identify the best predictor(s) from all the candidates are lacking in the literature. A novel combination of computational and statistical models based on serial magnetic resonance imaging (MRI) was introduced to quantify sensitivity and specificity of mechanical predictors to identify the best candidate for plaque rupture site prediction. Serial in vivo MRI data of carotid plaque from one patient was acquired with follow-up scan showing ulceration. 3D computational fluid-structure interaction (FSI) models using both baseline and follow-up data were constructed and plaque wall stress (PWS) and strain (PWSn) and flow maximum shear stress (FSS) were extracted from all 600 matched nodal points (100 points per matched slice, baseline matching follow-up) on the lumen surface for analysis. Each of the 600 points was marked "ulcer" or "nonulcer" using follow-up scan. Predictive statistical models for each of the seven combinations of PWS, PWSn, and FSS were trained using the follow-up data and applied to the baseline data to assess their sensitivity and specificity using the 600 data points for ulcer predictions. Sensitivity of prediction is defined as the proportion of the true positive outcomes that are predicted to be positive. Specificity of prediction is defined as the proportion of the true negative outcomes that are correctly predicted to be negative. Using probability 0.3 as a threshold to infer ulcer occurrence at the prediction stage, the combination of PWS and PWSn provided the best predictive accuracy with (sensitivity, specificity)?=?(0.97, 0.958). Sensitivity and specificity given by PWS, PWSn, and FSS individually were (0.788, 0.968), (0.515, 0.968), and (0.758, 0.928), respectively. The proposed computational-statistical process provides a novel method and a framework to assess the sensitivity and specificity of various risk indicators and offers the potential to identify the optimized predictor for plaque rupture using serial MRI with follow-up scan showing ulceration as the gold standard for method validation. While serial MRI data with actual rupture are hard to acquire, this single-case study suggests that combination of multiple predictors may provide potential improvement to existing plaque assessment schemes. With large-scale patient studies, this predictive modeling process may provide more solid ground for rupture predictor selection strategies and methods for image-based plaque vulnerability assessment.     

Dosimetric characterization of 3D printed bolus at different infill percentage for external photon beam radiotherapy

Background and purpose3D printing is rapidly evolving and further assessment of materials and technique is required for clinical applications. We evaluated 3D printed boluses with acrylonitrile butadiene styrene (ABS) and polylactide (PLA) at different infill percentage.Material and methodsA low-cost 3D printer was used. The influence of the air inclusion within the 3D printed boluses was assessed thoroughly both with treatment planning system (TPS) and with physical measurements. For each bolus, two treatment plans were calculated with Monte Carlo algorithm, considering the computed tomography (CT) scan of the 3D printed bolus or modelling the 3D printed bolus as a virtual bolus structure with a homogeneous density. Depth dose measurements were performed with Gafchromic films.ResultsHigh infill percentage corresponds to high density and high homogeneity within bolus material. The approximation of the bolus in the TPS as a homogeneous material is satisfying for infill percentages greater than 20%. Measurements performed with PLA boluses are more comparable to the TPS calculated profiles. For boluses printed at 40% and 60% infill, the discrepancies between calculated and measured dose distribution are within 5%.Conclusions3D printing technology allows modulating the shift of the build-up region by tuning the infill percentage of the 3D printed bolus in order to improve superficial target coverage.     

A customized fixation plate with novel structure designed by topological optimization for mandibular angle fracture based on finite element analysis

Background

The purpose of this study was to design a customized fixation plate for mandibular angle fracture using topological optimization based on the biomechanical properties of the two conventional fixation systems, and compare the results of stress, strain and displacement distributions calculated by finite element analysis (FEA).

Methods

A three-dimensional (3D) virtual mandible was reconstructed from CT images with a mimic angle fracture and a 1 mm gap between two bone segments, and then a FEA model, including volume mesh with inhomogeneous bone material properties, three loading conditions and constraints (muscles and condyles), was created to design a customized plate using topological optimization method, then the shape of the plate was referenced from the stress concentrated area on an initial part created from thickened bone surface for optimal calculation, and then the plate was formulated as “V” pattern according to dimensions of standard mini-plate finally. To compare the biomechanical behavior of the “V” plate and other conventional mini-plates for angle fracture fixation, two conventional fixation systems were used: type A, one standard mini-plate, and type B, two standard mini-plates, and the stress, strain and displacement distributions within the three fixation systems were compared and discussed.

Results

The stress, strain and displacement distributions to the angle fractured mandible with three different fixation modalities were collected, respectively, and the maximum stress for each model emerged at the mandibular ramus or screw holes. Under the same loading conditions, the maximum stress on the customized fixation system decreased 74.3, 75.6 and 70.6% compared to type A, and 34.9, 34.1, and 39.6% compared to type B. All maximum von Mises stresses of mandible were well below the allowable stress of human bone, as well as maximum principal strain. And the displacement diagram of bony segments indicated the effect of treatment with different fixation systems.

Conclusions

The customized fixation system with topological optimized structure has good biomechanical behavior for mandibular angle fracture because the stress, strain and displacement within the plate could be reduced significantly comparing to conventional “one mini-plate” or “two mini-plates” systems. The design methodology for customized fixation system could be used for other fractures in mandible or other bones to acquire better mechanical behavior of the system and improve stable environment for bone healing. And together with SLM, the customized plate with optimal structure could be designed and fabricated rapidly to satisfy the urgent time requirements for treatment.

     

Validation of a 3D computational fluid-structure interaction model simulating flow through an elastic aperture

This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two.