I already solved the stokes equation and was trying to generalize it to solve the steady NS equation.
I saw the step-57 tutorial but I have a mesh from file (while the tutorial create in dealii and refine, things that I'm not allowed to do) and also uses line search algorithm, which also can't perform. The general structure of the program is: setup vectors, sparsity and matrices (pressure and velocity), newton iteration where first solve the stokes problem to get the initial guess for newton (Re will always be < 200 in my test case, for now), and the start to iterate, at each iteration: - assemble the system - solve the linearized system I know the stokes solution is correct. But for the NS solution I have the residual decreasing at first 2 steps and then increasing, leading to divergence (around step 16 dealii fires the no convergence). The assembly is the same as the tutorial at step-57. It's becoming really frustrating, so if anyone can help I really appreciate. I give you the code in case my explanation wasn't clear. https://drive.google.com/drive/folders/1igXMgWq2g6UPRxaRE9Ot1jU5L_eUW9Wk?usp=sharing PS: the mesh is fine since already used for previous problems. -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/73ca42ee-e80b-4312-a8a4-9942b0785633n%40googlegroups.com.
