| Roesch |
February 20, 2020 04:55 |
FOAM FATAL ERROR: Continuity error cannot be removed by adjusting the outflow
I'm running a singlephase CFD-DEM-Simulation with openfoam coupled to liggghts.
It's the first time I'm trying to work with the coupld simulation.
I get the following error:
Quote:
Total # of neighbors = 31483
Ave neighs/atom = 3.25204
Neighbor list builds = 0
Dangerous builds = 0
LIGGGHTS finished
DICPCG: Solving for sSmoothField, Initial residual = 0.08496994, Final residual = 6.31568e-06, No Iterations 4
DICPCG: Solving for vSmoothFieldx, Initial residual = 5.387611e-06, Final residual = 5.387611e-06, No Iterations 0
DICPCG: Solving for vSmoothFieldy, Initial residual = 5.399834e-06, Final residual = 5.399834e-06, No Iterations 0
DICPCG: Solving for vSmoothFieldz, Initial residual = 6.051159e-06, Final residual = 6.051159e-06, No Iterations 0
timeStepFraction() = 1
DICPCG: Solving for vSmoothFieldx, Initial residual = 0.2320728, Final residual = 1.553256e-06, No Iterations 5
DICPCG: Solving for vSmoothFieldy, Initial residual = 0.2267464, Final residual = 1.473994e-06, No Iterations 5
DICPCG: Solving for vSmoothFieldz, Initial residual = 0.2813704, Final residual = 5.429863e-06, No Iterations 4
update Ksl.internalField()
TotalForceExp: (0 0 0)
TotalForceImp: (-0.2967192 -0.3316887 0.02308912)
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 7.638402e-05, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 7.389899e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 5.334098e-05, No Iterations 1
[0] [1]
[1]
[1] --> FOAM FATAL ERROR:
[1] Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux : 5.038763e-05
Specified mass inflow : 3.771036e-06
Specified mass outflow : 1.605388e-06
Adjustable mass outflow : 1.478737e-35
[1]
[1]
[1] From function bool Foam::adjustPhi(Foam::surfaceScalarField&, const volVectorField&, Foam::volScalarField&)
[1] in file cfdTools/general/adjustPhi/adjustPhi.C at line 107.
[1]
FOAM parallel run exiting
|
So I understand there is something wrogn with my boundaries. But bc I'm trying to wokring with a singlephase simulation there are no topics available with this problem. My p and u boundaries are listed below:
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
".*"
{
type fixedFluxPressure;
value uniform 0;
}
freesurface
{
type totalPressure;
p0 uniform 0;
}
}
// ************************************************** *********************** //
|
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
".*"
{
//type fixedValue;
//value uniform (0 0 0);
type noSlip;
//type zeroGradient;
}
freesurface
{
type slip;
|
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object Us;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
".*"
{
type zeroGradient;
}
freesurface
{
type slip;
}
}
|
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1;
boundaryField
{
".*"
{
type zeroGradient;
}
}
|
Thanks in advance for your help.
Regards
Oh and I'm using the PISO algorithm |
