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

Upwind scheme (in a cell with zero velocity)

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

Like Tree2Likes
  • 1 Post By FMDenaro
  • 1 Post By MollyTom

Reply
 
LinkBack Thread Tools Display Modes
Old   May 2, 2013, 12:13
Default Upwind scheme (in a cell with zero velocity)
  #1
New Member
 
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 4
MollyTom is on a distinguished road
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.
MollyTom is offline   Reply With Quote

Old   May 2, 2013, 12:19
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,604
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by MollyTom View Post
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....
FMDenaro is offline   Reply With Quote

Old   May 2, 2013, 13:17
Default
  #3
New Member
 
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 4
MollyTom is on a distinguished road
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?
MollyTom is offline   Reply With Quote

Old   May 2, 2013, 13:47
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,604
Rep Power: 23
FMDenaro will become famous soon enough
- 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.
FMDenaro is offline   Reply With Quote

Old   May 6, 2013, 14:04
Default
  #5
New Member
 
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 4
MollyTom is on a distinguished road
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.
MollyTom is offline   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
Difficulty in calculating angular velocity of Savonius turbine simulation alfaruk CFX 8 December 3, 2013 17:51
Are backward scheme and the upwind scheme the same thing in FDM? lnk Main CFD Forum 1 August 30, 2012 02:38
Velocity from an upwind cell vitke OpenFOAM Running, Solving & CFD 0 September 3, 2007 06:15
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15
DPM velocity at cell face Szabolcs Varga FLUENT 1 September 16, 2004 07:30


All times are GMT -4. The time now is 19:37.