[Ready33]

Hello,

I have a simulation with cylindrical inlet (rocket engine with solid propellant). Gas flow thru surface of cylinder, so I use surfaceNormalFixedValue

inlet

{

type surfaceNormalFixedValue;

refValue uniform -1.1269;

}

But there is a problem. Velocity of inlet gas depends on actual pressure

u = u0 . p^k

where u0 and k are constant. Do you have any idea how to do it?
Thaks.

[malv83]

Originally Posted by Ready33

Hello,

I have a simulation with cylindrical inlet (rocket engine with solid propellant). Gas flow thru surface of cylinder, so I use surfaceNormalFixedValue

inlet

{

type surfaceNormalFixedValue;

refValue uniform -1.1269;


}

But there is a problem. Velocity of inlet gas depends on actual pressure

u = u0 . p^k

where u0 and k are constant. Do you have any idea how to do it?
Thaks.

Hi,

Did you ever accomplish this?...Does any body knows how to do this?

The portion of the code in surfaceNormalFixedValue.C to modify would be:

Foam::newNormalVelFvPatchVectorField::
newNormalVelFvPatchVectorField
(
    const fvPatch& p,
    const DimensionedField<vector, volMesh>& iF,
    const dictionary& dict
)
:
    fixedValueFvPatchVectorField(p, iF, dict, false),
    refValue_("refValue", dict, p.size())
{
    fvPatchVectorField::operator=(refValue_*patch().nf());
}

I am Thinking that first we have to retrieve the pressure by (not sure):

const volScalarField& Pw = db().lookupObject<volScalarField>("p");

and then Multiply this by refValue_*patch().nf() which is a volVectorField. Is this possible? Any ideas?