# Formula for 2nd order upwind scheme for non-uniform grids?

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

 August 30, 2012, 03:57 Formula for 2nd order upwind scheme for non-uniform grids? #1 Senior Member   TWB Join Date: Mar 2009 Posts: 105 Rep Power: 8 Hi, I've seen the formula for 2nd order upwind scheme for uniform grids as: phi_e = 1.5*phi_P - 0.5*phi_W for u_e > 0 phi_e = 1.5*phi_E - 0.5*phi_EE for u_e < 0 Is there a formula meant for non-uniform grids as well? Thanks

 August 30, 2012, 04:57 #2 Senior Member   Join Date: Dec 2011 Location: Madrid, Spain Posts: 133 Rep Power: 6 Hi, you can try to derive it. Just take the Taylor expansions around the nodes E and EE, taking into account the different grid spacing for node i. The formulas you showed is a particular case where the grid spacing is constant. Under this assumption, the grid spacings appearing in numerator and denominator of your expressions when you solve for phi_e cancel out and disappears from the equation. However, the different will appear in the expression of phi_e if you consider grid non-uniformity. I hope it helps. Cheers, Michujo.

August 30, 2012, 17:25
#3
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,588
Rep Power: 20
Quote:
 Originally Posted by quarkz Hi, I've seen the formula for 2nd order upwind scheme for uniform grids as: phi_e = 1.5*phi_P - 0.5*phi_W for u_e > 0 phi_e = 1.5*phi_E - 0.5*phi_EE for u_e < 0 Is there a formula meant for non-uniform grids as well? Thanks
Be careful that upwind schemes are usually first, third, fifth ... order accurate... the reason is to have a local truncation error that has a dissipative behaviour ... the above formulas are based on a linear extrapolation of the value, that is highly unstable in general.

I suggest to construct a second degree polynomial on non-uniform stencil in such a way that if u_e>0 then Phi_W,Phi_P, Phi_E are involved (the counterpart for u_e<0 involves Phi_P,Phi_E, Phi_EE).
This way the use of Phi_P ensures better stability properties

 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 msarkar OpenFOAM Running, Solving & CFD 33 January 7, 2015 23:10 marine OpenFOAM 66 October 26, 2014 12:10 Sebaj OpenFOAM 7 October 31, 2012 15:16 mihaipruna OpenFOAM Running, Solving & CFD 5 April 24, 2012 17:18 Heinz Wilkening Main CFD Forum 2 November 3, 1998 15:33

All times are GMT -4. The time now is 20:03.