CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

sowfa single turbine

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

LinkBack Thread Tools 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: 5
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, 06:35
Default Re: SOWFA single turbine
New Member
Thomas Schulz
Join Date: Jul 2014
Posts: 17
Rep Power: 5
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
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Air Flow Analysis around a Wind Turbine NightHawkRaven FLUENT 4 August 12, 2014 17:18
FSI - Wind Turbine AUN CFX 13 August 29, 2012 16:44
Boundary condition setting regarding turbine simulation using CFX Lacerlacer CFX 11 March 12, 2012 10:32
Problem with torque output from a turbine simulation. pa-dundas FLUENT 3 August 30, 2011 12:25
P4 1.5 or Dual P3 800EB on Gibabyte board Danial FLUENT 4 September 12, 2001 11:44

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