[ErenC]

Greetings to all!

I am trying to calculate non-dimensional fields. Such as:

Tn
(
IOobject
(
"Tn",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
(T-Tc)/(Th-Tc)
);


Un
(
IOobject
(
"Un",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
U/U0
);



of course, I defined Th and Tc, U0. This should actually work but it doesnt. It correctly calculates initial values but it doesn't calculate it for all iterations. In last iteration I still have same Tn values of the first iteration. I should be missing something.

Of course, I can calculate it in paraview but I am trying to gain some time.

Any help would be great,
Eren.

[clapointe]

If you put those somewhere like createFields, they will just be initialized (not updated). To update them, you'll need to put

 Tn = (T - Tc)/(Th - Tc)

inside your run loop.

Caelan

[ErenC]

Originally Posted by clapointe

If you put those somewhere like createFields, they will just be initialized (not updated). To update them, you'll need to put

 Tn = (T - Tc)/(Th - Tc)

inside your run loop.

Caelan

Thank you sir! It is working. How about runtime calculation of gradients? I defined GUn in createFields and added following line to code:

GUn= con * fvc::grad(Un)

(where con is a constant to make dx,dy non dimensional)


But not working.
error: no match for 'operator='

[clapointe]

Well it is hard to diagnose why it might not be working without seeing the code you've written or the corresponding error, but I'd guess that you haven't specified GUn as a volVectorField. Another option is to take the magnitude of the gradient to make it a scalar field.

Caelan

[maHein]

Is GUn of type volTensorField? Calculating the gradient of a vector (here: velocity) results in a tensor. Maybe that's the reason for the error.

[Daniel_Khazaei]

Originally Posted by maHein

Is GUn of type volTensorField? Calculating the gradient of a vector (here: velocity) results in a tensor. Maybe that's the reason for the error.

yes, fvc::grad(volVectorField) is a volTensorField!

[ErenC]

And it quickly solves my error! I changed field and it compiled. Thank you all!
Hope this thread will help another fellas to gain some time too!

[HPE]

OpenFOAM-Doxygen is your friend to deduce the types of objects, defined operators and so on in OpenFOAM software:

OpenFOAM-Doxygen