|
[Sponsors] |
July 23, 2021, 05:38 |
Problem with spray combustion case
|
#1 |
New Member
Anjul Pandey
Join Date: Jul 2020
Posts: 16
Rep Power: 6 |
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. |
|
Tags |
injection, openfoam 5.x, spray |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |