[mturcios777]

From what I remember from the spray library, it uses the cell temperature. For each parcel, we are working under the assumption that the parcel is much smaller than the cell that contains it. Because the spray parcels will exchange momentum and energy with the ambient air, the temperature will change.

[conceptone]

Ok, I checked the temperature before it cracks(at time=0.00149s) and the temperature below is at 0.00144s, when the time sub-file output last time. There is no <=0 temperature in this file, they are nearly 900k, the ambient temperature:
...
899.999
899.999
899.999
899.999
899.998
899.998
899.998
899.998
899.997
899.997
899.997
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.996
899.997
899.997
899.997
899.997
899.998
899.998
899.998
899.998
899.999
899.999
899.999
899.999
900
900
900
900
900
900
900
900
)
;

boundaryField
{
walls
{
type zeroGradient;
}
}

and so is the pressure, they are nearly the ambient pressure 6e+6:

5.99976e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99978e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99977e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
5.99976e+06
)
;

boundaryField
{
walls
{
type zeroGradient;
}
}
Could you give me some suggestions?

[mturcios777]

At this point it would be easiest for you to send me a copy of your case, as there could be all kinds of things going on here, possibly in sprayProperties file, themophyiscalProperties, the mesh or the spray location. Post it here if you can, or email if it is too big.

[conceptone]

Originally Posted by mturcios777

At this point it would be easiest for you to send me a copy of your case, as there could be all kinds of things going on here, possibly in sprayProperties file, themophyiscalProperties, the mesh or the spray location. Post it here if you can, or email if it is too big.

I'd like to thank your warm help at first. I have attached the files. The combustion vessel geometry is 108mm cube and spray location (0.054,0.054,0.103). The input data provided from the website:http://www.sandia.gov/ecn/cvdata/dsearch/frameset.php and http://www.cmt.upv.es/ECN01.aspx

Regards

[achinta]

Hello Marco and conceptone,
Did you guys figure out what the error is? I am facing a similar error

Kind regards,
Achinta

[mturcios777]

Originally Posted by achinta

Hello Marco and conceptone,
Did you guys figure out what the error is? I am facing a similar error

Kind regards,
Achinta

Sorry Achinta, I've been busy with other work and have just started looking at the problem.

conceptone, I noticed that you haven't included your chemkin folder. If you could send it that would be great.

[conceptone]

Originally Posted by mturcios777

Sorry Achinta, I've been busy with other work and have just started looking at the problem.

conceptone, I noticed that you haven't included your chemkin folder. If you could send it that would be great.

Sorry about I forgetting to include that file. I'll attach again. I guess maybe I resolve the problem, which is not caused by the input parameters.

Another problem, actually more general problem, so how can I know the meaning of output files calculated by openfoam? I mean like in sprayFoam, in the some time file, like in 0.0012/lagrangian/sprayCloud/, what's the meaning of results in mass0 file, or positions file? I mean how can I know it by doing sth.?
Thanks

[conceptone]

Originally Posted by achinta

Hello Marco and conceptone,
Did you guys figure out what the error is? I am facing a similar error

Kind regards,
Achinta

Hi,actually, it's a little ... you should use the debug version to see what's wrong. If you can't build the debug version, you can post your error report and I'll check whether it's the same as mine.

You can just check the T file(temperature)in the time file which is close to the time collapse to see whether it exceeds the critical temperature of C12H26. My case is just so. If it's so, you got to modify the code of breakup model, it's finally demonstrated by this.