[rdbisme]

Originally Posted by hooman.4028

I appreciate the quick response. It is a simple fluidized bed case (3D). The files are in the following link:

https://www.dropbox.com/s/2lrcoyzymx...files.rar?dl=0

As a general advice, it's better to post directly only the parts of the code of the dictionaries that are of interest, mainly because people are lazy and downloading, extracting and opening all the files is somehow time consuming. But if you post directly the boundary conditions dictionaries we can read them directly also from the smartphone .

In addition to that if you post the code directly in your reply it stays for future reference, even if you delete the archive from Dropbox.

Since I'm from the phone at the moment I would kindly ask to edit your post with these modifications so I can read them and try to see if I can help...

Sent by my Honor 8 using Tapatalk

[hooman.4028]

Originally Posted by tidusuper91

As a general advice, it's better to post directly only the parts of the code of the dictionaries that are of interest, mainly because people are lazy and downloading, extracting and opening all the files is somehow time consuming. But if you post directly the boundary conditions dictionaries we can read them directly also from the smartphone .

In addition to that if you post the code directly in your reply it stays for future reference, even if you delete the archive from Dropbox.

Since I'm from the phone at the moment I would kindly ask to edit your post with these modifications so I can read them and try to see if I can help...

Sent by my Honor 8 using Tapatalk

Makes perfect sense :

U

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

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

internalField   uniform (-1.112096 0 0);

boundaryField
{
    inlet
    {
        type               interstitialInletVelocity;
        inletVelocity      uniform (-1.112096 0 0);
        alpha              alpha.water;
        value              $internalField;
    }

    outlet
    {
        type               pressureInletOutletVelocity;
        phi                phi.water;
        value              $internalField;
    }
 top
    {
        type               noSlip;

    }
   
}

p_rgh

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

internalField       uniform 1e5;

boundaryField
{
    inlet
    {
        type            fixedFluxPressure;
        value           $internalField;
    }

    outlet
    {
        type            prghPressure;
        p               $internalField;
        value           $internalField;
    }

    top
    {
        type               fixedFluxPressure;
        value              $internalField;
    }
  
}

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

solution output

Courant Number mean: 5.34725e-08 max: 0.9516
Max Ur Courant Number = 0.625617
deltaT = 1.9149e-16
Time = 8.00643e-05

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver:  Solving for alpha.particles, Initial residual = 6.52947e-18, Final residual = 1.56184e-18, No Iterations 1
alpha.particles volume fraction = 0.0618866  Min(alpha1) = -4.53324e-36  Max(alpha1) = 0.600345
Constructing momentum equations
particles min/max T 600.004 - 2.19393e+19
water min/max T 300 - 1.4779e+19
GAMG:  Solving for p_rgh, Initial residual = 0.1844, Final residual = 0.0011326, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.00100111, Final residual = 6.8761e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 6.74233e-06, Final residual = 5.48156e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 5.48074e-08, Final residual = 8.52694e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 8.52694e-09, Final residual = 3.51475e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.51475e-09, Final residual = 1.34518e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.34518e-09, Final residual = 5.5909e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.5909e-10, Final residual = 2.14042e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.14042e-10, Final residual = 8.94611e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.94611e-11, Final residual = 3.44695e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.44695e-11, Final residual = 1.44658e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.132845, Final residual = 0.000674916, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000617022, Final residual = 3.97858e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.78338e-06, Final residual = 3.46792e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.46763e-08, Final residual = 5.48935e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.48935e-09, Final residual = 2.10436e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.10436e-09, Final residual = 8.82825e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.82825e-10, Final residual = 3.38082e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.38082e-10, Final residual = 1.42845e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.42845e-10, Final residual = 5.4936e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.4936e-11, Final residual = 2.32655e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.32655e-11, Final residual = 9.07378e-12, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.124431, Final residual = 0.000599917, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000556152, Final residual = 4.06027e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.89279e-06, Final residual = 3.58228e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.58196e-08, Final residual = 5.67329e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.67328e-09, Final residual = 2.17595e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.17595e-09, Final residual = 9.11311e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 9.11311e-10, Final residual = 3.48985e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.48985e-10, Final residual = 1.47163e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.47163e-10, Final residual = 5.6625e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.6625e-11, Final residual = 2.3944e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.3944e-11, Final residual = 9.34572e-12, No Iterations 1
PIMPLE: iteration 2
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver:  Solving for alpha.particles, Initial residual = 6.55151e-18, Final residual = 1.56407e-18, No Iterations 1
alpha.particles volume fraction = 0.0618866  Min(alpha1) = -4.53324e-36  Max(alpha1) = 0.600345
Constructing momentum equations
particles min/max T 600.004 - 2.41872e+19
water min/max T 300 - 1.62932e+19
GAMG:  Solving for p_rgh, Initial residual = 0.160617, Final residual = 0.000597597, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000526756, Final residual = 2.94328e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 2.81125e-06, Final residual = 2.54894e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 2.54879e-08, Final residual = 9.69721e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 9.69721e-09, Final residual = 4.01961e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.01961e-09, Final residual = 1.53481e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.53481e-09, Final residual = 6.44996e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 6.44996e-10, Final residual = 2.46642e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.46642e-10, Final residual = 1.04385e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.04385e-10, Final residual = 4.0082e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.0082e-11, Final residual = 1.69983e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.125642, Final residual = 0.00063384, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000583315, Final residual = 3.90957e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.73053e-06, Final residual = 3.42426e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.42398e-08, Final residual = 5.42277e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.42277e-09, Final residual = 2.07962e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.07962e-09, Final residual = 8.71991e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.71991e-10, Final residual = 3.33975e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.33975e-10, Final residual = 1.41033e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.41033e-10, Final residual = 5.42492e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.42492e-11, Final residual = 2.29642e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.29642e-11, Final residual = 8.95838e-12, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.119471, Final residual = 0.000576083, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000535629, Final residual = 3.9422e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.7856e-06, Final residual = 3.48448e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.48419e-08, Final residual = 5.51853e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.51853e-09, Final residual = 2.11671e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.11671e-09, Final residual = 8.86369e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.86369e-10, Final residual = 3.39446e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.39446e-10, Final residual = 1.43116e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.43116e-10, Final residual = 5.50713e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.50713e-11, Final residual = 2.3285e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.3285e-11, Final residual = 9.08688e-12, No Iterations 1
PIMPLE: iteration 3
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver:  Solving for alpha.particles, Initial residual = 6.46559e-18, Final residual = 1.55648e-18, No Iterations 1
alpha.particles volume fraction = 0.0618866  Min(alpha1) = -4.53324e-36  Max(alpha1) = 0.600345
Constructing momentum equations
particles min/max T 600.004 - 2.76185e+19
water min/max T 300 - 1.86046e+19
GAMG:  Solving for p_rgh, Initial residual = 0.146758, Final residual = 0.000597587, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000532726, Final residual = 3.22656e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.08401e-06, Final residual = 2.81869e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 2.8185e-08, Final residual = 4.45865e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 4.45865e-09, Final residual = 1.70722e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.70722e-09, Final residual = 7.1702e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 7.1702e-10, Final residual = 2.7444e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.7444e-10, Final residual = 1.16091e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.16091e-10, Final residual = 4.46259e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.46259e-11, Final residual = 1.89177e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.89177e-11, Final residual = 7.37524e-12, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.120119, Final residual = 0.000604078, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.00055852, Final residual = 3.85352e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.68598e-06, Final residual = 3.38648e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.3862e-08, Final residual = 5.36435e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.36435e-09, Final residual = 2.05778e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.05778e-09, Final residual = 8.62458e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.62458e-10, Final residual = 3.30356e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.30356e-10, Final residual = 1.39434e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.39434e-10, Final residual = 5.36491e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.36491e-11, Final residual = 2.27011e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.27011e-11, Final residual = 8.859e-12, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.115279, Final residual = 0.000557476, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.000519602, Final residual = 3.86166e-06, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.713e-06, Final residual = 3.41831e-08, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 3.41803e-08, Final residual = 5.41388e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 5.41388e-09, Final residual = 2.07679e-09, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.07679e-09, Final residual = 8.69468e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.69468e-10, Final residual = 3.33002e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.33002e-10, Final residual = 1.40362e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.40362e-10, Final residual = 5.40171e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.40171e-11, Final residual = 2.28344e-11, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 2.28344e-11, Final residual = 8.91275e-12, No Iterations 1
smoothSolver:  Solving for Theta.particles, Initial residual = 0.354633, Final residual = 1.99366e-17, No Iterations 1
ExecutionTime = 352.74 s

fvSolution

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

solvers
{
    "alpha.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 2;

       //smoothLimiter   0.1;

        implicitPhasePressure yes;
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-9;
        relTol          0;
        minIter         1;
    }

    p_rgh
    {
        solver          GAMG;
        smoother        DIC;
        tolerance       1e-8;
        relTol          0.01;
        minIter         1;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "U.*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-5;
        relTol          0;
        minIter         1;
    }

    "(h|e).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-6;
        relTol          0;
        minIter         0;
		maxIter		0;
    }

    "Theta.*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-6;
        relTol          0;
        minIter         1;
    }

    "(k|epsilon).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-5;
        relTol          0;
        minIter         1;
    }
}

PIMPLE
{
    nOuterCorrectors 3;
    nCorrectors      3;
    nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
    equations
    {
        ".*"            0.5;
    }
}


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

fvSchemes

ddtSchemes
{
    default     Euler;
}

gradSchemes
{
    default     Gauss linear;
    /*div(phi,U)      bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;*/
}

divSchemes
{
    default                         none;

    "div\(phi,alpha.*\)"            Gauss vanLeer;
    "div\(phir,alpha.*\)"           Gauss vanLeer;

    "div\(alphaRhoPhi.*,U.*\)"      Gauss limitedLinearV 1;
    "div\(phi.*,U.*\)"              Gauss limitedLinearV 1;

    "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss limitedLinear 1;
    "div\(alphaRhoPhi.*,K.*\)"      Gauss limitedLinear 1;
    "div\(alphaPhi.*,p\)"           Gauss limitedLinear 1;

    div(alphaRhoPhi.particles,Theta.particles) Gauss limitedLinear 1;

    "div\(alphaRhoPhi.*,(k|epsilon).*\)"  Gauss limitedLinear 1;


    /*"div\(alphaRhoPhi.*,U.*\)"      bounded Gauss upwind;
    "div\(phi.*,U.*\)"              bounded Gauss upwind;

    "div\(alphaRhoPhi.*,(h|e).*\)"  bounded Gauss upwind;
    "div\(alphaRhoPhi.*,K.*\)"      bounded Gauss upwind;
    "div\(alphaPhi.*,p\)"           bounded Gauss upwind;

    div(alphaRhoPhi.particles,Theta.particles) bounded Gauss upwind;

    "div\(alphaRhoPhi.*,(k|epsilon).*\)"  bounded Gauss upwind;*/

    div((((alpha.water*thermo:rho.water)*nuEff.water)*dev2(T(grad(U.water))))) Gauss linear;

    div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I)))  Gauss linear;
}

laplacianSchemes
{
    default     Gauss linear uncorrected;
    bounded     Gauss linear uncorrected;


    /*default Gauss linear limited 0.5;
    bounded Gauss linear limited 0.5;*/


}

interpolationSchemes
{
    default     linear;
}

snGradSchemes
{
    default     uncorrected;
    bounded     uncorrected;
}

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

checkMesh

Checking geometry...
    Overall domain bounding box (-0.05560247 -0.05558248 0.0003519632) (0.09964492 0.05564084 0.2476475)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (5.655138e-15 -5.589869e-15 2.068406e-15) OK.
    Max cell openness = 3.333135e-16 OK.
    Max aspect ratio = 12.30445 OK.
    Minimum face area = 4.666168e-08. Maximum face area = 2.189186e-06.  Face area magnitudes OK.
    Min volume = 4.40432e-11. Max volume = 1.523349e-09.  Total volume = 0.002161009.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.98925 average: 2.613793
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.380155 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

Let me know if you need more information.

[hooman.4028]

Originally Posted by tidusuper91

As a general advice, it's better to post directly only the parts of the code of the dictionaries that are of interest, mainly because people are lazy and downloading, extracting and opening all the files is somehow time consuming. But if you post directly the boundary conditions dictionaries we can read them directly also from the smartphone .

In addition to that if you post the code directly in your reply it stays for future reference, even if you delete the archive from Dropbox.

Since I'm from the phone at the moment I would kindly ask to edit your post with these modifications so I can read them and try to see if I can help...

Sent by my Honor 8 using Tapatalk

The other files:

phase properties

type    heatAndMomentumTransferTwoPhaseSystem;

phases (particles water);

particles
{
    type          purePhaseModel;

    diameterModel constant;
    constantCoeffs
    {
        d               3.72e-6;
    }

    alphaMax        0.62;
    residualAlpha   1e-6;
}

water
{
    type          purePhaseModel;

    diameterModel constant;
    constantCoeffs
    {
        d               1;
    }

    residualAlpha   0;
}

blending
{
    default
    {
        type            none;
        continuousPhase water;
    }
}

surfaceTension
(
    (water and particles)
    {
        type            constant;
        sigma           0;
    }
);

aspectRatio
(
);

drag
(
    (particles in water)
    {
        type            GidaspowErgunWenYu;
        residualRe      1e-3;
        swarmCorrection
        {
            type        none;
        }
    }
);

virtualMass
(
);

heatTransfer
(
    (particles in water)
    {
        type            RanzMarshall;
        residualAlpha   1e-4;
    }
);

lift
(
);

wallLubrication
(
);

turbulentDispersion
(
);

// Minimum allowable pressure
pMin            10000;

alpha water

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }


}

[rdbisme]

Originally Posted by hooman.4028

The other files:

phase properties

type    heatAndMomentumTransferTwoPhaseSystem;

phases (particles water);

particles
{
    type          purePhaseModel;

    diameterModel constant;
    constantCoeffs
    {
        d               3.72e-6;
    }

    alphaMax        0.62;
    residualAlpha   1e-6;
}

water
{
    type          purePhaseModel;

    diameterModel constant;
    constantCoeffs
    {
        d               1;
    }

    residualAlpha   0;
}

blending
{
    default
    {
        type            none;
        continuousPhase water;
    }
}

surfaceTension
(
    (water and particles)
    {
        type            constant;
        sigma           0;
    }
);

aspectRatio
(
);

drag
(
    (particles in water)
    {
        type            GidaspowErgunWenYu;
        residualRe      1e-3;
        swarmCorrection
        {
            type        none;
        }
    }
);

virtualMass
(
);

heatTransfer
(
    (particles in water)
    {
        type            RanzMarshall;
        residualAlpha   1e-4;
    }
);

lift
(
);

wallLubrication
(
);

turbulentDispersion
(
);

// Minimum allowable pressure
pMin            10000;

alpha water

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }


}

For what I understand, you're using reactingTwoPhaseEulerFoam and not twoPhaseEulerFoam. In order to not solve the energy equation you need to use a isothermalPhaseModel as phase model (and not, instead, using the hacky workaround I suggested before for twoPhaseEulerFoam) like pureIsothermalPhaseModel.

[hooman.4028]

Originally Posted by tidusuper91

For what I understand, you're using reactingTwoPhaseEulerFoam and not twoPhaseEulerFoam. In order to not solve the energy equation you need to use a isothermalPhaseModel as phase model (and not, instead, using the hacky workaround I suggested before for twoPhaseEulerFoam) like pureIsothermalPhaseModel.

Thanks Ruben! The energy is not solved anymore, however, my solution blows up time step constantly decreases and then suddenly a jump in courant number.

I have tried everything... Low courant number ... Low relaxation factor... Upwind methods ... confusing ?

PS: The reason I switched to reactingTwoPhaseEulerFoam was its more stability based on some comments in the forum.

[us7]

Hi All,
I want your help in adding capillary pressure in the twoPhaseEulerFoam solver that I have modified for porous media (darcy source). May be if you can give some good suggestion.

Few things I want to check about this solver:

1: if I look into the solver (multiphase>twoPhaseEulerFoam) > PEqn.H : I can see its solving the pressure separately (each Phase) then its taking average of both pressures. am i correct?
2: If its solving two pressures in each domain if i want to add capillary pressure for porous media that would be addition of Pc in one equation only like Pc=Pa-Pb and Pa=Pb+Pc where Pc=Pdx(Saturation function model) then where should I add Pc in PEqn.H?

Can you please give me some suggestion on it?

Regards,
Umer

[rdbisme]

Originally Posted by us7

Hi All,
I want your help in adding capillary pressure in the twoPhaseEulerFoam solver that I have modified for porous media (darcy source). May be if you can give some good suggestion.

Few things I want to check about this solver:

1: if I look into the solver (multiphase>twoPhaseEulerFoam) > PEqn.H : I can see its solving the pressure separately (each Phase) then its taking average of both pressures. am i correct?
2: If its solving two pressures in each domain if i want to add capillary pressure for porous media that would be addition of Pc in one equation only like Pc=Pa-Pb and Pa=Pb+Pc where Pc=Pdx(Saturation function model) then where should I add Pc in PEqn.H?

Can you please give me some suggestion on it?

Regards,
Umer

Hello Umer, I highly suggest you to read the Rusche thesis that is linked in the header file of the solver. The scheme and the pressure equation are very well explained there, way better than I could do here
Moreover, I'm not very familiar with capillarity so, on the base of the informations you given, I do not feel I could really be of help. I'm sorry...

For the specific problem, moreover, I suggest you to open a new thread. This thread is for discussing how to suppress the energy equation!

[us7]

Originally Posted by tidusuper91

Hello Umer, I highly suggest you to read the Rusche thesis that is linked in the header file of the solver. The scheme and the pressure equation are very well explained there, way better than I could do here
Moreover, I'm not very familiar with capillarity so, on the base of the informations you given, I do not feel I could really be of help. I'm sorry...

For the specific problem, moreover, I suggest you to open a new thread. This thread is for discussing how to suppress the energy equation!

Hello Ruben Di Battista,
Thanks for your quick and nice reply. As i am not expert in CFD but just looking for few modifications in this solver. I will open new thread for capillarity as suggested but can you tell me how difficult is to add a new equation for mass transfer in it. M=KA(Cs-Ca). I am looking for gas dissolution rate of gas in water.

Regards,
Umer

[rdbisme]

Originally Posted by us7

Hello Ruben Di Battista,
Thanks for your quick and nice reply. As i am not expert in CFD but just looking for few modifications in this solver. I will open new thread for capillarity as suggested but can you tell me how difficult is to add a new equation for mass transfer in it. M=KA(Cs-Ca). I am looking for gas dissolution rate of gas in water.

Regards,
Umer

Well, for my limited experience I would say that adding a source term would be quite hard. (Obviously depending on your knowledge of C++) The solver is highly templated so you need to study a bit on it

[Vam]

Hi everyone,

I changed the geometry fo my simulation to a cilinder with 20 circles as the inlet using SALOME.
I had the following error while running my simulation:

--> FOAM FATAL ERROR: Negative initial temperature T0: -484.949

Does someone know how can I solve it?

[shanvach]

Originally Posted by tidusuper91

I also managed to avoid energy equation solving with following settings in fvSolutions:

"(h|e).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1;
        relTol          0;
        minIter         0;
    maxIter        0;
    }

It's somehow hacky but does the job.

Hi tidusuper91,

Just a couple of follow up questions.
So I am trying to simulate two phase flow with one phase dispersed in another. However, fluids are incompressible and I don't want to solve the energy equation. If you suppress the energy equation the way you did it, does it become an incompressible flow solver.
Your help in this matter is greatly appreciated.

Thanks and Regards,

[rdbisme]

Originally Posted by shanvach

Hi tidusuper91,

Just a couple of follow up questions.
So I am trying to simulate two phase flow with one phase dispersed in another. However, fluids are incompressible and I don't want to solve the energy equation. If you suppress the energy equation the way you did it, does it become an incompressible flow solver.
Your help in this matter is greatly appreciated.

Thanks and Regards,

Well, the only thing it does is not solving the energy equation. It does not change the equations to be solved so... What do you mean by incompressible? Somerhing like SimpleFoam but only for the liquid phase?

[shanvach]

Originally Posted by tidusuper91

Well, the only thing it does is not solving the energy equation. It does not change the equations to be solved so... What do you mean by incompressible? Somerhing like SimpleFoam but only for the liquid phase?

Since it is a compressible solver,I was under the impression that twophaseEulerFoam solves the energy equation so as to calculate the temperature and eventually the density using the equation of state. So if it does not solve the energy equation, that would mean that it does not calculate the temperature and hence density is not solved. So how does it give the solution if it does not solve the energy equation?.
Hope I was clear in my question. Also I just want to simulate two phase imcompressible flow without heat transfer. I am just interested in transport.

[rdbisme]

Originally Posted by shanvach

Since it is a compressible solver,I was under the impression that twophaseEulerFoam solves the energy equation so as to calculate the temperature and eventually the density using the equation of state. So if it does not solve the energy equation, that would mean that it does not calculate the temperature and hence density is not solved. So how does it give the solution if it does not solve the energy equation?.
Hope I was clear in my question. Also I just want to simulate two phase imcompressible flow without heat transfer. I am just interested in transport.

What I believe it does is then just keeping the density value you configured in the initial conditions (its long time I do not use it) except maybe for some drift due to numerical errors...

You could maybe tryout reactingTwoPhaseEulerFoam with an Isothermal phase model...