[Users] Timelevels on refined meshes

Michail Chabanov chabanov at th.physik.uni-frankfurt.de
Fri Oct 15 16:53:05 CDT 2021


Hello,

I have a rather general question on how previous timelevels are filled 
if one has more than one refinement level. So, assuming I have GF called 
RHO which has 3 timelevels and I have a mesh which consists of a coarse 
and a refined part, so 2 levels. On the coarse mesh I will have the 
timestep DT while on the fine mesh I will have dt = DT/2. The zeroth 
timelevel of RHO is the solution at time t on all levels.

My first question is whether RHO_p on the refined grid lives on t - dt 
or on t - DT ( t - 2*dt or t -2*DT for RHO_p_p ) before timelevel 
rotation takes place?

In the case that RHO_p lives on t - dt in the refined region, I would 
have another question:

Imagine I have the same setup as before, evolve for one coarse step such 
that the whole solution is at t + DT now, no timelevel rotation yet. If 
I now decide to regrid and coarsen the refined region such that I only 
have one refinement level in the whole domain, which values will be 
given in RHO_p and RHO_p_p in the region which was previously refined 
and is now coarse?

Do we have to extrapolate to get RHO_p_p at t - DT, because after the 
evolution has provided RHO on t + DT in the fine domain but before 
coarsening we will have RHO_p at t+dt and RHO_p_p at t if my 
understanding is correct. This means we would have to copy RHO_p_p to 
RHO_p and extrapolate RHO_p_p to t - DT. Otherwise, the three timelevels 
of RHO in this specific region wouldn't be equidistant anymore after the 
next evolution step of size DT ( RHO at t+ 2*DT, RHO_p at t + DT and 
RHO_p_p at t + dt).

Thanks,

Michail Chabanov



More information about the Users mailing list