CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Different one sided derivative approximation?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 5, 2006, 00:56
Default 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
  Reply With Quote

Old   January 5, 2006, 01:31
Default 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 Ch<sup>r</sup> and different schemes of the same order can have different values for C.

For non-uniforms grids you can fit a polynomial

p(y) = p<sub>s</sub> + p<sub>1</sub> y + p<sub>2</sub> y<sup>2</sup> + ...

and also use

dp/dy(y=0) = 0

to obtain the value of p<sub>s</sub>
  Reply With Quote

Old   January 11, 2006, 05:40
Default 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
  Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Derivative matrix jollage Main CFD Forum 0 September 21, 2011 18:16
can i get pressure derivative in UDF? ruirui389 Fluent UDF and Scheme Programming 4 April 8, 2011 04:12
scalar and its derivative sucker Main CFD Forum 1 April 15, 2008 13:28
Second derivative with cfx 5.7.1 Alessandro CFX 10 April 21, 2006 10:31
Boussinesq approximation again Gabriel Main CFD Forum 3 May 11, 2000 09:24


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