CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   oscillatingRotatingMotion: simulation time change with omega and amplitude values (https://www.cfd-online.com/Forums/openfoam-solving/188904-oscillatingrotatingmotion-simulation-time-change-omega-amplitude-values.html)

decibelle June 9, 2017 06:22

oscillatingRotatingMotion: simulation time change with omega and amplitude values
 
hello,

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

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:
Code:

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
Code:

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 June 10, 2017 15:19

Code:

#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 June 13, 2017 03:09

1 Attachment(s)
Hello,

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

Say me if you want more information.

Thanks

Attachment 56620

CHAYANIT November 14, 2018 06:28

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


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