|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 13 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
Quote:
where you have zero velocity there is no convective flux contributing to the balance.... |
||
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 13 ![]() |
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? |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
- 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 |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Molly Tom
Join Date: May 2013
Posts: 3
Rep Power: 13 ![]() |
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.
|
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Difficulty in calculating angular velocity of Savonius turbine simulation | alfaruk | CFX | 14 | March 17, 2017 06:08 |
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 |
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues | michele | OpenFOAM Meshing & Mesh Conversion | 2 | July 15, 2005 04:15 |
DPM velocity at cell face | Szabolcs Varga | FLUENT | 1 | September 16, 2004 07:30 |