
[Sponsors] 
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: 6 
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: 1,662
Rep Power: 23 
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: 6 
Quote:
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 1e6 or so. 

July 1, 2013, 05:42 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23 
Quote:
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: 6 

July 1, 2013, 05:59 

#6 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23 
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: 6 
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: 1,662
Rep Power: 23 
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: 6 
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: 1,662
Rep Power: 23 
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: 6 
Thanks, and if somebody already knows any refs, please let me know. Thanks


Thread Tools  
Display Modes  


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 