Different one sided derivative approximation?

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

 January 5, 2006, 00:56 Different one sided derivative approximation? #1 zonexo Guest   Posts: n/a In my FVM code, I am trying to approximate the face value of pressure at the boundary wall. Since it is a cell-center code, all pressure values are at the cell's center. For e.g. at the south/bottom wall, dp/dy=0. the simplest approx will be p_s=p(i,j) higher accuracy from Blazark's CFD book gives p_s=(1/8)*(15*p(i,j)-10*p(i,j+1)+3*p(i,j+2)). however, it is also possible to use dp/dy=(1/2h)*(3*p_s-4*p_n+p_nn) where h=delta, p_n=(p(i,j)+p(i,j+1))/2, p_nn=(p(i,j+1)+p(i,j+2))/2 simplifying, p_s=(2/3)*p(i,j)-(1/3)*p(i,j+1)-(1/6)*p(i,j+2) based on the no. of values used, they should be of the same accuracy, but why are their values different? Also, does anyone has similar expression for non-uniform grids? thanks

 January 5, 2006, 01:31 Re: Different one sided derivative approximation? #2 Praveen. C Guest   Posts: n/a They may be of the same order but the error can be different. Remember that the error is of the form Chr and different schemes of the same order can have different values for C. For non-uniforms grids you can fit a polynomial p(y) = ps + p1 y + p2 y2 + ... and also use dp/dy(y=0) = 0 to obtain the value of ps

 January 11, 2006, 05:40 Re: Different one sided derivative approximation? #3 sourabh Guest   Posts: n/a i think there is a mistake in ur simplification, it should come as Ps = (2/3)*Pi,j + (1/2)*Pi,j+1 - (1/6)*Pi,j+2

 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 jollage Main CFD Forum 0 September 21, 2011 18:16 ruirui389 Fluent UDF and Scheme Programming 4 April 8, 2011 04:12 sucker Main CFD Forum 1 April 15, 2008 13:28 Alessandro CFX 10 April 21, 2006 10:31 Gabriel Main CFD Forum 3 May 11, 2000 09:24

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