From dd22a8d8a0c9ba3f54dff586945170ea9c97341c Mon Sep 17 00:00:00 2001 From: henry Date: Wed, 28 May 2008 22:38:03 +0100 Subject: [PATCH] Implamented Rasmus' corrections to the particle force. --- .../multiphase/twoPhaseEulerFoam/alphaEqn.H | 16 ++-- .../bed/constant/polyMesh/boundary | 76 ++++++++----------- .../bed/constant/ppProperties | 2 +- .../twoPhaseEulerFoam/bed/system/fvSchemes | 2 +- .../bed2/constant/ppProperties | 2 +- .../twoPhaseEulerFoam/bed2/system/fvSchemes | 2 +- .../bubbleColumn/constant/ppProperties | 2 +- .../bubbleColumn/system/fvSchemes | 2 +- 8 files changed, 50 insertions(+), 54 deletions(-) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H index b889a6dc6b..e8b056ed7e 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H @@ -24,12 +24,18 @@ if (g0.value() > 0.0) { - surfaceScalarField alphaf = fvc::interpolate(alpha); + ppMagf = rUaAf*fvc::interpolate + ( + (1.0/(rhoa*(alpha + scalar(0.0001)))) + *g0*min(exp(preAlphaExp*(alpha - alphaMax)), expMax) + ); - // Correct the particle-particle force magnitude - ppMagf = g0*rUaAf*min(exp(preAlphaExp*(alphaf - alphaMax)), expMax); - - alphaEqn -= fvm::laplacian(ppMagf, alpha); + alphaEqn -= fvm::laplacian + ( + (fvc::interpolate(alpha) + scalar(0.0001))*ppMagf, + alpha, + "laplacian(alphaPpMag,alpha)" + ); } alphaEqn.relax(); diff --git a/tutorials/twoPhaseEulerFoam/bed/constant/polyMesh/boundary b/tutorials/twoPhaseEulerFoam/bed/constant/polyMesh/boundary index d5522e225c..490a269bd5 100644 --- a/tutorials/twoPhaseEulerFoam/bed/constant/polyMesh/boundary +++ b/tutorials/twoPhaseEulerFoam/bed/constant/polyMesh/boundary @@ -1,59 +1,49 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: 1.4.1 | +| \\ / O peration | Version: dev | | \\ / A nd | Web: http://www.openfoam.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ - FoamFile { - version 2.0; - format binary; - - root "/home/dm2/henry/OpenFOAM/OpenFOAM-1.4.1/tutorials/twoPhaseEulerFoam"; - case "bed"; - instance "constant"; - local "polyMesh"; - - class polyBoundaryMesh; - object boundary; + version 2.0; + format binary; + class polyBoundaryMesh; + object boundary; + location "constant/polyMesh"; } - // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 4 ( -top -{ - type patch; - physicalType inletOutlet; - nFaces 25; - startFace 6100; -} - -bottom -{ - type patch; - physicalType inlet; - nFaces 25; - startFace 6125; -} - -walls -{ - type wall; - physicalType wallFunctions; - nFaces 250; - startFace 6150; -} - -frontBack -{ - type empty; - nFaces 6250; - startFace 6400; -} + top + { + type patch; + physicalType inletOutlet; + nFaces 25; + startFace 6100; + } + bottom + { + type patch; + physicalType inlet; + nFaces 25; + startFace 6125; + } + walls + { + type wall; + physicalType wallFunctions; + nFaces 250; + startFace 6150; + } + frontBack + { + type empty; + nFaces 6250; + startFace 6400; + } ) // ************************************************************************* // diff --git a/tutorials/twoPhaseEulerFoam/bed/constant/ppProperties b/tutorials/twoPhaseEulerFoam/bed/constant/ppProperties index e8d0ed2c01..8d4525b777 100644 --- a/tutorials/twoPhaseEulerFoam/bed/constant/ppProperties +++ b/tutorials/twoPhaseEulerFoam/bed/constant/ppProperties @@ -30,7 +30,7 @@ expMax 1.0e+3; alphaMax 0.6; -g0 g0 [ 0 2 -2 0 0 0 0 ] 1000 ; +g0 g0 [1 -1 -2 0 0 0 0 ] 1000 ; packingLimiter no; diff --git a/tutorials/twoPhaseEulerFoam/bed/system/fvSchemes b/tutorials/twoPhaseEulerFoam/bed/system/fvSchemes index 503a4f8b30..abbb8bdf13 100644 --- a/tutorials/twoPhaseEulerFoam/bed/system/fvSchemes +++ b/tutorials/twoPhaseEulerFoam/bed/system/fvSchemes @@ -52,7 +52,7 @@ laplacianSchemes laplacian(nuEffa,Ua) Gauss linear corrected; laplacian(nuEffb,Ub) Gauss linear corrected; laplacian((rho*(1|A(U))),p) Gauss linear corrected; - laplacian(ppMagf,alpha) Gauss linear corrected; + laplacian(alphaPpMag,alpha) Gauss linear corrected; laplacian((alphaEps*nuEffb),epsilon) Gauss linear corrected; laplacian((alphak*nuEffb),k) Gauss linear corrected; } diff --git a/tutorials/twoPhaseEulerFoam/bed2/constant/ppProperties b/tutorials/twoPhaseEulerFoam/bed2/constant/ppProperties index d9c86b04ee..d3de84081b 100644 --- a/tutorials/twoPhaseEulerFoam/bed2/constant/ppProperties +++ b/tutorials/twoPhaseEulerFoam/bed2/constant/ppProperties @@ -30,7 +30,7 @@ expMax 1.0e+3; alphaMax 0.62; -g0 g0 [ 0 2 -2 0 0 0 0 ] 0 ; +g0 g0 [ 1 -1 -2 0 0 0 0 ] 0 ; packingLimiter no; diff --git a/tutorials/twoPhaseEulerFoam/bed2/system/fvSchemes b/tutorials/twoPhaseEulerFoam/bed2/system/fvSchemes index a0b372ff1f..4887cc3775 100644 --- a/tutorials/twoPhaseEulerFoam/bed2/system/fvSchemes +++ b/tutorials/twoPhaseEulerFoam/bed2/system/fvSchemes @@ -52,7 +52,7 @@ laplacianSchemes laplacian(nuEffa,Ua) Gauss linear corrected; laplacian(nuEffb,Ub) Gauss linear corrected; laplacian((rho*(1|A(U))),p) Gauss linear corrected; - laplacian(ppMagf,alpha) Gauss linear corrected; + laplacian(alphaPpMag,alpha) Gauss linear corrected; laplacian(Galphaf,alpha) Gauss linear corrected; laplacian((alphak*nuEffb),k) Gauss linear corrected; laplacian((alphaEps*nuEffb),epsilon) Gauss linear corrected; diff --git a/tutorials/twoPhaseEulerFoam/bubbleColumn/constant/ppProperties b/tutorials/twoPhaseEulerFoam/bubbleColumn/constant/ppProperties index c7569b3012..1778c4b4f7 100644 --- a/tutorials/twoPhaseEulerFoam/bubbleColumn/constant/ppProperties +++ b/tutorials/twoPhaseEulerFoam/bubbleColumn/constant/ppProperties @@ -30,7 +30,7 @@ expMax 10.0; alphaMax 0.6; -g0 g0 [ 0 2 -2 0 0 0 0 ] 0 ; +g0 g0 [ 1 -1 -2 0 0 0 0 ] 0 ; packingLimiter off; diff --git a/tutorials/twoPhaseEulerFoam/bubbleColumn/system/fvSchemes b/tutorials/twoPhaseEulerFoam/bubbleColumn/system/fvSchemes index 4f568ff98c..1bdee66fd5 100644 --- a/tutorials/twoPhaseEulerFoam/bubbleColumn/system/fvSchemes +++ b/tutorials/twoPhaseEulerFoam/bubbleColumn/system/fvSchemes @@ -52,7 +52,7 @@ laplacianSchemes laplacian(nuEffa,Ua) Gauss linear corrected; laplacian(nuEffb,Ub) Gauss linear corrected; laplacian((rho*(1|A(U))),p) Gauss linear corrected; - laplacian(ppMagf,alpha) Gauss linear corrected; + laplacian(alphaPpMag,alpha) Gauss linear corrected; laplacian((alphak*nuEffb),k) Gauss linear corrected; laplacian((alphaEps*nuEffb),epsilon) Gauss linear corrected; }