[kac24]

Hi all,

I am trying to use:
U
{
internalField uniform (0 0 2.5);

boundaryField
{
inlet
{
type codedFixedValue;
value uniform (0 0 1);

name parabolicVelocity;
code
#{
const vectorField& Cf = patch().Cf;
vectorField& field = *this;

const scalar R = 69
const scalar Vav = 2.5;

forAll(Cf, faceI)
{
const scalar x = Cf[faceI][0]
const scalar y = Cf[faceI][1]

field[faceI] = vector(0,0,2*Vav*(1-(pow(((pow(pow(x,2+pow(y,2),0.5))/R),2))
}

#};
}
but I am getting an error
Could not load "/alice-home/1/k/kac24/cavity_but_bendy/bubbler_18mar/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libparabolicVelocity_0dd38b0724eb3f538f626fd383ecd a3f8f71b9cd.so"





I could not see any documentation anywhere. The tutorial (yt) I based my code on is for version openfoam .org.




Does anyone have any working codedFixedValue script for version v2306 or similar that I could use as a template?

Thanks in advance

[fnohpat]

Originally Posted by kac24

Hi all,

I am trying to use:
U
{
internalField uniform (0 0 2.5);

boundaryField
{
inlet
{
type codedFixedValue;
value uniform (0 0 1);

name parabolicVelocity;
code
#{
const vectorField& Cf = patch().Cf;
vectorField& field = *this;

const scalar R = 69
const scalar Vav = 2.5;

forAll(Cf, faceI)
{
const scalar x = Cf[faceI][0]
const scalar y = Cf[faceI][1]

field[faceI] = vector(0,0,2*Vav*(1-(pow(((pow(pow(x,2+pow(y,2),0.5))/R),2))
}

#};
}
but I am getting an error
Could not load "/alice-home/1/k/kac24/cavity_but_bendy/bubbler_18mar/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libparabolicVelocity_0dd38b0724eb3f538f626fd383ecd a3f8f71b9cd.so"





I could not see any documentation anywhere. The tutorial (yt) I based my code on is for version openfoam .org.




Does anyone have any working codedFixedValue script for version v2306 or similar that I could use as a template?

Thanks in advance

Hi Karol!

This works perfectly in OpenFOAM 11.

inlet
{

type codedFixedValue;
value uniform (0 0 0);
name parabolicU; // user-defined name

code
#{
const fvPatch& boundaryPatch = this->patch();
const vectorField& Cf = boundaryPatch.Cf();

vectorField& v = *this;
scalar Uc = 0.264075; // centerline velocity, U_average
scalar R = 0.04; // inlet radius

forAll(Cf,faceI)
{
scalar z = Cf[faceI].z(), y = Cf[faceI].y(), rSq = z*z + y*y;
v[faceI] = vector(2*Uc*(1 - rSq/(R*R)),0, 0);
}
#};

I hope you find it useful.