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

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

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 7, 2009, 17:10
Default In ReactingParcelIO.C there is
  #1
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
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
gschaider is offline   Reply With Quote

Old   January 8, 2009, 21:27
Default Bernhard, The reacting parc
  #2
New Member
 
Greg Collecutt
Join Date: Mar 2009
Location: Brisbane, Queensland, Australia
Posts: 21
Rep Power: 8
gcollecutt is on a distinguished road
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.
gcollecutt is offline   Reply With Quote

Old   January 12, 2009, 12:13
Default Hi Greg! Yes, I'm working o
  #3
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
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
gschaider is offline   Reply With Quote

Old   January 12, 2009, 21:37
Default Bernhard, How have you stru
  #4
New Member
 
Greg Collecutt
Join Date: Mar 2009
Location: Brisbane, Queensland, Australia
Posts: 21
Rep Power: 8
gcollecutt is on a distinguished road
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.
gcollecutt is offline   Reply With Quote

Old   January 13, 2009, 08:17
Default Hi Greg! You're right. I re
  #5
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
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
gschaider is offline   Reply With Quote

Old   January 13, 2009, 15:29
Default Hi Berhhard, I have been wo
  #6
New Member
 
Greg Collecutt
Join Date: Mar 2009
Location: Brisbane, Queensland, Australia
Posts: 21
Rep Power: 8
gcollecutt is on a distinguished road
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.
gcollecutt is offline   Reply With Quote

Old   January 13, 2009, 17:34
Default Hi Greg! Got to check with
  #7
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
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
gschaider is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Volume fraction to mass fraction Torsaki CFX 9 February 21, 2015 11:46
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


All times are GMT -4. The time now is 14:15.