CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[SOWFA] sowfa single turbine

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   July 11, 2014, 05:18
Default sowfa single turbine
New Member
ruben p.
Join Date: Dec 2013
Location: Copenhagen, Denmark
Posts: 7
Rep Power: 12
ollebapur is on a distinguished road

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:

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 ??:?
 at ??:?
#9  __libc_start_main in "/lib/i386-linux-gnu/"
 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??

ollebapur is offline   Reply With Quote

Old   November 25, 2014, 05:35
Default Re: SOWFA single turbine
New Member
Thomas Schulz
Join Date: Jul 2014
Posts: 17
Rep Power: 12
cico0815 is on a distinguished road
Dear Ruben,

I recently started to work with the solver suite provided by the NREL and I ran into several problems with the usage. Most of them I already solved and I think your problem with those solvers is based on several modifications made to OpenFOAM during the version changes. It may compile correctly but that does not mean it is working correctly.

You may want to look into the source code of the pisoFoamTurbine solver included in the SOWFA suite and compare it to the original pisoFoam solver sources. You will find several changes (appart from the FAST modifications). If you take the original pisoFoam solver as a base and add the modifications according to the suggestions made by the SOWFA people you might be able to let it run correctly.

I also think that a good start might be the fastDuct tutorial included in the suite instead of the windPlant tutorial.

Was this helpful?


cico0815 is offline   Reply With Quote


sowfa, turbine

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Single moving reference frame vs Multiple reference frame for Wind Turbine CFD mrswordf1sh FLUENT 4 January 25, 2020 08:43
Single rotating frame wind turbine meshing wih ICEM xXAhmedXx CFD Freelancers 5 January 16, 2018 14:55
ANSYS Simulation of 3d single wind turbine blade jaye.wolf ANSYS 0 March 11, 2016 03:16
[SOWFA] Scale down on SOWFA Actuator Line Turbine Model korawin10287 OpenFOAM Community Contributions 0 August 9, 2014 18:56
Boundary condition setting regarding turbine simulation using CFX Lacerlacer CFX 11 March 12, 2012 09:32

All times are GMT -4. The time now is 18:45.