[Madeleine C]

I am attempting the write out the source terms in ther SRFPimpleFoam solver. Using a similar way than what is used to write out the Uabs:


createFields.H


// Create the Uabs
volVectorField Uabs
(
IOobject
(
"Uabs",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
SRF->Uabs()
);

// Create the Su
volVectorField Su
(
IOobject
(
"Su",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
SRF->Su()
);


This works fine for Uabs, but the problem comes in with Su and its separate terms Fcoriolis() + Fcentrifugal().


The difference is that Uabs is a is volVectorField while the source terms are internal fields, from SRFModel.H:



//- Return the coriolis force
tmp<volVectorField::Internal> Fcoriolis() const;

//- Return the centrifugal force
tmp<volVectorField::Internal> Fcentrifugal() const;

//- Source term component for momentum equation
tmp<volVectorField::Internal> Su() const;

//- Return velocity of SRF for complete mesh
tmp<volVectorField> U() const;

//- Return absolute velocity for complete mesh
tmp<volVectorField> Uabs() const;


How do I go about writing out internalFields?

[decanter]

Hello,

it has been a while but did you find any solution? I am looking for the exact same thing.

Cheers,
Ayca