Rereading of ReactingParcels fails if one of the Phases has a fraction of 0

gschaider · January 7, 2009, 16:10

In ReactingParcelIO.C there is an error that makes it fail if reading data from a previous run where some of the parcels have for instance a Y_liquid of 0 (division by zero). This patch fixes this:

ReactingParcelIO.patch

BTW: I haven't tested it, but it occurred to me during fixing this that if there was for instance a species H2O in the liquid and the solid phase of the parcel they both would get written to lagrangian/YH2O. I'd propose calling them YlH2O and YsH2O just to avoid such clashes (however unlikely they might be)

gcollecutt · January 8, 2009, 20:27

Bernhard,

The reacting parcel class was originally written for modelling combustion of coal particles. The class was intended to allow the parcels to have a liquid component but the calcMassTransfer() routine is not yet able to support this - hence the class was never tested with liquid components.

Have you included some form evaporation/condensation in a mass transfer model?

Greg.

gschaider · January 12, 2009, 11:13

Hi Greg!

Yes, I'm working on an evaporation model included. But the error also occurs for Y_gas and Y_solid equal to 0. I was only using Y_liquid as an example.

Bernhard

gcollecutt · January 12, 2009, 20:37

Bernhard,

How have you structured the inclusion of the evaporation model? The massTransfer() routine was designed to call run time selectable devolatilisation models. I think it would be unwise to inlcude the liquid evaporation physics into the devolatilisation model - it should go in as its own model. In this case the parcelProperties dictionary should specify a both a devolatilisation model (or NoDevolatilisation) and an evaporation model (or NoEvaporisation). The reacting parcel template should likewise be modified to call these to models seperately and get rid of the generic call to massTransfer() which already is misleading since the surface reaction model also causes mass transfer.

Greg.

gschaider · January 13, 2009, 07:17

Hi Greg!

You're right. I rewrote massTransfer for my parcels and introduced a separate EvaporationModel-hierarchy.

Bernhard

gcollecutt · January 13, 2009, 14:29

Hi Berhhard,

I have been working with one of the guys at OpenCFD to develop the reacting parcel class. Is there a chance we could have a look at your code?

Greg.

gschaider · January 13, 2009, 16:34

Hi Greg!

Got to check with our client whether it's OK if I removed the customer-specific stuff and gave it to you.

Bernhard

January 7, 2009, 16:10	In ReactingParcelIO.C there is	#1
gschaider Assistant Moderator Bernhard Gschaider Join Date: Mar 2009 Posts: 4,225 Rep Power: 51	In ReactingParcelIO.C there is an error that makes it fail if reading data from a previous run where some of the parcels have for instance a Y_liquid of 0 (division by zero). This patch fixes this: ReactingParcelIO.patch BTW: I haven't tested it, but it occurred to me during fixing this that if there was for instance a species H2O in the liquid and the solid phase of the parcel they both would get written to lagrangian/YH2O. I'd propose calling them YlH2O and YsH2O just to avoid such clashes (however unlikely they might be) __________________ Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

January 8, 2009, 20:27	Bernhard, The reacting parc	#2
gcollecutt New Member Greg Collecutt Join Date: Mar 2009 Location: Brisbane, Queensland, Australia Posts: 21 Rep Power: 17	Bernhard, The reacting parcel class was originally written for modelling combustion of coal particles. The class was intended to allow the parcels to have a liquid component but the calcMassTransfer() routine is not yet able to support this - hence the class was never tested with liquid components. Have you included some form evaporation/condensation in a mass transfer model? Greg.

January 12, 2009, 11:13	Hi Greg! Yes, I'm working o	#3
gschaider Assistant Moderator Bernhard Gschaider Join Date: Mar 2009 Posts: 4,225 Rep Power: 51	Hi Greg! Yes, I'm working on an evaporation model included. But the error also occurs for Y_gas and Y_solid equal to 0. I was only using Y_liquid as an example. Bernhard __________________ Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

January 12, 2009, 20:37	Bernhard, How have you stru	#4
gcollecutt New Member Greg Collecutt Join Date: Mar 2009 Location: Brisbane, Queensland, Australia Posts: 21 Rep Power: 17	Bernhard, How have you structured the inclusion of the evaporation model? The massTransfer() routine was designed to call run time selectable devolatilisation models. I think it would be unwise to inlcude the liquid evaporation physics into the devolatilisation model - it should go in as its own model. In this case the parcelProperties dictionary should specify a both a devolatilisation model (or NoDevolatilisation) and an evaporation model (or NoEvaporisation). The reacting parcel template should likewise be modified to call these to models seperately and get rid of the generic call to massTransfer() which already is misleading since the surface reaction model also causes mass transfer. Greg.

January 13, 2009, 07:17	Hi Greg! You're right. I re	#5
gschaider Assistant Moderator Bernhard Gschaider Join Date: Mar 2009 Posts: 4,225 Rep Power: 51	Hi Greg! You're right. I rewrote massTransfer for my parcels and introduced a separate EvaporationModel-hierarchy. Bernhard __________________ Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

January 13, 2009, 14:29	Hi Berhhard, I have been wo	#6
gcollecutt New Member Greg Collecutt Join Date: Mar 2009 Location: Brisbane, Queensland, Australia Posts: 21 Rep Power: 17	Hi Berhhard, I have been working with one of the guys at OpenCFD to develop the reacting parcel class. Is there a chance we could have a look at your code? Greg.

January 13, 2009, 16:34	Hi Greg! Got to check with	#7
gschaider Assistant Moderator Bernhard Gschaider Join Date: Mar 2009 Posts: 4,225 Rep Power: 51	Hi Greg! Got to check with our client whether it's OK if I removed the customer-specific stuff and gave it to you. Bernhard __________________ Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Volume fraction to mass fraction	Torsaki	CFX	10	March 10, 2018 03:41
equilibrium between phases	chouki	FLUENT	1	August 27, 2008 17:29
How set volume fraction and species mass fraction	Reginaldo Cotto	FLUENT	0	August 31, 2006 04:58
RunFoamX fails	jvn	OpenFOAM Installation	12	August 11, 2005 11:12
SOR Fails Near Walls	Chih-Wei Chiu	Main CFD Forum	0	May 9, 2002 01:05