Parametric uncertainty analysis of pulse wave propagation in a model of a human arterial network

Accepted versio

A new framework for large strain electromechanics based on convex multi-variable strain energies: Finite Element discretisation and computational implementation

In Gil and Ortigosa (2016), Gil and Ortigosa introduced a new convex multi-variable framework for the numerical simulation of Electro Active Polymers (EAPs) in the presence of extreme deformations and electric fields. This extends the concept of polyconvexity to strain energies which depend on non-strain based variables. The consideration of the new concept of multi-variable convexity guarantees the well posedness of generalised Gibbs’ energy density functionals and, hence, opens up the possibility of a new family of mixed variational principles. The aim of this paper is to present, as an example, the Finite Element implementation of two of these mixed variational principles. These types of enhanced methodologies are known to be necessary in scenarios in which the simpler displacement-potential based formulation yields non-physical results, such as volumetric locking, bending and shear locking, pressure oscillations and electro-mechanical locking, to name but a few. Crucially, the use of interpolation spaces in which some of the unknown fields are described as piecewise discontinuous across elements can be used in order to efficiently condense these fields out. This results in mixed formulations with a computational cost comparable to that of the displacement-potential based approach, yet far more accurate. Finally, a series of very challenging numerical examples are presented in order to demonstrate the accuracy, robustness and efficiency of the proposed methodology

Sparse Pseudospectral Approximation Method

Multivariate global polynomial approximations - such as polynomial chaos or stochastic collocation methods - are now in widespread use for sensitivity analysis and uncertainty quantification. The pseudospectral variety of these methods uses a numerical integration rule to approximate the Fourier-type coefficients of a truncated expansion in orthogonal polynomials. For problems in more than two or three dimensions, a sparse grid numerical integration rule offers accuracy with a smaller node set compared to tensor product approximation. However, when using a sparse rule to approximately integrate these coefficients, one often finds unacceptable errors in the coefficients associated with higher degree polynomials. By reexamining Smolyak's algorithm and exploiting the connections between interpolation and projection in tensor product spaces, we construct a sparse pseudospectral approximation method that accurately reproduces the coefficients of basis functions that naturally correspond to the sparse grid integration rule. The compelling numerical results show that this is the proper way to use sparse grid integration rules for pseudospectral approximation

A computational framework for large strain nearly and truly incompressible electromechanics based on convex multi-variable strain energies

The series of papers published by Gil and Ortigosa (Gil and Ortigosa, 2016; Ortigosa and Gil, 2016, 0000) introduced a new convex multi-variable variational and computational framework for the numerical simulation of Electro Active Polymers (EAPs) in scenarios characterised by extreme deformations and/or extreme electric fields. Building upon this body of work, five key novelties are incorporated in this paper. First, a generalisation of the concept of multi-variable convexity to energy functionals additively decomposed into isochoric and volumetric components. This decomposition is typical of nearly and truly incompressible materials, group which represents the majority of the most relevant EAPs. Second, convexification or regularisation strategies are applied to a priori non-convex multi-variable isochoric functionals to yield physically meaningful convex multi-variable functionals. Third, based on the mixed variational principles introduced in Gil and Ortigosa (2016) in the context of compressible electro-elasticity, a novel extended Hu–Washizu mixed variational principle for nearly and truly incompressible scenarios is presented. From the computational standpoint, a static condensation procedure is applied in order to condense out the element-wise extra fields, the resulting formulation having a comparable cost to the more standard three-field displacement-potential-pressure mixed formulation. Fourth, the computational framework for the three-field mixed variational principle in nearly and truly incompressible scenarios is also presented. In this case, the novelty resides in the consideration of convex multi-variable energy functionals. Ultimately, this leads to the definition of new tangent operators for the Helmholtz’s energy functional in the specific context of incompressible electro-elasticity. Fifth, a Petrov–Galerkin stabilisation technique is applied on the three-field formulation for the circumvention of the Ladyz˘enskaja–Babus˘ka–Brezzi (LBB) condition, enabling the use of linear tetrahedral finite elements for the interpolation of the unknowns of the problem. Finally, a series of challenging numerical examples is presented in order to provide an exhaustive comparison of the different variational formulations presented in this paper

Geometrical validity of curvilinear finite elements

In this paper, we describe a way to compute accurate bounds on Jacobian determinants of curvilinear polynomial finite elements. Our condition enables to guarantee that an element is geometrically valid, i.e., that its Jacobian determinant is strictly positive everywhere in its reference domain. It also provides an efficient way to measure the distortion of curvilinear elements. The key feature of the method is to expand the Jacobian determinant using a polynomial basis, built using Bézier functions, that has both properties of boundedness and positivity. Numerical results show the sharpness of our estimates