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

Rotation of wing with microtabs - PimpleDyMFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 7, 2011, 10:21
Default Rotation of wing with microtabs - PimpleDyMFoam
  #1
New Member
 
Emmanuel Haemers
Join Date: Nov 2010
Posts: 1
Rep Power: 0
s0170068 is on a distinguished road
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
        amplitude    -0.349066282;//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.

Thanks in advance!

Emmanuel
Attached Images
File Type: jpg Error with rotation of mesh.jpg (58.7 KB, 68 views)
File Type: jpg RotatingWing.0776.jpg (49.0 KB, 72 views)
Attached Files
File Type: pdf Error.pdf (88.0 KB, 28 views)
s0170068 is offline   Reply With Quote

Old   February 7, 2011, 11:21
Default
  #2
Senior Member
 
Join Date: Apr 2010
Posts: 151
Rep Power: 7
flowris is on a distinguished road
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.
flowris is offline   Reply With Quote

Old   December 17, 2012, 03:41
Default
  #3
New Member
 
Poria honarmand
Join Date: Jul 2012
Location: Iran
Posts: 3
Rep Power: 5
poriaFoam is on a distinguished road
Send a message via Yahoo to 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
poriaFoam is offline   Reply With Quote

Old   December 20, 2012, 05:16
Unhappy
  #4
Member
 
ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 4
ali jafari is on a distinguished road
Quote:
Originally Posted by poriaFoam View Post
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 ?

please help me about movingMesh !! (its too hard)
ali jafari is offline   Reply With Quote

Old   December 21, 2012, 06:09
Default
  #5
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
run
Code:
./Allrun
for more detailed take a look at Allrun file.
adambarfi is offline   Reply With Quote

Old   December 21, 2012, 14:24
Default
  #6
New Member
 
Poria honarmand
Join Date: Jul 2012
Location: Iran
Posts: 3
Rep Power: 5
poriaFoam is on a distinguished road
Send a message via Yahoo to poriaFoam
yes use Allrun
poriaFoam is offline   Reply With Quote

Reply

Tags
rotation pimpledym wing

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
3D Wing Analysis Ed FLUENT 4 August 3, 2012 05:57
Forward swept wing simulation abid CD-adapco 1 June 16, 2012 13:26
[Gridgen] Blocking topology for blunt wing siw Main CFD Forum 3 July 21, 2010 13:55
Isolating wing induced drag component siw CFX 2 June 22, 2010 10:21
How to set domain translation and rotation in CFX ricardo.halfeld CFX 7 September 8, 2009 18:22


All times are GMT -4. The time now is 20:27.