[PatrickCFD]

OF7 (no groovyBC)

What I have:
1.png
- Curved mesh that rotates around a fixed point outside of the mesh (via dynamicMesh)
- Inlet boundary condition is a polynomial profile, which means that the inlet flow varies across one axis (in the picture the red vertices depict the highest velocities). I have integrated the profile like this:
inlet
{
type fixedProfile;
profile polynomial
(
((-1.9 0 -1.0) (1 0 1))
);
direction (1 0 0);
origin 0;
}

The numbers in the double brackets determine the angle and the velocity of the inlet flow.

The Problem:
2.png
I have a transient simulation and the angle of the flow stays the same while the angle of the inlet changes during the rotation of the mesh.

What I would like to have:
I would like for the inlet flow to always be perpendicular to the inlet surface during rotation. I couldn`t find a suitable boundary condition that covers a mesh/inlet rotation plus a polynomial profile. I guess a fixedProfile type is the only one that allows polynomial profiles to be created?!

Thanks in advance.

[PatrickCFD]

I tried going down the "codedFixedValue" route. I use the following code:

code
#{
const fvPatch& boundaryPatch = patch();
const vectorField& Cf = boundaryPatch.Cf();
const vectorField& n = boundaryPatch.nf();
vectorField& field = *this;

const scalar vmin = 25.597;
const scalar vmax = 35;
const scalar dmin = 24.5;
const scalar h = 9;
forAll(Cf, faceI)
{
const scalar x = Cf[faceI].x();
const scalar z = Cf[faceI].z()-33.5;
const scalar distance = pow(x*x+z*z,0.5);
const scalar faceVel = vmin + (distance - dmin) * (vmax - vmin)/h;
field[faceI] = faceVel*n;
}
#}

Apart from the last line of code everything seems to work, at least run through. But the last line (57) messes up the whole thing.


The displayed errors are:

processor0/0/U/boundaryField/inlet:57:17: no match for »operator=« (operand types are »Foam::Vector<double>« and »Foam::tmp<Foam::Field<Foam::Vector<double> > >«)
- processor0/0/U/boundaryField/inlet:57:17: candidates are:
In file included from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/vector.H:39:0,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/point.H:35,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/pointField.H:35,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/edge.H:40,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/edgeList.H:32,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/PrimitivePatch.H:56,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/primitivePatch.H:35,
from /opt/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/polyPatch.H:43,
from /opt/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude/fvPatch.H:39,
from /opt/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude/fvPatchField.H:47,
from /opt/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude/fixedValueFvPatchField.H:54,
from /opt/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude/fixedValueFvPatchFields.H:29,
from fixedValueFvPatchFieldTemplate.H:37,
from fixedValueFvPatchFieldTemplate.C:25:

I guess OpenFOAM wants me to use a vector or anything. But by choosing a vector I determine a direction which I don't want since I rotate the whole mesh (including the inlet). Maybe someone can make something of this. I'm thankful for the tiniest of responses.