SetFields for LES simulation
 

Hello. I am running a LES simulation on a simple square pipe geometry using PisoFoam as solver. I would like to initialize some turbulence setting some different initial values of the velocity fields in different regions of my geometry using SetFieldsDict.
Without using SetFields, the simulation works properly. If I use SetFields, I get this error, and I do not understand what I should do to solve the problem. Thank you in advance.



/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : pisoFoam
Date : Apr 29 2019
Time : 11:55:41
Host : "lrspc49"
PID : 15432
I/O : uncollated
Case : /home/andrea/Desktop/pitzDaily
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

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
--> FOAM Warning :
From function void Foam::LESModels::cubeRootVolDelta::calcDelta()
in file LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C at line 55
Case is 2D, LES is not strictly applicable

bounding k, min: 0 max: 2e-05 average: 0
dynamicKEqnCoeffs
{
filter simple;
Ce 1.048;
}

No MRF models present

No finite volume options present
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::LESModels::dynamicKEqn<Foam::IncompressibleT urbulenceModel<Foam::transportModel> >::Ck(Foam::GeometricField<Foam::SymmTensor<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6 Foam::LESModels::dynamicKEqn<Foam::IncompressibleT urbulenceModel<Foam::transportModel> >::correctNut(Foam::GeometricField<Foam::SymmTenso r<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7 Foam::LESModels::dynamicKEqn<Foam::IncompressibleT urbulenceModel<Foam::transportModel> >::correctNut() at ??:?
#8 ? at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? at ??:?
Floating point exception (core dumped)

A Floating Point Exception often results from a division by zero. Without further information about your variables i'd guess that either the set velocity values is far to larger and instantly breaks the Courant number or one of the turbulence variables is far to small.