CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   BC for temperature added to interFoam (

dgadensg April 22, 2010 08:47

BC for temperature added to interFoam
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


#0  Foam::error::printStack(Foam::Ostream&) in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#1  Foam::sigSegv::sigSegvHandler(int) in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#2  ?? in "/lib/"
#3  Foam::IOobject::~IOobject() in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#4  Foam::regIOobject::~regIOobject() in "/home/dang/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#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/"
#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 April 22, 2010 10:27

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 April 22, 2010 15:16

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

All times are GMT -4. The time now is 23:53.