problem with wallShearstress at interFaom
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 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.
so what can i do?
can anybody help me?
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.
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.
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.
i think just multiply by nu_1 is not correct. In wallShearStress.C you have:
) & Reff.boundaryField()[patchi];and
Reff is defined, if for example you are running laminar, as:
In my opinion you have to multiply somewhere by alpha to eliminate the contribution from the zones where there is fluid 2. or not?
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?
do you define the transportModel in your constant/transportProperties-file? If not just add the line:
when you have a newtonian fluid.
|All times are GMT -4. The time now is 11:43.|