|
[Sponsors] |
December 10, 2008, 21:41 |
Hi there
I am trying to use
|
#1 |
New Member
ehsan vaghefi
Join Date: Mar 2009
Posts: 9
Rep Power: 17 |
Hi there
I am trying to use 'SimpleFoam' for a mesh where there is a known sink in each element. I have read the 'simplefoam' guide on line and it says : The first term of H{U} represents the matrix coefficients of the neighbouring cells multiplied by their velocity, while the second part contains the unsteady term and all the sources except the pressure gradient. So I have the idea of adding my sinks as VolScalarField to the Ueqn : tmp<fvvectormatrix> UEqn ( fvm::div(phi, U) + turbulence ->divDevReff(U) ); UEqn().relax(); eqnResidual = solve ( UEqn() == -fvc::grad(p) + SINK ).initialResidual(); maxResidual = max(eqnResidual, maxResidual); but it seems that no matter what I do, the answers don't change as if the SINK term is not there. How can I solve this problem ? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Define non uniform TKE boundary condition in simpleFOAM | qtian | OpenFOAM Running, Solving & CFD | 9 | July 30, 2007 14:54 |
mass sink UDF | H | FLUENT | 1 | November 8, 2004 22:44 |
define my own DEFINE macro? to model pithing naca | calogero | FLUENT | 0 | April 10, 2003 05:10 |
Element by element preconditioning | Jarmo Monttinen | Main CFD Forum | 0 | October 27, 2000 15:53 |
source or sink?? | frederic felten | Main CFD Forum | 2 | August 8, 2000 12:35 |