Dynamic finite-strain modelling of the human left ventricle in health and disease using an immersed boundary-finite element method

Detailed models of the biomechanics of the heart are important both for developing improved interventions for patients with heart disease and also for patient risk stratification and treatment planning. For instance, stress distributions in the heart affect cardiac remodelling, but such distributions are not presently accessible in patients. Biomechanical models of the heart offer detailed three-dimensional deformation, stress and strain fields that can supplement conventional clinical data. In this work, we introduce dynamic computational models of the human left ventricle (LV) that are derived from clinical imaging data obtained from a healthy subject and from a patient with a myocardial infarction (MI). Both models incorporate a detailed invariant-based orthotropic description of the passive elasticity of the ventricular myocardium along with a detailed biophysical model of active tension generation in the ventricular muscle. These constitutive models are employed within a dynamic simulation framework that accounts for the inertia of the ventricular muscle and the blood that is based on an immersed boundary (IB) method with a finite element description of the structural mechanics. The geometry of the models is based on data obtained non-invasively by cardiac magnetic resonance (CMR). CMR imaging data are also used to estimate the parameters of the passive and active constitutive models, which are determined so that the simulated end-diastolic and end-systolic volumes agree with the corresponding volumes determined from the CMR imaging studies. Using these models, we simulate LV dynamics from end-diastole to end-systole. The results of our simulations are shown to be in good agreement with subject-specific CMR-derived strain measurements and also with earlier clinical studies on human LV strain distributions

Particle Acceleration in Pulsar Wind Nebulae: PIC modelling

We discuss the role of particle-in-cell (PIC) simulations in unveiling the origin of the emitting particles in PWNe. After describing the basics of the PIC technique, we summarize its implications for the quiescent and the flaring emission of the Crab Nebula, as a prototype of PWNe. A consensus seems to be emerging that, in addition to the standard scenario of particle acceleration via the Fermi process at the termination shock of the pulsar wind, magnetic reconnection in the wind, at the termination shock and in the Nebula plays a major role in powering the multi-wavelength signatures of PWNe.Comment: 32 pages, 16 figures, to appear in the book "Modelling Nebulae" edited by D. Torres for Springer, based on the invited contributions to the workshop held in Sant Cugat (Barcelona), June 14-17, 201

Efficient fe strategies for springback prediction – material modelling and computational aspects

Blanks of sheet metal are characterized by an intrinsic plastic anisotropic behaviour resulting from the plastic deformation during the rolling of sheets. Another type of anisotropy is elastic anisotropy which might be essential especially during elastic recovery processes during unloading after forming and springback. Thus, this paper focuses on the study of the sensitivity of the amount of springback in unconstrained bending with respect to elastic anisotropy. A finite strain constitutive model for evolving elastic and plastic anisotropy combining nonlinear isotropic and kinematic hardening is discussed. The evolution of elastic anisotropy is described by representing the Helmholtz free energy as a function of a family of evolving structure tensors. In addition, plastic anisotropy is modelled via the dependence of the yield surface on the same family of structure tensors. The constitutive equations of the model are implemented as a user material subroutine UMAT in the commercial solver ABAQUS/Standard, which is then applied to the simulation of springback in unconstrained bending

Unsteady wake modelling for tidal current turbines

The authors present a numerical model for three-dimensional unsteady wake calculations for tidal turbines. Since wakes are characterised by the shedding of a vortex sheet from the rotor blades, the model is based on the vorticity transport equations. A vortex sheet may be considered a jump contact discontinuity in tangential velocity with, in inviscid hydrodynamic terms, certain kinematic and dynamic conditions across the sheet. The kinematic condition is that the sheet is a stream surface with zero normal fluid velocity; the dynamic condition is that the pressure is equal on either side of the sheet. The dynamic condition is explicitly satisfied at the trailing edge only, via an approximation of the Kutta condition. The shed vorticity is the span-wise derivative of bound circulation, and the trailed vorticity is the time derivative of bound circulation, and is convected downstream from the rotors using a finite-volume solution of vorticity transport equations thus satisfying the kinematic conditions. Owing to an absence in the literature of pressure data for marine turbines, results from the code are presented for the NREL-UAE Phase IV turbine. Axial flow cases show a close match in pressure coefficients at various spanwise stations; however, yawed flow cases demonstrate the shortcomings of a modelling strategy lacking viscosity

Effect of thermal expansion on the linear stability of planar premixed flames for a simple chain-branching model: The high activation energy asymptotic limit

The linear stability of freely propagating, adiabatic, planar premixed ames is investigated in the context of a simple chain-branching chemistry model consisting of a chain-branching reaction step and a completion reaction step. The role of chain-branching is governed by a crossover temperature. Hydrodynamic effects, induced by thermal expansion, are taken into account and the results compared and contrasted with those from a previous purely thermal-di�usive constant density linear stability study. It is shown that when thermal expansion is properly accounted for, a region of stable ames predicted by the constant density model disappears, and instead the ame is unstable to a long-wavelength cellular instability. For a pulsating mode, however, thermal expansion is shown to have only a weak e�ect on the critical fuel Lewis number required for instability. These e�ects of thermal expansion on the two-step chain-branching ame are shown to be qualitatively similar to those on the standard one-step reaction model. Indeed, as found by constant density studies, in the limit that the chain-branching crossover temperature tends to the adiabatic ame temperature, the two-step model can be described to leading order by the one-step model with a suitably de�ned e�ective activation energy

Quasi-static imaged-based immersed boundary-finite element model of human left ventricle in diastole

SUMMARY: Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid–structure interaction simulations of the heart

Multi-D magnetohydrodynamic modelling of pulsar wind nebulae: recent progress and open questions

In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting particles and the quest for the acceleration sites of particles of different energies along the termination shock, by using wisps motion as a diagnostic tool; the study of the magnetic dissipation process in high magnetization nebulae by means of new long-term three-dimensional simulations of the pulsar wind nebula evolution; the investigation of the relativistic tearing instability in thinning current sheets, leading to fast reconnection events that might be at the origin of the Crab nebula gamma-ray flares.Comment: 30 pages, 12 figure