[pg22]

I am using OF1.5 and rhopSonicFoam for a compressible external flow case. On the far field boundary I would like to apply a totalTemperature boundary condition:

topBoundary
{
type totalTemperature;
T0 uniform 300;
U U;
phi phi;
psi psi;
gamma 1.4;
value uniform 300;
}

I find that solver aborts with the following message:

request for volVectorField U from objectRegistry region0 failed
available objects of type volVectorField are

0
(
)

I have found that a totalPressure boundary condition (with fixedValue for T) runs for a short time before inflow leads to rapid divergence.

In the source files for totalPressure and totalTemperature I see that the vector field for velocity is accessed in a different way in each case:

For totalTemperatureFvPatchScalarField.C:

const fvPatchVectorField& Up =
patch().lookupPatchField<volVectorField, vector>(UName_);

For totalPressureFvPatchScalarField.C:

void totalPressureFvPatchScalarField::updateCoeffs()
{
updateCoeffs(patch().lookupPatchField<volVectorField, vector>(UName_));
}

Can anyone explain why the former BC fails and the latter
succeeds; and even suggest a solution?

Many thanks.