I wish to solve the two-dimensional Navier Stokes equations using multigrid method on a collocated grid using the Predictor-Corrector method mentioned below. But first, let me elaborate on what I had already did. Previously, I used a staggered grid and a artificial compressibility based solver. So my governing equations were:
Since my solver is explicit (the time term here is just to do the iterations). To implement multigrid in this kind of methodology, I defined the defect/residual on fine grid as:
This defect/residual and also the solution on the fine grid is then restricted to the coarse grid (Full Storage Approximation Scheme) using simple injection operator. Since the equations on the coarse grid remain almost the same in appearance, except that their RHS is not zero and instead is defined as:
where the dashed terms are m and not the corrections itself. So, after I solve the coarse mesh equations to find the correction, I interpolate them to fine mesh to correct the fine mesh solution.
This method works just fine.
Now, I want to implement multigrid on a colocated grid with the following predictor-corrector based discretization:
where the three steps are predictor, corrector and update steps respectively. Predictor is essentially a forward difference based on the values at old time level while corrector uses a backward difference using the values at predictor time level.
How will I define the defect/residual on the fine mesh for this scheme? Any help will be appreciated. And if there are any problems in my existing approach, please let me know.