Search

Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid-only to coupled fluid-structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

Models for cardiac mechanics require an activation mechanism properly representing the stress-strain relations in the contracting myocardium. In this paper, we propose a new activation model that accounts for the transmural heterogeneities observed in myocardial strain measurements. In order to take the anisotropy of the active mechanics into account, our model is based on an active strain formulation. Thanks to multiplicative decomposition of the deformation gradient tensor, in this formulation, the active strains orthogonal to the fibers can be naturally described. We compare the results of our novel formulation against different anisotropic models of the active contraction of the cardiac muscle, as well as against experimental data available in the literature. We show that with the currently available models, the strain distributions are not in agreement with the reported experimental measurements. Conversely, we show that our new transmurally heterogeneous orthotropic activation model improves the accuracy of shear strains related to in-plane rotations and torsion.

Biological tissues often experience drastic changes in their microstructure due to their pathophysiological conditions. Such microstructural changes could result in variations in mechanical properties, which can be used in diagnosing or monitoring a wide range of diseases, most notably cancer. This paves the avenue for non‐invasive diagnosis by instrumented palpation although challenges remain in quantitatively assessing the amount of diseased tissue by means of mechanical characterization. This paper presents a framework for tissue diagnosis using a quantitative and efficient estimation of the fractions of cancerous and non‐cancerous tissue without a priori knowledge of tissue microstructure. First, the sample is tested in a creep or stress relaxation experiment, and the behavior is characterized using a single term Prony series. A rule of mixtures, which relates tumor fraction to the apparent mechanical properties, is then obtained by minimizing the difference between strain energy of a heterogeneous system and an equivalent homogeneous one. Finally, the percentage of each tissue constituent is predicted by comparing the observed relaxation time with that calculated from the rule of mixtures. The proposed methodology is assessed using models reconstructed from histological samples and magnetic resonance imaging of prostate. Results show that estimation of cancerous tissue fraction can be obtained with a maximum error of 12% when samples of different sizes, geometries, and tumor fractions are presented. The proposed framework has the potential to be applied to a wide range of diseases such as rectal polyps, cirrhosis, or breast and prostate cancer whose current primary diagnosis remains qualitative.

The main objective of this paper is to propose a novel method that provides an opportunity to evaluate an orthodontic process at early phase of the treatment. This was accomplished by finding out a correlation between the applied orthodontic force and thermal variations in the tooth structure. To this end, geometry of the human tooth surrounded by the connective soft tissue called the periodontal ligament and the bone was constructed by employing dental CT scan images of a specific case. The periodontal ligament was modeled by finite strain viscoelastic model through a nonlinear stress-strain relation (hyperelasticity) and nonlinear stress-time relation (viscoelasticity). The tooth structure was loaded by a lateral force with 15 different quantities applied to 20 different locations, along the midedge of the tooth crown. The resultant compressive stress in the periodontal ligament was considered as the cause of elevated cell activity that was modeled by a transient heat flux in the thermal analysis. The heat flux value was estimated by conducting an experiment on a pair of rats. The numerical results showed that by applying an orthodontic force to the tooth structure, a significant temperature rise was observed. By measuring the temperature rise, the orthodontic process can be evaluated.

The membrane or interfacial viscosity is an important property in many multiphase and biofluidic situations, such as the red blood cell dynamics and emulsion stability. The immersed boundary method (IBM), which incorporates the dynamic flow-membrane interaction via force distribution and velocity interpolation, has been extensively employed in simulations of such systems. Unfortunately, direct implementation of membrane viscosity in IBM suffers severe numerical instability, which causes an IBM calculation to break down before generating any useful results. Few attempts have been recently reported; however, several concerns exist in these attempts, such as the inconsistency to the classical definition of membrane viscosity, the inability to model the shear and dilatational viscosities separately, the unjustified mathematical formulations, and the complicated algorithms and computation. To overcome these concerns, in this paper, we propose a finite difference approach for implementing membrane viscosity in immersed boundary simulations. The viscous stress is obtained via finite difference approximations to the differential strain-stress relationship, with the help of a subsampling scheme to reduce the numerical noise in the calculated strain rates. This simple method has also avoided the complicated matrix calculations in previous attempts, and hence, a better computational efficiency is expected. Detailed mathematical description of the method and key steps for its implementation in immersed boundary programs are provided. Validation and illustration calculations are performed, and our results are compared with analytical solutions and previous publications with satisfactory agreement. The influences of membrane mesh resolution and simulation time step are also examined; and the results show no indication that our finite difference method has downgraded the general IBM accuracy. Based on these simulations and analysis, we believe that our method would be a better choice for future IBM simulations of capsule dynamics with viscoelastic membranes.

This paper presents a numerical study in which several partitioned solution procedures for incompressible fluid-structure interaction are compared and validated against the results of an experimental fluid-structure interaction benchmark. The numerical methods discussed cover the three main families of coupling schemes: strongly coupled, semi-implicit, and loosely coupled. Very good agreement is observed between the numerical and experimental results. The comparisons confirm that strong coupling can be efficiently avoided, via semi-implicit and loosely coupled schemes, without compromising stability and accuracy. Copyright © 2016 John WileySons, Ltd.

In this paper, we present a novel rational macromodeling approach for the description of 1D blood flow in the human cardiovascular system, which is suitable for time-domain simulations. Using the analogy of the blood flow propagation problem with transmission lines and considering the hypothesis of linearized Navier-Stokes equations, a frequency-domain rational macromodel for each arterial segment has been built. The poles and the residues of each arterial segment macromodel have been calculated by means of the Vector Fitting technique. Finally, the rational macromodel of the whole cardiovascular system is obtained by properly combining the macromodels of the single arterial segments using an interconnect matrix. The rational form of the proposed cardiovascular model leads to a state-space or electrical circuit model suitable for time-domain analysis. The stability and passivity properties of the global cardiovascular model are discussed to guarantee stable time-domain simulations. The proposed macromodeling approach has been validated by pertinent numerical results. Copyright © 2015 John WileySons, Ltd.

A coupled computational model of the foot consisting of a three-dimensional soft tissue continuum and a one-dimensional (1D) transient blood flow network is presented in this article. The primary aim of the model is to investigate the blood flow in major arteries of the pathologic foot where the soft tissue stiffening occurs. It has been reported in the literature that there could be up to about five-fold increase in the mechanical stiffness of the plantar soft tissues in pathologic (e.g. diabetic) feet compared with healthy ones. The increased stiffness results in higher tissue hydrostatic pressure within the plantar area of the foot when loaded. The hydrostatic pressure acts on the external surface of blood vessels and tend to reduce the flow cross-section area and hence the blood supply. The soft tissue continuum model of the foot was modelled as a tricubic Hermite finite element mesh representing all the muscles, skin and fat of the foot and treated as incompressible with transversely isotropic properties. The details of the mechanical model of soft tissue are presented in the companion paper, Part 1. The deformed state of the soft tissue continuum because of the applied ground reaction force at three foot positions (heel-strike, midstance and toe-off) was obtained by solving the Cauchy equations based on the theory of finite elasticity using the Galerkin finite element method. The geometry of the main arterial network in the foot was represented using a 1D Hermite cubic finite element mesh. The flow model consists of 1D Navier-Stokes equations and a nonlinear constitutive equation to describe vessel radius-transmural pressure relation. The latter was defined as the difference between the fluid and soft tissue hydrostatic pressure. Transient flow governing equations were numerically solved using the two-step Lax-Wendroff finite difference method. The geometry of both the soft tissue continuum and arterial network is anatomically-based and was developed using the data derived from visible human images and magnetic resonance images of a healthy male volunteer. Simulation results reveal that a two-fold increase in tissue stiffness leads to about 28% reduction in blood flow to the affected region.

This paper deals with the modeling of guided waves propagation in in vivo cortical long bone, which is known to be anisotropic medium with functionally graded porosity. The bone is modeled as an anisotropic poroelastic material by using Biot's theory formulated in high frequency domain. A hybrid spectral/finite element formulation has been developed to find the time-domain solution of ultrasonic waves propagating in a poroelastic plate immersed in two fluid halfspaces. The numerical technique is based on a combined Laplace-Fourier transform, which allows to obtain a reduced dimension problem in the frequency-wavenumber domain. In the spectral domain, as radiation conditions representing infinite fluid halfspaces may be exactly introduced, only the heterogeneous solid layer needs to be analyzed by using finite element method. Several numerical tests are presented showing very good performance of the proposed procedure. A preliminary study on the first arrived signal velocities computed by using equivalent elastic and poroelastic models will be presented.