[daven]

Hello together,


we recently use OpenFOAM in coupled multiphysical simulations with other solvers and need to recompute certain timesteps in each solver in order to use implicit coupling algorithms. Therefore, we reload in OpenFOAM internal fields (e.g. p, U) and reevaluate the boundaries by using the OpenFOAM function correctBoundaryConditions().

Reevaluating boundaries works fine for most of the boundary condition types, but it fails for some wall function types (e.g. epsilonWallFunction), where OpenFOAM raises the following error message:

--> FOAM FATAL ERROR: 

    request for volScalarField::Internal kEpsilon:G from objectRegistry region0 failed
    available objects of type volScalarField::Internal are



23
(
thermo:mu
thermo:psi
nut
epsilon_0
K
K_0
h
p_rgh_0
rho
k
h_0
p_rgh
dpdt
gh
alphat
rho_0
p
T
k_0
p_0
thermo:rho
epsilon
thermo:alpha
)

We couldn't find a workaround, since the wallFunction is only applied for specific patches and we would like to reload all other patches, where boundaries have other (unproblematic) boundary condition types.


If you are interested in a more extensive problem description/history or like to contribute on github: This is currently an open issue in precice/openfoam-adapter#4.


Any idea or hint is very welocme!


Best regards,
David

[mAlletto]

The Problem in your case is that the field G is only known inside the turbulence model class and non in the global object registry. In order to evaluate the wallfunvions from outside the turbulence model you have to make the required fields visible globally by inserting them into the mesh database

[daven]

Thanks for your reply.


Do you know how to access turbulence models in general? Inserting it in the base class should not be a problem. But accessing it results in my case in an -fpermissive compiler error.

[mAlletto]

You can get a constant reference to most fields see https://www.openfoam.com/documentati...enceModel.html. I believe also to the G field

[mAlletto]

A and in almost all solvers there is a reference to the turbulence model. You can find an example here https://github.com/OpenFOAM/OpenFOAM...createFields.H