CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Rotation of wing with microtabs - PimpleDyMFoam (http://www.cfd-online.com/Forums/openfoam-solving/83683-rotation-wing-microtabs-pimpledymfoam.html)

s0170068 January 7, 2011 10:21

Rotation of wing with microtabs - PimpleDyMFoam
 
3 Attachment(s)
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

flowris February 7, 2011 11:21

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.

poriaFoam December 17, 2012 03:41

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

ali jafari December 20, 2012 05:16

Quote:

Originally Posted by poriaFoam (Post 397849)
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)

adambarfi December 21, 2012 06:09

run
Code:

./Allrun
for more detailed take a look at Allrun file.

poriaFoam December 21, 2012 14:24

yes use Allrun


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