CFD Online Discussion Forums - 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:

http://www.cfd-online.com/Forums/mem...rash-point.jpg

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:

Code:

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