CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Bug in interDyMFoam 6DoF (http://www.cfd-online.com/Forums/openfoam-bugs/62334-bug-interdymfoam-6dof.html)

hsieh January 24, 2009 13:21

latest OF-1.5.x fixed the prob
 
latest OF-1.5.x fixed the problem when the domain rotating to 180 degree.

I just ran a case where the domain moves from 0 degree to 180 degree, then, rotates back from 180 to 0 degree. As time reached 0 degree, solver crashed. If I changed from 0 degree to 0.0001 degree, then, run continues without problem.

Pei

henry January 24, 2009 14:25

Could you post the motion desc
 
Could you post the motion description so we can reproduce the problem here?

Thanks

H

hsieh January 25, 2009 16:32

Hi, Henry, Below is the err
 
Hi, Henry,

Below is the error message.

I will see if I can post the motion profile.

Pei
---------------------
solidBodyMotionFunctions::SKA::transformation(): Time = 3.13106 transformation: ((0 0 0) (1 (0 0 9.15736e-06)))
Execution time for mesh.update() = 0.07 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.336652 Min(gamma) = 5.68349e-06 Max(gamma) = 0.998867
MULES: Solving for gamma
Liquid phase volume fraction = 0.33665 Min(gamma) = 5.6871e-06 Max(gamma) = 0.998867
MULES: Solving for gamma
Liquid phase volume fraction = 0.336648 Min(gamma) = 5.6871e-06 Max(gamma) = 0.998867
GAMG: Solving for pd, Initial residual = 0.0441962, Final residual = 0.000147613, No Iterations 2
time step continuity errors : sum local = 4.67942e-06, global = -5.79765e-17, cumulative = 5.09779e-16
GAMGPCG: Solving for pd, Initial residual = 0.00885975, Final residual = 4.89186e-10, No Iterations 6
time step continuity errors : sum local = 1.55695e-11, global = -5.80112e-17, cumulative = 4.51768e-16
ExecutionTime = 823.46 s ClockTime = 828 s

Courant Number mean: 0.0425424 max: 0.446283
deltaT = 0.00378788
Time = 3.13485

solidBodyMotionFunctions::SKA::transformation(): Time = 3.13485 transformation: ((0 0 0) (1 (0 0 0)))
Execution time for mesh.update() = 0.07 s
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::surfaceInterpolation::makeWeights() const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::surfaceInterpolation::weights() const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::linear<double>::weights(Foam::GeometricField <double,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<doubl e,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libinterfaceProperties.s o"
#7 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::tmp<foam::geometricfield<foam::outerproduct< foam::vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e,> const&, Foam::word const&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libinterfaceProperties.s o"
#9 Foam::tmp<foam::geometricfield<foam::outerproduct< foam::vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e,> const&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libinterfaceProperties.s o"
#10 Foam::LimitedScheme<double,>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#11 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<do uble,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libinterfaceProperties.s o"
#12 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<doubl e,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libinterfaceProperties.s o"
#13 Foam::fv::gaussConvectionScheme<double>::interpola te(Foam::GeometricField<double ,> const&, Foam::GeometricField<double,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#14 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double,> const&, Foam::GeometricField<double,> const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#15 Foam::tmp<foam::geometricfield<double,> > Foam::fvc::flux<double>(Foam::GeometricField<doubl e,> const&, Foam::GeometricField<double,> const&, Foam::word const&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/interDyMFoa m"
#16 Foam::tmp<foam::geometricfield<double,> > Foam::fvc::flux<double>(Foam::tmp<foam::geometricf ield<double,> > const&, Foam::tmp<foam::geometricfield<double,> > const&, Foam::word const&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/interDyMFoa m"
#17 main in "/home/phsieh/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/interDyMFoa m"
#18 __libc_start_main in "/lib64/libc.so.6"
#19 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception
phsieh@huwei:~/OpenFOAM/phsieh-1.5/run/sloshingTank3D6DoF>
---------------------

hsieh January 25, 2009 16:38

Hi, Henry, Let's see if the
 
Hi, Henry,

Let's see if the data file is small enough to upload here.

http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif 6DoF.dat

Pei

henry January 26, 2009 08:17

Thanks for the bug report and
 
Thanks for the bug report and the case setup. We have now fixed the problem and will push the fix to 1.5.x shortly.

H


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