CFD Online Discussion Forums - [SOWFA] sowfa single turbine

Hi,

I am going to use SOWFA to simulate a windfarm and compare the results against real measurements. To start, I try to simulate a single wind turbine using the solver "windPlantPisoSolver", and after a few seconds, the execution stops with this error:

Code:

Courant Number mean: 2.0237e+06 max: 1.4461e+08

Time = 0.15    Time Step = 3

C_s = 0.168

u* = 15871 m/s    |<devR>|^(1/2) = 15871 m/s    <U>_1/2 = (-1.7322e+05 -1.3412e+05 0) m/s

Qs = 0 K*(m/s)    |<q>| = 0 K*(m/s)    L = 1e+30 m    phi_m = 1    <T>_1/2 = 532.83 K

Thrust from Body Force = 1.1134e+10    Thrust from Act. Line = 5.5954e+10    Ratio = 0.19898

Torque from Body Force = -2.3774e+09    Torque from Act. Line = -1.1848e+10    Ratio = 0.20066

DILUPBiCG:  Solving for Ux, Initial residual = 0.99957, Final residual = 4.5524e+10, No Iterations 1001

DILUPBiCG:  Solving for Uy, Initial residual = 0.99775, Final residual = 12931, No Iterations 1001

DILUPBiCG:  Solving for Uz, Initial residual = 0.99917, Final residual = 32041, No Iterations 1001

#0  Foam::error::printStack(Foam::Ostream&) at ??:?

#1  Foam::sigFpe::sigHandler(int) at ??:?

#2  Uninterpreted: 

#3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const at ??:?

#4  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?

#5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?

#6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?

#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?

#8  

 at ??:?

#9  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"

#10  

 at ??:?

Floating point exception (core dumped)

If I run the exact same case, with a way larger volume defined in blockMesh, it executes without any problem (but the results are incorrect). From now, even with the 0.15 sec of simulation I can see that the results make more sense, but I don't understand what's making the courant number explode...

Any ideas??

Thanks,
Ruben

Attachment 32266

Attachment 32267