CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Calculation of integral boundary values (https://www.cfd-online.com/Forums/openfoam-solving/58109-calculation-integral-boundary-values.html)

mattijs March 17, 2005 09:10

This is the fix. http://ww
 
This is the fix.

http://www.cfd-online.com/OpenFOAM_D...s/mime_txt.gif SortableList.C

ali April 19, 2005 12:25

The algorithm posted by Jarrod
 
The algorithm posted by Jarrod a few threads above (for evaluating flux at the patches such as inlet etc) works fine. But what if I have a channel of say length L=1 (from x=0 to x=1) and I want to calculate flux through the middle point (x=0.5). How do I introduce a face at x=0.5 which is not a boundary patch? I found the velocities at the middle to be more accurate, so I thought calculating flux at the middle would be better. inlet flux changes its sign while it shouldn't be like that. I guess something may be wrong there.

cedric_duprat January 14, 2009 11:29

Dear all, Keeping the idea
 
Dear all,

Keeping the idea of Ali, I need to calculate a mass flow average throught a plane in my geometry.
I'll like also to define this plane (coordinate, ect...) after my calculation because I can't guess the result.
I use these few lines to create a plane inside the geometry where I can get every value I need.
// Plane 1
point pnt(0,0,0);
vector direction(0,0,1);
plane pl1(pnt,direction);
cuttingPlane cutPlane1(mesh,pl1);
const labelList& cutCells1 = cutPlane1.cells();
word setName("someCells");
cellSet currentSet1(mesh, setName, cutCells1);

// Create mesh subsetting engine
fvMeshSubset subsetter1(mesh);
label patchI = -1;
subsetter1.setLargeCellSubset(currentSet1, patchI, true);

But, how can I use them to define a patch like boundary patch which wouldn't be at the boundary of my geometry ?
So I could use previous piece of code to get different value like mass flow average or pressure average, ect ...

Thank you for your help,

Cedric

mabinty September 9, 2009 07:44

dear all!!

added the posted codes to chtMultiRegionFoam to get the net mass flux in the fluid region(s). it compiles but giving me the following warnings:

************************************************** *******************************
fluid/setRegionFluidFields.H: In function ‘int main(int, char**)’:
fluid/setRegionFluidFields.H:1: warning: unused variable ‘mesh’
fluid/setRegionFluidFields.H:5: warning: unused variable ‘K’
fluid/setRegionFluidFields.H:6: warning: unused variable ‘U’
fluid/setRegionFluidFields.H:7: warning: unused variable ‘phi’
fluid/setRegionFluidFields.H:8: warning: unused variable ‘g’
fluid/setRegionFluidFields.H:10: warning: unused variable ‘turb’
fluid/setRegionFluidFields.H:11: warning: unused variable ‘DpDt’
fluid/setRegionFluidFields.H:14: warning: unused variable ‘psi’
fluid/setRegionFluidFields.H:15: warning: unused variable ‘h
************************************************** *******************************

but the code is running. should I stil care about the warnings??

greatly appreciate your comments!

aram

sankarv April 21, 2010 14:10

Is this mass flux or mass flow rate ?
 
Hello,

I have a question about the code you have posted here.
Does it calculate Mass Flux or Mass flow rate ?

Thanks
Vaidya


Quote:

Originally Posted by mattijs (Post 202068)
Hello Jarrod,



we located found the problem in SortableList. It is because setSize when it grows the list does not adapt the list of indices. Until we fix it you can work around it by creating your list big enough and shrinking it (i.e. setSize with something smaller than the current size)

Mattijs


immortality March 30, 2013 14:11

hi
I have a question about about where should we add the includes and their related headers?

haakon March 31, 2013 06:12

immortality: Have you heard about the phrase thread necromancy before?

I really think it is bad practice to post the same question twice on the same forum at the same time. Post it once, and wait for answer.

immortality March 31, 2013 07:06

i didn't received a related answer yet.

immortality March 31, 2013 07:12

i didn't received a related answer yet.


All times are GMT -4. The time now is 21:00.