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

pimpleDymFoam for tidal turbine

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

Like Tree1Likes
  • 1 Post By andre.weiner

Reply
 
LinkBack Thread Tools Display Modes
Old   August 18, 2014, 05:19
Default pimpleDymFoam for tidal turbine
  #1
Member
 
chen long
Join Date: Dec 2012
Posts: 32
Rep Power: 4
Jackie Chen is on a distinguished road
Hi

I try to simulate the flow around tidal turbines by using pimpleDyMFoam. however, the simulation stopped right after i executing the "Allrun" command. The msg in the log file are shown as below, anyone know the reason?

Thanks

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

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone innercylinder
Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 2660 source faces and 2660 target faces
AMI: Patch source weights min/max/average = 4.83927e-09, 1.00063, 0.834913
AMI: Patch target weights min/max/average = 0, 1.20756, 0.834843
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}

Reading field rAU if present

No field sources present


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


Starting time loop

Courant Number mean: 9.59912e-05 max: 0.0060191
deltaT = 1.1999e-05
Time = 1.1999e-05

solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 1.1999e-05 transformation: ((0 0 0) (1 (0 -5.02615e-05 0)))
AMI: Creating addressing and weights between 2660 source faces and 2660 target faces
AMI: Patch source weights min/max/average = 1.91434e-05, 1.21386, 0.961821
AMI: Patch target weights min/max/average = 0, 1.24578, 0.960982
PIMPLE: iteration 1
[1] #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[1] #5
[1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[1] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7
[1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[ubuntu:05018] *** Process received signal ***
[ubuntu:05018] Signal: Floating point exception (8)
[ubuntu:05018] Signal code: (-6)
[ubuntu:05018] Failing at address: 0x3e80000139a
[ubuntu:05018] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f15e79524a0]
[ubuntu:05018] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f15e7952425]
[ubuntu:05018] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f15e79524a0]
[ubuntu:05018] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0x24) [0x7f15e8aa2754]
[ubuntu:05018] [ 4] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7vol MeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_ 11dimensionedIdEERKS8_+0x2a0) [0x43e180]
[ubuntu:05018] [ 5] pimpleDyMFoam() [0x4186cb]
[ubuntu:05018] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f15e793d76d]
[ubuntu:05018] [ 7] pimpleDyMFoam() [0x41c38d]
[ubuntu:05018] *** End of error message ***
Jackie Chen is offline   Reply With Quote

Old   August 18, 2014, 06:30
Default
  #2
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 27
Rep Power: 4
andre.weiner is on a distinguished road
The solver crashes because some weights are zero.

Code:
AMI: Patch source weights min/max/average = 1.91434e-05, 1.21386, 0.961821
AMI: Patch target weights min/max/average = 0, 1.24578, 0.960982
You can try the low weight correction available in OF 2.3

Code:
AMIPatch 
{ 
    type            cyclicAMI; 
    ... 
    lowWeightCorrection 0.2; 
}
Also your mesh seems to be poor at the AMI-patches (average should be close to one)
so maybe you should improve the mesh quality first.

Best regards!
Jackie Chen likes this.
andre.weiner is offline   Reply With Quote

Old   August 18, 2014, 08:13
Default
  #3
Member
 
chen long
Join Date: Dec 2012
Posts: 32
Rep Power: 4
Jackie Chen is on a distinguished road
Andre

Thanks a lot. Normally how you improve the quality of mesh?
Jackie Chen is offline   Reply With Quote

Old   August 18, 2014, 08:33
Default
  #4
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 27
Rep Power: 4
andre.weiner is on a distinguished road
A better mesh quality regarding to AMI faces in OF you achieve e.g. with :

- a better approximation of your geometry
- AMI cell faces of the same shape and size

How you accomplish this depends on your meshing tool.
andre.weiner is offline   Reply With Quote

Old   August 18, 2014, 08:41
Default
  #5
Member
 
chen long
Join Date: Dec 2012
Posts: 32
Rep Power: 4
Jackie Chen is on a distinguished road
Andre

what meshing tool are you using?

Best regards
Jackie Chen is offline   Reply With Quote

Old   August 18, 2014, 08:47
Default
  #6
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 27
Rep Power: 4
andre.weiner is on a distinguished road
That pretty much depends on the available resources :-)
You can use snappyHexMesh like in the propeller tutorial

cd $FOAM_TUTORIALS/incompressible/pimpleDyMFoam/propeller/
andre.weiner is offline   Reply With Quote

Old   August 18, 2014, 11:09
Default
  #7
Member
 
chen long
Join Date: Dec 2012
Posts: 32
Rep Power: 4
Jackie Chen is on a distinguished road
Andre

My case is based on thsi tutorial. M using snappyhexmesh utility, but i only simple execute the command dunno how to set the details to control the quality of the mesh.
May be i need to understand the snappyhexMeshDict so that i can produce more quality mesh of my turbine.
Jackie Chen 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
Wind Turbine from SRFSimpleFoam to pimpleDymFoam be_inspired OpenFOAM 3 March 11, 2015 11:58
Air Flow Analysis around a Wind Turbine NightHawkRaven FLUENT 4 August 12, 2014 17:18
simulating a wind turbine s.q FLUENT 1 October 17, 2013 08:11
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


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