# Revision of Solving incompressible finite elasticity without tears from Thu, 2019-02-07 15:57

Solving incompressible elasticity has been quite challenging numerically. The conventional approach for handling incompressibility is the so-called penalty method. One straightforward problem is that the penalty parameter goes directly into the tangent matrix. The bigger the penalty parameter, the worse the condition number of the matrix. When the problem size is not big, direct solvers can be used and we probably do not realize the ill-conditioning. This is not the case when the problem size gets bigger and gets coupled in the multiphysics setting (e.g. fluid-structure interaction). This drives us to think of things in a different way. Mathematically, the penalty method can be transformed into a saddle-point problem, just like what has been done in fluid mechanics. In thermodynamics, this is achieved through the Legendre transformation between potentials. For a saddle-point problem, there are many ways of solving the equations [1]. The trick is to devise a good approximation for the Schur complement. In a recent JCP paper [2], we developed an effective way of solving the finite-strain elasticity, using an iterative algorithm for sparse matrices. The algorithm can drive the residual to machine tolerance within a few iterations; it works well for compressible and incompressible materials; it is scalable on parallel machines. For the thermomechanical foundation of this framework, readers are also referred to [3].

[1] Michele Benzi, Gene H. Golub, and Jörg Liesen. "Numerical solution of saddle point problems." Acta Numerica 14 (2005): 1-137.

[2] Ju Liu and Alison L. Marsden. "A robust and efficient iterative method for hyper-elastodynamics with nested blockpreconditioning." Journal of Computational Physics 383 (2019): 72-93.

https://authors.elsevier.com/a/1YXBH508HiFTO

https://www.sciencedirect.com/science/article/pii/S0021999119300440

[3] Ju Liu and Alison L. Marsden. "A unified continuum and variational multiscale formulation for fluids, solids, and fluid-structure interaction." Computer Methods in Applied Mechanics and Engineering 337 (2018): 549-597.

https://www.sciencedirect.com/science/article/pii/S0045782518301701