We present a numerical study of the sensitivity of the grad-div stabilization parameter for mixed finite element discretizations of incompressible flow problems. For incompressible isothermal and nonisothermal Stokes equations and Navier-Stokes equations, we develop the associated sensitivity equations for changes in the grad-div parameter. Finite element schemes are devised for computing solutions to the sensitivity systems, analyzed for stability and accuracy, and finally tested on several benchmark problems. Our results reveal that solutions are most sensitive for small values of the parameter, near obstacles and corners, when the pressure is large, and when the viscosity is small.
All Science Journal Classification (ASJC) codes
- Computational Mathematics