pbalz |
February 10, 2015 08:28 |
interDyMFoam with VOF + 6DOF instable
1 Attachment(s)
Dear Foamers,
I try to calculate the pitch angle and draft of a floating body in water, using a VOF approach with dynamic remeshing. The DTCHull tutorial for the interDyMFoam solver is basically exactly what I need, but of course with another geometry.
So, after testing some different parameters, I took the DTCHull tutorial and adapted it to my geometry and boundary conditions (i.e. u=3.7 m/s instead of 1.6 m/s). I have not made any change to the major options regarding discretization or solver setup. By the way: I'm running both OF 2.3.0 and 2.3.x on Ubuntu 12.04.
This is where the problem all starts.
OF starts to calculate the solution and after a few time steps the pressure forces start to become unstable.
Afterwards the linear and angular velocity of my solid also start to oscillate - with a rising amplitude... A few timesteps later the velocity has reached really high values, which leads to a high courant number and therefore to a decreasing delta t. Then alpha becomes unbounded and the solution quits with a floating point exception (obviously because of the decreasing time step size).
Here are the first two time steps of my logfile. You can clearly see that the pressure forces are already unstable:
Code:
Starting time loop
Interface Courant Number mean: 0 max: 0
Courant Number mean: 8.0273e-05 max: 0.0100382
deltaT = 2.39992e-05
Time = 2.39992e-05
Restraint translationDamper: force (-0 -0 -0)
Restraint rotationDamper: moment (-0 -0 -0)
6-DoF rigid body motion
Centre of rotation: (0 0 0)
Centre of mass: (0 0 0)
Orientation: (1 0 0 0 1 0 0 0 1)
Linear velocity: (0 0 1.23235e-05)
Angular velocity: (0 2.32008e-09 0)
Execution time for mesh.update() = 1.2 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0828371, No Iterations 4
time step continuity errors : sum local = 9.49826e-10, global = -3.30579e-11, cumulative = -4.13543e-11
smoothSolver: Solving for alpha.water, Initial residual = 5.63182e-09, Final residual = 1.19013e-13, No Iterations 1
Phase-1 volume fraction = 0.790864 Min(alpha1) = 0 Max(alpha1) = 1.00001
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.790864 Min(alpha1) = -8.01534e-16 Max(alpha1) = 1.00001
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.000695745, No Iterations 6
time step continuity errors : sum local = 3.84999e-10, global = -3.62078e-11, cumulative = -7.75621e-11
GAMG: Solving for p_rgh, Initial residual = 0.0246502, Final residual = 1.58645e-05, No Iterations 5
time step continuity errors : sum local = 3.40771e-11, global = -3.0063e-12, cumulative = -8.05684e-11
GAMG: Solving for p_rgh, Initial residual = 0.00308546, Final residual = 4.63604e-08, No Iterations 11
time step continuity errors : sum local = 9.57511e-14, global = -7.38217e-15, cumulative = -8.05757e-11
smoothSolver: Solving for omega, Initial residual = 2.37895e-06, Final residual = 2.26823e-12, No Iterations 1
smoothSolver: Solving for k, Initial residual = 1, Final residual = 3.01431e-11, No Iterations 2
ExecutionTime = 35.37 s ClockTime = 36 s
forces forces output:
sum of forces:
pressure : (-2.45527e+07 958.752 -128685)
viscous : (-49.4466 0.000300508 -0.00483469)
porous : (0 0 0)
sum of moments:
pressure : (226245 -2.52705e+08 -4.63274e+07)
viscous : (0.00899411 17.5397 -104.706)
porous : (0 0 0)
Interface Courant Number mean: 0 max: 0
Courant Number mean: 9.63109e-05 max: 0.0328231
deltaT = 2.87988e-05
Time = 5.2798e-05
Restraint translationDamper: force (-0 -0 -0.211866)
Restraint rotationDamper: moment (-0 -5.37608e-05 -0)
6-DoF rigid body motion
Centre of rotation: (0 0 7.09807e-10)
Centre of mass: (0 0 7.09807e-10)
Orientation: (1 0 1.33631e-13 0 1 0 -1.33631e-13 0 1)
Linear velocity: (0 0 -8.34663e-05)
Angular velocity: (0 -0.0432865 0)
Execution time for mesh.update() = 1.39 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0872902, No Iterations 2
time step continuity errors : sum local = 2.05229e-12, global = -1.18994e-13, cumulative = -8.06947e-11
smoothSolver: Solving for alpha.water, Initial residual = 3.72614e-09, Final residual = 1.68466e-13, No Iterations 1
Phase-1 volume fraction = 0.790864 Min(alpha1) = -2.42957e-16 Max(alpha1) = 1.00001
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.790864 Min(alpha1) = -9.02248e-10 Max(alpha1) = 1.00001
GAMG: Solving for p_rgh, Initial residual = 0.00749169, Final residual = 4.46451e-06, No Iterations 9
time step continuity errors : sum local = 1.29741e-11, global = 8.75379e-13, cumulative = -7.98194e-11
GAMG: Solving for p_rgh, Initial residual = 0.0401748, Final residual = 1.36224e-05, No Iterations 5
time step continuity errors : sum local = 2.42329e-11, global = 9.33856e-13, cumulative = -7.88855e-11
GAMG: Solving for p_rgh, Initial residual = 0.00439205, Final residual = 2.06961e-08, No Iterations 10
time step continuity errors : sum local = 3.57504e-14, global = 2.9144e-15, cumulative = -7.88826e-11
smoothSolver: Solving for omega, Initial residual = 2.85482e-06, Final residual = 3.16525e-12, No Iterations 1
smoothSolver: Solving for k, Initial residual = 0.0587749, Final residual = 3.24357e-08, No Iterations 1
ExecutionTime = 57.99 s ClockTime = 59 s
forces forces output:
sum of forces:
pressure : (-1.27512e+06 1470.39 126192)
viscous : (-49.1517 0.000384695 -0.0032971)
porous : (0 0 0)
sum of moments:
pressure : (-241862 -8.11241e+06 -2.46272e+06)
viscous : (0.00708633 18.9842 -104.16)
porous : (0 0 0)
Interface Courant Number mean: 1.17239e-11 max: 0.0258687
Courant Number mean: 0.000115569 max: 0.0277234
deltaT = 3.45578e-05
Time = 8.73558e-05
And finally the last time steps before everything explodes:
Code:
Restraint translationDamper: force (-0 -0 -6.62811e+10)
Restraint rotationDamper: moment (-0 1.08837e+08 -0)
6-DoF rigid body motion
Centre of rotation: (0 0 -0.153159)
Centre of mass: (0 0 -0.153159)
Orientation: (0.999959 0 -0.00901383 0 1 0 0.00901383 0 0.999959)
Linear velocity: (0 0 -2.27725e+07)
Angular velocity: (0 -2.02891e+06 0)
Execution time for mesh.update() = 1.35 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0625312, No Iterations 4
time step continuity errors : sum local = 0.000358987, global = -3.31254e-05, cumulative = -3.35953e-05
smoothSolver: Solving for alpha.water, Initial residual = 0.000479791, Final residual = 4.8796e-11, No Iterations 3
Phase-1 volume fraction = 0.790864 Min(alpha1) = -1.29611e-10 Max(alpha1) = 1.04721
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.790864 Min(alpha1) = -1.29447e-10 Max(alpha1) = 1.04721
GAMG: Solving for p_rgh, Initial residual = 0.801051, Final residual = 0.000572986, No Iterations 9
time step continuity errors : sum local = 6.1879e-06, global = -6.77958e-07, cumulative = -3.42732e-05
GAMG: Solving for p_rgh, Initial residual = 0.000450071, Final residual = 2.31368e-07, No Iterations 7
time step continuity errors : sum local = 1.55703e-06, global = 6.38958e-07, cumulative = -3.36343e-05
GAMG: Solving for p_rgh, Initial residual = 0.999972, Final residual = 3.22908e-08, No Iterations 20
time step continuity errors : sum local = 0.00749391, global = -0.000607174, cumulative = -0.000640809
smoothSolver: Solving for omega, Initial residual = 1, Final residual = 9.37478e-08, No Iterations 121
smoothSolver: Solving for k, Initial residual = 0.00732862, Final residual = 9.28197e-08, No Iterations 140
ExecutionTime = 2566.2 s ClockTime = 2594 s
forces forces output:
sum of forces:
pressure : (1.37784e+34 2.66482e+34 1.18088e+36)
viscous : (1.90643e+18 -1.43781e+19 2.28481e+19)
porous : (0 0 0)
sum of moments:
pressure : (-2.14866e+36 2.87143e+35 1.80333e+34)
viscous : (-7.09308e+19 1.06297e+19 6.54609e+18)
porous : (0 0 0)
Interface Courant Number mean: 113849 max: 3.45359e+08
Courant Number mean: 1.26635e+07 max: 6.24547e+08
deltaT = 2.00261e-21
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 10 to 11 to distinguish between timeNames at time 0.0823347
Time = 0.082334668882
Restraint translationDamper: force (0 0 4.57787e+11)
Restraint rotationDamper: moment (-0 4.6905e+10 -0)
6-DoF rigid body motion
Centre of rotation: (0 0 -0.153159)
Centre of mass: (0 0 -0.153159)
Orientation: (0.999959 0 -0.00901383 0 1 0 0.00901383 0 0.999959)
Linear velocity: (0 0 4.84996e+10)
Angular velocity: (0 3.44134e+09 0)
Execution time for mesh.update() = 1.35 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0574656, No Iterations 4
time step continuity errors : sum local = 2.40604e+13, global = 1.49371e+12, cumulative = 1.49371e+12
smoothSolver: Solving for alpha.water, Initial residual = 0.391807, Final residual = 8.20113e-11, No Iterations 58
Phase-1 volume fraction = 0.848726 Min(alpha1) = 0 Max(alpha1) = 6.78522
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 2.79008e+76 Min(alpha1) = -1.87585e+97 Max(alpha1) = 2.50798e+97
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00714166, No Iterations 1000
time step continuity errors : sum local = 5.71911e+32, global = -1.41715e+31, cumulative = -1.41715e+31
GAMG: Solving for p_rgh, Initial residual = 3.65482e-37, Final residual = 3.65482e-37, No Iterations 0
time step continuity errors : sum local = 1.2657e+55, global = -5.14288e+54, cumulative = -5.14288e+54
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 4.96727e-08, No Iterations 246
time step continuity errors : sum local = 4.77892e+102, global = -3.54492e+102, cumulative = -3.54492e+102
[1] #0 [3] Foam::error::printStack(Foam::Ostream&)#0 Foam::error::printStack(Foam::Ostream&) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #8 Foam::fvMatrix<double>::solve() in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[3] #8 Foam::fvMatrix<double>::solve() in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[3] #9 Foam::incompressible::RASModels::kOmegaSST::correct() in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #9 Foam::incompressible::RASModels::kOmegaSST::correct() in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[3] #10 in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #10
[1] in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[1] #11 __libc_start_main[3] in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[3] #11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #12 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #12
[3] in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[balz-VirtualBox:03920] *** Process received signal ***
[balz-VirtualBox:03920] Signal: Floating point exception (8)
[balz-VirtualBox:03920] Signal code: (-6)
[balz-VirtualBox:03920] Failing at address: 0x3e800000f50
[balz-VirtualBox:03920] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36150) [0x7f2b5747b150]
[balz-VirtualBox:03920] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f2b5747b0d5]
[balz-VirtualBox:03920] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36150) [0x7f2b5747b150]
[balz-VirtualBox:03920] [ 3] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x43c) [0x7f2b58584bcc]
[balz-VirtualBox:03920] [ 4] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x2d) [0x7f2b5858505d]
[balz-VirtualBox:03920] [ 5] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x480) [0x7f2b5857dbd0]
[balz-VirtualBox:03920] [ 6] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7f2b5ad9eff7]
[balz-VirtualBox:03920] [ 7] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x11c) [0x7f2b5bece1ec]
[balz-VirtualBox:03920] [ 8] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam8fvMatrixIdE5solveEv+0xca) [0x7f2b5bece54a]
[balz-VirtualBox:03920] [ 9] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0xc5b) [0x7f2b5bf144cb]
[balz-VirtualBox:03920] [10] interDyMFoam() [0x42dbc9]
[balz-VirtualBox:03920] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f2b5746676d]
[balz-VirtualBox:03920] [12] interDyMFoam() [0x437abd]
[balz-VirtualBox:03920] *** End of error message ***
[1] in "/home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[balz-VirtualBox:03918] *** Process received signal ***
[balz-VirtualBox:03918] Signal: Floating point exception (8)
[balz-VirtualBox:03918] Signal code: (-6)
[balz-VirtualBox:03918] Failing at address: 0x3e800000f4e
[balz-VirtualBox:03918] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36150) [0x7f7632c6c150]
[balz-VirtualBox:03918] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f7632c6c0d5]
[balz-VirtualBox:03918] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36150) [0x7f7632c6c150]
[balz-VirtualBox:03918] [ 3] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x43c) [0x7f7633d75bcc]
[balz-VirtualBox:03918] [ 4] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x2d) [0x7f7633d7605d]
[balz-VirtualBox:03918] [ 5] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x480) [0x7f7633d6ebd0]
[balz-VirtualBox:03918] [ 6] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7f763658fff7]
[balz-VirtualBox:03918] [ 7] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x11c) [0x7f76376bf1ec]
[balz-VirtualBox:03918] [ 8] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam8fvMatrixIdE5solveEv+0xca) [0x7f76376bf54a]
[balz-VirtualBox:03918] [ 9] /home/balz/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0xc5b) [0x7f76377054cb]
[balz-VirtualBox:03918] [10] interDyMFoam() [0x42dbc9]
[balz-VirtualBox:03918] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f7632c5776d]
[balz-VirtualBox:03918] [12] interDyMFoam() [0x437abd]
[balz-VirtualBox:03918] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 3920 on node balz-VirtualBox exited on signal 8 (Floating point exception).
I really don't know how to get this case running. After a few disappointing days of trying I ran out of ideas.
Things I have tried, but without any success:
1. setting the time stepping to fixed and decrease the time step size
2. playing around with the boundary conditions (mainly pressure and velocity)
3. changing the mesh size, varying from 500k cells to 3 mio cells
4. changed some SixDofRigidBodyMotionCoeffs; inner+outer Distance, damping factors, more constraints
5. serial processing
6. increasing the volume of my flow region
I really want to know how to get a stable solution without any fluctuations, since this worked perfectly with the dtchull geometry. First I thought it could be some ill conditioned boundary conditions... But at the moment it's more likely that it has something to do with the geometry, because I couldn't achieve a solution with Fluent either.
For those who are interested or have an idea on how to solve this problem, please find attached my case directory, including a (simplified) stl.
Please, I really need some help on this topic! Any idea on how to reach a solution is highly appreciated!
Regards,
Pascal
|