[faiazk]

Hi All,

I am trying to implement turbulent wind flow over a small building in a domain. At first I used a special technique to generate turbulent flow varying at each time step. I managed to make openfoam read those velocity files at each time step. At first I tried to make the flow develop over the entire domain. Right now, I have my target velocity profiles and turbulence in the entire domain.

Up to this point I did not have to use snappyHexMesh as I did not use the building inside the domain. Everything ran proper in empty domain both in single and parallel processors. But, when I use snappyHexMesh before pisoFoam I get errors. Snappy runs properly though both in single and parallel processors but after using snappy if I use pisoFoam I am receiving the following errors. The courant number gets excessively high when I use snappy before pisoFoam. I ran this case with snappy before with uniform flow without turbulence. But, may be something is going wrong when I use snappy with turbulent inflow.

Please shed some light on this. I would really appreciate any comments and suggestions to make this run properly.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.0-dbb428a3a855
Exec : pisoFoam
Date : Oct 24 2018
Time : 13:30:26
Host : "faiaz-Lenovo-ideapad-320-15IKB"
PID : 5595
I/O : uncollated
Case : /home/faiaz/pracCDRFG/parallel/withsnappy/5origdomdralysug
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PISO: Operating solver in PISO mode

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type LES
Selecting LES turbulence model dynamicKEqn
Selecting LES delta type cubeRootVol
dynamicKEqnCoeffs
{
filter simple;
Ce 1.048;
}

No MRF models present

No finite volume options present


Starting time loop

fieldAverage fieldAverage1:
Starting averaging at time 0

Reading surface description:
nearWall

Time = 0.02

Courant Number mean: 0.00100704711813 max: 0.2285478
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 9.54297156366e-07, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 8.79318426854e-07, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 9.00098713879e-07, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0739342676304, No Iterations 6
time step continuity errors : sum local = 0.00014907945604, global = -7.73755610755e-06, cumulative = -7.73755610755e-06
GAMG: Solving for p, Initial residual = 0.00165381437962, Final residual = 5.17384700419e-07, No Iterations 8
time step continuity errors : sum local = 3.79416905931e-07, global = -8.60204715924e-08, cumulative = -7.82357657914e-06
smoothSolver: Solving for k, Initial residual = 0.999999999935, Final residual = 2.13607265526e-07, No Iterations 2
ExecutionTime = 42.99 s ClockTime = 43 s

Reading/initialising field UMean
Reading/initialising field pMean
Reading/initialising field UPrime2Mean
Reading/initialising field pPrime2Mean
fieldAverage fieldAverage1 write:
Calculating averages

scalarTransport write:
smoothSolver: Solving for s, Initial residual = 1, Final residual = 3.64629548373e-06, No Iterations 4

Time = 0.04

Courant Number mean: 0.325402878711 max: 10.3144914466
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 ? at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? at ??:?
Floating point exception

Thanks in advance.

Faiaz