CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Failed to simulate wind turbine based on wingmotion tutorial with pimpleDyMFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 8, 2014, 04:12
Default Failed to simulate wind turbine based on wingmotion tutorial with pimpleDyMFoam
  #1
New Member
 
Carl Liu
Join Date: Apr 2014
Posts: 3
Rep Power: 3
gohome is on a distinguished road
Hi, everyone
I am quite new to OF and trying to simulate a horizontal axis wind turbine with wingmotion tutorial in pimpleDyMFoam.
Any help or suggestion would be greatly appreciated!
PS: I am using OF230.
PPS: The cases are over 100k and can be downloaded from the link.

My thought is to treat the wind turbine as a rotating rotor which is driven by the inlet flow.
This rotor is constrained by a axis and a fixed point (= its geometry center), and its rotation is restrained by a damper.

But the integration crashes after several time steps.
You can check this by running the "Case1" here.
My guess is that the meshes rotate with the wind turbine and some are squeezed too much as seen in "screenshot1.png"
Attachment 32057

So I turned to AMI which would work for rotation objects.
I modified the propeller tutorial based on my limited understanding
and it is the 'Case2' here.
The meshes are generated rightly but it seems the initial conditions are not right as the flow is nearly zero during the whole integration.
And the wind turbine is not rotating any more...

Any suggestions?
Thanks in advance!

------------------------

The error mesg is similar, I dont know how to fix it.
Any advice? thanks a lot!

Code:
 [power:…ingMotion2D_pimpleDyMFoam]$ cat log 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : pimpleDyMFoam
Date   : Jul 08 2014
Time   : 15:57:24
Host   : "power"
PID    : 28690
Case   : /home/gohome/Study/openfoam230/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Translational constraint tensor (0 0 0 0 0 0 0 0 0)
Rotational constraint tensor (1 0 0 0 0 0 0 0 0)

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

AMI: Creating addressing and weights between 15678 source faces and 15822 target faces
AMI: Patch source sum(weights) min/max/average = 0.5010934349, 1.474859027, 1.000194821
AMI: Patch target sum(weights) min/max/average = 0.09389027004, 1.168376359, 0.9990085188
Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85034;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 0.0007921648591 max: 0.1782758803

Starting time loop

Courant Number mean: 0.0007921648591 max: 0.1782758803
deltaT = 1e-05
Time = 1e-05


Restraint rotationDamper:  moment (-0 -0 -0)
Centre of mass: (0 -0.053 -0.089)
Linear velocity: (0 0 0)
Angular velocity: (-1.276190562e-05 0 0)
AMI: Creating addressing and weights between 15678 source faces and 15822 target faces
AMI: Patch source sum(weights) min/max/average = 0.5010934349, 1.474859027, 1.000194821
AMI: Patch target sum(weights) min/max/average = 0.09389027004, 1.168376359, 0.9990085188
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.01332130797, No Iterations 6
time step continuity errors : sum local = 2.281601699e-09, global = -7.114131082e-11, cumulative = -7.114131082e-11
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.9999999971, Final residual = 0.01170120474, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0150662549, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.01508156985, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.006603088823, No Iterations 7
time step continuity errors : sum local = 1.132741528e-09, global = -7.336139937e-11, cumulative = -1.445027102e-10
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.1250681082, Final residual = 3.712558956e-07, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.2256651087, Final residual = 6.60532103e-07, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.2263188243, Final residual = 6.689522455e-07, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.4776872574, Final residual = 5.675631116e-08, No Iterations 20
time step continuity errors : sum local = 1.899086223e-14, global = -3.752338358e-15, cumulative = -1.445064625e-10
smoothSolver:  Solving for omega, Initial residual = 6.276897656e-05, Final residual = 6.053036459e-07, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 7.811010559e-07, No Iterations 3
ExecutionTime = 12.28 s  ClockTime = 12 s

Courant Number mean: 0.0007924403218 max: 0.1849597784
deltaT = 1e-05
Time = 2e-05


Restraint rotationDamper:  moment (2.552381125e-08 -0 -0)
Centre of mass: (0 -0.053 -0.089)
Linear velocity: (0 0 0)
Angular velocity: (-1.235007842 0 0)
AMI: Creating addressing and weights between 15678 source faces and 15822 target faces
AMI: Patch source sum(weights) min/max/average = 0.5010934349, 1.474859027, 1.000194821
AMI: Patch target sum(weights) min/max/average = 0.09389027004, 1.168376359, 0.9990085188
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.01120405059, No Iterations 2
time step continuity errors : sum local = 3.999152841e-14, global = -3.703542962e-15, cumulative = -1.445101661e-10
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.08762828166, Final residual = 0.0002195620438, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.2149869513, Final residual = 0.0005942728825, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.2522776861, Final residual = 0.0007084008448, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.4194730027, Final residual = 0.002164392235, No Iterations 6
time step continuity errors : sum local = 1.333936083e-09, global = 1.287276558e-10, cumulative = -1.578251025e-11
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.01067637411, Final residual = 1.592974448e-07, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.04016122951, Final residual = 9.17971348e-07, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.04743291374, Final residual = 1.178594246e-07, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.4885660006, Final residual = 6.938117308e-08, No Iterations 19
time step continuity errors : sum local = 3.235710057e-14, global = 3.14924399e-15, cumulative = -1.5779361e-11
smoothSolver:  Solving for omega, Initial residual = 6.627998761e-06, Final residual = 4.274861074e-10, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.9795699155, Final residual = 7.435090749e-07, No Iterations 3
ExecutionTime = 20.74 s  ClockTime = 21 s

Courant Number mean: 0.0007924656356 max: 0.2009499482
deltaT = 1e-05
Time = 3e-05


Restraint rotationDamper:  moment (0.002469990161 -0 -0)
Centre of mass: (0 -0.053 -0.089)
Linear velocity: (0 0 0)
Angular velocity: (-2.244016598 0 0)
AMI: Creating addressing and weights between 15678 source faces and 15822 target faces
AMI: Patch source sum(weights) min/max/average = 0.5010851984, 1.474838948, 1.000167963
AMI: Patch target sum(weights) min/max/average = 0.09385795596, 1.168373659, 0.9990083961
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.008189369104, No Iterations 6
time step continuity errors : sum local = 5.587740237e-11, global = 6.714291038e-13, cumulative = -1.51079319e-11
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.01169287521, Final residual = 3.289517688e-05, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.04629253508, Final residual = 0.0001506767796, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.06114240829, Final residual = 0.0002021080614, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.6402490417, Final residual = 0.004007032311, No Iterations 5
time step continuity errors : sum local = 1.103559705e-09, global = -7.064827283e-11, cumulative = -8.575620473e-11
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.002485534071, Final residual = 4.017024319e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.01193583705, Final residual = 2.926351742e-07, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.01718045466, Final residual = 4.354358227e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.4945616154, Final residual = 6.032139435e-08, No Iterations 21
time step continuity errors : sum local = 1.158522461e-14, global = -2.282932827e-15, cumulative = -8.575848766e-11
smoothSolver:  Solving for omega, Initial residual = 6.431649157e-06, Final residual = 3.787475548e-10, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.3114458663, Final residual = 2.366442508e-07, No Iterations 3
ExecutionTime = 29.57 s  ClockTime = 30 s

Courant Number mean: 0.0007964539046 max: 0.1759835946
deltaT = 1e-05
Time = 4e-05


Restraint rotationDamper:  moment (0.002018043034 -0 -0)
Centre of mass: (0 -0.053 -0.089)
Linear velocity: (0 0 0)
Angular velocity: (-2.030410901 0 0)
AMI: Creating addressing and weights between 15678 source faces and 15822 target faces
AMI: Patch source sum(weights) min/max/average = 0, 1.474822556, 1.000108614
AMI: Patch target sum(weights) min/max/average = 0.09383137474, 1.168371466, 0.9990082437
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.01490658934, No Iterations 5
time step continuity errors : sum local = 1.862936667e-11, global = -2.374184112e-13, cumulative = -8.599590607e-11
PIMPLE: iteration 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#5  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"

Last edited by wyldckat; October 19, 2014 at 05:59. Reason: Reconstructed accidentally deleted post and merged with related post. attached image was lost
gohome is offline   Reply With Quote

Reply

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
[GAMBIT] Gambit Start Up error .... HELP PLZ! Abhya ANSYS Meshing & Geometry 2 July 12, 2012 00:40
Reading LESModel & Thermophysical Properties In Similar Manner to Function Object deji OpenFOAM 0 July 27, 2010 10:54
plz help its urgent garima chaudhary FLUENT 0 September 25, 2007 05:28
plz help again! sameera Main CFD Forum 9 August 16, 2007 08:54
plz ......help me out.....************************ krishna FLUENT 2 January 10, 2007 03:53


All times are GMT -4. The time now is 15:26.