[anjul]

Hello everyone,
I am trying to simulate a spray flame. I completed the case setup and when I ran it, it crashed.
tauRelaxTarget = 10
tauRelaxDelta = 1
tauRealax = 4.83408618186e-05
DILUPBiCGStab: Solving for f, Initial residual = 6.45931836899e-07, Final residual = 1.49710628766e-09, No Iterations 1
Execution Time for dens coupling= 0 s
DILUPBiCGStab: Solving for O2, Initial residual = 2.82605568944e-05, Final residual = 2.78121903272e-10, No Iterations 1
DILUPBiCGStab: Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CO, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CH2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for C2H4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for C2H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CH3CHO, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for C2H5OH, Initial residual = 2.61917655922e-06, Final residual = 4.45483070883e-09, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 3.23813837811e-06, Final residual = 5.8483957574e-09, No Iterations 1
GAMG: Solving for p, Initial residual = 4.4176940908e-05, Final residual = 4.34979798893e-07, No Iterations 3
GAMG: Solving for p, Initial residual = 1.27090201634e-06, Final residual = 9.64216730093e-09, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.20757251109e-13, global = 4.54238314028e-14, cumulative = -6.63505497962e-11
bounding Cs, min: -23.4296830462 max: 30.8757728882 average: 0.0927099029716
Dsgs max/min :0.000613419518875 0
bounding k, min: -2.61631518873 max: 5.30127090894 average: 0.0171170682025
bounding Cs, min: -23.4296830462 max: 30.8757728882 average: 0.0927099029716
LES ExecutionTime = 1235.72 s
Moving Particles
SP -> Inflow
3 particles inserted at inflow boundaries.
SP -> Mixing
SP -> Solve ODE
447266 particles moved. 0 wall hits. 54 rebounds off noOutflow patches. 0 rebounds off noRandomOutflow patches. 3 particles left the domain. 0 particles hit periodic BCs. 2573 particles passed between processors.
0 wedge hits
0 symmetry Plane hits
Particle number & weight control
Mass in fv domain is 0.00750683569203. Mass in particle field is 0.00747006759734
using finiteRate Chemistry
SP -> Update composition
Moving Fuel Particles from solver

Solving 3-D cloud fuelCloud1

Cloud: fuelCloud1 injector: model1
Added 9 new parcels

[1] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 ? in "/lib64/libc.so.6"
[1] #3 ? at ??:?
[1] #4 ? at ??:?
[1] #5 ? at ??:?
[1] #6 ? at ??:?
[1] #7 ? at ??:?
[1] #8 ? at ??:?
[1] #9 ? at ??:?
[1] #10 main at ??:?
[1] #11 __libc_start_main in "/lib64/libc.so.6"
[1] #12 ? at ??:?

I am using a modified form of an already available solver in OpenFOAM which was developed by my seniors. It is already tested so I don't think it has any issues. The details of the injection model which I gave are:
injectionModels
{
model1
{
type ConditionalDropletInjection;
SOI 0;
massTotal 4.72e-04;//[kg/s]
parcelBasisType fixed;
injectionMethod disc;
flowType ConditionalDropletVelocity;
outerDiameter 2.0e-03;
innerDiameter 0;
duration 1;
position (1.0e-06 0 0);
direction (1 0 0);
parcelsPerSecond 1900000;
nParticle 1;
flowRateProfile constant 1;
Cd constant 5.3649e-03;
thetaInner constant 0.0;
thetaOuter constant 0.0;
U_bulk 35.69;//Calculated

sizeDistribution
{
type RosinRammler;

RosinRammlerDistribution
{
minValue 5e-06;
maxValue 80e-06;
d 40e-06;
n 3;
}
}
}
}
I have also modified the injection model in OpenFOAM which I named as ConditionalDropletInjection and associated velocity as ConditionalDropletVelocity. These are the modifications I did in ConeNozzleInjection model in setProperties part:
//Edited
parcel.d() = sizeDistribution_->sample();
Info << "Parcel Diameter = " << parcel.d() << nl;
scalar lam = rndGen.sample01<scalar>();

scalar theta_minp = 20 + 15 * (parcel.d() - 60*10e-06)/(40e-06);
scalar theta_maxp = 50 - 5 * (parcel.d() - 60*10e-06)/(40e-06);
scalar coneAngle;

if (lam >= 0 && lam < 0.5)
coneAngle = theta_minp + 2*lam * (40 - theta_minp);
else
coneAngle = 40 + (2*lam - 1) * (theta_maxp - 40);
//Editing upto here

coneAngle *= deg2Rad;
Info << "ConeAngle = " << coneAngle << nl;
scalar alpha = sin(coneAngle);
scalar dcorr = cos(coneAngle);

vector normal = alpha*normal_;
vector dirVec = dcorr*direction_;
dirVec += normal;
dirVec /= mag(dirVec);
Velocity of droplets is calculated as:
case ftConditionalDropletVelocity:
{
scalar eta_p = mag(1 - 2*lam);
scalar Umag = (1 - 2 * 0.05 * eta_p + 0.05) * U_bulk;
parcel.U() = Umag*dirVec;
Info << "Parcel velocity = " << parcel.U() << nl;
break;
}

Any idea what can be the issue. Any help would be appreciated.