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

turbulent jet simulation

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

Like Tree1Likes
  • 1 Post By msarkar

Reply
 
LinkBack Thread Tools Display Modes
Old   October 19, 2009, 07:27
Default turbulent jet simulation
  #1
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 9
antonio_ing is on a distinguished road
hi Foamers.

I'm trying to simulate the coaxial jet flow field for my PhD thesis. I tried to create the mesh and set up the case but there is some floating point error, maybe due to wrong boundary conditions. Does someone have some simulation of an axisymmetric round jet or even suggest me the boundary conditions?

thanks in advance
antonio_ing is offline   Reply With Quote

Old   April 1, 2010, 11:58
Default
  #2
Senior Member
 
Join Date: Dec 2009
Posts: 112
Rep Power: 8
heavy_user is on a distinguished road
Hi Antonio,

did you get any further with your problems?
I am having a similar issue right now..

regards!
heavy_user is offline   Reply With Quote

Old   April 1, 2010, 13:45
Default
  #3
New Member
 
Nat Trask
Join Date: Mar 2009
Location: Providence, RI
Posts: 18
Rep Power: 9
natrask is on a distinguished road
Hi Antonio,

I'm actually just finishing my thesis work writing a solver to predict turbulent atomization of a similar type of injector (gas-centered swirl coaxial). You can find a paper describing our work here (http://oai.dtic.mil/oai/oai?verb=get...fier=ADA513613). We did a number of validation cases, one of which was a coaxial fuel injector, and another being a diesel plain orifice injection.

The problem with predicting these guys in OpenFOAM is that the only appropriate solver that's released is interFoam (or one of the other interface tracking variations) which require mesh resolution of the scale of the droplet size (which is prohibitively small for the Weber numbers we're usually looking at). You could try a Lagrangian method (like the one that's set up in dieselFoam) but that won't account for the volume occupied by the liquid phase. For coax injectors, the liquid core profile can significantly alter the recirculation zone that forms between the two phases.

I'd be happy to share my code once I get it polished up (maybe by the end of April?).
natrask is offline   Reply With Quote

Old   April 6, 2010, 04:57
Default
  #4
Senior Member
 
Join Date: Dec 2009
Posts: 112
Rep Power: 8
heavy_user is on a distinguished road
Hi natrask,

nice work! Your results seem to be pretty close to the "real world" !

Would you mind posting a setup of your case?
With BC and stuff?

regards!
heavy_user is offline   Reply With Quote

Old   May 10, 2010, 08:07
Default
  #5
JMP
New Member
 
Josť M. Pastor
Join Date: May 2009
Location: Valencia (Spain)
Posts: 3
Rep Power: 9
JMP is on a distinguished road
Hi Nat,

Your approach and results seems quite interesting. We've got some experience with lagrangain methods, both on StarCD and OpenFOAM, but your approach could solve some issues when applying to dense sprays. Indeed I'm quite interested on coding this type (sigma-Y or 'ELSA') spray models into OpenFOAM. If you could share some details about your implementation or any additional publication it would be quite helpful.
Thanks

Regards

Jose
JMP is offline   Reply With Quote

Old   September 16, 2010, 02:31
Default
  #6
Member
 
MSarkar
Join Date: Dec 2009
Posts: 99
Rep Power: 8
msarkar is on a distinguished road
Hi All,

I am trying to simulate a heated air jet dispersion in an open area. The heated air is coming out from a pipe leak ( leak diameter 1cm) with a velocity of 135 m/s and a temperature of 100 C. In this case, the Mach-number is around 0.4 which is mildly compressible flow. So, I was using rhoSimpleFoam and rhoPImpleFoam. I tried several boundary conditions and mesh sizes but with out any success. The results do not look reasonable. As initial case I am considering a 2D rectangular geometry. Could you please suggest me which solver I should use?

The mesh and boundary file and one set of used boundary conditions are attached here. If you have any suggestions regarding boundary conditions, Please let me know. Any help would be very much appreciated.

top:
alphat: type calculated

epsilon:
type inletOutlet;
inletValue uniform 0.01;
value uniform 0.01;
k:
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}

mut:
type calculated;

p:
{
type fixedValue;
value uniform 100000;
}

T:
{
type zeroGradient;
value uniform 333;
}

U:
{
type pressureNormalInletOutletVelocity;
//type slip;
phi phi;
rho rho;
value uniform (0 0 0);
}

bottom:
alphat: type calculated

epsilon:
type inletOutlet;
inletValue uniform 0.01;
value uniform 0.01;

k:
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}

mut:
type calculated;

p:
{
type fixedValue;
value uniform 100000;
}

T:
{
type zeroGradient;
value uniform 333;
}

U:
{
type pressureNormalInletOutletVelocity;
//type slip;
phi phi;
rho rho;
value uniform (0 0 0);
}


inlet:
alphat: type calculated

epsilon:
{
type compressible::turbulentMixingLengthDissipationRate Inlet;
mixingLength 0.0005; // 0.01*0.05
value uniform 0.01;
}

k:
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.1;
}

mut:
type calculated;

p:
{
type zeroGradient;
value uniform 100000;
}

T:
{
type fixedValue;
value uniform 373;
}

U:
{
type fixedValue;
value uniform (135.74 0 0);
}


outlet:
alphat: type calculated

epsilon:
{
type inletOutlet;
inletValue uniform 0.01;
value uniform 0.01;
}

k:
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}

mut:
type calculated;

p:
{
type fixedValue;
value uniform 100000;
}

T:
{
type zeroGradient;
value uniform 333;
}

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

lowerleft:
alphat:
type alphatWallFunction;
Prt 0.85;
value uniform 0;

epsilon:
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}

k:
{
type compressible::kqRWallFunction;
value uniform 0.1;
}

mut:
{
type mutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}

p:
{
type fixedValue;
value uniform 100000;
}

T:
{
type zeroGradient;
value uniform 333;
}

U:
{
type fixedValue;
value uniform (0 0 0);
}


upperleft:
alphat:
type alphatWallFunction;
Prt 0.85;
value uniform 0;

epsilon:
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}

k:
{
type compressible::kqRWallFunction;
value uniform 0.1;
}

mut:
{
type mutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}

p:
{
type fixedValue;
value uniform 100000;
}

T:
{
type zeroGradient;
value uniform 333;
}

U:
{
type fixedValue;
value uniform (0 0 0);
}


frontAndBack:

For all variables it is empty.
Attached Images
File Type: jpg mesh.jpg (51.8 KB, 129 views)
Attached Files
File Type: txt boundary.txt (1.6 KB, 66 views)
samurai_01 likes this.
msarkar 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
LES simulation of sonic jet Jianping CFX 4 January 24, 2005 23:01
Modelling a turbulent jet and k-epsilon constants Ant CD-adapco 3 January 24, 2005 16:56
round jet simulation Robert Havercroft FLUENT 2 April 16, 2004 21:51
Turbulent Jet Ulf E CD-adapco 9 February 13, 2001 09:25
air jet simulation Yue zou Main CFD Forum 0 September 28, 1998 08:30


All times are GMT -4. The time now is 21:17.