[jinyaqing]

Hi, everyone, I am a new user in openfoam. Currently I am running a case of dynamic mesh and I hope to run LES rather than RANS. It seems that I could not find LES in all kinds of dynamic models (interDyMFoam, pimpleDyMFoam..) Could anyone tell me if it is possible to run that case? Or could I download Openfoam-extended to do that?
Thank you very much.

[jinyaqing]

Hi, everyone. I am stilling considering the possibility of running Les in dynamic mesh. I saw in tutorial there are dynamics mesh cases with RANS. So I am trying to copy all the LES files from other cases and run dynamic mesh. It seems that the program does not report error, but I am not sure whether the results are correct...Could that be a possible approach to run the Les cases?
Sincerely Yaqing

[Jetfire]

Hi ,

I dont think you can solve LES with dynamicMesh
There is no solver that can solve this as of now.
which solver did you use?

[jinyaqing]

Hi, thank you for reply. In fact, at the beginning, I also thought that there is no way to do that as no tutorial file shows LES in dynamic mesh. However, I found several paper that apply LES in pimpleDyMfoam, and I am trying LES in interdymfoam. It seems that the program is running and gives me some right simulation results compared to experiments, although it sometimes collapses somewhere due to the sudden increase of Co number. But I am not sure if my LES is right as I just copied all Les files from other place and force it to run LES.

[wyldckat]

Greetings to all!

I'm using OpenFOAM 2.3 as a reference here. If you look at the tutorial "multiphase/interFoam/les/nozzleFlow2D", you'll see how a case can be configured for interFoam to use a LES turbulence model. The same set-up will then apply to the interDyMFoam case, with the small adjustment for the dynamic mesh details.

As for the broader question: almost any turbulence model can be chosen in most solvers. As for RAS or LES, it then really depends on whether it's applicable to your case or not. And study the tutorials for each solver in the "tutorials" folder, to see the differences between each situation, so that you can have a better perspective on this.
In addition, study the OpenFOAM User Guide: http://www.openfoam.org/docs/user/turbulence.php

Best regards,
Bruno

[jinyaqing]

Hi, thank you for your reply. Just as I expected, I copied the LES file to interDyMFoam and I see the simulation is running. I compared the results with experiments and it looks reasonable. However, for cases just like floating object (in fact I am doing falling object, but they should be very similar), I always get bounding k problems and simulation crushes there. I tried to use upwind scheme, limitedlinear and they are not working. I tried to change the tol residual and it also failed. Is there any better method to solve the bounding k problem. Thank you.

smoothSolver:  Solving for k, Initial residual = 0.99999979199, Final residual = 9.72474117538e-08, No Iterations 197
bounding k, min: -19455911245.9 max: 7.55487814355e+12 average: 14005690.2459
ExecutionTime = 40057.28 s  ClockTime = 40237 s

Interface Courant Number mean: 3.2670186072e-05 max: 0.0179877082776
Courant Number mean: 0.124426369256 max: 146743.539065
Time = 0.05041

PIMPLE: iteration 1

Centre of mass: (0.0499955921349 0.0500036220831 0.0832505228839)
Linear velocity: (0.0427515936267 -0.00935890129684 0.0682692127567)
Angular velocity: (-0.0493536141679 2.68556573617 0.0469246636873)
Execution time for mesh.update() = 0.540000000001 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 8.492982505e-07, No Iterations 48
time step continuity errors : sum local = 7.37054124342e-11, global = -6.40814789378e-11, cumulative = 1.20026128904e-06
smoothSolver:  Solving for alpha.water, Initial residual = 1.56439755255e-05, Final residual = 2.95605164287e-11, No Iterations 1
Phase-1 volume fraction = 0.740064988695  Min(alpha1) = 0  Max(alpha1) = 1.0000000001
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.740064988695  Min(alpha1) = 0  Max(alpha1) = 1.0000000001
GAMG:  Solving for p_rgh, Initial residual = 0.00474665629495, Final residual = 7.06910656876e-06, No Iterations 8
time step continuity errors : sum local = 0.000203369786927, global = -1.84153106314e-06, cumulative = -6.41269774101e-07
PIMPLE: iteration 2

Centre of mass: (0.0499957222774 0.0500035924211 0.0832507728155)
Linear velocity: (0.00260703226214 0.00702417422465 -0.0416072201563)
Angular velocity: (-0.688305748981 0.0605522618482 0.0578033106325)
Execution time for mesh.update() = 0.520000000004 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 8.42479300602e-07, No Iterations 47
time step continuity errors : sum local = 3.01678472303e-10, global = -2.60500843379e-10, cumulative = -6.41530274944e-07
smoothSolver:  Solving for alpha.water, Initial residual = 7.88100158276e-05, Final residual = 1.26086360883e-09, No Iterations 1
Phase-1 volume fraction = 0.740064988196  Min(alpha1) = 0  Max(alpha1) = 1.00000000026
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.740064988196  Min(alpha1) = 0  Max(alpha1) = 1.00000000026
GAMG:  Solving for p_rgh, Initial residual = 7.44569890105e-05, Final residual = 7.91054415732e-06, No Iterations 1
time step continuity errors : sum local = 0.65147750279, global = 0.000154140055944, cumulative = 0.000153498525669
PIMPLE: iteration 3

Centre of mass: (0.0499951906893 0.0500037859139 0.0832494241196)
Linear velocity: (0.521457059417 -0.381389915571 5.58662713013)
Angular velocity: (5.16393367775 -7.75232159694 2.57187237157)
Execution time for mesh.update() = 0.539999999994 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 9.25190165543e-07, No Iterations 34
time step continuity errors : sum local = 5.97806145586e-07, global = -5.17292350383e-07, cumulative = 0.000152981233318
smoothSolver:  Solving for alpha.water, Initial residual = 0.00328601719605, Final residual = 1.04371685399e-09, No Iterations 4
Phase-1 volume fraction = 0.740064988581  Min(alpha1) = 0  Max(alpha1) = 1.00000049184
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.740064988581  Min(alpha1) = -1.27800220843e-07  Max(alpha1) = 1.00000141625
GAMG:  Solving for p_rgh, Initial residual = 2.13289570658e-05, Final residual = 1.76935908146e-06, No Iterations 1
time step continuity errors : sum local = 5.29352545666, global = -0.000154374298359, cumulative = -1.39306504068e-06
PIMPLE: iteration 4

Centre of mass: (0.0500009107777 0.0499997082802 0.083307055159)
Linear velocity: (-0.59781948808 -2.96023303452 7.91018762314)
Angular velocity: (219.887434029 -42.8696415905 -15.5762561118)
Execution time for mesh.update() = 0.540000000001 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 7.79460508711e-07, No Iterations 36
time step continuity errors : sum local = 4.46906201536e-05, global = 3.86545640612e-05, cumulative = 3.72614990206e-05
smoothSolver:  Solving for alpha.water, Initial residual = 0.313181986022, Final residual = 0.000315153687459, No Iterations 1000
Phase-1 volume fraction = 0.638437340966  Min(alpha1) = 0  Max(alpha1) = 1.00000237053
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.638437340966  Min(alpha1) = -49802.7733395  Max(alpha1) = 10131.8287269
GAMG:  Solving for p_rgh, Initial residual = 5.49507292549e-05, Final residual = 4.43802042657e-06, No Iterations 2
time step continuity errors : sum local = 63.9433429684, global = 0.61171820884, cumulative = 0.611755470339
PIMPLE: iteration 5

Centre of mass: (0.0499839979238 0.0499779974827 0.0832726597245)
Linear velocity: (154970.308008 -121662.100577 164345.814281)
Angular velocity: (4688049.2481 4911152.6324 918300.367423)
Execution time for mesh.update() = 0.519999999997 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 9.19928549276e-07, No Iterations 30
time step continuity errors : sum local = 0.000262045146201, global = 0.000140923248851, cumulative = 0.611896393588
smoothSolver:  Solving for alpha.water, Initial residual = 0.99999970292, Final residual = 3.49114126684e-05, No Iterations 1000
Phase-1 volume fraction = -72.5156235818  Min(alpha1) = -106.716133272  Max(alpha1) = 1
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = -72.5156235819  Min(alpha1) = -5.96508281878e+12  Max(alpha1) = 315538337970
[0] #0  [1] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigHandler(int) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2   in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #5  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #6  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #7  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #7  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #8  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #8  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platfor in ms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #"/h9  ome/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #9  

[1]  at interDyMFoam.C:0
[1] #10  [0]  at interDyMFoam.C:0
[0] #10  

[0]  in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #12  [1]  in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/interDyMFoam"
[1] #11  __libc_start_main
 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #12  
[0]  in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/interDyMFoam"
[yaqing-HP-Z230-Tower-Workstation:10840] *** Process received signal ***
[yaqing-HP-Z230-Tower-Workstation:10840] Signal: Floating point exception (8)
[yaqing-HP-Z230-Tower-Workstation:10840] Signal code:  (-6)
[yaqing-HP-Z230-Tower-Workstation:10840] Failing at address: 0x3e800002a58
[yaqing-HP-Z230-Tower-Workstation:10840] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fd9c7ec04a0]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fd9c7ec0425]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fd9c7ec04a0]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 3] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d) [0x7fd9c915409d]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 4] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x71e) [0x7fd9c8fc167e]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 5] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarsestLevelERNS_5FieldIdEERKS2_+0x503) [0x7fd9c8fde653]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 6] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0xded) [0x7fd9c8fe051d]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 7] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x583) [0x7fd9c8fe2383]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 8] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7fd9cb7a3197]
[yaqing-HP-Z230-Tower-Workstation:10840] [ 9] interDyMFoam() [0x43acb2]
[yaqing-HP-Z230-Tower-Workstation:10840] [10] interDyMFoam() [0x42ec61]
[yaqing-HP-Z230-Tower-Workstation:10840] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fd9c7eab76d]
[yaqing-HP-Z230-Tower-Workstation:10840] [12] interDyMFoam() [0x43689d]
[yaqing-HP-Z230-Tower-Workstation:10840] *** End of error message ***
[1]  in "/home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/interDyMFoam"
[yaqing-HP-Z230-Tower-Workstation:10841] *** Process received signal ***
[yaqing-HP-Z230-Tower-Workstation:10841] Signal: Floating point exception (8)
[yaqing-HP-Z230-Tower-Workstation:10841] Signal code:  (-6)
[yaqing-HP-Z230-Tower-Workstation:10841] Failing at address: 0x3e800002a59
[yaqing-HP-Z230-Tower-Workstation:10841] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f6f48d194a0]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f6f48d19425]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f6f48d194a0]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 3] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d) [0x7f6f49fad09d]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 4] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x71e) [0x7f6f49e1a67e]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 5] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarsestLevelERNS_5FieldIdEERKS2_+0x503) [0x7f6f49e37653]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 6] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0xded) [0x7f6f49e3951d]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 7] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x583) [0x7f6f49e3b383]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 8] /home/yaqing/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7f6f4c5fc197]
[yaqing-HP-Z230-Tower-Workstation:10841] [ 9] interDyMFoam() [0x43acb2]
[yaqing-HP-Z230-Tower-Workstation:10841] [10] interDyMFoam() [0x42ec61]
[yaqing-HP-Z230-Tower-Workstation:10841] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f6f48d0476d]
[yaqing-HP-Z230-Tower-Workstation:10841] [12] interDyMFoam() [0x43689d]
[yaqing-HP-Z230-Tower-Workstation:10841] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 10841 on node yaqing-HP-Z230-Tower-Workstation exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
2 total processes killed (some possibly by mpirun during cleanup)

2 total processes killed (some possibly by mpirun during cleanup)

[jinyaqing]

Hi, everyone. I am currently trying to run a free falling object in openfoam. I basically referred to floating object tutorial and it worked for a while. However, it seems that after the object falls a while, the simulation crashed. I checked the reason and I think it is the mesh deformation. Some mesh at bottom is compressed too much and the thickness is almost 0. Therefore, I am looking for how can I use sixDoFRigidBodyMotion and also improve the mesh deformation.
I know that in moving cone tutorial the dynamic mesh is moving properly. but I just have no idea how to apply that scheme to sixDoFRigidBodyMotion case. How could I include something like movingconeTopoFvMesh and how could I edit that?
Cheers in advance and I appreciate any reply.
Yaqing

[jinyaqing]

Hi, everyone. I am trying to use layeradditionremoval for dynamic mesh. The simulation is free falling object in water. I have found some paper on internet writing that we can apply layeradditionremoval as:
right
{
type layerAdditionRemoval;
faceZoneName rightExtrusionFaces;
minLayerThickness 2e-4;
maxLayerThickness 5e-4;
oldLayerThickness -1;
active on;
}
But where shall I put this commend? I tried to include the layeradditionremoval.C file in dynamicmeshDict but it does't work...Anyone can provide me some example about how to use it? Thank you very much and cheer in advance.
Yaqing

[wyldckat]

Hi Yaqing,

Sorry for the late reply, but only now was I able to come back to your questions. I've also moved the latest two posts from their own threads, into this thread, because they are related to this one.

As far as I can figure out, the "layerAdditionRemoval" feature is something that is meant to be used with "engineFoam" or "engineDyMFoam", which are some somewhat old solvers.

I've done a bit more researching and this seems to be something that is only working in the foam-extend fork. Have a look at foam-extend 3.1 and the tutorial case "incompressible/icoDyMFoam/movingConeTopo".

I don't have enough time to look deeper into this. It might be possible to adapt the "movingBodyTopoFvMesh" class from foam-extend 3.1 to OpenFOAM 2.3.

Best regards,
Bruno