|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
ON
Join Date: Jun 2021
Posts: 6
Rep Power: 5 ![]() |
Greetings,
I am using multiphaseEulerFoam to simulate a micro-bioreactor (31mm*18mm*87mm) with a tiny mixer (d=11.4mm) at 300 rpm. It contains water and gas phases. The sparger pumps air into the mixer at -0.00103745 m/s and the gas leaves the micro-bioreactor from the outlet. I got bounding epsilonm as shown below, sometimes bounding km. And soon the equation will diverge and obtain huge negative pressure or temperature. Code:
bounding epsilonm, min: -0.00814733 max: 16.3522 average: 0.0880985 My Code:
ddtSchemes { default Euler; } gradSchemes { default cellLimited leastSquares 0.5; nHat Gauss linear; //Face unit interface normal } divSchemes { default none; div(phi,alpha.air) Gauss vanLeer; div(phi,alpha.water) Gauss vanLeer; div(phir,alpha.water,alpha.air) Gauss vanLeer; div(phir,alpha.air,alpha.water) Gauss vanLeer; "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinear 1; "div\(phi.*,U.*\)" Gauss limitedLinear 1; "div\(alphaRhoPhi.*,(h|e).*\)" Gauss limitedLinear 1; "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; "div\(alphaRhoPhi.*,\(p\|thermo:rho.*\)\)" Gauss limitedLinear 1; "div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1; "div\(phim,(k|epsilon)m\)" Gauss limitedLinear 1; "div\(\(\(\(alpha.**thermo:rho.*\)*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; } laplacianSchemes { default Gauss linear limited 1; } interpolationSchemes { default linear; } snGradSchemes { default limited 1; } Code:
solvers { "alpha.*" { nAlphaCorr 2; nAlphaSubCycles 2; nLimiterIter 10; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-8; relTol 0; minIter 1; } p_rgh { solver GAMG; smoother DIC; tolerance 1e-8; relTol 0.01; minIter 3; } p_rghFinal { $p_rgh; relTol 0; minIter 3; } "U.*" { solver PBiCGStab; preconditioner DILU; tolerance 1e-08; relTol 0; minIter 3; } "e.*" { solver PBiCGStab; preconditioner DILU; tolerance 1e-08; relTol 0; minIter 0; } "(k|epsilon).*" { solver PBiCGStab; preconditioner DILU; tolerance 1e-08; relTol 0; minIter 3; } } PIMPLE { consistent yes; nOuterCorrectors 3; //1-3 nCorrectors 2; //1-2 nNonOrthogonalCorrectors 1; } relaxationFactors { fields { ".*" 0.7; } equations { ".*" 0.7; } } |
|
![]() |
![]() |
![]() |
Tags |
bounding epsilon, bounding k |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[isoAdvector] Is it possible to add isoAdvector capabilities to interphaseChangeFoam? | chaohui | OpenFOAM Community Contributions | 9 | June 21, 2024 12:21 |
Converging Diverging Nozzle with dbnsTurbFoam | Saleh Abuhanieh | OpenFOAM Running, Solving & CFD | 4 | December 13, 2019 10:26 |
[solids4Foam] How to calculate drag coeff when using solids4Foam | amuzeshi | OpenFOAM CC Toolkits for Fluid-Structure Interaction | 15 | November 7, 2019 12:50 |
Internal Flow Instability with layers | Zephiro88 | OpenFOAM Running, Solving & CFD | 1 | June 20, 2019 04:20 |
Bounding k and epsilon & coldEngineFoam | sasanghomi | OpenFOAM | 1 | September 13, 2013 12:12 |