[Lewis Liang]

Dear Foamers,

I modify the governing equation by adding a source term in a solver, and the source term is a list of scalar which is solved by another solver. But the mesh geometry is kept the same by employing the same blockMeshDict. The source term here call sigma0 as shown below:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0.4";
    object      sigma0;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   nonuniform List<scalar> 
1125000
(
187900
187900
187900
187900
187900
187900
187900
187901
187901
187902
187902
187903
187903
187904
187905
187906
187907
187908
187909
187910
187911
187912
187913
187914
187916
187917
187918
187920
....
)
;
    }
}

And I try to add this list of data into the governing equation of new solver, the governing equation is shown:

fvScalarMatrix p2Eqn
(
fvm::ddt(p2) == fvm::laplacian(Cv, p2) - sigma0
);

And I found solution online, which is applying a dictionary to achieve this purpose, then I make a header file to read this list of data, but I got a error after compiling the solver.

    
 Info<< "Reading sigma0 field\n" << endl;

 IOdictionary sigmaD0
    (
        IOobject
        (
            "sigmaD0",    // name of the file
            runTime.constant(),
            mesh, 
            IOobject::MUST_READ,
            IOobject::NO_WRITE
        )
    );

    volScalarField  sigma0
   (
        sigmaD0.lookup("sigma0")
   );

Does anyone meet this kind of problem? Thanks in advance!

Cheers!

Lewis

[zhangyan]

IOobject sigma0_Header
(
    "sigma0",
    runTime.constant(),
    mesh,
    IOobject::MUST_READ
);          

volScalarField sigma0(sigma0_Header, mesh);

[Lewis Liang]

Hi!@zhangyan Thanks for your sharing. The solver is compiled successfully, and the list of data seem to be read in the governing equation so far. For this code:

volScalarField sigma0(sigma0_Header, mesh);

I have tried to figure out the original syntax in OpenFOAM class, but I did't find them from the volScalarField class. Could you please tell the location of the original OpenFOAM code?

Best regards,

Lewis