[Hale]

Hi Bernhard,

Thanks for your response. I am actually using the code calcMassFlow developed by Philippose, which is also mentioned earlier in this thread. I downloaded the file from http://www.cfd-online.com/Forums/ope...-boundary.html

This code seems to be the same as the one you suggest: http://openfoamwiki.net/index.php/Contrib_calcMassFlow

The thing is that I am very new to C++ and OpenFOAM and I don't know how to extract the information. I tried http://openfoamwiki.net/index.php/Co...e_calcMassFlow

functions
{
patchMassFlow
{
type patchExpression;
accumulations (
sum
);
patches (
inlet
outlet
);
expression "phi";
verbose true;
}
}

but it gives error for the patchExpression

--> FOAM FATAL ERROR:
Unknown function type patchExpression
Table of functionObjects is empty

From function functionObject::New(const word& name, const Time&, const dictionary&)
in file db/functionObjects/functionObject/functionObject.C at line 74.
FOAM exiting


What should I replace with the patchExpression?

And my last question is that will this function written above give me a file with the mass fluxes for the desired patches? Or should I add some sprict to it?

Thank you for your time

[Hale]

Hi again,

I am now trying to compile the swak4Foam given here http://openfoamwiki.net/index.php/Contrib/swak4Foam

but I cannot execute it by typing wmake all!! It gives the following error:

-bash-4.1$ wmake all
/appl/OpenFOAM/OpenFOAM-2.0.1/wmake/wmake: ./Allwmake: /bin/bash^M: bad interpreter: No such file or directory

What do I do wrong?

[gschaider]

Originally Posted by Hale

Hi again,

I am now trying to compile the swak4Foam given here http://openfoamwiki.net/index.php/Contrib/swak4Foam

but I cannot execute it by typing wmake all!! It gives the following error:

-bash-4.1$ wmake all
/appl/OpenFOAM/OpenFOAM-2.0.1/wmake/wmake: ./Allwmake: /bin/bash^M: bad interpreter: No such file or directory

What do I do wrong?

Nothing. Your Linux distro is "broken" (not really) ... it doesn't have a bash ... what kind of distro is this?

You can try replacing in the first line of Allwmake /bin/bash witg /bin/sh but in the past this introduced other problems. So the best fix (IMHO) is to install bash on that machine (I'm sure there is a package for compatibility)

[baila12]

Originally Posted by Hale

Hi,

Is there any body who knows how to write the fluxes given by calcMassFlow in a text file and not on the terminal where it normally does?

I want to plot the values on a graph but I don't know how to write them in a file!

NB. calcMassFlow is a ultility that gives the mass flux through patches.

Thank a lot

Hale,
If you want to write the output given by calcMassFlow to a text file, you can us the command calcMassFlow | tee log. Your output will then appear in the log file located in the case directory. This same command can be used for other solvers and utilities if you just replace calcMassFlow with the name of whichever solver you want information from.

[Hale]

Originally Posted by gschaider

Nothing. Your Linux distro is "broken" (not really) ... it doesn't have a bash ... what kind of distro is this?

You can try replacing in the first line of Allwmake /bin/bash witg /bin/sh but in the past this introduced other problems. So the best fix (IMHO) is to install bash on that machine (I'm sure there is a package for compatibility)

Thanks a lot for your response. As you said there is no bash installed on the distro. I actually don't know what kind of distro it is. I use linux through my university network, so I have no idea about their distro. I will contact them to hear more about it.

[Hale]

Originally Posted by baila12

Hale,
If you want to write the output given by calcMassFlow to a text file, you can us the command calcMassFlow | tee log. Your output will then appear in the log file located in the case directory. This same command can be used for other solvers and utilities if you just replace calcMassFlow with the name of whichever solver you want information from.

Thank you for the response baila12. I actually tried with calcMassFlow > log and it worked well. Now I'm having problem reading that file in matlab. This is another issue

[Hale]

Originally Posted by philippose

Hello Gabriela,

Good evening to you!

Sooo.... here are the steps you need to follow to compile and use calcMassFlow...

I am assuming that you have a working OpenFOAM installation, and that you are able to compile things using wmake....

1. Unzip the file "calcMassFlow.zip" to the folder: [home]/OpenFOAM/[username]-[OFVersion]/applications

2. Go into the created folder "calcMassFlow"

3. type "wmake"

4. The utility compiles, and installs the executable file in the folder:
[home]/OpenFOAM/[username]-[OFVersion]/applications/bin/linuxGCC4DOpt

The above four steps basically complete the compilation / installation of the utility.

To use it... you need to do the following:

1. In the folder "calcMassFlow" mentioned earlier, there is a sample "calcMassFlowDict" file.... copy this file into the "system" folder of your case.

2. Modify the entries under "patchNames" to the names of the patches through which you want to calculate the massflow.

3. Optionally, you can uncomment the "rho" line, and specify the density of the medium you simulated, in case you want the massflow in units "kg/s" instead of "m^3/s"

4. Go to the root folder of your case, and as usual, type:
calcMassFlow <root> <case> [-time t]

or

calcMassFlow <root> <case> [-latestTime]

or

calcMassFlow <root> <case>


That should be it....

Hope everything works as expected :-)!

Enjoy!

Philippose

Hi,

I have used this utility to calculate the mass flow at the boundaries and it works very well. The only problem I'm facing is that I can't find a way to sort all the information generated by calcMassFlow at each time step. The goal is to plot the fluxes at each time step and monitor the variation as function of time. I tried matlab without any success. Is there any smart way to do this?

Thanks a lot
Hale

[nostra]

Hi guys,

I try to do it with OF 2.2 and it works well. Thanks a lot

[Hale]

I'm interested in getting the flux at each cell of my boundaries (at each time step) and not the sum of the fluxes at each boundary, is there any solution for my problem? I know that the data for each cell is given in the phi file but since I'm running in parallel it is not possible to find the data for a specific boundary in the that file!

[gschaider]

Originally Posted by Hale

I'm interested in getting the flux at each cell of my boundaries (at each time step) and not the sum of the fluxes at each boundary, is there any solution for my problem? I know that the data for each cell is given in the phi file but since I'm running in parallel it is not possible to find the data for a specific boundary in the that file!

Yes it is. The thing is: if the boundary doesn't have faces on processorX then the file on processorX won't have that info. You've got to check the other processors. Or reconstruct.

There are ways of getting that data out during the run but you've got to be more specific on how you want to postprocess that (usually I'd say "single values are only interesting during development because they characterize the discretization but have got nothing to do with the physical system")

[Hale]

Thanks for your reply Bernhard.

I thought that the phi file only gives the flux of each cell without stating which boundary they belong to, but after I read your post I looked carefully and I found what I was looking for

But now the problem is that I only need the values that are in the middle of the boundary. Is there any way to check where each mesh cell is located?

[skyinventorbt]

Dear Foamers,

If you are simulating incompressible flows,

in Paraview there is an FILTER called "SURFACE FLOW"

Usage:
Create slice and apply surface flow.
Open Spread sheet view and find the value.

(Note: This is only velocity integrated over area, so multiply by rho to get mass flux)

-
KANNAN

[zqlhzx]

Hi philippose,
Thank you for upload your attach.It works well.It is easy to get the mass flow rate at the knowed boundry,e.g,inlet,outlet.But if I want to get the mass flow rate on any patch not only the boundry,how can I do?
For example,Can I give the command like "basePoint(0 0 0);normalVetor(0 0 1)"to define the patch what I want to calculate mass flow rate?
Thanks in advance!

[hy1112006]

Very useful utility, thank you very much!!

[hy1112006]

Originally Posted by zqlhzx

Hi philippose,
Thank you for upload your attach.It works well.It is easy to get the mass flow rate at the knowed boundry,e.g,inlet,outlet.But if I want to get the mass flow rate on any patch not only the boundry,how can I do?
For example,Can I give the command like "basePoint(0 0 0);normalVetor(0 0 1)"to define the patch what I want to calculate mass flow rate?
Thanks in advance!

Hello zqlhzx,

Have you make that work ? I also have same situation as yours.

If you did, could you please tell me how to do it?

Thanks,

Yi Han

[Adlak]

Hi everyone,
I want to calculate mixing efficiency through the formulae attached in a file.

Can anyone help me how to do that ???

[gschaider]

Originally Posted by Adlak

Hi everyone,
I want to calculate mixing efficiency through the formulae attached in a file.

Can anyone help me how to do that ???

A formula is useless if you don't explain what the symbols mean. I think this can be done with swakfoam but I refuse to guess things like where "mass flow of H2" is defined

[Adlak]

Originally Posted by gschaider

A formula is useless if you don't explain what the symbols mean. I think this can be done with swakfoam but I refuse to guess things like where "mass flow of H2" is defined

Hi,
mass flow of H2 is constant at fuel inlet. Mostly I want to know about the numerator part where I have to do integration. I have to integrate alpha*(density of mixture)*(mass fraction of h2)*(axial velocity) over a cross section in the combustor. I will do the same at many cross sections and plot efficiency graph along length of combustor.

[gschaider]

Originally Posted by Adlak

Hi,
mass flow of H2 is constant at fuel inlet. Mostly I want to know about the numerator part where I have to do integration. I have to integrate alpha*(density of mixture)*(mass fraction of h2)*(axial velocity) over a cross section in the combustor. I will do the same at many cross sections and plot efficiency graph along length of combustor.

Are the cross-sections aligned with cell-faces (and have you got faceSets or faceZones for that) or are you planing to use sampledSurfaces.

And what is the phi in the definition of alpha? Surly not the phi-field?

Anyway: is doable in swak4foam

[Adlak]

Originally Posted by gschaider

Are the cross-sections aligned with cell-faces (and have you got faceSets or faceZones for that) or are you planing to use sampledSurfaces.

And what is the phi in the definition of alpha? Surly not the phi-field?

Anyway: is doable in swak4foam

Hi,

Phi is local equivalence ratio. Yes, cross sections are aligned with cell faces. I dont have any other things.