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
|