# A 2D test case for Linear reconstruction FV scheme?

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

 July 1, 2013, 04:10 A 2D test case for Linear reconstruction FV scheme? #1 New Member   Faraz Join Date: Jan 2011 Location: Netherlands Posts: 10 Rep Power: 8 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

July 1, 2013, 04:45
#2
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,564
Rep Power: 40
Quote:
 Originally Posted by Faraz 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

July 1, 2013, 05:02
#3
New Member

Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 8
Quote:
 Originally Posted by FMDenaro 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.

July 1, 2013, 05:42
#4
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,564
Rep Power: 40
Quote:
 Originally Posted by Faraz 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 ...

July 1, 2013, 05:55
#5
New Member

Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 8
Quote:
 Originally Posted by FMDenaro but u=v=constant is a uniform translation velocity field ...
That's true, but would it cause a problem?

 July 1, 2013, 05:59 #6 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,564 Rep Power: 40 what about the initial condition and BCs ?

July 1, 2013, 06:48
#7
New Member

Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 8
Quote:
 Originally Posted by FMDenaro 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 !

 July 1, 2013, 07:07 #8 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,564 Rep Power: 40 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.

 July 1, 2013, 07:49 #9 New Member   Faraz Join Date: Jan 2011 Location: Netherlands Posts: 10 Rep Power: 8 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.

 July 1, 2013, 08:16 #10 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,564 Rep Power: 40 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

 July 2, 2013, 05:15 #11 New Member   Faraz Join Date: Jan 2011 Location: Netherlands Posts: 10 Rep Power: 8 Thanks, and if somebody already knows any refs, please let me know. Thanks

 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 atareen64 OpenFOAM Running, Solving & CFD 5 May 1, 2011 17:17 Fabio88 OpenFOAM Installation 21 June 2, 2010 03:01 chiven OpenFOAM Installation 3 December 1, 2009 08:21 cfdproject OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 0 April 14, 2009 15:45 Alex FLUENT 0 April 9, 2006 08:23

All times are GMT -4. The time now is 18:15.