Trying to create a new solver: phaseChangeHeatDyMFoam
 

I am trying to incorporate dynamic meshing with phaseChangeHeatFoam solver. The new solver is compiled by copying and imitating the files of interDyMFoam and interFoam. I made a new solver called phaseChangeHeatDyMFoam.

But, while running a test code i get this error:

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vect or<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#4 Foam::smoothInterfaceProperties::calculateK() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception (core dumped)


What can i to solve the problem.

Here is my boundary conditions folder:
https://drive.google.com/drive/folde...gX?usp=sharing

Quick question: Have you confirmed if the case runs with interDyMFoam or pimpleDyMFoam?

Because that's the first thing you should do, namely to have a case with which you can compare apples with apples, otherwise you are risking trying to solve 2 to 10 problems at the same time. The only guaranteed way to figure out how to solve each problem is if you are able to isolate each problem. And in this case, you need to first test if the case even works with the other solvers.


On the other hand, if you have already tested the case with the two other solvers and it worked, then you should have written indicating that as well. ;)

Hi,

I have tested a case by running it with interDyMFoam, where it works fine. But, It do not work in phaseChangeHeatDyMFoam.

Debug using gdb
 

Go to your solver's directory, type this in the terminal:

Then when you are running your test case, type
If you get an error, type
This will show you exactly where the problem is.