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

A 2D test case for Linear reconstruction FV scheme?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 1, 2013, 04:10
Default A 2D test case for Linear reconstruction FV scheme?
  #1
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Hello,

Does anybody know a 2D test case for accuracy study of a linear reconstruction Finite Volume or maybe the MUSCL scheme?
I would prefer a test case without periodic boundary conditions.

Many thanks,
Faraz
Faraz is offline   Reply With Quote

Old   July 1, 2013, 04:45
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Faraz View Post
Hello,

Does anybody know a 2D test case for accuracy study of a linear reconstruction Finite Volume or maybe the MUSCL scheme?
I would prefer a test case without periodic boundary conditions.

Many thanks,
Faraz
You could test a linear solution of Df/Dt=0 in a 2D velocity filed u=-omega*y, v= omega*x which is a rigid rotation. Set some initial field f(x,y,0), after one complete rotation the solution should be equal to the initial state
FMDenaro is offline   Reply With Quote

Old   July 1, 2013, 05:02
Default
  #3
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
You could test a linear solution of Df/Dt=0 in a 2D velocity filed u=-omega*y, v= omega*x which is a rigid rotation. Set some initial field f(x,y,0), after one complete rotation the solution should be equal to the initial state

Thanks,

That is a very interesting example. Actually I did something similar to that, but within the domain x = [-1,1] & y = [0,1] for which there seems to be a singularity of velocity at (x,y) = (0,0).

However, I did another example for an angular convection:

x = [0,1], y = [0,1],

and veloctiy: u = Omega, v = Omega

So, it converges with correct rate for N = 8, 16, 32, ..., 512 (conv rate = 2 for linear)

But after N = 512, 1024 the convergence rate degrades for some reason (conv rate = 1.7 e.g. for linear), which is quite strange ! and it is not due to zero error space because the error is still 1e-6 or so.
Faraz is offline   Reply With Quote

Old   July 1, 2013, 05:42
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Faraz View Post
Thanks,

That is a very interesting example. Actually I did something similar to that, but within the domain x = [-1,1] & y = [0,1] for which there seems to be a singularity of velocity at (x,y) = (0,0).

However, I did another example for an angular convection:

x = [0,1], y = [0,1],

and veloctiy: u = Omega, v = Omega

So, it converges with correct rate for N = 8, 16, 32, ..., 512 (conv rate = 2 for linear)

But after N = 512, 1024 the convergence rate degrades for some reason (conv rate = 1.7 e.g. for linear), which is quite strange ! and it is not due to zero error space because the error is still 1e-6 or so.

but u=v=constant is a uniform translation velocity field ...
FMDenaro is offline   Reply With Quote

Old   July 1, 2013, 05:55
Default
  #5
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
but u=v=constant is a uniform translation velocity field ...
That's true, but would it cause a problem?
Faraz is offline   Reply With Quote

Old   July 1, 2013, 05:59
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
what about the initial condition and BCs ?
FMDenaro is offline   Reply With Quote

Old   July 1, 2013, 06:48
Default
  #7
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
what about the initial condition and BCs ?
Well, for angular convection, the initial profile at inflow ghost cell is either a Gaussian or a sine/cosine distribution on ghost cells adjacent to either x = 0 and y =0, and zero everywhere else.

Regarding BC, for left state of the face at the inflow BC, an average of the ghost and the adjacent internal cell is performed. Because with the constant reconstruction at the inflow BC, the convergence gets close to 1st order other than 2nd order.

Is there a better boundary treatment for the inflow?
I also noticed that the convergence rate is a bit dependent on initial profile !
Faraz is offline   Reply With Quote

Old   July 1, 2013, 07:07
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
x=0 and y=0 are inflow boundaries where you prescribe Dirichlet BC, x=1 and y=1 are outflow boundaries where you can prescribe homogeneous Neumann BC. What you can do is to run in time until a steady state (for steady inflow) is reached. Then compute the errors and repeat for via via refined grids.
Use the same dt for all run, take value that allow a stable computation on the most refined grid.
FMDenaro is offline   Reply With Quote

Old   July 1, 2013, 07:49
Default
  #9
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Yes indeed, but unfortunately the same degradation of order in higher number of cells remains. But, as I mentioned, it is also a bit dependent on the initial profile I define on the Dirichlet BC.

Apart from this problem, do you think if with the rigid rotation example I should obtain the right order of convergence (in spite of that singularity in (0,0) )? Do you, by any chance, know any references to check with that?

I really can't find any relevant examples in the literature or papers.

Thanks again.
Faraz is offline   Reply With Quote

Old   July 1, 2013, 08:16
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
the initial and boundary conditions must be smooth to check the correct order of accuracy. I remember several tests in journals, check in J.Comp.Phys as well as in Computer & FLuids, Int.J. Num. Meth.FLuids.

Have also a look at the book of LeVeque
FMDenaro is offline   Reply With Quote

Old   July 2, 2013, 05:15
Default
  #11
New Member
 
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 6
Faraz is on a distinguished road
Thanks, and if somebody already knows any refs, please let me know. Thanks
Faraz 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
rhoSimpleFoam - test case atareen64 OpenFOAM Running, Solving & CFD 5 May 1, 2011 17:17
critical error during installation of openfoam Fabio88 OpenFOAM Installation 21 June 2, 2010 03:01
Problems in compiling paraview in Suse 10.3 platform chiven OpenFOAM Installation 3 December 1, 2009 08:21
Need help with the error cfdproject OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 0 April 14, 2009 15:45
Porous Media test case Alex FLUENT 0 April 9, 2006 08:23


All times are GMT -4. The time now is 11:37.