[Tscar]

Hi,

I am trying to make velocity profile on inlet according to the followind expression:
u_z = -0.25 * (1/*mu) * dpl * [ ( ln(r/r1) / ln(r2/r1) ) * (r2*r2 - r1*r1) - (r*r - r1*r1) ]

r is radial distance from the center of global coordinate system and it is r*r = x*x + y*y

inlet
{
type codedFixedValue;
value $internalField;
redirectType velocityProfile;

code
#{
scalar mu=1.87e-5;
scalar dpl=50;
scalar r1=2e-3;
scalar r2=4e-3;
fixedValueFvPatchVectorField myPatch(*this);
forAll(this->patch().Cf(),i)
{
myPatch[i]=vector(0,0, -0.25*(1/mu)*dpl* ((((Foam::log(Foam::sqrt(Foam:ow(this->patch().Cf()[i].x(), 2) + Foam:ow(this->patch().Cf()[i].y(), 2)) / r1)) / Foam::log(r2/r1))) * (r2*r2 - r1*r1) - (Foam:ow(this->patch().Cf()[i].x(), 2) + Foam:ow(this->patch().Cf()[i].x(), 2) - r1*r1) ) );
}
operator==(myPatch);
#};
}

I checked the coded expression many times and I did not find any mistake. Yet the computed velocity profile is not right. Can you please help me?

[Tscar]

And now I am embarased. There should be

(Foamw(this->patch().Cf()[i].x(), 2) + Foamw(this->patch().Cf()[i].y(), 2) - r1*r1) ) );

Problem solved.

[gouravjee]

Originally Posted by Tscar

Hi,

I am trying to make velocity profile on inlet according to the followind expression:
u_z = -0.25 * (1/*mu) * dpl * [ ( ln(r/r1) / ln(r2/r1) ) * (r2*r2 - r1*r1) - (r*r - r1*r1) ]

r is radial distance from the center of global coordinate system and it is r*r = x*x + y*y

inlet
{
type codedFixedValue;
value $internalField;
redirectType velocityProfile;

code
#{
scalar mu=1.87e-5;
scalar dpl=50;
scalar r1=2e-3;
scalar r2=4e-3;
fixedValueFvPatchVectorField myPatch(*this);
forAll(this->patch().Cf(),i)
{
myPatch[i]=vector(0,0, -0.25*(1/mu)*dpl* ((((Foam::log(Foam::sqrt(Foam:ow(this->patch().Cf()[i].x(), 2) + Foam:ow(this->patch().Cf()[i].y(), 2)) / r1)) / Foam::log(r2/r1))) * (r2*r2 - r1*r1) - (Foam:ow(this->patch().Cf()[i].x(), 2) + Foam:ow(this->patch().Cf()[i].x(), 2) - r1*r1) ) );
}
operator==(myPatch);
#};
}

I checked the coded expression many times and I did not find any mistake. Yet the computed velocity profile is not right. Can you please help me?

hello,
i have been working on a solidification/melting problem in which i am supposed to give a gaussian profile as a input heat flux i have written a code for that.but it is not giving me desired profile.
can you tell me where is the problem??
[CODE]
input
{
type codedFixedValue;
value uniform 573;
redirectType temperatureProfile;

code
#{

fixedValueFvPatchScalarField myPatch(*this);
forAll(this->patch().Cf(),i)
{
myPatch[i]=2076000*exp(-Foam:ow((this->patch().Cf()[i].x()-122.5)/6.067,2));
}
operator==(myPatch);
#};
}

[\CODE]