Free Tools and Strategies for the Generation of 3D Finite Element Meshes: Modeling of the Cardiac Structures.
Ontology highlight
ABSTRACT:
SUBMITTER: Pavarino E
PROVIDER: S-EPMC3670575 | biostudies-literature | 2013
REPOSITORIES: biostudies-literature
Publications
Free Tools and Strategies for the Generation of 3D Finite Element Meshes: Modeling of the Cardiac Structures.
Pavarino E E Neves L A LA Machado J M JM de Godoy M F MF Shiyou Y Y Momente J C JC Zafalon G F D GF Pinto A R AR Valêncio C R CR
International journal of biomedical imaging 20130516
The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocess ...[more]
Similar Datasets
Project description:IntroductionMitral Valve (MV) 3D structural data can be easily obtained using standard transesophageal echocardiography (TEE) devices but quantitative pre- and intraoperative volume analysis of the MV is presently not feasible in the cardiac operation room (OR). Finite element method (FEM) modelling is necessary to carry out precise and individual volume analysis and in the future will form the basis for simulation of cardiac interventions.MethodWith the present retrospective pilot study we describe a method to transfer MV geometric data to 3D Slicer 2 software, an open-source medical visualization and analysis software package. A newly developed software program (ROIExtract) allowed selection of a region-of-interest (ROI) from the TEE data and data transformation for use in 3D Slicer. FEM models for quantitative volumetric studies were generated.ResultsROI selection permitted the visualization and calculations required to create a sequence of volume rendered models of the MV allowing time-based visualization of regional deformation. Quantitation of tissue volume, especially important in myxomatous degeneration can be carried out. Rendered volumes are shown in 3D as well as in time-resolved 4D animations.ConclusionThe visualization of the segmented MV may significantly enhance clinical interpretation. This method provides an infrastructure for the study of image guided assessment of clinical findings and surgical planning. For complete pre- and intraoperative 3D MV FEM analysis, three input elements are necessary: 1. time-gated, reality-based structural information, 2. continuous MV pressure and 3. instantaneous tissue elastance. The present process makes the first of these elements available. Volume defect analysis is essential to fully understand functional and geometrical dysfunction of but not limited to the valve. 3D Slicer was used for semi-automatic valve border detection and volume-rendering of clinical 3D echocardiographic data. FEM based models were also calculated.MethodA Philips/HP Sonos 5500 ultrasound device stores volume data as time-resolved 4D volume data sets. Data sets for three subjects were used. Since 3D Slicer does not process time-resolved data sets, we employed a standard movie maker to animate the individual time-based models and visualizations. Calculation time and model size were minimized. Pressures were also easily available. We speculate that calculation of instantaneous elastance may be possible using instantaneous pressure values and tissue deformation data derived from the animated FEM.
| S-EPMC1421418 | biostudies-literature
Project description:Masticatory overload on dental implants is one of the causes of marginal bone resorption. The implant-abutment connection (IAC) design plays a critical role in the quality of the stress distribution, and, over the years, different designs were proposed. This study aimed to assess the mechanical behavior of three different types of IAC using a finite element model (FEM) analysis. Three types of two-piece implants were designed: two internal conical connection designs (models A and B) and one internal flat-to-flat connection design (model C). This three-dimensional analysis evaluated the response to static forces on the three models. The strain map, stress analysis, and safety factor were assessed by means of the FEM examination. The FEM analysis indicated that forces are transmitted on the abutment and implant's neck in model B. In models A and C, forces were distributed along the internal screw, abutment areas, and implant's neck. The stress distribution in model B showed a more homogeneous pattern, such that the peak forces were reduced. The conical shape of the head of the internal screw in model B seems to have a keystone role in transferring the forces at the surrounding structures. Further experiments should be carried out in order to confirm the present suppositions.
| S-EPMC6803952 | biostudies-literature
Project description:Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace's equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions.
| S-EPMC5100939 | biostudies-literature
Project description:Silicon nanomembranes are ultrathin, highly permeable, optically transparent and biocompatible substrates for the construction of barrier tissue models. Trans-epithelial/endothelial electrical resistance (TEER) is often used as a non-invasive, sensitive and quantitative technique to assess barrier function. The current study characterizes the electrical behavior of devices featuring silicon nanomembranes to facilitate their application in TEER studies. In conventional practice with commercial systems, raw resistance values are multiplied by the area of the membrane supporting cell growth to normalize TEER measurements. We demonstrate that under most circumstances, this multiplication does not 'normalize' TEER values as is assumed, and that the assumption is worse if applied to nanomembrane chips with a limited active area. To compare the TEER values from nanomembrane devices to those obtained from conventional polymer track-etched (TE) membranes, we develop finite element models (FEM) of the electrical behavior of the two membrane systems. Using FEM and parallel cell-culture experiments on both types of membranes, we successfully model the evolution of resistance values during the growth of endothelial monolayers. Further, by exploring the relationship between the models we develop a 'correction' function, which when applied to nanomembrane TEER, maps to experiments on conventional TE membranes. In summary, our work advances the the utility of silicon nanomembranes as substrates for barrier tissue models by developing an interpretation of TEER values compatible with conventional systems.
| S-EPMC5770343 | biostudies-literature
Project description:We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L(?) estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme are demonstrated with FETK through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.
| S-EPMC3178275 | biostudies-literature
Project description:Endovascular interventions using coil embolization and flow diversion are becoming the mainstream treatment for intracranial aneurysms (IAs). To help assess the effect of intervention strategies on aneurysm hemodynamics and treatment outcome, we have developed a finite-element-method (FEM)-based technique for coil deployment along with our HiFiVS technique for flow diverter (FD) deployment in patient-specific IAs. We tested four clinical intervention strategies: coiling (1-8 coils), single FD, FD with adjunctive coils (1-8 coils), and overlapping FDs. By evaluating post-treatment hemodynamics using computational fluid dynamics (CFD), we compared the flow-modification performance of these strategies. Results show that a single FD provides more reduction in inflow rate than low packing density (PD) coiling, but less reduction in average velocity inside the aneurysm. Adjunctive coils add no additional reduction of inflow rate beyond a single FD until coil PD exceeds 11%. This suggests that the main role of FDs is to divert inflow, while that of coils is to create stasis in the aneurysm. Overlapping FDs decreases inflow rate, average velocity, and average wall shear stress (WSS) in the aneurysm sac, but adding a third FD produces minimal additional reduction. In conclusion, our FEM-based techniques for virtual coiling and flow diversion enable recapitulation of complex endovascular intervention strategies and detailed hemodynamics to identify hemodynamic factors that affect treatment outcome.
| S-EPMC4801175 | biostudies-literature
Project description:High-order cubic Hermite finite elements have been valuable in modeling cardiac geometry, fiber orientations, biomechanics, and electrophysiology, but their use in solving three-dimensional problems has been limited to ventricular models with simple topologies. Here, we utilized a subdivision surface scheme and derived a generalization of the "local-to-global" derivative mapping scheme of cubic Hermite finite elements to construct bicubic and tricubic Hermite models of the human atria with extraordinary vertices from computed tomography images of a patient with atrial fibrillation. To an accuracy of 0.6 mm, we were able to capture the left atrial geometry with only 142 bicubic Hermite finite elements, and the right atrial geometry with only 90. The left and right atrial bicubic Hermite meshes were G1 continuous everywhere except in the one-neighborhood of extraordinary vertices, where the mean dot products of normals at adjacent elements were 0.928 and 0.925. We also constructed two biatrial tricubic Hermite models and defined fiber orientation fields in agreement with diagrammatic data from the literature using only 42 angle parameters. The meshes all have good quality metrics, uniform element sizes, and elements with aspect ratios near unity, and are shared with the public. These new methods will allow for more compact and efficient patient-specific models of human atrial and whole heart physiology.
| S-EPMC3660421 | biostudies-literature
Project description:Photoemission electron microscopy and imaging X-ray photoelectron spectroscopy are today frequently used to obtain chemical and electronic states, chemical shifts, work function profiles within the fields of surface- and material sciences. Lately, because of recent technological advances, these tools have also been valuable within life sciences. In this study, we have investigated the power of photoemission electron microscopy and imaging X-ray photoelectron spectroscopy for visualization of human neutrophil granulocytes. These cells, commonly called neutrophils, are essential for our innate immune system. We hereby investigate the structure and morphology of neutrophils when adhered to gold and silicon surfaces. Energy-filtered imaging of single cells are acquired. The characteristic polymorphonuclear cellular nuclei divided into 2-5 lobes is visualized. Element-specific imaging is achieved based on O 1s, P 2p, C 1s, Si 2p, and N 1s core level spectra, delivering elemental distribution with submicrometer resolution, illustrating the strength of this type of cellular morphological studies.
| S-EPMC7809688 | biostudies-literature
Project description:Finite element analysis has been an increasingly widely applied biomechanical modeling method in many different science and engineering fields over the last decade. In the biological sciences, there are many examples of FEA in areas such as paleontology and functional morphology. Despite this common use, the modeling of trabecular bone remains a key issue because their highly complex and porous geometries are difficult to replicate in the solid mesh format required for many simulations. A common practice is to assign uniform model material properties to whole or portions of models that represent trabecular bone. In this study we aimed to demonstrate that a physical, element reduction approach constitutes a valid protocol for addressing this problem in addition to the wholesale mathematical approach. We tested a customized script for element reduction modeling on five exemplar trabecular geometry models of carnivoran temporomandibular joints, and compared stress and strain energy results of both physical and mathematical trabecular modeling to models incorporating actual trabecular geometry. Simulation results indicate that that the physical, element reduction approach generally outperformed the mathematical approach: physical changes in the internal structure of experimental cylindrical models had a major influence on the recorded stress values throughout the model, and more closely approximates values obtained in models containing actual trabecular geometry than solid models with modified trabecular material properties. In models with both physical and mathematical adjustments for bone porosity, the physical changes exhibit more weight than material properties changes in approximating values of control models. Therefore, we conclude that maintaining or mimicking the internal porosity of a trabecular structure is a more effective method of approximating trabecular bone behavior in finite element models than modifying material properties.
| S-EPMC6925947 | biostudies-literature
Project description:Concussion is a significant public health problem affecting 1.6-2.4 million Americans annually. An alternative to reducing the burden of concussion is to reduce its incidence with improved protective equipment and injury mitigation systems. Finite element (FE) models of the brain response to blunt trauma are often used to estimate injury potential and can lead to improved helmet designs. However, these models have yet to incorporate how the patterns of brain connectivity disruption after impact affects the relay of information in the injured brain. Furthermore, FE brain models typically do not consider the differences in individual brain structural connectivities and their purported role in concussion risk. Here, we use graph theory techniques to integrate brain deformations predicted from FE modeling with measurements of network efficiency to identify brain regions whose connectivity characteristics may influence concussion risk. We computed maximum principal strain in 129 brain regions using head kinematics measured from 53 professional football impact reconstructions that included concussive and non-concussive cases. In parallel, using diffusion spectrum imaging data from 30 healthy subjects, we simulated structural lesioning of each of the same 129 brain regions. We simulated lesioning by removing each region one at a time along with all its connections. In turn, we computed the resultant change in global efficiency to identify regions important for network communication. We found that brain regions that deformed the most during an impact did not overlap with regions most important for network communication (Pearson's correlation, ? = 0.07; p = 0.45). Despite this dissimilarity, we found that predicting concussion incidence was equally accurate when considering either areas of high strain or of high importance to global efficiency. Interestingly, accuracy for concussion prediction varied considerably across the 30 healthy connectomes. These results suggest that individual network structure is an important confounding variable in concussion prediction and that further investigation of its role may improve concussion prediction and lead to the development of more effective protective equipment.
| S-EPMC7174699 | biostudies-literature