CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   problem with wallShearstress at interFaom (https://www.cfd-online.com/Forums/openfoam/98813-problem-wallshearstress-interfaom.html)

jignesh_thaker2007 March 20, 2012 02:45

problem with wallShearstress at interFaom
 
Hi everyone

Can you tel me how to write the value of wallShearstress in laminar pipe flow.

Actually i gave the command of wallShearstress but it cant write the value of shearstress but it gave me some error like
Create time

Create mesh for time = 0

Time = 0
Reading field U

Reading/calculating face flux field phi



--> FOAM FATAL IO ERROR:
keyword transportModel is undefined in dictionary "/home/turbo3/Desktop/jignesh/1inchpipe/pref0sample/constant/transportProperties"

file: /home/turbo3/Desktop/jignesh/1inchpipe/pref0sample/constant/transportProperties from line 20 to line 34.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 400.

FOAM exiting


so what can i do?
can anybody help me?

Jignesh

Andrea_85 March 20, 2012 04:51

Hi Jignesh,

i think wallShearStress utility is not suitable for multiphase. I have a similar problem with force calculation (have a look here http://www.cfd-online.com/Forums/ope...tml#post350266). So if you want to use it in your multiphase cases you have to write your own code and read the correct "transportProperties" dict.

best

andrea

Arnoldinho March 20, 2012 09:15

In case you are only interested in the stresses e.g. for the water-phase, you could simply copy one of the one-phase dictionaries (with values for water) into your case and temporarily use that one for the utility.

Arne

jignesh_thaker2007 March 21, 2012 02:38

Quote:

Originally Posted by Arnoldinho (Post 350417)
In case you are only interested in the stresses e.g. for the water-phase, you could simply copy one of the one-phase dictionaries (with values for water) into your case and temporarily use that one for the utility.

Arne

thanx for ur reply but when i will used only water phase, is it gave me correct value? because in my pipe water and air both are flowing at the same time and my flow pattern is stratified. so is it possible to use this in two phase flow?

thanx
jignesh

Arnoldinho March 21, 2012 03:55

That depends on what you are interested in.
If its only the water, then yes, if its both, then no. In the latter case you indeed have to modify the utility to account for both phases. I thought you could only be interested in one phase, therefore my hint.

Andrea_85 March 21, 2012 04:32

Hi all.
i think just multiply by nu_1 is not correct. In wallShearStress.C you have:

wallShearStress.boundaryField()[patchi] =
(
-mesh.Sf().boundaryField()[patchi]
/mesh.magSf().boundaryField()[patchi]
) & Reff.boundaryField()[patchi];
and

Reff is defined, if for example you are running laminar, as:
-nu()*dev(twoSymm(fvc::grad(U_))).


In my opinion you have to multiply somewhere by alpha to eliminate the contribution from the zones where there is fluid 2. or not?

best

andrea



jignesh_thaker2007 March 21, 2012 04:54

Quote:

Originally Posted by Andrea_85 (Post 350585)
Hi all.
i think just multiply by nu_1 is not correct. In wallShearStress.C you have:

wallShearStress.boundaryField()[patchi] =
(
-mesh.Sf().boundaryField()[patchi]
/mesh.magSf().boundaryField()[patchi]
) & Reff.boundaryField()[patchi];
and

Reff is defined, if for example you are running laminar, as:
-nu()*dev(twoSymm(fvc::grad(U_))).


In my opinion you have to multiply somewhere by alpha to eliminate the contribution from the zones where there is fluid 2. or not?

best

andrea




Thanx sir

but i am solving two phase laminar flow pattern.

so if possible then please tel me how i modified wallShearstress.C?

i used two fluid air and water.

so how can i calculate wallShearstress.C?
Thanx
jignesh

miael March 30, 2012 09:43

Hello jignesh,

do you define the transportModel in your constant/transportProperties-file? If not just add the line:

transportModel Newtonian;

when you have a newtonian fluid.

Greetings


All times are GMT -4. The time now is 03:24.