CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

SimpleFoam Define a Sink on every element

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 10, 2008, 21:41
Default Hi there I am trying to use
  #1
New Member
 
ehsan vaghefi
Join Date: Mar 2009
Posts: 9
Rep Power: 7
ehsan_vaghefi is on a distinguished road
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 ?
ehsan_vaghefi is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 21:41.