| gohome |
July 8, 2014 04:12 |
Failed to simulate wind turbine based on wingmotion tutorial with pimpleDyMFoam
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"
|
