Jet simulation
 

Greetings everyone !

I'm trying to simulate 2D air jet flow with OpenFoam 1.7 , and i've got some troubles. To solve this problem i picked sonicFoam. Initial geometry are attached. Some clarifications: there is a tank with compressed air under pressure 0f 150 atm, a pipe with diameter 0.0015 m, and reservoir. Pressure in pipe and in reservoir equals 1 atm.

And there is a problem: results obtained by sonicFoam are to smoothed. It looks like there is no shockwave in reservoir, and Mach disk is moving, while it should be fixed.

There is my BC for P, U, and T. Bounds of reservoir are under outlet condition and bounds of pipe and tank are walls.

P : outlet
{
type waveTransmissive;
value uniform 100000;
field p;
gamma 1.3;
phi phi;
rho rho;
psi psi;
fieldInf 100000;
lInf 1;
}
wall
{
type zeroGradient;
}

U : outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}

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

}
T : outflow
{
type zeroGradient;
}

wall
{
type zeroGradient;
}

Content of my fvSchemes file:

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phid,p) Gauss limitedLinear 0;
div(phi,e) Gauss limitedLinear 0;
div(phiU,p) Gauss limitedLinear 0;
div((muEff*dev2(grad(U).T()))) Gauss linear 0;
}

laplacianSchemes
{
default none;
laplacian((rho*(1|A(U))),p) Gauss linear uncorrected;
laplacian(muEff,U) Gauss linear uncorrected;
laplacian(alphaEff,e) Gauss linear uncorrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default uncorrected;
}

fluxRequired
{
default no;
p ;
}

i have red on the forum that rhoSonicFoam is better than you deal with large Mach numbers, but rhoSonicFoam using phiv instead phi for mass flux, and i failed to adopt my BC for phiv. Although i changed entries for phi like that:
phi phiv;
and added similar entry for inletOutlet U condition, it gives me an
Error on start solving:

request for surfaceScalarField UBlendingFactor from objectRegistry region0 failed
available objects of type surfaceScalarField are ...

I will be grateful for any help.

Best regards !