I am trying to solve an almost incompressible problem with heterogeneous properties by domain decomposition. Solution with CG converges slowly or divergerces completely. My problem becomes ill-conditioned (condition number in the order of 10^5) due to the poissons ratio approaching to 0.5 and stiffness differences in between materials. I detected that in element stiffness matrix, several additional zero energy modes (hourglass modes etc) are formed in addition to standard rigid body modes.
I am planning to prevent these modes by MatNullSpace object in Petsc. Since the computation of nullspaces with SVD or QR is expensive, I am trying to find a cheaper way of predicting the nullspace of floating subdomains.
Is it possible to predict nullspaces for subdomain from contributing finite element's null spaces?
I read a broad sentence in a publication says that "The rigid body modes of a collection of elements is equal to the assembly of the rigid body modes of the individual elements taking into account the multiplicity of those degrees of freedom that lie in multiple neighboring elements."
does that really make sense?