[openfoammaofnepo]

Dear All,

I used the following to write out a new variable:

U_new=U;
U_new.write();

However, I found that in the U_new file, for the boundary conditions, their types are all "calculated", not the original ones, like:

wall
{
    type calculated;
    value uniform (0 0 0); 
}

Originally, it should be like:

wall
{
    type fixedValue;
    value uniform (0 0 0); 
}

Does anybody know how to write out the boundary type information correctly? Thank you in advance.

OFFO

[adhiraj]

Does this work?

volVectorField U_new
                    (
                     IOobject(
                              "U_new",
                              mesh.time().timeName(),
                              mesh,
                              IOobject::NO_READ,
                              IOobject::NO_WRITE
                             ),
                      mesh,   
                      dimensionedVector("zero",dimensionSet(0,1,-1,0,0,0,0),Foam::vector(0,0,0) ),
                      U.boundaryField().types()                       
                    )

Then you can assign like this:

U_new=U;
U_new.write();

[openfoammaofnepo]

It works well. Thank you very much!

[openfoammaofnepo]

Dear Adhiraj,

I found that when I use the above lines to write out the U_new, the values at the boundaries are still not correct, although the boundary type is correct now. for example:

about the turbulentInlet BC:

type turbulentInlet
value nonuniform<List>
(100e+100 100e+200 100d+200
...
...)

Here the values of the BCs are not reasonable, they are correctly given from U. Do you have any comments about this?Thank you.
OFFO

[openfoammaofnepo]

Hello,

it seems I still have some difficulty in writing out the new fields. The task is like: I have a old field, called U, and then did some treatment in the interalFields of U. I will try to write out the new internalField but old boundaryField into a new field, called U_c.

The code is as follows:

            volVectorField U_c
            (
                IOobject
                (
                    "U_c",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                mesh,
                dimensionedVector
                (
                 "U_c",
                 U.dimensions(),
                 vector(0,0,0)
                ),
                U.boundaryField().types()
            );

            vectorField& Ucell = U.internalField();
            vectorField& U_ccell = U_c.internalField();
            forAll(cellcentre,celli)
            {
                U_ccell[celli] = Ucell[celli];
                //.................some operations here. 
            }

            forAll(U.boundaryField(), patchi)
            {
                fvPatchVectorField& pU   = U.boundaryField()[patchi];
                fvPatchVectorField& pU_c = U_c.boundaryField()[patchi];

                forAll(pU, facei)
                {
                    pU_c[facei] = pU[facei];
                }
            }

            U_c.write();

The problem is:
The boundary type is correct written out. However, the values on the patches are not correctly output. The values on some patches are random numbers-very small or large numbers.

Could anyone tell me what I am missing in the above code?

[wyldckat]

Greetings to all!

@OFFO:

Originally Posted by openfoammaofnepo

Could anyone tell me what I am missing in the above code?

Have you been able to solve this issue?

If you haven't, are you certain you're providing the complete code that you are using for manipulating "U_c"?
I ask this because "U_c" is created with the option "IOobject::AUTO_WRITE", which means that at the end of an iteration, it will also write the field once again. This could mean that the "U_c" field is being altered between the time it was written with the explicit call to ".write()" and the end of the run time loop.

Best regards,
Bruno