|
[Sponsors] |
September 4, 2015, 07:10 |
phi in scalarTransportFoam, why not U
|
#1 |
New Member
novo
Join Date: Jun 2015
Posts: 11
Rep Power: 10 |
Dear Openfoamers,
sorry for the ignorant question but I would appreciate if anybody could revert me to the proper info since I have been looking around and am still quite puzzled. In the scalarTransportFoam the equation to be solved includes the scalar "T" which is transported both by diffusion constant ("DT") and convection vector ("U"). But in the scalarTransportFoam the equation solved includes "phi" instead of "U". Code:
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) { // Building and solving the scalar tranport equation solve ( fvm::ddt(T) // Unsteady term + fvm::div(phi, T) // Convective term - fvm::laplacian(DT, T) // Diffusive term ); } |
|
December 18, 2015, 17:14 |
|
#2 |
Senior Member
Thomas Oliveira
Join Date: Apr 2015
Posts: 114
Rep Power: 12 |
"phi" is the scalar product of the of Uf and the area vector. Uf is the velocity at the face, while U is the velocity at the the center of the cell. If you used U at the discretized form of the divergence, you would have the divergence equals to zero no matter how your velocity field were.
|
|
Tags |
openfaom, scalartransportfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM.org] OpenFOAM build on Intel Xeon Phi | asaijo | OpenFOAM Installation | 31 | July 26, 2017 10:35 |
mixerVesselAMI2D's mass is not balancing | sharonyue | OpenFOAM Running, Solving & CFD | 6 | June 10, 2013 09:34 |
Summing phi over all faces in patch | NewtoFOAM | OpenFOAM Programming & Development | 2 | May 8, 2013 05:37 |
Recreation of phi file, icoFoam | Madeleine P. Vincent | OpenFOAM Running, Solving & CFD | 0 | May 7, 2012 11:17 |
Turbulence Model phi vs phi_ | doug | OpenFOAM Running, Solving & CFD | 4 | November 10, 2009 04:33 |