[Rudiy]

Hallo, i am working at a chtMultiRegionSimpleFoam Solver to calculate the jouleHeating at a boundary with contact resistance.

For the calculations i would need to access the refValue of a certain boundary in the boundaryField.

I managed to get the average Value at the boundary field with:
const fvMesh &PoP1 = V_.db().parent().objectRegistry::lookupObject<fvMe sh>(region1);
const volScalarField &solidT1 = PoP1.thisDb().objectRegistry::lookupObject<volScal arField>(field);
double Pot1 = gAverage(solidT1.boundaryField()[patchID]);

but i would need the refValue of
solidT1.boundaryField()[patchID]

does someone know how to do this?

[olesen]

Originally Posted by Rudiy

Hallo, i am working at a chtMultiRegionSimpleFoam Solver to calculate the jouleHeating at a boundary with contact resistance.

For the calculations i would need to access the refValue of a certain boundary in the boundaryField.

I managed to get the average Value at the boundary field with:
const fvMesh &PoP1 = V_.db().parent().objectRegistry::lookupObject<fvMe sh>(region1);
const volScalarField &solidT1 = PoP1.thisDb().objectRegistry::lookupObject<volScal arField>(field);
double Pot1 = gAverage(solidT1.boundaryField()[patchID]);

but i would need the refValue of
solidT1.boundaryField()[patchID]

does someone know how to do this?

I would typically code something like this:

const auto& mesh = V_.db().time().lookupObject<fvMesh>(region1);
auto* solidT1 = mesh.getObjectPtr<volScalarField>(field);

double Pot1 = 0;


if (solidT1)
 {
    auto& bfld = solidT1->boundaryFieldRef()[patchID]);
    Pot1 = gAverage(bfld); 

    bfld = Zero;  // Or whatever ...
 }

The getObjectPtr method is very convenient since it combines a test-if-exists and a const_cast.

If you stick with most of your own code, you will need to do a const_cast on the volScalarField in order to use the boundaryFieldRef() access.

[Rudiy]

Thank you very much for your reply. I would have an other question.

If i use your code i get for the bfld variable with:

Info << bfld << " = bfld" << endl;

the following result:

type compressible::PotentialRadCoupledMixed;
refValue uniform 0.571153;
refGradient uniform 0;
valueFraction uniform 0.0161289;
source uniform 0;
value uniform 0.923078;
Tnbr myjouleHeatingSource:V;
sigmaLayer constant 66666;
sigma constant 200000;
kappaMethod solidThermo;
= bfld

now i would need the values of
refValue uniform 0.571153; and
sigmaLayer constant 66666;

is there a way to get these values ? (instead of using gAverage(bfld))

[olesen]

Originally Posted by Rudiy

Thank you very much for your reply. I would have an other question.

If i use your code i get for the bfld variable with:

Info << bfld << " = bfld" << endl;

...
now i would need the values of
refValue uniform 0.571153; and
sigmaLayer constant 66666;

is there a way to get these values ? (instead of using gAverage(bfld))

You will definitely need to read some documents there. The general boundary field is generally derived from a Field of Type, so that part is fairly straight forward to get at. The other bits such as "refValue uniform ..." and "sigmaLayer constant ..." are either still available as plain dictionary entries that you can scrape back together (if you have a coded access to a dictionary, for example), or if you know what the underlying boundary type actually is, you can dynamic_cast (or refCast for OpenFOAM) into the corresponding type. However, I don't know offhand if these elements are directly exposed as methods for any particular subclass. In either case, they are a PatchFunction1 or Function1, which means you will need to evaluate them with a given time value.