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

incompressible multiphase - jet simulation

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By smog

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 10, 2015, 11:09
Default incompressible multiphase - jet simulation
  #1
Member
 
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11
smog is on a distinguished road
Hello all, I am doing an incompressible, multiphase simulation of a water jet shooting into a cylinder containing stagnant air. I am using interFoam solver. I have refined the grid to have the minimum base cell dimension of 0.2 mm. I am having a really hard time :
1) Getting the solution to stay stable (there are severe spurious velocities in the internal field), whenever I use any 2nd order divergence schemes. I used Gauss linearUpwind grad (U), Gauss vanLeerV 1, Gauss limitedLinearV 1. None of them stop my simulation from blowing up.

2) On a coarser mesh, I can get some simulations to run without blowing up, but the water jet shoots through the domain without showing any signs of break-up. Just a straight beam-like behaviour.

I have come to believe I am doing something wrong with the IC/BCs. Here are my files, I have changed them recently, let me know if they make sense in the context of the problem.



k:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 1.5;

boundaryField
{
inflow
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 1.5;
//type fixedValue;
//value uniform 1.5; // calculated as per [1.5*pow(U*I,2)], U = mean flow velocity, I = turbulence intensity (assumed 5%)
}
outflow
{
type inletOutlet;
inletValue uniform 1.5;
value uniform 1.5;
//type zeroGradient;
}
wall_imported
{
type kqRWallFunction; // Average y+ = 30
value uniform 1.5;
}
inoutflow
{
type fixedValue;
value uniform 0; // No BL on this patch
}
"proc.*"
{
type processor;
value uniform 0;
}
}


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


nut:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0; // can be calculated = 1.4e-04, as per [sqrt(3/2)*(U*I*l)], I = turb intensity (assumed 10%), l = turb length scale

boundaryField
{
inflow
{
type calculated;
value uniform 0;
}
outflow
{
type calculated;
value uniform 0;
}
wall_imported
{
type nutkWallFunction; // Average y+ = 30
value uniform 0;
}
inoutflow
{
type fixedValue;
value uniform 0; // No BL on this patch
}
"proc.*"
{
type processor;
value uniform 0;
}
}


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

Omega:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 1e04;

boundaryField
{
inflow
{
type fixedValue;
value uniform 1e04; // clculated as per [sqrt(k)/l], k = turbulence kinetic energy, l = turbulence length scale = (0.038*D) for pipe flows
}
outflow
{
type inletOutlet;
inletValue uniform 1e04;
value uniform 1e04;
//type zeroGradient;
}
wall_imported
{
type omegaWallFunction; // Average y+ = 30
value uniform 1e04;
}
inoutflow
{
type fixedValue;
value uniform 0; // No BL on this patch
}
"proc.*"
{
type processor;
value uniform 0;
}
}


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

p_org:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p_org;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
inflow
{
type zeroGradient;
}
outflow
{
type fixedValue;
value uniform 0;
}
wall_imported
{
type zeroGradient;
}
inoutflow
{
type zeroGradient;
}
"proc.*"
{
type processor;
value uniform 0;
}
}


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

U:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 0 0);

boundaryField
{
inflow
{
type flowRateInletVelocity;
massFlowRate constant 0.1;
value uniform (0 0 0);
//type fixedValue;
//value uniform (-20 0 0);
}
outflow
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
inletValue uniform (0 0 0);
//type zeroGradient;
}
wall_imported
{
type fixedValue;
value uniform (0 0 0);
}
inoutflow
{
type fixedValue;
value uniform (0 0 0);
}
"proc.*"
{
type processor;
value uniform 0;
}
}


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






Sorry for the long post. Any comments will be highly appreciated.
smog is offline   Reply With Quote

Old   September 14, 2015, 11:06
Default
  #2
Member
 
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11
smog is on a distinguished road
No takers?
smog is offline   Reply With Quote

Old   October 18, 2015, 20:22
Default
  #3
New Member
 
Join Date: Sep 2015
Location: Germany
Posts: 11
Rep Power: 11
Frenk_T is on a distinguished road
Hello Shashank, I am working on a similar case; did you solve the problem?
Frenk_T is offline   Reply With Quote

Old   October 19, 2015, 11:16
Default
  #4
Member
 
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11
smog is on a distinguished road
Which one? The stability problem or the jet break up one?
smog is offline   Reply With Quote

Old   October 19, 2015, 15:13
Default
  #5
New Member
 
Join Date: Sep 2015
Location: Germany
Posts: 11
Rep Power: 11
Frenk_T is on a distinguished road
Actually both. I am using interFoam to simulate the formation of a Taylor bubble in a riser with an initial condition of riser full of (static) liquid.
You are using wall function: how did you calculate y+ for multiphase flow?
Thanks
Frenk_T is offline   Reply With Quote

Old   October 19, 2015, 15:28
Default
  #6
Member
 
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11
smog is on a distinguished road
That approach was not working out for me. I should have known better, I was using RANS modeling to simulate secondary break-up. So I started with LES, and with a finer mesh, had the jet to break. I tried on a 2D mesh first (though not applicable to 2D, LES gave me very promising first results). I cannot keep up with the size of the mesh anymore in 3D, so am looking at the DES (SA) models OpenFOAM has. That should let me run coarser for the bulk flow without needing me to heavily resolve the near wall region. With the DES models, I will start using the wall functions as the near wall treatment is RANS.
Frenk_T likes this.
smog is offline   Reply With Quote

Old   May 16, 2022, 11:07
Default Question
  #7
New Member
 
PRASANNA M
Join Date: Jan 2021
Posts: 1
Rep Power: 0
Prassy is on a distinguished road
Hey Shashank,

I know its been a while since you posted this, but I need some help regarding setting up a similar simulation. Do you use the setregions to inject the jet or an inlet, please let me know.

Thanks.
Prassy is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
simple Liquid jet simulation in an open domain canev.civelek FLUENT 3 January 23, 2024 00:20
Unbalanced energies in jet formation simulation npuh002 FLUENT 0 September 1, 2015 00:02
pressure equation in Jet simulation kumar OpenFOAM Running, Solving & CFD 0 June 9, 2010 13:18
Incompressible multiphase solver with laplaceelectrostatic solution coastal593 OpenFOAM Running, Solving & CFD 2 January 9, 2008 13:57
New to multiphase fluent simulation, any tutorials Camps FLUENT 0 September 9, 2006 12:41


All times are GMT -4. The time now is 22:53.