Thank you for your advice. I would make the problem simpler. Do you get the right solution if you > have a constant viscosity? If you iterate, do you get the solution after > one iteration? >
Yes, I checked this. After one iteration, I get the solution after one iteration when constant viscosity case. > > I will also note that your manufactured solution is symmetric, but your > computational solution is not. This provides you with a powerful way > because you don't actually need to look at the convergence, it's enough > to check that on a *coarse* mesh the solution is *symmetric*. It may be > that solution is already non-symmetric on 1 cell, or on a 2x2 mesh -- in > which case you already know that something is wrong. The question then > is why that is so -- is your rhs correct, for example? > After modified some of mistakes, now I have following result for pressure field. With same manufactured pressure solution... <https://lh3.googleusercontent.com/-f-k63HKKEis/WJxNxEUa3lI/AAAAAAAAA_U/xiH9SS-bToAom3KyMF4eq2k9hna-6TzawCLcB/s1600/manufacturedp.png> <Manufactured Solution> <https://lh3.googleusercontent.com/-8Y_ZBl2EiNk/WJxOjKwv54I/AAAAAAAAA_c/pkkQ4nfpGfoFxRJY42_j7osOrMFAw2uNQCLcB/s1600/mesh1.png> <2x2 mesh> <https://lh3.googleusercontent.com/-1jqshnea3YQ/WJxOl6Hfm0I/AAAAAAAAA_g/mMUyP0zHc_siKD_p92HrIBekhV3p0DJWQCLcB/s1600/mesh2.png> <4x4 mesh> <https://lh3.googleusercontent.com/-ijDs6iJvzyo/WJxOoe_5o9I/AAAAAAAAA_k/wKnfLtF5O1kskPAH6hI3_3GeUwgdZQdzACLcB/s1600/mesh4.png> <more refined mesh> For velocity field, I still have at least qualitatively good result as I posted previous comment. Tho it does not converge with mesh refinement when it is compared to manufactured solution in L2_norm base. > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > <javascript:> > www: http://www.math.colostate.edu/~bangerth/ > -- 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. For more options, visit https://groups.google.com/d/optout.
