# Manufactured solutions for incompressible Navier-Stokes

 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 24, 2013, 10:56 Manufactured solutions for incompressible Navier-Stokes #1 New Member   Italo Join Date: Sep 2009 Posts: 4 Rep Power: 9 I am testing a solver for the incompressible Navier-Stokes equations using manufactured solutions. The code uses a 1st order projection method on staggered grid, very similar to the MAC method. I use 1st order time stepping and second order spatial derivatives. If I choose a manufactured solution that is divergence-free, all works fine. I get second order convergence for a stationary problem. But if the manufactured solution is not divergence-free, the method drops to 1st order, even though I added the source term for the Poisson equation. I found some publications that also use divergence-free manufactured solutions, but none explain why. Must the manufactured solution be divergence-free? Isn't it possible to test the code with a more generic solution? It seems intuitive that it should work. When solving the Poisson equation on a square surrounded by Neumann boundary conditions, I must prescribe the pressure at one point. And that is exactly on that point that the error is big. And the error only shows up when using non-divergence-free manufactured solutions.

March 25, 2013, 10:14
#2
New Member

Mlohry
Join Date: Mar 2013
Posts: 2
Rep Power: 0
Quote:
 Originally Posted by Italo I am testing a solver for the incompressible Navier-Stokes equations using manufactured solutions. The code uses a 1st order projection method on staggered grid, very similar to the MAC method. I use 1st order time stepping and second order spatial derivatives. If I choose a manufactured solution that is divergence-free, all works fine. I get second order convergence for a stationary problem. But if the manufactured solution is not divergence-free, the method drops to 1st order, even though I added the source term for the Poisson equation.
If the solution you're using is not divergence free, then it wouldn't satisfy incompressible NS, and an incompressible solver couldn't accurately recreate the solution. Or is there something else I'm missing?

Similarly, your numerical method would have to drop to first order if your solution contains discontinuities, so I wouldn't be surprised to see only first order convergence for any such solution.

March 25, 2013, 10:20
#3
New Member

Italo
Join Date: Sep 2009
Posts: 4
Rep Power: 9
Quote:
 Originally Posted by mlohry If the solution you're using is not divergence free, then it wouldn't satisfy incompressible NS, and an incompressible solver couldn't accurately recreate the solution. Or is there something else I'm missing?
The manufactured solution is not a solution for the original NS equations. A source term must be added to the equations so that it becomes a solution.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Palani Velladurai Main CFD Forum 7 September 6, 2013 02:51 Raju Main CFD Forum 5 July 29, 2006 14:15 Jose Choy Main CFD Forum 3 October 24, 2003 02:28 Rajil Saraswat Main CFD Forum 2 June 9, 2003 07:21 Jose Choy Main CFD Forum 2 May 18, 2000 05:45

All times are GMT -4. The time now is 19:54.