[dgadensg]

Hi all.

I have, almost succesfully, implemented the temperature equation in the interFoam solver for simulationg the filling stage of injection moulding.

My problem is that the solver crashes when the mold is almost filled. as shown on the picture below:



The end patches left and right are defined as zeroGradient for both velocity U, temperature T, and volumefraction alpha. The pressure p is defined as fixedvalue 0.

The error message is shown below:

Courant Number mean: 0.00623793 max: 0.495492
deltaT = 9.80392e-05
Time = 0.6

MULES: Solving for alpha1

Opnået 0.962465/98 pct. fyldning.

DILUPBiCG:  Solving for T, Initial residual = 0.00029924, Final residual = 2.10118e-08, No Iterations 1
DICPCG:  Solving for p, Initial residual = 2.94435e-05, Final residual = 8.68614e-08, No Iterations 56
time step continuity errors : sum local = 1.39854e-09, global = -1.18444e-12, cumulative = -0.00169153
ExecutionTime = 709.3 s  ClockTime = 721 s

End

#0  Foam::error::printStack(Foam::Ostream&) in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigSegv::sigSegvHandler(int) in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib/libc.so.6"
#3  Foam::IOobject::~IOobject() in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4  Foam::regIOobject::~regIOobject() in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5  Foam::twoPhaseMixture::~twoPhaseMixture() in "/home/dang/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/interTempFoam"
#6  main in "/home/dang/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/interTempFoam"
#7  __libc_start_main in "/lib/libc.so.6"
#8  _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/x86_64/elf/start.S:116
Segmentation fault

Hope someone can help me with an answer.

// Dan G

[herbert]

Hi Dan,

I'm a little bit confused about your error massage, because the error seems to happen after the simulation is finished (just because there is displayed "End" in front of the error message). Or did you add some code under this command?

[dgadensg]

Hi herbert.

I do not know why it prints out "end", because the simulation is programmed to stop as soon as the mold is 98 percent filled.

I have tried to run the same case with the standard interFoam solver, and it works fine.

the problem is also occuring for other cases where as soon as the polymer phase reaches an outlet, the solver stops.

// Dan