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

TwoPhaseEulerFoam alpha "ripples"

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 28, 2012, 17:45
Default TwoPhaseEulerFoam alpha "ripples"
  #1
New Member
 
Michael Blatnik
Join Date: Jan 2012
Posts: 2
Rep Power: 0
MichaelB88 is on a distinguished road
Hello,

I am trying to simulate a biomass fluidized bed with Geldart type A particles (zeolite particles, d=100e-06m, rho=1750 kg/m^3). I am using an axi-symmetric assumption. For now, I am first doing the model with kinetic theory off, to avoid complications. With kinetic theory on, I was getting the same problems, plus more.

I am using a fluidizing velocity of 1.5 cm/s. The actual minimum fluidization velocity for the experimental bed I'm modeling is around 0.6 cm/s. There is also a center feeding tube, from which a gas velocity will come out (a drop tube for cellulose), but for now I've set that to 0 velocity.

My problem is ripples or waves in alpha where the fluidizing gas is expelled from the bottom. I believe it is the result of numerical errors.

I have tried several different types of numerical schemes, starting with the default limitedLinear, then trying vanLeer (crashed, made no difference), QUICK (crashes and gives relatively large negative alpha values), and upwind (1st order, solved the ripple-problem but gave unrealistic answers, no bubbles).

Any ideas? Please take a look at the photos for a better look at what I mean.

Here are the fvSchemes and fvSolutions files:
fvSchemes:
Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phia,Ua)    Gauss limitedLinearV 1;
    div(phib,Ub)    Gauss limitedLinearV 1;
    div(phib,k)     Gauss limitedLinear 1;
    div(phib,epsilon) Gauss limitedLinear 1;
    div(phi,alpha)  Gauss limitedLinear01 1;
    div(phir,alpha) Gauss limitedLinear01 1;
    div(phi,Theta)  Gauss limitedLinear 1;
    div(Rca)        Gauss linear;
    div(Rcb)        Gauss linear;
}

laplacianSchemes
{
    default         none;
    laplacian(nuEffa,Ua) Gauss linear corrected;
    laplacian(nuEffb,Ub) Gauss linear corrected;
    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
    laplacian(alphaPpMag,alpha) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}
fvSolution:
Code:
solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0;
        smoother        DIC;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    pFinal
    {
        $p;
        tolerance       1e-08;
        relTol          0;
    }

    "(k|epsilon)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-05;
        relTol          0.1;
    }

    "(k|epsilon)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-05;
        relTol          0;
    }

    alpha
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-10;
        relTol          0;
    }

    alphaFinal
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-10;
        relTol          0;
    }
}

PIMPLE
{
    nCorrectors     4;
    nNonOrthogonalCorrectors 0;
    nAlphaCorr      2;
    correctAlpha    yes;
    pRefCell        0;
    pRefValue       0;
}
If you need anymore information, let me know.
Some additional information:
Turbulence off, laminar
alphaCorrector yes

Thanks!

Michael
Attached Images
File Type: jpg screenshot.0005.jpg (15.7 KB, 76 views)
File Type: jpg screenshot.0011.jpg (18.8 KB, 68 views)
File Type: jpg screenshot.0015.jpg (21.8 KB, 73 views)
MichaelB88 is offline   Reply With Quote

Old   July 9, 2015, 04:20
Default Hi!
  #2
New Member
 
Su Choi
Join Date: Jul 2015
Posts: 3
Rep Power: 11
S choi is on a distinguished road
Hi.
I'm conducting a simulation similar with yours and I find your post written a few years ago.
But this post has no answer.

I'm wondering whether you solve this problem and the way you solve.
S choi is offline   Reply With Quote

Reply

Tags
fluidized bed, numerical error, numerical scheme, openfoam, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
TwoPhaseEulerFoam and Boundary conditions raagh77 OpenFOAM Running, Solving & CFD 99 February 6, 2018 18:31
Error, running a modified version of twoPhaseEulerFoam foamer OpenFOAM Running, Solving & CFD 0 March 11, 2012 16:17
about bounded alpha equation for bubbleFoam and twoPhaseEulerFoam kaifu OpenFOAM 2 May 6, 2011 05:33
mpi problem with twoPhaseEulerFoam piterw82 OpenFOAM Bugs 0 June 21, 2010 08:45
Determining alpha and beta for porous baffle Liaqat Khan Siemens 1 October 27, 2000 04:44


All times are GMT -4. The time now is 11:43.