|
[Sponsors] |
October 1, 2012, 13:26 |
Call face value of a parameter
|
#1 |
Member
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 14 |
Hi FOAMERS
I am trying to solve a set of equations in the following form -------- (Density)*fvm::ddt(U) + (Density/gL)*fvm::div(phi, U) - (Visc)*(fvm::laplacian(U)) - GravityVector*g*Density*((BetaT*(T-TNot))+BetaC*(C-CNot)) + fvm::Sp(gL*Visc/Perm,U) -------- Perm is the permeability which is calculated by the following relation Perm[celli] = (pow(SDAS.value(),2))*pow((gL[celli]),3)/(180*(pow((1-gL[celli]),2))); Is there any way to make sure that the permeability at the cell face is calculated by mean harmonic interpolation? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Netgen] Import netgen mesh to OpenFOAM | hsieh | OpenFOAM Meshing & Mesh Conversion | 32 | September 13, 2011 05:50 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 07:36 |
fluent add additional zones for the mesh file | SSL | FLUENT | 2 | January 26, 2008 11:55 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 14:00 |
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues | michele | OpenFOAM Meshing & Mesh Conversion | 2 | July 15, 2005 04:15 |