On 1/12/20 9:17 PM, David Eaton wrote: > My inflow condition is uniform. This formulation and mesh is tested in a > simple C++ code without library. TheĀ large mesh near the inflow does not > give > this problem. > Yes. I am using C0 element. I did calculation using tecplot. However, the > results from a my C++ code does not give this issue either. Just now, I check > the formulation again. Although I use Q2Q1 Taylor-Hood element without any > stabilization, these issues are still happening.
David -- we don't really know what formulation you are using, how you are implementing it, what you are comparing against, and a number of other factors. If you have a formulation that computes u,p, and you are plotting the vorticity, you need to expect that the isocontours are discontinuous for the reasons Praveen already stated. If you are getting results that make no sense to you, then the first step would be to ensure that your program is converging as expected. To do this, choose a solution that you know and compute the error norm; then ensure that the program yields error norms that decrease as expected with mesh refinement. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/e1b55e95-a139-d62e-e786-61393cf84dc7%40colostate.edu.