# Rotation of wing with microtabs - PimpleDyMFoam

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

January 7, 2011, 09:21
Rotation of wing with microtabs - PimpleDyMFoam
#1
New Member

Emmanuel Haemers
Join Date: Nov 2010
Posts: 1
Rep Power: 0
Hi all,

It try to solve a case with PimpleDyMFoam;

I use a self-made grid from a wing with snappyHexMesh, made it semi-3D so that case is only 2D (via extrudeMesh, autopatch 90 etc.)

First I used a wing (exported via Autodesk Inventor in STL) and let it rotate to obtain aerodynamic coeff.
Afterwards I substituted my adapted wing with a microtab and the solver stops after a few seconds with the following error:

Code:
```Time = 0.0571638

GAMG:  Solving for cellDisplacementx, Initial residual = 0.00128750962120888624677650380334625879186205565929412841796875, Final residual = 9.103979218607043519785539283095676665880091604776680469512939453125e-07, No Iterations 1
GAMG:  Solving for cellDisplacementy, Initial residual = 0.0011809155049652248457847658613673047511838376522064208984375, Final residual = 6.615464252064893674386107658147437859952333383262157440185546875e-07, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#2   in "/lib/libc.so.6"
#3  void Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::Field<Foam::Vector<double> >&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#4  Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#5  Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#6  Foam::fv::gaussDivScheme<Foam::Tensor<double> >::fvcDiv(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#7  Foam::tmp<Foam::GeometricField<Foam::innerProduct<Foam::Vector<double>, Foam::Tensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Tensor<double> >(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleTurbulenceModel.so"
#8  Foam::tmp<Foam::GeometricField<Foam::innerProduct<Foam::Vector<double>, Foam::Tensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Tensor<double> >(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleTurbulenceModel.so"
#9  Foam::incompressible::RASModels::kOmegaSST::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#10
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/pimpleDyMFoam"
#11  __libc_start_main in "/lib/libc.so.6"
#12
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/pimpleDyMFoam"
Floating point exception```
After investigation the boundaries of my mesh are rotating in a wrong way (attachments). One is with refinment box, other without.

This is strange because the same code for a normal wing works good ( picture rotatingwing) and everything else (except the polymesh map) remains the same:

Code:
```boundaryField
{
wing
{
type        angularOscillatingDisplacement;
axis        (0 0 1);
origin        (5 0 0); // 1/2-chord rotation
angle0        0;//-0.235619;//13.5 deg = 0.235619 rad
omega        6.283;//1 Hz * 2 * Pi
value        uniform (0 0 0);
}

"(front|back)"
{
type            empty;
}

"(inlet|topAndBottom|outlet)"
{
type     fixedValue;
value   uniform (0 0 0);

}

".*"
{
type            fixedValue;
value           uniform (0 0 0);
}```
I used the tutorial wingmotion from OF 1.7 to start from. This tutorial also runs for a limited time.

I'm a beginner in OF and I have no idea how to make this work.

I hope somebody can help.

Emmanuel
Attached Images
 Error with rotation of mesh.jpg (58.7 KB, 72 views) RotatingWing.0776.jpg (49.0 KB, 78 views)
Attached Files
 Error.pdf (88.0 KB, 30 views)

 February 7, 2011, 10:21 #2 Senior Member   Join Date: Apr 2010 Posts: 151 Rep Power: 15 Hello, Maybe try OpenFOAM-1.7.1. The same tutorial is right now running for me. I am at time 0.29 seconds and everything seems stable.

 December 17, 2012, 02:41 #3 New Member   Poria honarmand Join Date: Jul 2012 Location: Iran Posts: 3 Rep Power: 13 hello Dear I am working on wing motion in OF 2.1.1, I am using the pimpleDyMFoam solver, and i cant put frequency and amplitude in point displacement. and what is the meaning of ".*" in file 0, point displacement. thank you

December 20, 2012, 04:16
#4
Member

ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 13
Quote:
 Originally Posted by poriaFoam hello Dear I am working on wing motion in OF 2.1.1, I am using the pimpleDyMFoam solver, and i cant put frequency and amplitude in point displacement. and what is the meaning of ".*" in file 0, point displacement. thank you

hi

I want to work on wing motion OF 2.1.1 . but I could not run tutorial. there wasnt blockmeshdic !

did my problem related to movingmesh ?

 December 21, 2012, 05:09 #5 Senior Member     Mostafa Mahmoudi Join Date: Jan 2012 Posts: 322 Rep Power: 14 run Code: `./Allrun` for more detailed take a look at Allrun file.

 December 21, 2012, 13:24 #6 New Member   Poria honarmand Join Date: Jul 2012 Location: Iran Posts: 3 Rep Power: 13 yes use Allrun

 Tags rotation pimpledym wing

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Ed FLUENT 5 April 13, 2019 13:07 abid Siemens 1 June 16, 2012 13:26 siw Main CFD Forum 3 July 21, 2010 13:55 siw CFX 2 June 22, 2010 10:21 ricardo.halfeld CFX 7 September 8, 2009 18:22

All times are GMT -4. The time now is 23:55.