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

temperature blow-up: laser melting with icoReactingMultiPhaseInterFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2022, 07:49
Default temperature blow-up: laser melting with icoReactingMultiPhaseInterFoam
  #1
New Member
 
Jon
Join Date: Mar 2013
Posts: 13
Rep Power: 11
CHUIKOV is on a distinguished road
Dear all,

I have a divergence problem when modelling a laser melting case using icoReactingMultiPhaseInterFoam. OF version is ESI 2006.

It is a multiphase problem with solid, liquid and gas. The geometry is simple, it's a metal block with air (or argon gas) above. At time t=0 the bottom of the domain is 100% solid and the top is air. At t=0 a laser source is activated, which moves along from an end to the other of the block.

Physics includes a Voller Prakash type porosity term, Lee-type interphase mass transfer and a surface tension term. I use funkySetFields to set the inital alphas. As advised elsewhere in the Forum I'm using localDensityAbsorptionEmissionCoeffs so that the gas phase is transparent to the laser.

I'm interested in the temperature and shape once the metal has melted and solidified.

The problem I have is that the temperature explodes a few moments after the melting point of the metal is reached, it begins to grow uncontrollably until the calculation diverges.

I have tried changing discretization schemes, several boundary conditions, finer meshes but the temperature always explodes when I activate the term massTransfer in phaseProperties. With this term deactivated the problem converges.

I attach some of the case files and a paraFoam slice to show the problematic area (temp of order(1e5))) in case anybody can spot the problem. Link to the log file (be aware its 160Mb):

https://www.dropbox.com/s/2aen93q5u52hiv9/log.run?dl=0

Many thanks to all sharing their knowledge and experience on this topic.

cheers!

Jon


edit to add the last chunk of the log info:


Code:


fieldMinMax minMax write:
    min(T) = 25.46526211 in cell 394 at location (-0.0079 0 0.002) on processor 19
    max(T) = 177401.7077 in cell 242 at location (-0.0081 0 0.0005) on processor 18

limitFields limitLiquid:
: min(-0.09999716743): max(0.1357963488)
limitFields limitGas:
: min(1.188988269e-103): max(1.000000832)
Courant Number mean: 3.950596932e-06 max: 0.1010351509
Interface Courant Number mean: 6.219196043e-13 max: 5.850225857e-08
Maximum ddtAlpha : 0
Maximum DiffNum : 6.788873975e-108
deltaT = 3.56069572e-111
Time = 0.0436389805

smoothSolver:  Solving for alpha.gas, Initial residual = 3.097301169e-12, Final residual = 3.097301169e-12, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 2.423811728e-08, Final residual = 2.159409001e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 2.132342514e-12, Final residual = 2.132342514e-12, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.547563538e-13, Final residual = 6.547563538e-13, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 1.704959661e-08, Final residual = 4.132314019e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 1.031058202e-16, Final residual = 1.031058202e-16, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.426438812e-16, Final residual = 6.426438812e-16, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 4.149905504e-15, Final residual = 4.149905504e-15, No Iterations 0
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 7.03327967e-17, Final residual = 7.03327967e-17, No Iterations 0
MULES: Solving for alpha.solid
Phase-sum volume fraction, min, max = 1 0.9999999987 1.000000001
alpha.gas volume fraction = 0.2000019159  Min(alpha) = 0  Max(alpha) = 1.00000083
alpha.liquid volume fraction = 3.182570312e-06  Min(alpha) = -0.09999716682  Max(alpha) = 0.1357963487
alpha.solid volume fraction = 0.7999949016  Min(alpha) = 0  Max(alpha) = 1
PIMPLE: iteration 1
Move particles...1500
DILUPBiCG:  Solving for T, Initial residual = 4.029891791e-11, Final residual = 4.029891791e-11, No Iterations 0
min/max(T) = 32.13018728, 177401.7077
DICPCG:  Solving for p_rgh, Initial residual = 0.2135799472, Final residual = 9.54001272e-07, No Iterations 78
DICPCG:  Solving for p_rgh, Initial residual = 8.435986918e-07, Final residual = 6.205931562e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.1142211298, Final residual = 9.667296335e-07, No Iterations 75
DICPCG:  Solving for p_rgh, Initial residual = 8.943219162e-07, Final residual = 5.844954791e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.01102589883, Final residual = 9.61143165e-07, No Iterations 60
DICPCG:  Solving for p_rgh, Initial residual = 9.502350641e-07, Final residual = 6.188257697e-07, No Iterations 3
ExecutionTime = 11270.59 s  ClockTime = 13856 s

fieldMinMax minMax write:
    min(T) = 32.13018728 in cell 1417 at location (-0.0085 -0.0006 0.002) on processor 18
    max(T) = 177401.7077 in cell 242 at location (-0.0081 0 0.0005) on processor 18

limitFields limitLiquid:
: min(-0.09999716682): max(0.1357963487)
limitFields limitGas:
: min(1.188988269e-103): max(1.00000083)
Courant Number mean: 4.14147262e-06 max: 0.1085182211
Interface Courant Number mean: 6.456123611e-13 max: 6.072139692e-08
Maximum ddtAlpha : 0
Maximum DiffNum : 6.719318882e-108
deltaT = 3.281196176e-111
Time = 0.0436389805

smoothSolver:  Solving for alpha.gas, Initial residual = 3.062627957e-12, Final residual = 3.062627957e-12, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 2.34558473e-08, Final residual = 2.113657392e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 2.062629214e-12, Final residual = 2.062629214e-12, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.3312189e-13, Final residual = 6.3312189e-13, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 1.648793026e-08, Final residual = 4.003088614e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 1.194041634e-16, Final residual = 1.194041634e-16, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 5.918849759e-16, Final residual = 5.918849759e-16, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 5.889882691e-15, Final residual = 5.889882691e-15, No Iterations 0
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 8.865686765e-17, Final residual = 8.865686765e-17, No Iterations 0
MULES: Solving for alpha.solid
Phase-sum volume fraction, min, max = 1 0.9999999975 1.000000002
alpha.gas volume fraction = 0.2000019159  Min(alpha) = 0  Max(alpha) = 1.000000828
alpha.liquid volume fraction = 3.182570254e-06  Min(alpha) = -0.09999716624  Max(alpha) = 0.1357963486
alpha.solid volume fraction = 0.7999949016  Min(alpha) = 0  Max(alpha) = 1
PIMPLE: iteration 1
Move particles...1500
DILUPBiCG:  Solving for T, Initial residual = 3.932466412e-11, Final residual = 3.932466412e-11, No Iterations 0
min/max(T) = 25.46527599, 177401.7077
DICPCG:  Solving for p_rgh, Initial residual = 0.4000196739, Final residual = 9.924328698e-07, No Iterations 54
DICPCG:  Solving for p_rgh, Initial residual = 1.049551454e-06, Final residual = 5.188065951e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.05500427564, Final residual = 8.83073489e-07, No Iterations 48
DICPCG:  Solving for p_rgh, Initial residual = 8.551950879e-07, Final residual = 5.20441678e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.04381858212, Final residual = 9.920301417e-07, No Iterations 43
DICPCG:  Solving for p_rgh, Initial residual = 9.65572392e-07, Final residual = 5.328944926e-07, No Iterations 3
ExecutionTime = 11270.85 s  ClockTime = 13857 s

fieldMinMax minMax write:
    min(T) = 25.46527599 in cell 394 at location (-0.0079 0 0.002) on processor 19
    max(T) = 177401.7077 in cell 242 at location (-0.0081 0 0.0005) on processor 18

limitFields limitLiquid:
: min(-0.09999716624): max(0.1357963486)
limitFields limitGas:
: min(1.188988269e-103): max(1.000000828)
Courant Number mean: 3.950519745e-06 max: 0.1010351484
Interface Courant Number mean: 6.220016912e-13 max: 5.851177863e-08
Maximum ddtAlpha : 0
Maximum DiffNum : 6.191880788e-108
deltaT = 3.247578916e-111
Time = 0.0436389805

smoothSolver:  Solving for alpha.gas, Initial residual = 3.134744233e-12, Final residual = 3.134744233e-12, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 2.423335043e-08, Final residual = 2.15862556e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 2.132453075e-12, Final residual = 2.132453075e-12, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.545856046e-13, Final residual = 6.545856046e-13, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 1.704388079e-08, Final residual = 4.084817164e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 1.183278073e-16, Final residual = 1.183278073e-16, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.489081603e-16, Final residual = 6.489081603e-16, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 4.058992165e-15, Final residual = 4.058992165e-15, No Iterations 0
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 8.557479418e-17, Final residual = 8.557479418e-17, No Iterations 0
MULES: Solving for alpha.solid
Phase-sum volume fraction, min, max = 1 0.9999999991 1.000000002
alpha.gas volume fraction = 0.2000019159  Min(alpha) = 0  Max(alpha) = 1.000000826
alpha.liquid volume fraction = 3.182570192e-06  Min(alpha) = -0.09999716563  Max(alpha) = 0.1357963484
alpha.solid volume fraction = 0.7999949016  Min(alpha) = 0  Max(alpha) = 1
PIMPLE: iteration 1
Move particles...1500
DILUPBiCG:  Solving for T, Initial residual = 4.030030412e-11, Final residual = 4.030030412e-11, No Iterations 0
min/max(T) = 32.13022215, 177401.7077
DICPCG:  Solving for p_rgh, Initial residual = 0.2135857511, Final residual = 9.817938124e-07, No Iterations 78
DICPCG:  Solving for p_rgh, Initial residual = 8.681679619e-07, Final residual = 6.098755124e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.1142241114, Final residual = 9.311502075e-07, No Iterations 75
DICPCG:  Solving for p_rgh, Initial residual = 8.614046166e-07, Final residual = 6.011152458e-07, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.01102597145, Final residual = 9.635440573e-07, No Iterations 60
DICPCG:  Solving for p_rgh, Initial residual = 9.526068307e-07, Final residual = 6.098392712e-07, No Iterations 3
ExecutionTime = 11271.08 s  ClockTime = 13857 s

fieldMinMax minMax write:
    min(T) = 32.13022215 in cell 1417 at location (-0.0085 -0.0006 0.002) on processor 18
    max(T) = 177401.7077 in cell 242 at location (-0.0081 0 0.0005) on processor 18

limitFields limitLiquid:
: min(-0.09999716563): max(0.1357963484)
limitFields limitGas:
: min(1.188988269e-103): max(1.000000826)
Courant Number mean: 4.141393774e-06 max: 0.1085182239
Interface Courant Number mean: 6.456967843e-13 max: 6.073118882e-08
Maximum ddtAlpha : 0
Maximum DiffNum : 6.128442318e-108
deltaT = 2.992657639e-111
Time = 0.0436389805

smoothSolver:  Solving for alpha.gas, Initial residual = 3.060185192e-12, Final residual = 3.060185192e-12, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 2.345123846e-08, Final residual = 2.1129505e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 2.062700743e-12, Final residual = 2.062700743e-12, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 6.329373282e-13, Final residual = 6.329373282e-13, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 1.648253937e-08, Final residual = 4.011439035e-13, No Iterations 1
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 1.157772044e-16, Final residual = 1.157772044e-16, No Iterations 0
MULES: Solving for alpha.solid
smoothSolver:  Solving for alpha.gas, Initial residual = 5.904856583e-16, Final residual = 5.904856583e-16, No Iterations 0
MULES: Solving for alpha.gas
smoothSolver:  Solving for alpha.liquid, Initial residual = 5.905009211e-15, Final residual = 5.905009211e-15, No Iterations 0
MULES: Solving for alpha.liquid
smoothSolver:  Solving for alpha.solid, Initial residual = 8.504773186e-17, Final residual = 8.504773186e-17, No Iterations 0
MULES: Solving for alpha.solid
Phase-sum volume fraction, min, max = 1 0.9999999976 1.000000002
alpha.gas volume fraction = 0.2000019159  Min(alpha) = 0  Max(alpha) = 1.000000824
alpha.liquid volume fraction = 3.182570134e-06  Min(alpha) = -0.09999716504  Max(alpha) = 0.1357963483
alpha.solid volume fraction = 0.7999949016  Min(alpha) = 0  Max(alpha) = 1
PIMPLE: iteration 1
Move particles...1500
DILUPBiCG:  Solving for T, Initial residual = 3.932596598e-11, Final residual = 3.932596598e-11, No Iterations 0
[17] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[17] #1  Foam::sigFpe::sigHandler(int) at ??:?
[17] #2  ? in /lib/x86_64-linux-gnu/libc.so.6
[17] #3  Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[17] #4  Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[17] #5  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[17] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[17] #7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
[17] #8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[17] #9  ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/icoReactingMultiphaseInterFoam
[17] #10  ? in /lib/x86_64-linux-gnu/libc.so.6
[17] #11  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[17] #12  ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/icoReactingMultiphaseInterFoam
[jon-X11DPi-N-T:99970] *** Process received signal ***
[jon-X11DPi-N-T:99970] Signal: Floating point exception (8)
[jon-X11DPi-N-T:99970] Signal code:  (-6)
[jon-X11DPi-N-T:99970] Failing at address: 0x3e800018682
[jon-X11DPi-N-T:99970] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fe90d55c520]
[jon-X11DPi-N-T:99970] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7fe90d5b0a7c]
[jon-X11DPi-N-T:99970] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7fe90d55c476]
[jon-X11DPi-N-T:99970] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fe90d55c520]
[jon-X11DPi-N-T:99970] [ 4] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEENS_13scalarProductIT_S2_E4typeERKNS_5UListIS2_EES8_+0x46)[0x7fe90e2e7996]
[jon-X11DPi-N-T:99970] [ 5] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG11scalarSolveERNS_5FieldIdEERKS2_h+0x6ac)[0x7fe90e0cda4c]
[jon-X11DPi-N-T:99970] [ 6] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x38)[0x7fe90e0ce8e8]
[jon-X11DPi-N-T:99970] [ 7] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x18f)[0x7fe911f44b6f]
[jon-X11DPi-N-T:99970] [ 8] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x44b)[0x7fe9113d271b]
[jon-X11DPi-N-T:99970] [ 9] /home/jon/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x28)[0x7fe911378f08]
[jon-X11DPi-N-T:99970] [10] icoReactingMultiphaseInterFoam(+0x3afd6)[0x55cd2a129fd6]
[jon-X11DPi-N-T:99970] [11] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7fe90d543d90]
[jon-X11DPi-N-T:99970] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7fe90d543e40]
[jon-X11DPi-N-T:99970] [13] icoReactingMultiphaseInterFoam(+0x3de65)[0x55cd2a12ce65]
[jon-X11DPi-N-T:99970] *** End of error message ***
--------------------------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 17 with PID 0 on node jon-X11DPi-N-T exited on signal 8 (Floating point exception).
Attached Images
File Type: jpg Temperature_Blowup.jpg (25.8 KB, 15 views)
Attached Files
File Type: zip case.zip (12.7 KB, 5 views)

Last edited by CHUIKOV; August 22, 2022 at 09:18. Reason: add information
CHUIKOV is offline   Reply With Quote

Old   August 22, 2022, 08:43
Default
  #2
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 10
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Post the error message
Kummi is offline   Reply With Quote

Old   August 22, 2022, 09:18
Default
  #3
New Member
 
Jon
Join Date: Mar 2013
Posts: 13
Rep Power: 11
CHUIKOV is on a distinguished road
edited in the original post, many thanks!
CHUIKOV is offline   Reply With Quote

Old   August 22, 2022, 13:55
Default
  #4
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 10
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hi CHUIKOV,
Basic errors: It seems you messed up with some libraries.
As you can see from the errors in [jon-X11DPi-N-T:99970] [4] - [9]
I believe you are creating your customized solver in default OpenFOAM-v2006. If that's true, from now on copy the default solver in your user directory and make further modifications.
I couldn't able to judge other errors for now. If basic errors are solved, there might be a chance of disappearance of other errors.

So what to do now:
As you mentioned activating mass transfer results in complications, try rechecking the entire solver from the start where you modified. Modify your solver step by step, run the case, track the errors promptly.
Feel free to contact further.
Thank you
Kummi is offline   Reply With Quote

Old   August 23, 2022, 07:37
Default
  #5
New Member
 
Jon
Join Date: Mar 2013
Posts: 13
Rep Power: 11
CHUIKOV is on a distinguished road
hi again and thanks for looking into it!



Indeed I was trying to add some things to the solver but this is not the problem. I reverted back to the original solver and it still crashes, even tried the newest ESI version and same thing.


By using extremely small time steps O(1e-10) the solvers advances, but in this way it would take too long to get a result on my target domain.



I was wondering if any expert in multiphase problems can provide some insight?
CHUIKOV is offline   Reply With Quote

Old   August 24, 2022, 01:46
Default
  #6
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 10
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
If you try running the same case in debug mode, the errors will be located with line numbers, where it's easy for you to track down the mistakes. Try the link below:*
https://sites.google.com/site/foamgu...-debug-version
Secondly, use VS Code for operations like debugging and task running. This helps you navigate the location of errors while you run the case.

I'm not an expert, although let's arrange for a meeting and discuss it if nothing goes your way.

Thank you
Kummi is offline   Reply With Quote

Old   August 27, 2022, 08:47
Default
  #7
New Member
 
Jon
Join Date: Mar 2013
Posts: 13
Rep Power: 11
CHUIKOV is on a distinguished road
thanks again!


here an extract of the error after compiling OF in debug mode:


Code:

[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Tensor is not invertible due to the zero determinant:det(Tensor) = 1.129723342e-23
[1] 
[1]     From Foam::Tensor<Cmpt> Foam::inv(const Foam::Tensor<Cmpt>&, Cmpt) [with Cmpt = double]
[1]     in file /home/jon/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/lnInclude/TensorI.H at line 738.
[1] 
FOAM parallel run aborting
[1] 



15] #0  [29] #0  [13] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[23] #0  [31] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)
[11] #0  [6] #0  [2] --> FOAM FATAL ERROR: Foam::error::printStack(Foam::Ostream&)[3] #0  
Foam::error::printStack(Foam::Ostream&)[2] [14] #0  [1] #0  Foam::error::printStack(Foam::Ostream&)Tensor is not invertible due to the zero determinant:det(Tensor) = 1.8e-23
Foam::error::printStack(Foam::Ostream&)[2] 
Foam::error::printStack(Foam::Ostream&)[18] #0  [22] #0  [2] [7] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[30] #0      From Foam::error::printStack(Foam::Ostream&)[19] #0  [24] #0  [27] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[0] #0  [8] #0  Foam::Tensor<Cmpt> Foam::inv(const Foam::Tensor<Cmpt>&, Cmpt) [with Cmpt = double]Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)
[2] [16] #0  [9] #0  [17] #0  [25] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)    in file /home/jon/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/lnInclude/TensorI.H at line [26] #0  738Foam::error::printStack(Foam::Ostream&).
[2] 



Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::inv<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[13] #11   in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::reconstruct<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&)[23] #10  [jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
[jon-Super-Server:321145] PMIX ERROR: UNREACHABLE in file ../../../src/server/pmix_server.c at line 2193
 at ~/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/lnInclude/TensorI.H:741
[1] #5  Foam::Tensor<double> Foam::inv<double>(Foam::Tensor<double> const&) at ~/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/lnInclude/TensorI.H:754
[1] #6  Foam::inv(Foam::Field<Foam::Tensor<double> >&, Foam::UList<Foam::Tensor<double> > const&) at ~/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/fields/Fields/tensorField/tensorField.C:106 (discriminator 3)
[1] #7  void Foam::inv<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[1] #8  Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::inv<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[1] #9  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::reconstruct<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[1] #10  ? in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[1] #11  ? in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
[1] #12  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[1] #13  ? in ~/OpenFOAM/OpenFOAM-v2006-debug/platforms/linux64GccDPInt32Debug/bin/icoReactingMultiphaseInterFoam
 at ~/OpenFOAM/OpenFOAM-v2006-debug/src/OpenFOAM/lnInclude/TensorI.H:741
full error log attached.



any idea?



still have to learn how to properly debug the case.

cheers
Attached Files
File Type: txt error.txt (99.8 KB, 8 views)
CHUIKOV is offline   Reply With Quote

Old   August 28, 2022, 09:36
Default
  #8
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 10
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
A little help with compressibleInterFoam on OpenFoam2012
Have you checked this thread ?
This should help you to remove tensor computation errors.
Kummi is offline   Reply With Quote

Old   August 28, 2022, 18:58
Default
  #9
New Member
 
Jon
Join Date: Mar 2013
Posts: 13
Rep Power: 11
CHUIKOV is on a distinguished road
yes, already tried that suggestion, unfortunately the error remains
CHUIKOV is offline   Reply With Quote

Old   August 29, 2022, 13:04
Default
  #10
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 10
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Let's fix a meeting and discuss one day.
Kummi is offline   Reply With Quote

Reply

Tags
icoreactinmultiphaseif, laser melting

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
icoReactingMultiphaseInterFoam Laser DTRM OPFO OpenFOAM Running, Solving & CFD 42 August 19, 2022 02:29
Selective Laser Melting in Fluent Cooper24 FLUENT 0 June 23, 2021 06:38
Mapping bulk temperature of channel and moutend cylinder flow cases into field focus111 OpenFOAM Post-Processing 0 June 20, 2019 14:05
Melting too fast and not isothermal? Kronecker delta FLUENT 3 July 21, 2018 08:55
outlet temperature jigneshrohit99 FLUENT 1 March 25, 2016 14:26


All times are GMT -4. The time now is 04:08.