After the changes in #1102 where the architecture of the Newton solver was changed to be able to reuse the preconditioner, the initial conditions related to viscosity were not working as expected in the matrix-free application when using multigrid preconditioners.
Description of the solution
Now the kinematic viscosity of all the operators of the different multigrid levels is set to the temporary viscosity when needed to ensure the correct preconditioning.
How Has This Been Tested?
The flow around a sphere benchmark, which uses the viscosity ramp as initial condition, obtains now the expected results.
Description of the problem
After the changes in #1102 where the architecture of the Newton solver was changed to be able to reuse the preconditioner, the initial conditions related to viscosity were not working as expected in the matrix-free application when using multigrid preconditioners.
Description of the solution
Now the kinematic viscosity of all the operators of the different multigrid levels is set to the temporary viscosity when needed to ensure the correct preconditioning.
How Has This Been Tested?
The flow around a sphere benchmark, which uses the viscosity ramp as initial condition, obtains now the expected results.