CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Manufactured solutions for incompressible Navier-Stokes

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 24, 2013, 10:56
Default Manufactured solutions for incompressible Navier-Stokes
  #1
New Member
 
Italo
Join Date: Sep 2009
Posts: 4
Rep Power: 8
Italo is on a distinguished road
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.
Italo is offline   Reply With Quote

Old   March 25, 2013, 10:14
Default
  #2
New Member
 
Mlohry
Join Date: Mar 2013
Posts: 2
Rep Power: 0
mlohry is on a distinguished road
Quote:
Originally Posted by Italo View Post
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.
mlohry is offline   Reply With Quote

Old   March 25, 2013, 10:20
Default
  #3
New Member
 
Italo
Join Date: Sep 2009
Posts: 4
Rep Power: 8
Italo is on a distinguished road
Quote:
Originally Posted by mlohry View Post
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.
Italo is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Filtered navier stokes equation..LES:: Palani Velladurai Main CFD Forum 7 September 6, 2013 02:51
ILU for Navier stokes problems Raju Main CFD Forum 5 July 29, 2006 14:15
Navier stokes compresible viscid flow fea, somebody can help? Jose Choy Main CFD Forum 3 October 24, 2003 02:28
Newbie:Viscoelasticity and Navier stokes equation Rajil Saraswat Main CFD Forum 2 June 9, 2003 07:21
help: I am trying to solve Navier Stokes compressible and viscid flow Jose Choy Main CFD Forum 2 May 18, 2000 05:45


All times are GMT -4. The time now is 09:25.