Dear all, I cannot find anything about velocity correction in the code of Step-35.
I read the references listed there and some other papers about projection methods, and know that after obtaining the intermediate values for velocity, say v_n, we need to correct the velocity by u_n = v_n + grad(2dt/3 phi_n), which is mentioned in the introduction of Step-35, but it seems that Step-35 only corrects the pressure, and I cannot find anything in the code that implements the above equation. However, when I correct the velocity at every time step, the algorithm turns out to be divergent (iterating more than 1000 times at some time step). I am very confused, did I misunderstand the algorithm or the code? Thanks in advance. Best, Jiaqi -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
