A Parallel High-Order Fictitious Domain Approach for Biomechanical Applications

The focus of this contribution is on the parallelization of the Finite Cell Method (FCM) applied for biomechanical simulations of human femur bones. The FCM is a high-order fictitious domain method that combines the simplicity of Cartesian grids with the beneficial properties of hierarchical approximation bases of higher order for an increased accuracy and reliablility of the simulation model. A pre-computation scheme for the numerically expensive parts of the finite cell model is presented that shifts a significant part of the analysis update to a setup phase of the simulation, thus increasing the update rate of linear analyses with time-varying geometry properties to a range that even allows user interactive simulations of high quality. Paralellization of both parts, the pre-computation of the model stiffness and the update phase of the simulation is simplified due to a simple and undeformed cell structure of the computation domain. A shared memory parallelized implementation of the method is presented and its performance is tested for a biomedical application of clinical relevance to demonstrate the applicability of the presented method

Modelling mitral valvular dynamics–current trend and future directions

Dysfunction of mitral valve causes morbidity and premature mortality and remains a leading medical problem worldwide. Computational modelling aims to understand the biomechanics of human mitral valve and could lead to the development of new treatment, prevention and diagnosis of mitral valve diseases. Compared with the aortic valve, the mitral valve has been much less studied owing to its highly complex structure and strong interaction with the blood flow and the ventricles. However, the interest in mitral valve modelling is growing, and the sophistication level is increasing with the advanced development of computational technology and imaging tools. This review summarises the state-of-the-art modelling of the mitral valve, including static and dynamics models, models with fluid-structure interaction, and models with the left ventricle interaction. Challenges and future directions are also discussed

Multi-level hp-finite cell method for embedded interface problems with application in biomechanics

This work presents a numerical discretization technique for solving three-dimensional material interface problems involving complex geometry without conforming mesh generation. The finite cell method (FCM), which is a high-order fictitious domain approach, is used for the numerical approximation of the solution without a boundary-conforming mesh. Weak discontinuities at material interfaces are resolved by using separate FCM meshes for each material sub-domain, and weakly enforcing the interface conditions between the different meshes. Additionally, a recently developed hierarchical hp-refinement scheme is employed to locally refine the FCM meshes in order to resolve singularities and local solution features at the interfaces. Thereby, higher convergence rates are achievable for non-smooth problems. A series of numerical experiments with two- and three-dimensional benchmark problems is presented, showing that the proposed hp-refinement scheme in conjunction with the weak enforcement of the interface conditions leads to a significant improvement of the convergence rates, even in the presence of singularities. Finally, the proposed technique is applied to simulate a vertebra-implant model. The application showcases the method's potential as an accurate simulation tool for biomechanical problems involving complex geometry, and it demonstrates its flexibility in dealing with different types of geometric description

SOLID-SHELL FINITE ELEMENT MODELS FOR EXPLICIT SIMULATIONS OF CRACK PROPAGATION IN THIN STRUCTURES

Crack propagation in thin shell structures due to cutting is conveniently simulated using explicit finite element approaches, in view of the high nonlinearity of the problem. Solidshell elements are usually preferred for the discretization in the presence of complex material behavior and degradation phenomena such as delamination, since they allow for a correct representation of the thickness geometry. However, in solid-shell elements the small thickness leads to a very high maximum eigenfrequency, which imply very small stable time-steps. A new selective mass scaling technique is proposed to increase the time-step size without affecting accuracy. New ”directional” cohesive interface elements are used in conjunction with selective mass scaling to account for the interaction with a sharp blade in cutting processes of thin ductile shells

A finite strain nonlinear human mitral valve model with fluid structure interaction

A simulated human mitral valve under a physiological pressure loading is developed using a hybrid finite element immersed boundary method, which incorporates experimentally based constitutive laws in a three-dimensional fluid-structure interaction framework. A transversely isotropic material constitutive model is used for characterizing the mechanical behaviour of the mitral valve tissue based on recent mechanical tests of healthy human mitral leaflets. Our results show good agreement, in terms of the flow rate and the closing and opening configurations, with the measurements from the magnetic resonance images. The stresses in the anterior leaflet are found to be higher than those in the posterior leaflet, and concentrated around the annulus trigons and free edges of the valve leaflets. Those areas are located where the leaflet has the highest curvature. Effects of the chordae tendineae in the material model are studied and the results show that these chordae play an important role in providing a secondary orifice for the flow when valve opens. Although there are some discrepancies to be overcome in future works, our simulations show that the developed computational model is promising in mimicking the in vivo mitral valve dynamics and providing important information that are not obtainable by in vivo measurements. This article is protected by copyright. All rights reserved

A voxelized immersed boundary (VIB) finite element method for accurate and efficient blood flow simulation

We present an efficient and accurate immersed boundary (IB) finite element (FE) method for internal flow problems with complex geometries (e.g., blood flow in the vascular system). In this study, we use a voxelized flow domain (discretized with hexahedral and tetrahedral elements) instead of a box domain, which is frequently used in IB methods. The proposed method utilizes the well-established incremental pressure correction scheme (IPCS) FE solver, and the boundary condition-enforced IB (BCE-IB) method to numerically solve the transient, incompressible Navier--Stokes flow equations. We verify the accuracy of our numerical method using the analytical solution for the Poiseuille flow in a cylinder, and the available experimental data (laser Doppler velocimetry) for the flow in a three-dimensional 90{\deg} angle tube bend. We further examine the accuracy and applicability of the proposed method by considering flow within complex geometries, such as blood flow in aneurysmal vessels and the aorta, flow configurations that would otherwise be difficult to solve by most IB methods. Our method offers high accuracy, as demonstrated by the verification examples, and high applicability, as demonstrated through the solution of blood flow within complex geometry. The proposed method is efficient, since it is as fast as the traditional finite element method used to solve the Navier--Stokes flow equations, with a small overhead (not more than 5$\%$) due to the numerical solution of a linear system formulated for the IB method.Comment: arXiv admin note: substantial text overlap with arXiv:2007.0208