CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Advective BC: local velocity?

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

LinkBack Thread Tools Display Modes
Old   July 3, 2009, 08:48
Default Advective BC: local velocity?
New Member
Daniel Heidorn
Join Date: Jul 2009
Location: Graz,Austria
Posts: 3
Rep Power: 8
dany is on a distinguished road
Dear OpenFoamers,

I want to simulate a high-speed atomizing liquid jet with a high density ratio using the lesInterFoam code (OpenFoam 1.5) setting the LESModel to laminar (= "DNS"). The outlet BC gives me a headache. Unfortunately the advective BC is, as the simulation shows, not the right condition to use.

What I understand of the advective boundary condition (/OpenFOAM/OpenFOAM-1.5/src/finiteVolume/fields/fvPatchFields/derived/advective) is that (without defining fieldInf and lInf and lets say using the Euler ddtScheme) the refValue is the old velocity field of the BC.
This refValue is used in the mixedFvPatchField.C file (in /OpenFOAM/OpenFOAM-1.5/src/finiteVolume/fields/fvPatchFields/basic/mixed/) where the new LOCAL velocity is calculated switching between fixedValue and zeroGradient BC (or a blend between this two BC's).

As I see it, I must admit to be a noob in OF, this LOCAL velocity is used in the advective term of the equation:

delta U / delta t = - U_LOCAL * (delta U / delta x)

which tends to be more unstable as to use the MEAN velocity:

delta U / delta t = - U_MEAN * (delta U / delta x)

.) Am I wrong with my assumption that the code uses a LOCAL velocity in the advective term?
.) If not, is it the easiest way just to overwrite the solution of the LOCAL velocity in the mixedFvPatchField.C with a new calculated MEAN velocity?

Thx for your help.

With best regards
dany is offline   Reply With Quote


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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07
Variables Definition in CFX Solver 5.6 R P CFX 2 October 26, 2004 02:13

All times are GMT -4. The time now is 08:16.