CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Problem with spray combustion case

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 23, 2021, 05:38
Default Problem with spray combustion case
  #1
New Member
 
Anjul Pandey
Join Date: Jul 2020
Posts: 16
Rep Power: 6
anjul is on a distinguished road
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.
anjul is offline   Reply With Quote

Reply

Tags
injection, openfoam 5.x, spray

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Restarting case problem - transient introduced Legwan FLUENT 2 January 1, 2021 03:58
ECFM 3Z combustion problem Bücherknopf STAR-CD 12 September 15, 2020 04:39
Coal Combustion Problem svsadatia CFX 0 March 7, 2016 04:33
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
Problem w/ liquid spray combustion(no reaction) faudiS8 FLUENT 1 June 7, 2012 13:20


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