# Upwind scheme (in a cell with zero velocity)

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 2, 2013, 12:13 Upwind scheme (in a cell with zero velocity) #1 New Member   Molly Tom Join Date: May 2013 Posts: 3 Rep Power: 5 I am somewhat new to CFD and am wondering how the first order upwind scheme would be applied to a cell with a zero velocity. It is my understanding that, in the upwind scheme, the sign of the cell velocity is used to determine which cell velocities are applied at the cell faces. But what happens when the cell velocity is zero (how is the upwind direction determined)? This is probably not important in a large region of zero velocity cells, but I'm thinking about a zero velocity cell that is adjacent to a non-zero velocity cell.

May 2, 2013, 12:19
#2
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,694
Rep Power: 33
Quote:
 Originally Posted by MollyTom I am somewhat new to CFD and am wondering how the first order upwind scheme would be applied to a cell with a zero velocity. It is my understanding that, in the upwind scheme, the sign of the cell velocity is used to determine which cell velocities are applied at the cell faces. But what happens when the cell velocity is zero (how is the upwind direction determined)? This is probably not important in a large region of zero velocity cells, but I'm thinking about a zero velocity cell that is adjacent to a non-zero velocity cell.

where you have zero velocity there is no convective flux contributing to the balance....

 May 2, 2013, 13:17 #3 New Member   Molly Tom Join Date: May 2013 Posts: 3 Rep Power: 5 But wouldn't it matter in a case like this? |------------------|-------------------|-----------------| |-------i-1---------|---------i---------|------- i+1-------| |-----u(i-1)=1----- |-------u(i)=0------|----- u(i+1)=0---| |------------------ |------------------|------------------| ---------------u(i-1/2)=???---------u(i+1/2)=0 So, for cell (i), shouldn't the first order upwind scheme assign a value of u(i-1/2)=1? But how would the upwind scheme determine this, because the velocity at cell (i) is zero?

 May 2, 2013, 13:47 #4 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 2,694 Rep Power: 33 - In FD scheme, you evaluate u*du/dx which is zero at node i as u(i)=0. - In FV scheme you need to evaluate the flux difference f(i+1/2) - f(i-1/2). Now you can evaluate the velocity at the interface by means of linear interpolation: u(i-1/2) = 0.5*(u(i)+u(i-1)); u(i+1/2) = 0.5*(u(i)+u(i+1)) Then the sign will be the direction for evaluating the upwinded flux. Note that your example is a discontinuity travelling immortality likes this.

 May 6, 2013, 14:04 #5 New Member   Molly Tom Join Date: May 2013 Posts: 3 Rep Power: 5 Yes. Thank you very much, I think I understand now. For computing the upwind flux direction, the velocity at the cell interface is computed by averaging the velocities at the cell centers that bound the interface. But, for the first-order upwind scheme (once the upwind direction has been determined), the interface fluxes are evaluated using the appropriate cell center velocity. Anne Lincke likes this.

 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 alfaruk CFX 8 December 3, 2013 17:51 lnk Main CFD Forum 1 August 30, 2012 02:38 vitke OpenFOAM Running, Solving & CFD 0 September 3, 2007 06:15 michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15 Szabolcs Varga FLUENT 1 September 16, 2004 07:30

All times are GMT -4. The time now is 22:41.

 Contact Us - CFD Online - Top