|
[Sponsors] |
September 4, 2015, 08:10 |
phi in scalarTransportFoam, why not U
|
#1 |
New Member
novo
Join Date: Jun 2015
Posts: 11
Rep Power: 11 |
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, 18: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 |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM.org] OpenFOAM build on Intel Xeon Phi | asaijo | OpenFOAM Installation | 31 | July 26, 2017 11:35 |
mixerVesselAMI2D's mass is not balancing | sharonyue | OpenFOAM Running, Solving & CFD | 6 | June 10, 2013 10:34 |
Summing phi over all faces in patch | NewtoFOAM | OpenFOAM Programming & Development | 2 | May 8, 2013 06:37 |
Recreation of phi file, icoFoam | Madeleine P. Vincent | OpenFOAM Running, Solving & CFD | 0 | May 7, 2012 12:17 |
Turbulence Model phi vs phi_ | doug | OpenFOAM Running, Solving & CFD | 4 | November 10, 2009 05:33 |