[decibelle]

hello,

I work on a simulation with pimpleDyMFoam for an oscillating and rotating motion.
I have choose these parameters in the controlDict

application     pimpleDyMFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         1;

deltaT          1e-2;
writeControl    adjustableRunTime;
writeInterval   0.01;

////- For testing with moveDynamicMesh
//deltaT          0.01;
//writeControl    timeStep;
//writeInterval   1;

purgeWrite      0;

writeFormat     binary;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable false;

adjustTimeStep  no;

maxCo           2;

when I run pimpleDyMFoam, the run is interrupted with this kind of error:

Courant Number mean: 616728 max: 3.53165e+09
Time = 0.07                                 

AMI: Creating addressing and weights between 19452 source faces and 19984 target faces
AMI: Patch source sum(weights) min/max/average = 0.50888, 1.3735, 1.00185             
AMI: Patch target sum(weights) min/max/average = 0, 1.38576, 0.995279                 
PIMPLE: iteration 1                                                                   
smoothSolver:  Solving for Ux, Initial residual = 0.61257, Final residual = 0.0523612, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.944341, Final residual = 0.0355865, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.781761, Final residual = 0.0242732, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"                                                                                                                   
#1  Foam::sigFpe::sigHandler(int) in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" 
#2  ? in "/lib64/libc.so.6"                                                                                                 
#3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"                                                                             
#4  ? in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"                              
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"                                   
#6  ? in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"                              
#7  __libc_start_main in "/lib64/libc.so.6"                                                                                 
#8  ? in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"                              
Floating point exception (core dumped)

my first question is: someone know what's does this error means?

and after if I change the omega and the amplitude values in the dynamicMeshDict

dynamicFvMesh   solidBodyMotionFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solidBodyMotionFvMeshCoeffs
{
    cellZone        capsule_zone;

    solidBodyMotionFunction  oscillatingRotatingMotion;
    oscillatingRotatingMotionCoeffs
    {
        origin      (0 0 0);
        amplitude        (1.5 0 0);
        omega       2; // rad/s
    }
}

this error appears at a different simulation time. For example for amplitude = 1.5 and omega =2 the error appears at 0.07 time step. For amplitude = 3 and omega =6 the error appears at 0.03.

so my second question is: someone know why the omega and amplitude values have an impact on the time simulation?

Thanks in advance for your help

[khedar]

#3  Foam::divide(Foam::Field<double>&, double const&, F

This is one of the most common error and normally occurs when there is a division by zero somewhere. Usual solution is to check your boundary conditions, initial conditions, mesh again. if the case is small enough upload it here.

[decibelle]

Hello,

I joined at this reply my case without polymesh bescause it's too big.

Say me if you want more information.

Thanks

capsule_prop.zip

[CHAYANIT]

Can you tell me what does amplitude here refers to?I mean we already stated this parameter in the pointDisplacement file in the oscillating foil case.
I think Amplitude (1.5 0 0) means it will rotate about x-axis but what does the value 1.5 refers to?
Any help would be appreciated.
Regards
Chayanit