An elliptically shaped mesh in rectangular Cartesian coordinates is inflated.

Instead of enforceing incompressibility through the hydro-static pressure constraint a penalty function is added in the pole-zero cellML code: J^2-2J+1. This results in an extra penalty term for each stress-component.

This example is the same problem as example 5g3, except it is using this nearly incompressible approach and not the incompressibility choice in the ipequa file. The output should be relatively similar to 5g3.

Created by Espen Remme May 2004.