|July 3, 2009, 08:48||
Advective BC: local velocity?
Join Date: Jul 2009
Posts: 3Rep Power: 8
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
|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|