A virtual simulation approach to assess the effect of trocar-site placement and scar characteristics on the abdominal wall biomechanics.

Analyses of registries and medical imaging suggest that laparoscopic surgery may be penalized with a high incidence of trocar-site hernias (TSH). In addition to trocar diameter, the location of the surgical wound (SW) may affect TSH incidence. The intra-abdominal pressure (IAP) exerted on the abdominal wall (AW) might also influence the appearance of TSH. In the present study, we used finite element (FE) simulations to predict the influence of trocar location and SW characteristics (stiffness) on the mechanical behavior of the AW subject to an IAP. Two models of laparoscopy patterns on the AW, with trocars in the 5-12 mm range, were generated. FE simulations for IAP values within the 4 kPa-20 kPa range were carried out using the Code Aster open-source software. Different stiffness levels of the SW tissue were considered. We found that midline-located surgical wounds barely deformed, even though they moved outwards along with the regular LA tissue. Laterally located SWs hardly changed their location but they experienced significant variations in their volume and shape. The amount of deformation of lateral SWs was found to strongly depend on their stiffness. Trocar incisions placed in a LA with non-diastatic dimensions do not compromise its mechanical integrity. The more lateral the trocars are placed, the greater is their deformation, regardless of their size. Thus, to prevent TSH it might be advisable to close lateral trocars with a suture, or even use a prosthetic reinforcement depending on the patient's risk factors (e.g., obesity).

The effect of thoracic dimensions on compression depth during cardiopulmonary resuscitation.

The effect of the dimensions of the thoracic cage on the resuscitation outcome of cardiopulmonary resuscitation (CPR) maneuvers has long been debated. In this study, the effect of changes in the rib cage dimensions on the achieved compression depth was investigated using finite element simulations. A total of 216 different rib cage geometry models were considered and, in each case, the result of applying different levels of compression force up to 600 N were simulated. The Haller Index of the rib cage is defined as the ratio of the transverse diameter and the antero-posterior diameter. Our results suggest that, with a fixed level of compression force, performing CPR on rib cages having a low Haller Index and/or a larger height leads to compression depths below the average. Alternatively, if a target compression depth is set for CPR, in general a lower compression force would be required for individuals with higher Haller Index and/or lower chest height. In addition, present results indicate that wider chested individuals will experience lower stress levels on their ribs to achieve the required CPR target depth. Moreover, in the present study we propose predictive models, based on anthropometric parameters, for compression depth and rib stress during chest compressions. In particular, the model suggests that in future correlations of empirical CPR data the patients' Haller index and vertical (sagittal) cross-area are the best parameters to be used as independent variables in a fit.

Virtual simulation of the biomechanics of the abdominal wall with different stoma locations.

An ostomy is a surgical procedure by which an artificial opening in the abdominal wall, known as a stoma, is created. We assess the effects of stoma location on the abdominal wall mechanics. We perform three-dimensional finite element simulations on an anatomy model which was generated on the basis of medical images. Our simulation methodology is entirely based on open source software. We consider seventeen different locations for the stoma incision (trephine) and we simulate the mechanical response of the abdominal wall when an intraabdominal pressure as high as 20 kPa is applied. We focus on factors related to the risk of parastomal hernia development such as the deformation experienced by the abdominal wall, the stress levels supported by its tissues and the corresponding level of trephine enlargement. No significant dependence was found between stoma location and the levels of abdominal wall deformations or stress supported by tissues, except for the case with a stoma located on the linea alba. Trephine perimeter and area respectively increased by as much as [Formula: see text] and [Formula: see text]. The level of trephine deformation depends on stoma location with considerably higher trephine enlargements found in stomas laterally located with respect to the rectus abdominis muscle.

Biomechanical modelling of the pelvic system: improving the accuracy of the location of neoplasms in MRI-TRUS fusion prostate biopsy.

BACKGROUND: An accurate knowledge of the relocation of prostate neoplasms during biopsy is of great importance to reduce the number of false negative results. Prostate neoplasms are visible in magnetic resonance images (MRI) but it is difficult for the practitioner to locate them at the time of performing a transrectal ultrasound (TRUS) guided biopsy. In this study, we present a new methodology, based on simulation, that predicts both prostate deformation and lesion migration during the biopsy. METHODS: A three-dimensional (3-D) anatomy model of the pelvic region, based on medical images, is constructed. A finite element (FE) numerical simulation of the organs motion and deformation as a result of the pressure exerted by the TRUS probe is carried out using the Code-Aster open-source computer software. Initial positions of potential prostate lesions prior to biopsy are taken into consideration and the final location of each lesion is targeted in the FE simulation output. RESULTS: Our 3-D FE simulations show that the effect of the pressure exerted by the TRUS probe is twofold as the prostate experiences both a motion and a deformation of its original shape. We targeted the relocation of five small prostate lesions when the TRUS probe exerts a force of 30 N on the rectum inner wall. The distance travelled by these lesions ranged between 5.6 and 13.9 mm. CONCLUSIONS: Our new methodology can help to predict the location of neoplasms during a prostate biopsy but further studies are needed to validate our results. Moreover, the new methodology is completely developed on open-source software, which means that its implementation would be affordable to all healthcare providers.

Biomechanical response of human rib cage to cardiopulmonary resuscitation maneuvers: Effects of the compression location.

The biomechanical response of a human rib cage to cardiopulmonary resuscitation maneuvers was investigated by means of finite element simulations. We analyzed the effect of the location where the force was applied on the achieved compression depths and stress levels experienced by the breastbone and ribs. For compression locations on the breastbone, a caudal shift of the application area toward the breastbone tip resulted in a 17% reduction of the force required to achieve a target 5 cm compression depth. We found that the use of compression regions located on the costal cartilages would involve higher risk of rib fractures.

Hemodynamic effects of blood clots trapped by an inferior vena cava filter.

The alteration of blood flow around an OPTEASE inferior vena cava filter with one or two blood clots attached was investigated by means of computational fluid dynamics. We used a patient-specific vein wall geometry, and we generated different clot models with shapes adapted to the filter and vein wall geometries. A total of eight geometries, with one or two clots and a total clot volume of 0.5 or 1 cm3 , were considered. A non-Newtonian model for blood viscosity was adopted and the possible development of turbulence was accounted for by means of a three-equation model. Two blood flow rates were considered for each case, representative for rest and exercise conditions. In exercise conditions, flow unsteadiness and even turbulence was detected in some cases. Pressure and wall shear stress (WSS) distributions were modified in all cases. Clots attached to the filter downstream basket considerably increased averaged WSS values by up to almost 50%. In all the cases a flow recirculation region appeared downstream of the clot. The degree of flow stagnation in these regions, an indicator of propensity to thrombogenesis, was estimated in terms of mean residence times and mean blood viscosity. High levels of flow stagnation were detected in rest conditions in the wake of those clots that were placed upstream from the filter. Our results suggest that one downstream placed big clot, showing a higher tendency to induce flow instabilities and turbulence, might be more harmful than two small clots placed in tandem.

Implementation of a new constitutive model for abdominal muscles.

BACKGROUND AND OBJECTIVE: Abdominal hernia repair is one of the most often performed surgical procedures worldwide. Numerical simulations of the abdominal wall mechanics can be a valuable tool to devise actions aimed at preventing hernia formation. A first step towards this goal is the development of consistent constitutive models for the tissues that form the human abdominal wall. In this study we propose, for each of the tissues involved, a new formulation of the so-called transversely isotropic hyperelastic model (TIHM). METHODS: We propose a new TIHM for the human abdominal wall tissues and we present a systemic view of the methodology that we have implemented in the present study. First we consider the mathematical background of the TIHM. The novelty of our formulation is that both the isotropic and the fiber contributions to the strain energy function are characterized exclusively by polynomial convex functions of certain invariant quantities. Then, we provide a detailed description on how the constitutive model is implemented into an open source finite element (FE) software. In our approach we use the specific interface provided by the MFront software to incorporate our TIHM formulation into the Code Aster FE solver. For each of the tissues considered, the values of the TIHM constants are adjusted by means of a numerical simulation of previous experimental data from tensile tests. RESULTS: We studied the following abdominal wall tissues: linea alba, rectus sheath, external oblique muscle, internal oblique muscle, transversus abdominis muscle and rectus abdominis muscle. Our formulation closely reproduces tensile test data for each tissue in the corresponding FE numerical simulation. CONCLUSIONS: The new TIHM formulation is suitable for a future numerical investigation of the abdominal wall, which will in turn help us to assess the best zone to practice a colostomy. The methodology implemented in the present study can be easily extended in the future to develop and implement a TIHM for active muscles and/or a different type of constitutive model which might be suitable to characterize other tissues of biomedical interest.

A comparative CFD study of four inferior vena cava filters.

Computational fluid dynamics was used to simulate the flow of blood within an inferior vena cava (IVC) geometry model that was reconstructed from computed tomography images obtained from a real patient. The main novelty of the present work is that we simulated the implantation of 4 different filter models in this realistic IVC geometry. We considered different blood flow rates in the range between Vin =20 and Vin =80 cm3 /s, and all simulations were performed with both the Newtonian and a non-Newtonian model for the blood viscosity. We compared the hemodynamics performance of the different filter models, and we paid a special attention to the total drag force, Fd , exerted by the blood flow on the filter surface. This force is the sum of 2 contributions: the viscous skin friction force, which was found to be roughly proportional to the filter surface area, and the pressure force, which depended on the particular filter geometry design. The Fd force is relevant because it must be balanced by the total force exerted by the filter hooks/struts on the IVC wall at the attachment locations. For the highest Vin value investigated, the variation in Fd among filters was from 116 to 308 dyne. We also showed how the present results can be extrapolated to obtain good estimates of the drag forces if the blood viscosity levels change, ie, if the patient with a filter implanted is treated with anticoagulant therapy.

Effects of walking in deep venous thrombosis: a new integrated solid and fluid mechanics model.

Deep venous thrombosis (DVT) is a common disease. Large thrombi in venous vessels cause bad blood circulation and pain; and when a blood clot detaches from a vein wall, it causes an embolism whose consequences range from mild to fatal. Walking is recommended to DVT patients as a therapeutical complement. In this study the mechanical effects of walking on a specific patient of DVT were simulated by means of an unprecedented integration of 3 elements: a real geometry, a biomechanical model of body tissues, and a computational fluid dynamics study. A set of computed tomography images of a patient's leg with a thrombus in the popliteal vein was employed to reconstruct a geometry model. Then a biomechanical model was used to compute the new deformed geometry of the vein as a function of the fiber stretch level of the semimembranosus muscle. Finally, a computational fluid dynamics study was performed to compute the blood flow and the wall shear stress (WSS) at the vein and thrombus walls. Calculations showed that either a lengthening or shortening of the semimembranosus muscle led to a decrease of WSS levels up to 10%. Notwithstanding, changes in blood viscosity properties or blood flow rate may easily have a greater impact in WSS.

Effect of anticoagulant treatment in deep vein thrombosis: A patient-specific computational fluid dynamics study.

A methodology that might help physicians to establish a diagnostic and treatment tailored for each specific patient with a pathological thrombus is presented. A realistic model for the geometry of a popliteal vein with a thrombus just above the knee was reconstructed from in vivo computed tomography images acquired from one specific patient and then it was used to perform computational fluid dynamics (CFD) simulations. The wall shear stress (WSS) response to the administration of anticoagulant drugs and the influence of viscosity on the shape of the velocity distribution were investigated. Both a Newtonian and a non-Newtonian viscosity model were implemented for different blood flow rates in the range 3-7 cm(3)/s. The effect of anticoagulants on the blood was simulated by setting three different levels of viscosity in the Newtonian model (µ/µ∞=0.60, 0.80 and 1 with µ∞=3.45×10(-3) Pas). A reduction of µ by a given amount always led to a more modest reduction, typically by a factor of two, of the resulting WSS levels. Moreover, for a given flow rate the calculation with the non-Newtonian viscosity model yielded WSS levels between 20% and 40% larger than those obtained in the corresponding Newtonian fluid simulation. It was also found that blood moves slowly in the region between the thrombus and the vein wall, a fact that will favor the growth of the thrombotic mass. Both the mean WSS levels and the degree of sluggishness of the blood flow can be described by functions of the Reynolds number.