CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

oscillatingRotatingMotion: simulation time change with omega and amplitude values

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 9, 2017, 07:22
Default oscillatingRotatingMotion: simulation time change with omega and amplitude values
  #1
Member
 
Join Date: May 2017
Posts: 38
Rep Power: 8
decibelle is on a distinguished road
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
decibelle is offline   Reply With Quote

Old   June 10, 2017, 16:19
Default
  #2
Senior Member
 
khedar
Join Date: Oct 2016
Posts: 111
Rep Power: 9
khedar is on a distinguished road
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.
khedar is offline   Reply With Quote

Old   June 13, 2017, 04:09
Default
  #3
Member
 
Join Date: May 2017
Posts: 38
Rep Power: 8
decibelle is on a distinguished road
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
decibelle is offline   Reply With Quote

Old   November 14, 2018, 07:28
Default
  #4
New Member
 
Chayanit Nigaltia
Join Date: Jan 2018
Posts: 29
Rep Power: 8
CHAYANIT is on a distinguished road
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
CHAYANIT is offline   Reply With Quote

Reply

Tags
motion, oscillating, pimpledymfoam, rotating, simulation time

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
interDyMFoam with VOF + 6DOF instable pbalz OpenFOAM Running, Solving & CFD 11 October 9, 2020 06:19
How to calculate oscillating roll of 2D body (with variable amplitude and omega) lavrov OpenFOAM Running, Solving & CFD 0 March 19, 2015 07:31


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