CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (https://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Error in channelOodles tutorial (https://www.cfd-online.com/Forums/openfoam-bugs/62407-error-channeloodles-tutorial.html)

torvic July 24, 2008 20:11

Hello Hope you're doing fin
 
Hello

Hope you're doing fine.
I have a question about "fieldAverage"

I saw its use in the controlDict of oodles, but where can I find the code for the "fieldAverage" function. I checked out "sampling" directory but only found "probes".

I want to take the average of other quantities, say a chemical specie.

Does "fieldAverage" allow for this ?

thanks in advance

Best

V

lakeat July 24, 2008 20:53

1. It's in src/postProcessing.
 
1. It's in src/postProcessing.
2. I think so.

lakeat July 26, 2008 08:44

Well, I know where the problem
 
Well, I know where the problem is.
According to the numbering of cells, postChannel's channelIndexing.C should have a few modifications.

Best,
Daniel

torvic July 28, 2008 17:56

Hi Daniel sorry for the lat
 
Hi Daniel

sorry for the late reply.

you're right, fieldAverage it's already there. I will see how to add the chemical specie field

Thanks so much for your help !

best

V

nishant_hull September 3, 2008 12:23

Hi foamers... I can see
 
Hi foamers...

I can see that the postChannel utility is not working as the structure of new channelOodles solver has been changed. Now the R field is not being written.
however I want to compare the stress components across the channel and therefore i need the urms and vrms value of the the velocity. ( I took these urms and vrms terms from postchannel utility declaration). How can we obtain these stress values now?

I again have some doubt about the UPrime2Mean. is this UPrime2Mean is stress-components along xx, yy and zz?? If yes, then what about the stress component in XY direction? (which was termed as "uv" in postChannel utility)

your comments plase!

cheers!
Nishant

lakeat September 3, 2008 13:21

1. Easy, just copy some lines
 
1. Easy, just copy some lines from OF-1.4.1
2. Please read the ProgrammerGuide P-16, item "Symmetric rank 2".

Daniel

nishant_hull September 3, 2008 14:47

Thanks Daniel 2) I unders
 
Thanks Daniel

2) I understand that the 6 components will give me the all stresses and I can see those 6 stresses in the file UPrime2Mean. but How to extract/sample the individual component of the stress from the file UPrime2Mean? I tried sample utility but they do not sample UPrime2Mean field. Doing it in excel is another option but isn't there some available method?

1) Do you mean, I should copy those postchannel-utility files and recompile my openfoam solver? Could you please tell me which files you are referring to?

Thanks once again.
Nishant

lakeat September 3, 2008 21:04

1) Compare version 1.4.1 and 1
 
1) Compare version 1.4.1 and 1.5, note their differences, and you will find a file called calculateAverages.H, try to use it in OF-1.5, and recompile the channeloodles solver.

2) For me, I use postChannel Utility as a template. Note how the utility treat R field and get u'u', u'v', etc., you will know what I mean.

Regards, Daniel

nishant_hull September 4, 2008 09:39

Thanks Daniel Its been grea
 
Thanks Daniel
Its been great help for me.
I have done some editing in postChannel utility and I have considered UPrime2Mean and R to be same symmmTensor and replaced R with UPrime2Mean and recompiled the utility. Is it a right thing to do?
see the readField.c and collapse.H file
http://www.cfd-online.com/OpenFOAM_D...s/mime_txt.gif readFields.H
http://www.cfd-online.com/OpenFOAM_D...s/mime_txt.gif collapse.H

nishant_hull September 4, 2008 09:49

Another question is: Is the
 
Another question is:
Is there any recent modification in channelIndex.C file? I can see that you suggested some modification in that file.
Actually, I am getting some ambiguous result when I compare the results for channel flow using postChannel utility (modified by me) and sample utility available with foam. See the attached picture.. the Umean's x-component is plotted where green line is plotted using postchannel utility and red line is plotted with sample utility.
http://www.cfd-online.com/OpenFOAM_D...s/126/9051.png

My channelIndex.C file is as follows:
http://www.cfd-online.com/OpenFOAM_D...s/mime_txt.gif channelIndex.C

lakeat September 4, 2008 12:38

Hi Nishant! 1) Compare vers
 
Hi Nishant!

1) Compare version 1.4.1 and 1.5, note their differences, then you will see what R means
2) Don't use UPrime2Mean, but use R to get what you want.
3) Include "calculateAverages.H" from OF-1.4.1 so that you can get R field.

Cheers,
Daniel

nishant_hull September 4, 2008 12:54

hi Daniel There is define
 
hi Daniel

There is definetely some problem here. because I have a openfoam-1.4.1 installed on other system and when I plot the Uf with R field then I found that the same difference in Umean's x-component and Uf (which is also umean's x -component defined in postChannel/collapse.H).
What could be the reason for this difference? Why there is different plot for same variable plot? I might be wrong somewhere in my assumption.

Can anybody explain?

Nishant

lakeat September 4, 2008 21:32

Sorry for my last post, I exam
 
Sorry for my last post, I examined the case and found (please correct me if I am wrong),

1st, Huh... Don't know how did you get that graph!

2nd, the following are their relationship.
1.4.1 - R field - uniform/nAveragingSteps.raw
1.5 - UPrime2Mean field - uniform/fieldAveragingProperties

R=UPrime2Mean
They are both running average. So you don't nned to copy line from 1.4.1.

Remember everytime you do the reaveraging, don't forget to delete uniform dir.

Daniel

nishant_hull September 5, 2008 08:36

Hi Daniel, Thnaks for repl
 
Hi Daniel,
Thnaks for reply. Could you please try to compare your Umean.component(0) profile and Uf profile obtained from postchannel utility?? As per the code, they should be the same.
In case you are busy, then could you please send the tar dir of the channel395 case which provides you good Uf profile?

If any one can please see the problem mention on this forum dated September 04, 2008 and give me some suggestion then it would be of great help.

Nishant

nishant_hull September 5, 2008 09:18

I got the solution to my probl
 
I got the solution to my problem.
its here..
http://www.cfd-online.com/OpenFOAM_D...es/1/3532.html

lakeat September 5, 2008 11:41

Great! http://www.cfd-online.c
 
Great! http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

mattijs September 29, 2008 12:44

Hi Nishant, Daniel. I wrote
 
Hi Nishant, Daniel.

I wrote a postChannel that works on any layered and/or decomposed mesh and is simpler to use. Could you test it and let me know if it is correct? The sample dictionary should be self explanatory.

http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif postChannel.tgz

lakeat October 1, 2008 22:05

Thank you, Mattijs, I will try
 
Thank you, Mattijs, I will try. (I have just returned from my National Day Holiday).
Regards, Daniel

lakeat October 7, 2008 00:25

Hi, A nice day to you. I
 
Hi,

A nice day to you.

I am thinking how to plot a E(k)-k curve? Could anyone give me some hints?
Is E(k) the same as file "k" in each time directory? I found k value is very small.
And if k=2*pi*f/U, then how can I get f?

The following are k versus 2*Pi/mag(U). I guess it is sgs turbulent energy.
http://www.cfd-online.com/OpenFOAM_D...s/126/9472.jpg

How did you guys do that? In matlab or what?

Thank you.

Ps, the new postChannel works well, thank you Mattijs.

Regards, Daniel

zebu83 July 16, 2009 05:30

Hello everybody!

I download the new postChannel directory and i modify the postChannelDict, but it doesn't work!

old postChannelDict:
Quote:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Nx 40;
Ny (50);
Nz 30;

// ************************************************** *********************** //
new postChannelDict:
Quote:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Nx 40;
Ny
(
25
25
);
Nz 30;

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
I have to modify anithing else?

thank you

MT

zebu83 July 17, 2009 04:00

Sorry I'm new in OF... yesterday, after I post the comment, I download the
patched sources for OF-1.5, I compiled it and now "postChannel" works fine!

Now I have another problem: I can't visualise the result.
Usually I type
Quote:

foamToVTK
(because paraFoam doesn't works) then I visualise everythings with paraview.
I'm going to search into others post in the forum.

Best regard

MT

zebu83 July 22, 2009 09:20

I solved also this problem!

But now I'm going to ask everybody: how can I calculate Re_tau in Channel?

Thank you

MT

Fransje June 8, 2010 07:13

Large Eddy Simulation of turbulent channel flows.
 
Dear all,

I was searching for the paper "Large Eddy Simulation of Turbulent Channel Flows" by Furby, Gosman, Tabor and Weller mentioned various times in this post, and got slightly confused. Google Scholar seems to think that there is a 1997 version and 2007 version, and although I found the 1997 version, I only found a broken link for the 2007 paper. Would any body happen to know whether there is a difference between them?

All help is appreciated!

Regards,

Francois.


All times are GMT -4. The time now is 13:31.