[petsc-users] multi-grid for anisotropy

[Justin Chang]

Hi all,

I think I have asked a previous question related to this, and have seen other threads about related issues, but say I have this anisotropic dispersion tensor:

D = (D_m + alpha_L*norm(v))*Identity + (alpha_L - alpha_T)*outer(v,v)/norm(v)

The anisotropy of this tensor comes from the ratio of alpha_L/alpha_T and velocity v, and is mesh independent. That said, can tuning the multi-grid preconditioners parameters of GAMG/HYPRE/etc help prevent things like negative concentrations and violations of the discrete maximum principles?

Thanks,