[petrus]

https://www.dropbox.com/home?preview...%3A24%3A15.png

https://www.dropbox.com/home?preview...%3A29%3A07.png


Hi.

My first post, I hope give all necessary information to get a pleasant, good and quickly answer.
Thanks in advance.

So, trying simulate a jet injection inside a combustion chamber (without combustion of course). Nitrogen gas.
The mesh is axisymmetric, he nozzle of injector and the combustion chamber.
Split in 3 blocks.
blocks
(
hex ( 0 2 2 0 1 3 9 8) (20 1 8) simpleGrading (1.00 1.00 0.3) // block number0
hex ( 2 5 5 2 3 6 11 9) (100 1 8) simpleGrading (10.00 1.00 0.3) // block number1
hex ( 3 6 11 9 4 7 12 10) (100 1 25) simpleGrading (10.00 1.00 35.00) // block number2
);

It run correctly, but the animation is not as expected, better said it is very strange.

p;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 2.14e6;

boundaryField
{

inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 2.14e6;
}

wall
{
type zeroGradient;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}

}
---------------------------------------------------------------------------------------------
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 298;

boundaryField
{

inlet
{
type fixedValue;
value uniform 298;
}

outlet
{
type zeroGradient;
}

wall
{
type zeroGradient;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * /

object U;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{

inlet
{
type fixedValue;
value uniform (400 0 0);
}

outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}


wall
{
type fixedValue;
value uniform (0 0 0);
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

object sig;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 -1 0 0 0 0 0];

internalField uniform 0;

boundaryField
{

inlet
{
type fixedValue;
value uniform 0;
}

outlet
{
type zeroGradient;
}

wall
{
type fixedValue;
value uniform 0;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 14903;

boundaryField
{
inlet
{
type fixedValue;
value uniform 14903;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type compressible::epsilonWallFunction;
value uniform 14903;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 14903;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 37.5;

boundaryField
{
inlet
{
type fixedValue;
value uniform 37.5;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type compressible::kqRWallFunction;
value uniform 37.5;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object mut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type mutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 -1 0 0 0 0];

internalField uniform 1118000;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1118000;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type omegaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
value uniform 1118000;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

blockMesh is ok, checkMesh also is ok.

I think the problem is in the boundary conditions, but I dunno where.



Thanking again any suggestion.