[gbope7]

Hello,



I am using the timeVaryingMappedFixedValue (tVMFV) boundary condition to map heat fluxes onto a surface of an object. From there I want to observe the temperature evolution in that object.

I define the points and corresponding values of heat flux for each timestep in constant/boundaryData/<patchname>/points and constant/boundaryData/<patchname>/<time>/HF respectively.


Note that I define the volScalarField 'HF' in the .C file of a custom function similar to laplacianFoam. Beyond the definition

        Info<< "Reading Heat Flux HF from fcn file\n" << endl;
        volScalarField HF
        (
            IOobject
            (
                "HF",
                runTime.timeName(),
                mesh,
                IOobject::MUST_READ,
                IOobject::AUTO_WRITE
            ),
            mesh
        );

no operations are being performed on the volScalarField 'HF'.



Whenever I run my custom function I get the following error:

--> FOAM FATAL ERROR: (openfoam-2406)
cannot find file "/Users/promisea/OpenFOAM/ibukunOluwa/testCases/blockTest/2D_Trials/nstx132460/test2D-NSTX-tdep_timeVaryBC_10nodes/0.000159/HF"

    From virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::readStream(regIOobject &, const fileName &, const word &, const bool) const
    in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 629.

FOAM exiting

where test2D_NSTX.../ is case directory and 0.000159/ is first timestep after 0/.


Oddly enough, when I adjust the function so that instead of solving for temperature in the heat equation, it solves for heat flux

        while (simple.correctNonOrthogonal())
        {
            fvScalarMatrix TEqn
            (
             //    fvm::ddt(T) - fvm::laplacian(Diff, T)
             // ==
             //    fvOptions(T)
            // fvm::ddt(T) - Diff*fvm::laplacian(T)
            //  ==
                fvOptions(T)
                    fvm::ddt(HF) - fvm::laplacian(Diff, HF)
                 ==
                    fvOptions(HF)
            );

            fvOptions.constrain(TEqn);
            TEqn.solve();
            fvOptions.correct(T);
        }

        #include "write.H"

        runTime.printExecutionTime(Info);

the values from constant/boundaryData/<patchname>/ are mapped correctly. Can someone advise on why this is happening?



Thank you!