AttouFerschneider: Correction to cube-root of volume fraction ratio

applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/AttouFerschneider/AttouFerschneider.C

@@ -52,7 +52,10 @@ Foam::dragModels::AttouFerschneider::KGasLiquid
     const phaseModel& solid = gas.fluid().phases()[solidName_];
 
     const volScalarField oneMinusGas(max(1 - gas, liquid.residualAlpha()));
-    const volScalarField cbrtR(solid/oneMinusGas);
+    const volScalarField cbrtR
+    (
+        cbrt(max(solid, solid.residualAlpha())/oneMinusGas)
+    );
     const volScalarField magURel(mag(gas.U() - liquid.U()));
 
     return
@@ -70,7 +73,10 @@ Foam::dragModels::AttouFerschneider::KGasSolid
 ) const
 {
     const volScalarField oneMinusGas(max(1 - gas, solid.residualAlpha()));
-    const volScalarField cbrtR(solid/oneMinusGas);
+    const volScalarField cbrtR
+    (
+        cbrt(max(solid, solid.residualAlpha())/oneMinusGas)
+    );
 
     return
         E1_*gas.mu()*sqr(oneMinusGas/solid.d())*sqr(cbrtR)

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/T.air

@@ -20,14 +20,7 @@ internalField       uniform 300;
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type               inletOutlet;
-        phi                phi.air;
-        inletValue         $internalField;
-        value              $internalField;
-    }
-    "(top|walls)"
     {
         type               zeroGradient;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/T.solid

@@ -10,7 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       volScalarField;
-    object      T.solid;
+    object      T.air;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/T.water

@@ -20,14 +20,7 @@ internalField       uniform 300;
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type               inletOutlet;
-        phi                phi.water;
-        inletValue         $internalField;
-        value              $internalField;
-    }
-    "(top|walls)"
     {
         type               zeroGradient;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/U.air

@@ -20,13 +20,7 @@ internalField       uniform (0 0 0);
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type                pressureInletOutletVelocity;
-        phi                 phi.air;
-        value               $internalField;
-    }
-    "(top|walls)"
     {
         type               noSlip;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/U.water

@@ -20,13 +20,7 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type                pressureInletOutletVelocity;
-        phi                 phi.water;
-        value               $internalField;
-    }
-    "(top|walls)"
     {
         type               noSlip;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/alpha.air.orig

@@ -20,14 +20,7 @@ internalField   uniform 0.5;
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type            inletOutlet;
-        phi             phi.air;
-        inletValue      $internalField;
-        value           $internalField;
-    }
-    "(top|walls)"
     {
         type            zeroGradient;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/alpha.water.orig

@@ -20,14 +20,7 @@ internalField   uniform 0;
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type            inletOutlet;
-        phi             phi.water;
-        inletValue      $internalField;
-        value           $internalField;
-    }
-    "(top|walls)"
     {
         type            zeroGradient;
     }

tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/trickleBed/0/p_rgh

@@ -20,16 +20,7 @@ internalField       uniform 1e5;
 
 boundaryField
 {
-    bottom
+    "(bottom|top|walls)"
-    {
-        type            prghTotalPressure;
-        p0              $internalField;
-        U               U.air;
-        phi             phi.air;
-        rho             thermo:rho.air;
-        value           $internalField;
-    }
-    "(top|walls)"
     {
         type            fixedFluxPressure;
         value           $internalField;