[Roesch]

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:

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:

/*--------------------------------*- 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;

}


}

// ************************************************** *********************** //

/*--------------------------------*- 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;

/*--------------------------------*- 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;
}

}

/*--------------------------------*- 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

[mabrouka]

Hi there,
Could you please share the coupling properties and the FvSolution files ?