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

Suppress twoPhaseEulerFoam energy

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

Like Tree27Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 11, 2017, 16:20
Default
  #21
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by hooman.4028 View Post
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
rdbisme is offline   Reply With Quote

Old   July 11, 2017, 16:30
Default
  #22
Member
 
Hooman
Join Date: Apr 2011
Posts: 35
Rep Power: 15
hooman.4028 is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
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

Code:
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

Code:
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

Code:
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

Code:
/*--------------------------------*- 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

Code:
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

Code:
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 is offline   Reply With Quote

Old   July 11, 2017, 17:16
Default
  #23
Member
 
Hooman
Join Date: Apr 2011
Posts: 35
Rep Power: 15
hooman.4028 is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
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

Code:

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

Code:
dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

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


}
hooman.4028 is offline   Reply With Quote

Old   July 12, 2017, 06:51
Default
  #24
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by hooman.4028 View Post
The other files:

phase properties

Code:

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

Code:
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.
rarnaunot likes this.
rdbisme is offline   Reply With Quote

Old   July 12, 2017, 17:20
Default
  #25
Member
 
Hooman
Join Date: Apr 2011
Posts: 35
Rep Power: 15
hooman.4028 is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
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.

Last edited by hooman.4028; July 12, 2017 at 19:11.
hooman.4028 is offline   Reply With Quote

Old   November 8, 2017, 12:55
Default
  #26
us7
New Member
 
Umer
Join Date: Aug 2016
Posts: 29
Rep Power: 10
us7 is on a distinguished road
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
us7 is offline   Reply With Quote

Old   November 9, 2017, 09:52
Default
  #27
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by us7 View Post
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!
rdbisme is offline   Reply With Quote

Old   November 10, 2017, 06:09
Default
  #28
us7
New Member
 
Umer
Join Date: Aug 2016
Posts: 29
Rep Power: 10
us7 is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
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
us7 is offline   Reply With Quote

Old   November 10, 2017, 16:45
Default
  #29
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by us7 View Post
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
rdbisme is offline   Reply With Quote

Old   May 24, 2018, 01:51
Default twophaseeulerfoam bubble column negative temperature
  #30
Vam
New Member
 
Join Date: Apr 2018
Posts: 2
Rep Power: 0
Vam is on a distinguished road
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?
Vam is offline   Reply With Quote

Old   September 25, 2018, 16:43
Default
  #31
Member
 
Join Date: Apr 2016
Posts: 30
Rep Power: 10
shanvach is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
I also managed to avoid energy equation solving with following settings in fvSolutions:

Code:
"(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,
shanvach is offline   Reply With Quote

Old   September 25, 2018, 18:21
Default
  #32
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by shanvach View Post
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?
rdbisme is offline   Reply With Quote

Old   September 25, 2018, 18:38
Default
  #33
Member
 
Join Date: Apr 2016
Posts: 30
Rep Power: 10
shanvach is on a distinguished road
Quote:
Originally Posted by tidusuper91 View Post
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.
shanvach is offline   Reply With Quote

Old   September 25, 2018, 18:45
Default
  #34
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by shanvach View Post
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...
rdbisme is offline   Reply With Quote

Reply

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
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 1 September 24, 2019 01:15
energy spectra from LES simulations kumar OpenFOAM Post-Processing 2 February 23, 2016 13:25
Natural convection, problems with total energy balance MdoNascimento STAR-CCM+ 0 March 13, 2014 13:15
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
SIMPLE and energy equation convergence Fabio Main CFD Forum 0 June 1, 2007 07:06


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