|
[Sponsors] |
September 10, 2015, 11:09 |
incompressible multiphase - jet simulation
|
#1 |
Member
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11 |
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. |
|
September 14, 2015, 11:06 |
|
#2 |
Member
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11 |
No takers?
|
|
October 18, 2015, 20:22 |
|
#3 |
New Member
Join Date: Sep 2015
Location: Germany
Posts: 11
Rep Power: 11 |
Hello Shashank, I am working on a similar case; did you solve the problem?
|
|
October 19, 2015, 11:16 |
|
#4 |
Member
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11 |
Which one? The stability problem or the jet break up one?
|
|
October 19, 2015, 15:13 |
|
#5 |
New Member
Join Date: Sep 2015
Location: Germany
Posts: 11
Rep Power: 11 |
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 |
|
October 19, 2015, 15:28 |
|
#6 |
Member
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11 |
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.
|
|
May 16, 2022, 11:07 |
Question
|
#7 |
New Member
PRASANNA M
Join Date: Jan 2021
Posts: 1
Rep Power: 0 |
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. |
|
|
|
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 |