| s0170068 |
January 7, 2011 09: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 |
