[Msoltanh]

Hi Foamer,


I am running OpenFoam-5x on Ubuntu.
I am trying to simulate reacting spray of acetone and I have already succeeded in simulating non-reacting acetone spray using sprayFoam solver.


So I modified sprayFoam/aachenBomb tutorial to apply my case,

First I applied detailed reaction mechanism for acetone in chemkin file as single reaction mechanism gives non convenient Temperature profile as compared with the experimental data.


Second I applied TDAC Model in chemistryProperties in order to have efficient simulation & save Cpu timing, so I have copied the TDAC Model from chemistryProperties of reactingFoam. And here the problem starts the simulation crashes after two time steps of injecting acetone.


a)I tried to solve this FoamError by changing the odecoeffs solver from seulex to rosenbrock34 as I thought applying 3rd and 4th order as this will help with stiff ordinary differential equations and also gave error.
b)I tried also to apply my case using chemFoam by applying just the reduction part , it gave also error so I have changed the initial set of species from CO & HO2 to CO2 &H2O and it worked. but when I tried it in sprayFoam solver gave the same error
c) I have also compiled the modified TDAC model for multiphase simulation which is uploaded in this bug and also the same error as I thought that the problem starts with the evaporation of acetone. https://bugs.openfoam.org/view.php?id=2753
d) I tried also to increase the total mass of acetone injection and increaing the duration of injection in order to increase the reaction of acetone and avoid any cooling down , but also chemFoam proved to me that the mechanism and amount of injection is convenient


after all my trials to solve the error the same error occurs.


I am afraid that TDAC coded in openfoam is not applied for lagrangian Model which leads to this error.
P.S I am using EDC combustion model & I have compiled acetone thermodynamic new liquid enthalpy in my solver (right enthalpy)



Courant Number mean: 0.02124973285 max: 0.1996819226
deltaT = 8.100804782e-06
Time = 4.02251051750598165


Solving 3-D cloud sprayCloud

Cloud: sprayCloud injector: model1
Added 6 new parcels

Cloud: sprayCloud
Current number of parcels = 7847
Current mass in system = 3.647952155e-06
Linear momentum = (9.694114795e-07 2.607167937e-08 9.903571468e-05)
|Linear momentum| = 9.904046255e-05
Linear kinetic energy = 0.001563244103
model1:
number of parcels added = 10993
mass introduced = 7.068583663e-06
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 258.3393594, 325.7480859
Mass transfer phase change = 3.420621223e-06
D10, D32, Dmax (mu) = 14.65702415, 19.9563246, 68.88376756
Liquid penetration 95% mass (m) = 0.389303597

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 0.0002886408108, Final residual = 9.344992496e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.0003361847375, Final residual = 6.681054853e-08, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1.986470185e-05, Final residual = 6.055225738e-08, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib64/libc.so.6"
#3 ? in "/lib64/libm.so.6"
#4 Foam::NonEquilibriumReversibleReaction<Foam::React ion, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >, Foam::ArrheniusReactionRate>::kf(double, double, Foam::Field<double> const&) const at ??:?
#5 Foam::TDACChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >:mega(Foam::Reaction<Foam::sutherlandTransport<Fo am::species::thermo<Foam::janafThermo<Foam:erfectG as<Foam::specie> >, Foam::sensibleEnthalpy> > > const&, Foam::Field<double> const&, double, double, double&, double&, int&, double&, double&, int&) const at ??:?
#6 Foam::TDACChemistryModel<Foam: siChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >:
mega(Foam::Field<double> const&, double, double, Foam::Field<double>&) const at ??:?
#7 Foam::TDACChemistryModel<Foam: siChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::derivatives(double, Foam::Field<double> const&, Foam::Field<double>&) const at ??:?
#8 Foam::Rosenbrock34::solve(double, Foam::Field<double> const&, Foam::Field<double> const&, double, Foam::Field<double>&) const at ??:?
#9 Foam::adaptiveSolver::solve(Foam::ODESystem const&, double&, Foam::Field<double>&, double&) const at ??:?
#10 Foam::ODESolver::solve(double&, Foam::Field<double>&, Foam::ODESolver::stepState&) const at ??:?
#11 Foam::ODESolver::solve(double, double, Foam::Field<double>&, double&) const at ??:?
#12 Foam:de<Foam::TDACChemistryModel<Foam:siChemistryM odel, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::solve(Foam::Field<double>&, double&, double&, double&, double&) const at ??:?
#13 double Foam::TDACChemistryModel<Foam: siChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:
erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
#14 Foam::combustionModels::EDC<Foam::combustionModels :
siChemistryCombustion>::correct() at ??:?
#15 ? at ??:?
#16 __libc_start_main in "/lib64/libc.so.6"
#17 ? at ??:?