CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Jet simulation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 2, 2011, 06:39
Default Jet simulation
  #1
New Member
 
Ivan
Join Date: Sep 2010
Location: Russia , Moscow.
Posts: 13
Rep Power: 7
Zato_Ichi is on a distinguished road
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 !


Attached Images
File Type: png init.png (8.0 KB, 42 views)
Zato_Ichi is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Free jet simulation msarkar OpenFOAM Running, Solving & CFD 34 November 3, 2015 11:12
Jet fan and Tunnel simulation ahlo7 CFX 7 March 7, 2012 11:25
Help!simulation of free water jet Annie FLUENT 3 June 29, 2009 00:51
grid stretching function for jet simulation Bobby Main CFD Forum 0 May 29, 2005 13:46
air jet simulation Yue zou Main CFD Forum 0 September 28, 1998 08:30


All times are GMT -4. The time now is 15:59.