[Maxxxx]

Hey,


I have coded a modified version of icoFoam with another equation. Now I would like to restart the calculation after specific time meaning I would like to reset U and p to initial Values of time 0 (sounds confusing but is required because my implemented equation lead to a huge change in local U field... thus using the actual flow field lead to a divergence... that´s why I would like to start from scratch there).


I tried to do this by resetting phi, U and p field to initial values (which actually works) in each time step meaning I should get the same result on each step. But that doesn´t work. Hence I tried to figure out the line which destroys my simulation (printing phi, U and p after several lines of code).


I found out that the mistake is done with " fvc::grad(p) " which is not calculated correctly. Now my Question: What exactly happens with fvc::grad(p)? Is there anything I need to reset to let it calculate correctly? For time step 1 it is zero in the internal field (defined by my init p field which is 0 at the whole mesh). At the second time step the gradient is not zero for each cell (but p is same as initial step).


Thanks.