[jgs55]

Hi everyone,

I want to implement a new discretization scheme for the convective term in OpenFOAM. At every timestep, I have a coefficient file generated by python with the format:

... header

internalField   nonuniform List<vector> 
16384
(
(1.03544e-05 -1.06153e-05 1.03544e-05 1.03544e-05 ) // Cell 0
(3.06175e-05 -1.01041e-05 -1.01041e-05 -1.01041e-05) // Cell 1
(5.11404e-05 -9.46919e-06 -9.46919e-06 -9.46919e-06) // Cell 2
(7.20913e-05 -9.58538e-06 -9.58538e-06 -9.58538e-06 ) // ... 
(9.31207e-05 -9.87433e-06 -9.87433e-06 -9.87433e-06 )
(0.000115088 -1.02901e-05 -1.02901e-05 -1.02901e-05 )

I need to read the file "coefficients" like any other field (e.g. u,p,k,etc) but at every timestep and load it to OpenFOAM as a volumeField.

I guess I should implement this on the runTime.loop(), since createFields.H is only executed once at the beginning of the main loop, but I do not have a clear idea of how to do that.

Thank you in advance!

[wangsen992]

Have you tried including this as a volVectorField (just like U) but in the arguments for IOobject, use MUST_READ_IF_MODIFIED?

[jgs55]

Thank you Mr Wang for your reply, it was what I was looking for. For anyone with the same doubt, the solution was to modify the createFields.H of your solver, and add the following lines:

volVectorField coefficients
(
    IOobject
    (
        "coefficients",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED, //If object is registered its timestamp will be checked every timestep and possibly re-read.
        IOobject::NO_WRITE  //No automatic writing, but can be written explicitly
    ),
    mesh
);