|
[Sponsors] |
twoPhaseEulerFoam confined plunging jet simulation failed |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 17, 2020, 05:53 |
twoPhaseEulerFoam confined plunging jet simulation failed
|
#1 |
Member
Veskov Eugene
Join Date: Feb 2011
Posts: 31
Rep Power: 15 |
Hello, Foamers!
My confined plunging jet simulation in twoPhaseEulerFoam failed after 3-rd PIMPLE iteration on the first time step: Boundary condition FoamFile { version 2.0; format ascii; class volScalarField; object epsilonm; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 1.5e-4; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phim; inletValue $internalField; value $internalField; } wall { type zeroGradient; value $internalField; } wallinner { type zeroGradient; value $internalField; } } dimensions [0 2 -2 0 0 0 0]; internalField uniform 3.75e-5; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phim; inletValue $internalField; value $internalField; } wall { type zeroGradient; value $internalField; } wallinner { type zeroGradient; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 1e5; boundaryField { opening { type calculated; value $internalField; } inlet { type calculated; value $internalField; } outlet { type calculated; value $internalField; } wall { type calculated; value $internalField; } wallinner { type calculated; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { opening { type totalPressure; p0 uniform 0; U U.air; phi phi.air; value uniform 0; } inlet { type fixedFluxPressure; value uniform 0; } outlet { type fixedFluxPressure; value prghPressure; } wall { type fixedFluxPressure; value uniform 0; } wallinner { type fixedFluxPressure; value uniform 0; } } FoamFile { version 2.0; format ascii; class volScalarField; object Theta; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.0; boundaryField { inlet { type fixedValue; value uniform 1.0e-7; } outlet { type inletOutlet; inletValue uniform 1.0e-7; value uniform 1.0e-7; } wall { type zeroGradient; } wallinner { type zeroGradient; } opening { type fixedValue; value uniform 1.0e-7; } } FoamFile { version 2.0; format ascii; class volScalarField; object alphat.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { opening { type calculated; value $internalField; } inlet { type calculated; value $internalField; } outlet { type calculated; value $internalField; } wall { type compressible::alphatWallFunction; Prt 0.85; value $internalField; } wallinner { type compressible::alphatWallFunction; Prt 0.85; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object epsilon.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 1.5e-4; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value uniform 0.000371; } outlet { type inletOutlet; phi phi.air; inletValue $internalField; value $internalField; } wall { type epsilonWallFunction; value $internalField; } wallinner { type epsilonWallFunction; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object k.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 3.75e-5; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phi.air; inletValue $internalField; value $internalField; } wall { type kqRWallFunction; value $internalField; } wallinner { type kqRWallFunction; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object kappai.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 -1 0 0 0 0 0]; internalField uniform 2000; boundaryField { wall { type zeroGradient; } wallinner { type zeroGradient; } opening { type inletOutlet; phi phi.air; inletValue $internalField; value $internalField; } inlet { type zeroGradient; } outlet { type zeroGradient; } } FoamFile { version 2.0; format ascii; class volScalarField; object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 1e-8; boundaryField { opening { type calculated; value uniform 0; } inlet { type fixedValue; value uniform 1.3416e-05; } outlet { type calculated; value $internalField; } wall { type nutkWallFunction; value uniform 0; } wallinner { type nutkWallFunction; value uniform 0; } } FoamFile { version 2.0; format ascii; class volScalarField; object T.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phi.air; inletValue $internalField; value $internalField; } wall { type zeroGradient; } wallinner { type zeroGradient; } } FoamFile { version 2.0; format binary; class volVectorField; object U.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { opening { //type fluxCorrectedVelocity; type pressureInletOutletVelocity; value uniform (0 0 0); } inlet { type fixedValue; value uniform (0 0 0); } outlet { type pressureInletOutletVelocity; phi phi.air; value $internalField; } wall { type noSlip; } wallinner { type noSlip; } } FoamFile { version 2.0; format ascii; class volScalarField; object alphat.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { opening { type calculated; value $internalField; } inlet { type calculated; value $internalField; } outlet { type calculated; value $internalField; } wall { type compressible::alphatWallFunction; Prt 0.85; value $internalField; } wallinner { type compressible::alphatWallFunction; Prt 0.85; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object epsilon.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 1.5e-4; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value uniform 0.000371; } outlet { type inletOutlet; phi phi.water; inletValue $internalField; value $internalField; } wall { type epsilonWallFunction; value $internalField; } wallinner { type epsilonWallFunction; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object k.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 3.75e-5; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phi.water; inletValue $internalField; value $internalField; } wall { type kqRWallFunction; value $internalField; } wallinner { type kqRWallFunction; value $internalField; } } FoamFile { version 2.0; format ascii; class volScalarField; object nut.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 1e-8; boundaryField { opening { type calculated; value uniform 0; } inlet { type fixedValue; value uniform 1.3416e-05; } outlet { type calculated; value $internalField; } wall { type nutkWallFunction; value uniform 0; } wallinner { type nutkWallFunction; value uniform 0; } } FoamFile { version 2.0; format ascii; class volScalarField; object T.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { opening { type fixedValue; value $internalField; } inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; phi phi.air; inletValue $internalField; value $internalField; } wall { type zeroGradient; } wallinner { type zeroGradient; } } FoamFile { version 2.0; format binary; class volVectorField; object U.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { opening { //type fluxCorrectedVelocity; type pressureInletOutletVelocity; value uniform (0 0 0); } inlet { type flowRateInletVelocity; massFlowRate 0.205; extrapolateProfile yes; rho rho; rhoInlet 1.0; value uniform (0 0 0); } outlet { type flowRateInletVelocity; massFlowRate 0.205; extrapolateProfile yes; rho rho; rhoInlet 1.0; value uniform (0 0 0); } wall { type noSlip; } wallinner { type noSlip; } } FoamFile { version 2.0; format ascii; class dictionary; location "system"; object setFieldsDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // defaultFieldValues ( volScalarFieldValue alpha.air 1 volScalarFieldValue alpha.water 0 ); regions ( boxToCell { box (-0.04 0 -0.04) (0.04 1.0 0.04); fieldValues ( volScalarFieldValue alpha.air 0 volScalarFieldValue alpha.water 1 ); } ); |
|
June 17, 2020, 05:55 |
solver settings and error
|
#2 |
Member
Veskov Eugene
Join Date: Feb 2011
Posts: 31
Rep Power: 15 |
FoamFile
{ version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { alpha.air { nAlphaCorr 1; nAlphaSubCycles 2; MULESCorr yes; // Switches on semi-implicit MULES nLimiterIter 8; // Number of MULES iterations over the limiter alphaApplyPrevCorr true; } p_rgh { solver GAMG; smoother DIC; tolerance 1e-8; relTol 0; } p_rghFinal { $p_rgh; relTol 0; } "U.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; minIter 1; } "e.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; minIter 1; } "(k|epsilon|Theta).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; minIter 1; } } PIMPLE { nOuterCorrectors 3; nCorrectors 1; nNonOrthogonalCorrectors 1; } relaxationFactors { equations { ".*" 0.1; } } Create time Create mesh for time = 0 PIMPLE: no residual control data found. Calculations will employ 3 corrector loops Reading g Reading hRef Creating twoPhaseSystem Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.air Selecting diameterModel for phase air: IATE Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState perfectFluid; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.water Selecting diameterModel for phase water: constant Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } Selecting default blending method: hyperbolic Selecting aspectRatioModel for (air in water): constant Selecting aspectRatioModel for (water in air): constant Selecting dragModel for (air and water): segregated Selecting swarmCorrection for (air and water): none Selecting dragModel for (air in water): SchillerNaumann Selecting swarmCorrection for (air in water): Tomiyama Selecting dragModel for (water in air): SchillerNaumann Selecting swarmCorrection for (water in air): Tomiyama Selecting virtualMassModel for (air in water): constantCoefficient Selecting virtualMassModel for (water in air): constantCoefficient Selecting heatTransferModel for (air in water): RanzMarshall Selecting heatTransferModel for (water in air): RanzMarshall Selecting liftModel for (air in water): LegendreMagnaudet Selecting liftModel for (water in air): LegendreMagnaudet Selecting wallLubricationModel for (air in water): Frank Selecting wallLubricationModel for (water in air): Frank Selecting turbulentDispersionModel for (air in water): constantCoefficient Selecting turbulentDispersionModel for (water in air): constantCoefficient Calculating field g.h Reading field p_rgh Creating field dpdt Creating field kinetic energy K No MRF models present No finite volume options present Courant Number mean: 0 max: 0 Max Ur Courant Number = 0 Calculating field DDtU1 and DDtU2 Starting time loop fieldAverage fieldAverage1: Restarting averaging for fields: U.air: starting averaging at time 0 U.water: starting averaging at time 0 alpha.air: starting averaging at time 0 p: starting averaging at time 0 Courant Number mean: 0 max: 0 Max Ur Courant Number = 0 deltaT = 6e-08 Time = 6e-08 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.09078137948 Min(alpha.air) = 0 Max(alpha.air) = 1 smoothSolver: Solving for kappai.air, Initial residual = 1, Final residual = 0.001264025234, No Iterations 1000 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 4.792804359e-05, Final residual = 4.211162828e-05, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.4524787668, Final residual = 0.4138812702, No Iterations 1000 min T.air 300 min T.water 300 Selecting patchDistMethod meshWave GAMG: Solving for p_rgh, Initial residual = 0.001385065939, Final residual = 2.795003837e-15, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 9.083520552e-16, Final residual = 9.083520552e-16, No Iterations 0 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.09077876671 Min(alpha.air) = -11.41765647 Max(alpha.air) = 1 smoothSolver: Solving for kappai.air, Initial residual = 0.9999999892, Final residual = 3.872285928e-08, No Iterations 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.9999978023, Final residual = 1.312056595e-10, No Iterations 2 smoothSolver: Solving for e.water, Initial residual = 0.9999999985, Final residual = 7.404774966e-10, No Iterations 2 min T.air 299.2787486 min T.water 274.0068754 GAMG: Solving for p_rgh, Initial residual = 0.0001938205006, Final residual = 8.062402264e-14, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1.396477599e-16, Final residual = 1.396477599e-16, No Iterations 0 PIMPLE: iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.09077953577 Min(alpha.air) = -2.486427675 Max(alpha.air) = 1 |
|
June 17, 2020, 06:08 |
solvers setting from cfd-online and errors
|
#3 |
Member
Veskov Eugene
Join Date: Feb 2011
Posts: 31
Rep Power: 15 |
with this settings Multiphase
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-5; relTol 0; minIter 1; } p_rgh { solver GAMG; smoother DIC; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 500; agglomerator faceAreaPair; mergeLevels 1; tolerance 1e-5; relTol 0; } 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-5; relTol 0; minIter 1; } "Theta.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-5; relTol 0; minIter 1; } "(k|epsilon).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-5; relTol 0; minIter 1; } "epsilonFinal" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-08; relTol 0; } } PIMPLE { nOuterCorrectors 3; nCorrectors 1; nNonOrthogonalCorrectors 1; pRefCell 0; pRefValue 0; residualControl { U { tolerance 1e-5; relTol 0; } p { tolerance 5e-4; relTol 0; } } } relaxationFactors { equations { ".*" 1; } } Courant Number mean: 0 max: 0 Max Ur Courant Number = 0 deltaT = 6e-08 Time = 6e-08 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.09078137948 Min(alpha.air) = 0 Max(alpha.air) = 1 smoothSolver: Solving for kappai.air, Initial residual = 1, Final residual = 0.0009141209896, No Iterations 1000 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 4.819361816e-05, Final residual = 4.231010327e-05, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.4591847899, Final residual = 0.4243305836, No Iterations 1000 min T.air 300 min T.water 300 Selecting patchDistMethod meshWave GAMG: Solving for p_rgh, Initial residual = 0.00148599746, Final residual = 6.291118185e-13, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 2.053537248e-13, Final residual = 2.053537248e-13, No Iterations 0 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air --> FOAM FATAL ERROR: cannot be called for a calculatedFvPatchField on patch outlet of field alpha.air in file "/mnt/c/Users/admin/tutorial/confjet/confjettfef/0/alpha.air" You are probably trying to solve for a field with a default boundary condition. From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientIntern alCoeffs() const [with Type = double] in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 188. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::calculatedFvPatchField<double>::gradientInte rnalCoeffs() const at ??:? #3 Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #4 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #5 Foam::twoPhaseSystem::solve() at ??:? #6 ? at ??:? #7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #8 ? at ??:? Aborted (core dumped) |
|
Tags |
confined plunging jet, twophaseeulerfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Executing simulation in parallel. orte_init failed | vidyadhar | OpenFOAM Running, Solving & CFD | 2 | June 30, 2017 01:13 |
3D DES simulation of turbulent jet | Kostas_K | FLUENT | 9 | April 14, 2017 10:19 |
Initial conditions for uniform flow | andreas | OpenFOAM | 5 | November 16, 2012 16:00 |
DeardorffDiffStress turbulence model | Zuixy | OpenFOAM Running, Solving & CFD | 0 | November 14, 2011 09:44 |
[OpenFOAM] ParaView/Parafoam error when making animation | Disco_Caine | ParaView | 6 | September 28, 2010 10:54 |