[Yongqiang66]

Hi Foamers,

I got an output array in Openfoam7. The values in the array correspond to each cell of the mesh.

I want to output the values in the array to a volscalarfield and write to files for visualization.

I have defined the field in "createfields.H". But I don't know how to assign values between array and volscalarfields.

Regards,
Yongqiang

[Yongqiang66]

In other words, is it allowed to assign values cell by cell in OpenFOAM?

Originally Posted by Yongqiang66

Hi Foamers,

I got an output array in Openfoam7. The values in the array correspond to each cell of the mesh.

I want to output the values in the array to a volscalarfield and write to files for visualization.

I have defined the field in "createfields.H". But I don't know how to assign values between array and volscalarfields.

Regards,
Yongqiang

[PositronCascade]

Are you trying to do something like this?

volScalarField p
(
    IOobject
    (
        "p",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    p_rgh + rho*gh
);

For instance, here the solver iterates p_rgh and then, using the expression

p = p_rgh + rho*gh

it converts p_rgh to p. So, you can define an expression like this at that line, if I understood your aim correctly.

[Yongqiang66]

Hi Positron,

Thank you for your reply.
I coupled OpenFOAM with another chemical simulator. They chemical simulator outputs an array. However, I aim to transfer the data in an array to an volscalarfield so that I can do visualization in ParaView.

It is to assign value for each of volscalarfield from an array. It seems that the following code is not what intends to do.

Thanks,
Yongqiang

Originally Posted by PositronCascade

Are you trying to do something like this?

volScalarField p
(
    IOobject
    (
        "p",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    p_rgh + rho*gh
);

For instance, here the solver iterates p_rgh and then, using the expression

p = p_rgh + rho*gh

it converts p_rgh to p. So, you can define an expression like this at that line, if I understood your aim correctly.

[ybapat]

You can use forAll loop.

[Yongqiang66]

Hi ybapat,

Thanks for your reply. I have solved this problem. The code can be run successfully.

Regards,
Yongqiang

Originally Posted by ybapat

You can use forAll loop.

[PositronCascade]

Originally Posted by Yongqiang66

Hi ybapat,

Thanks for your reply. I have solved this problem. The code can be run successfully.

Regards,
Yongqiang

Would you share the solution for the future possible visitors to this thread, so they can find a solution as well?

[Yongqiang66]

It's just as you suggested. Use forAll to loop over all the volscalarfield and = to directly assign values from the array . It is straightforward.

Originally Posted by PositronCascade

Would you share the solution for the future possible visitors to this thread, so they can find a solution as well?