
[Sponsors] 
uniform discretization scheme on nonuniform grid 

LinkBack  Thread Tools  Search this Thread  Display Modes 
March 27, 2019, 12:02 
uniform discretization scheme on nonuniform grid

#1 
New Member
Theo
Join Date: Mar 2009
Posts: 26
Rep Power: 17 
imagine using a discretization scheme for a uniform grid, e.g.
(x_(i+1) + x_(i1)  2*x_i)/h² for a second derivative, even if the grid is not uniform, i.e. h_(i+1), h(i1) and h_i are not equal. Does that introduce a significant error? or is there a certain grid stretching below which it is acceptable? Reason for this question is that uniform schemes are obviously much easier to implement and faster to compute. 

March 27, 2019, 12:06 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,831
Rep Power: 73 
Quote:
Of course the answer is that you have a wrong formula on non uniform grid. To check this try to expand the local truncation error. The correct formula is very easy to write, you can write the lagrangian polynomial of second degree on three non equidistant nodes. The second derivative is constant in such interval. 

March 27, 2019, 12:49 

#3 
New Member
Theo
Join Date: Mar 2009
Posts: 26
Rep Power: 17 
sure, but i wonder if there is a limit up to which this error is acceptable?


March 27, 2019, 12:53 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,831
Rep Power: 73 
Quote:
Compute the local truncation and compare the two formula. Check what happens for vanishing mesh size. Remember that to hope in a physically relevant solution you need consistence and stabilty when you perform a discretization of an equation. Have also a look to the book of Peric and Ferziger about the discretization on non uniform grids 

March 27, 2019, 13:39 

#5 
Senior Member

Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14 
For constant problems the error is acceptable. For anything else, dont do it.


Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
dsmcFoam setup  hherbol  OpenFOAM PreProcessing  1  November 19, 2021 01:52 
rhoPimpleFoam hardship  petrus  OpenFOAM Running, Solving & CFD  0  October 7, 2016 02:41 
Time continuity error & FAN patch  Zephiro88  OpenFOAM Running, Solving & CFD  4  April 22, 2015 12:39 
simpleFoam Floating point error (Core Dumped)  sam.ho  OpenFOAM Running, Solving & CFD  2  December 31, 2013 06:57 
Grid Independence with discretization schemes  saisanthoshm88  CFX  19  May 22, 2013 10:10 