[fatirishman53]

I am in the process of developing a custom solver based off of the icoFoam solver. I would like to output a value, after each time step, based off the following:



Currently, I am judging quasi-steadiness of turbulent flow based off of the velocity of a probed location at the center of the mesh. However, being able to output this value would give me more information before having to do any post processing.

If anyone could help, I would appreciate it. Thanks.

[agustinvo]

Hi,

I think you can use a "probe function", and each time you write the data, a file will be created.

If you want to see it in the terminal, you just include a line like:

Info << "my variable" << my_variable << endl;

But it will give you all the values. Maybe with the probe function you can see the value also in the terminal

[fatirishman53]

Originally Posted by agustinvo

Hi,

I think you can use a "probe function", and each time you write the data, a file will be created.

If you want to see it in the terminal, you just include a line like:

Info << "my variable" << my_variable << endl;

But it will give you all the values. Maybe with the probe function you can see the value also in the terminal

Thanks for your reply! However, the equation does not represent a standard variable which can be output with the use of probes. I figured it would be easiest to add it to the solver output messages since I have already been adding features to it anyway.

The crux of my problem is that I am unsure of how to obtain the velocity field values, within the solver, that are required to solve this equation. I think that if I can figure that out, I will have my problem whipped .

Specifically, I need to be able to (if I am interpreting the equation correctly)
- obtain all instantaneous u and v velocity components, multiply them, and average the products for the first term. (Reynolds stress)

- obtain an average of the u component of velocity across each mesh xz-plane and find the derivative with respect to y, for the last term.

- Calculate the Reynolds number, which should be an easy task once I figure out how to do the above two tasks.

Sorry, I probably should have included these details before.

[agustinvo]

Originally Posted by fatirishman53

- obtain all instantaneous u and v velocity components, multiply them, and average the products for the first term. (Reynolds stress)

You could define a variable UU to to that, and by using "fieldAverage" you can get the mean value

Originally Posted by fatirishman53

- obtain an average of the u component of velocity across each mesh xz-plane and find the derivative with respect to y, for the last term.

As previously, but first you have to define the gradient of U

In all these cases, you have to select the correct term in the tensor you obtain.

[fatirishman53]

Originally Posted by agustinvo

In all these cases, you have to select the correct term in the tensor you obtain.

Thanks again! Apparently I have a lot of research to do

If you have any tips for implementation, I certainly welcome them!

[fatirishman53]

Ok, so I think I can get all of the information I need from

U.mesh.C()

by looping through all of the cell centers by index. I just need to figure out how to get access to the cell center coordinates so that I can average along horizontal planes. Could someone point me in the direction of the methods/members within U.mesh.C(), or tell me how I can return cell center coordinates?

EDIT:
Does "U.mesh.C()" return cell center coordinates or velocities? Perhaps "U.mesh.C().component()" returns cell center coordinate in a particular direction and "U.mesh.C().internalField()" returns the velocity values at the cell centers?

[fatirishman53]

Here is what I have figured out so far:

For cell center positions:

const vectorField centers = u.mesh().C().internalField();

For cell center velocities:

scalarField velocityX = U.internalField()

I still can't seem to figure out how to correlate the values in my "scalarField" with the y-value of the "vectorField". I know that "centers.component(1)" will return the y-position of all cell centers, but how do I apply that to "velocityX" so that I can retrieve all velocities at a particular y-position?

BTW, I am completely aware that I may be going about this completely wrong. Any help would be much appreciated!