[sega]

Hello.

I want to add a small procedure into the interFoam.
Its basically a solution for calculating the pressure inside a spherical bubble, but also limited to this case.

The "code" would read like this:

for all cells
********if rho > 0.99 * rho_phase1
****************nCell = nCell + 1
****************sumPressure = sumPressure + Pressure_cell
********end if
********Pressure = sumPressure / nCell
end for

How can I implement this code into C++?
I'm still not familar enought with C++ to do this on my own.

[sega]

Meanwhile I stumbled along the codes of interFoam and some postprocessing tools.

I tried to include some field-reading like in the interFoam createFields.H and put them inside the C-file of my little code:

surfaceTensionPressureDrop.C

The tool is compiling, but when I am running it there is some output from the inlcuded createPhi.H and from the selected transport model:


Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Time = 24.9
Number of cells = 316
Summarized pressure = 296.902
Pressure jump = 0.939563

Is there a way of suppressing these outputs?