[s0170068]

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:

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:

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]

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]

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]

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 ?

please help me about movingMesh !! (its too hard)

[adambarfi]

run

./Allrun

for more detailed take a look at Allrun file.

[poriaFoam]

yes use Allrun