Hi there
I am trying to use
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 ? |
All times are GMT -4. The time now is 02:29. |