[rmz]

hello,
I am working on a simulation of wind on buildings with a complex Mesh.
I am using a RASModel kEspilon with the simpleFoam solver.
I am applying ABL conditions (atmospheric boundary layer).

I am facing problems with the solution convergence. the following is the begining of the output of simpleFoam:

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

Create time

Create mesh for time = 0

SIMPLE: convergence criteria
field p tolerance 0.001
field U tolerance 0.0001
field "(k|epsilon)" tolerance 0.0001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon

bounding k, min: 0 max: 46.975528397 average: 1.29999999999
bounding epsilon, min: 0 max: 173.377103808 average: 0.0100000000001
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.11;
C3 -0.33;
sigmak 1;
}

No MRF models present

No finite volume options present

and this is the end

Time = 1401

--> FOAM Warning :
From function Foam::fv::gaussConvectionScheme<Type>::gaussConvec tionScheme(const Foam::fvMesh&, const surfaceScalarField&, Foam::Istream&) [with Type = Foam::Vector<double>; Foam::surfaceScalarField = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 124
Reading "/home/ingerop/OpenFOAM/ingerop-4.1/run_Dell/PAP_run/PAP_case_06_30/system/fvSchemes.divSchemes.div(phi,U)" at line 32
Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
To remove this warning switch off 'warnUnboundedGauss' in "/opt/openfoam4/etc/controlDict"
smoothSolver: Solving for Ux, Initial residual = 1.73686761346e-09, Final residual = 1.73686761346e-09, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 8.49670942582e-09, Final residual = 8.49670942582e-09, No Iterations 0
smoothSolver: Solving for Uz, Initial residual = 7.97749175341e-10, Final residual = 7.97749175341e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 3.26294168382e-05, Final residual = 2.04903628538e-08, No Iterations 5
GAMG: Solving for p, Initial residual = 1.80088423459e-12, Final residual = 1.80088423459e-12, No Iterations 0
GAMG: Solving for p, Initial residual = 1.80088423459e-12, Final residual = 1.80088423459e-12, No Iterations 0
time step continuity errors : sum local = 1.92360928162e+41, global = -3.25979779589e+25, cumulative = -3.25979779589e+25
smoothSolver: Solving for epsilon, Initial residual = 2.08398503235e-17, Final residual = 2.08398503235e-17, No Iterations 0
smoothSolver: Solving for k, Initial residual = 8.88019064951e-08, Final residual = 7.39368688472e-11, No Iterations 5
bounding k, min: -6.88777988968e+34 max: 3.77003683949e+47 average: 1.86101835569e+41
ExecutionTime = 29170.4 s ClockTime = 29707 s

SIMPLE solution converged in 1401 iterations

End

Finalising parallel run

so my SIMPLE solution converged but I I have strange time step continuity error:
sum local = 1.92360928162e+41, global = -3.25979779589e+25, cumulative = -3.25979779589e+25

I am trynig to understand the meaning of "time step continuity errors". I searches the forum but didn't find a clear description.

What are the possible solutions to this problem.

thank you