[fisch]

Hi,

In my code to build a scalarTransportDyMFoam the mesh is changed with mesh.update() corretely.
Now i need a adjusted velocity field to calcuate phi out of it and to make it relative for my calculation of the PDE (scalarTransport Eqn.).

Do you know how to adjust the U field?

Thanks for any advice.
Rupert

[tomislav_maric]

This is from interDyMFoam.C:

        {
            // Calculate the relative velocity used to map the relative flux phi
            volVectorField Urel("Urel", U);

            if (mesh.moving())
            {
                Urel -= fvc::reconstruct(fvc::meshPhi(U));
            }

            // Do any mesh changes
            mesh.update();
        }

Does this help? You also have a correctPhi.H to include later. Check out the source.

[fisch]

Hi Tomislav,
thanks for the post.

In your given code you first calculate the U_relative and then you do the meshMotion, right?

What i want is to move the mesh and then i want to get the corresponding U_rel and U_abs.

Is this possible?
Thanks
Rupert

[tomislav_maric]

Originally Posted by fisch

Hi Tomislav,
thanks for the post.

In your given code you first calculate the U_relative and then you do the meshMotion, right?

What i want is to move the mesh and then i want to get the corresponding U_rel and U_abs.

Is this possible?
Thanks
Rupert

Hi Rupert,

why do you need to calculate the relative flux in this order, when it is actually prescribed with a given velocity field and the existing, non-moved mesh?

Tomislav

[fisch]

Ok, the whole story:

I want to build a scalarTransportFoam solver with moving meshes.

I think it works the following way:

After mesh.update() i have to figure out my current U_rel or U_abs. And calculate the corresponding phi field.
Then i have to make the phi relative as an explicit input of the equation.
Finally i want to store my U_abs field in my output files.

Right now the velocity field is not updated due to the mesh motion. It stays constant on the initial field.

Do you get my concern?
Thanks

[Thangam]

Hi,

I know this is a old thread. Just wondering if you managed to solve the issue.

http://www.cfd-online.com/Forums/ope...ic-meshes.html


Thanks
Thangam.