CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Having problems setting up surface reaction, catalyst, adsorption and desorption

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

Like Tree10Likes
  • 3 Post By Juun
  • 5 Post By villager
  • 1 Post By villager
  • 1 Post By grapheeene

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 31, 2014, 13:14
Question Having problems setting up surface reaction, catalyst, adsorption and desorption
  #1
New Member
 
Juun
Join Date: Aug 2013
Location: UK
Posts: 12
Rep Power: 13
Juun is on a distinguished road
(sorry, this is long, just tried to describe my problem clearer)

Hi, everyone

I am doing a simulation on catalyst reaction with ammonia (NH3). The adsorption and desorption rate is important.

I stuck at setting up surface adsorption and desorption reactions. I am not sure about the meaning of ‘selected site species’, so I just copy 2nd NH3 gas from fluent database and add to this site species (with name NH3-S), (the 1st NH3 gas is added to ‘selected species’ and used in inlet mixture).

Then, I set 2 reactions (backward reaction would do it, but I have different Arrhenius rate for each)

Adsorption NH3 = NH3(S)
Desorption NH3(S)=NH3

My first question(s):
Is this the way to set up ad/desorption on the surface in Fluent?
How to use the site-species correctly? (Because I not sure about my understanding, maybe it should be the surface material of the catalyst??)

Then, I want to include surface coverage dependent in desorption reaction, so I tick the option (coverage dependent reaction) at the bottom.

To test it, (only desorption reaction is selected), inlet NH3 is set to 0, and in the reaction mechanisms-site density-initial coverage NH3 is set to 1. I expect to see a small amount of NH3 is released from the surface, but the reaction did not happen during calculation.

However, when I turn off coverage dependent reaction; set initial coverage NH3 to 0, the result shows me NH3 is released from surface, which should not happen.

It seems like the option of coverage dependent reaction controls the ‘on’ and ‘off’ of the desorption reaction, no matter what values is given in initial coverage.

My expectation for this desorption reaction is: surface releases NH3 only when surface coverage > 0, and it stop when surface coverage = 0. The coverage-dependent reaction should do it (by looking at the equation), but I don’t know where is wrong in my model.

The second question(s): why is this happen? Do I have to use UDF to set up anything? Or any other ways can get to the expectation?

Thank you very much for your time and help
silverra1n, jumaah and granzer like this.
Juun is offline   Reply With Quote

Old   April 1, 2014, 09:52
Smile
  #2
Member
 
vlg
Join Date: Jul 2011
Location: My home :)
Posts: 81
Rep Power: 18
villager is on a distinguished road
Good day!

First, when you turn on the Species Model, FLUENT automatically creates the mixture. Mixture is comprised of several species (gas/fluid phase) and site species. Selected (that means, you select only that you need, you can add a lot of species to materials) site species are the species that are present on catalyst surface. All that you should specify is species molar mass (also, entropy and enthalpy of site species could be specified, but they are, generally, not known). When copying from existing material molar mass is duplicated. So, adsorption/desorption mass conservation seems to be satisfied.
But the second thing to account for is the conservation of number of sites. In your reaction NH3(s) suddenly "appears" from gas. How much NH3(s) can be adsorbed on surface?
FLUENT provides a model (see Theory Guide or "Surface Chemkin manual" in the web for explanation) that treats catalyst like a surface with constant site density - a number of catalytic (active) sites (in kmol) per m2 of surface (let it be active surface).
To specify this, you go into Materials -> your mixture -> Mechanism -> number of sites (1) -> site density (kmol/m2).
Only monolayer adsorption (true for chemisorption, chemical bonding occurs) is under consideration (see "g").
Then for the conservation of site species number, we must write (e.g. take Platinum):
NH3 + Pt(s) => NH3(s)
Molar mass of NH3(s) now should be the sum of NH3 m.m. and Pt m.m. (I think, Pt molar mass can be set arbutrary unless you simulate an etching)
Desorption is described as the reverse reaction.

Second, set a plausible rate parameters (analyze literature, try this article at first).

At third, do you simulate wall surface reaction (catalytic wall - e.g. tube with coating)? Then you should specify coefficient (on wall b.c.) - m2 of active surface per m2 of you wall (look for Surface Area Washcoat Factor).

At fourth, turning on "coverage dependent reaction" must be assisted with additional parameters input ("specify" button will be active).

And at last, when we do a steady-state (I suppose, you do it) calculation, we want to find a steady-state solution. Initial conditions then (in "mechanism" - initial site coverage) are only initial guess for the solution. You can specify 0 for some species and have e.g. 0.7 on your catalyst after simulation. The behaviour of the model, though, can be dependent on initial guess quality.
Note: Check the sum of initial coverages is 1.

Best wishes, Villager.
silverra1n, tmu, Juun and 2 others like this.
villager is offline   Reply With Quote

Old   April 1, 2014, 14:09
Default
  #3
New Member
 
Juun
Join Date: Aug 2013
Location: UK
Posts: 12
Rep Power: 13
Juun is on a distinguished road
Hi, Villager, Thank you so much for your help. The answer is very informative and detail.


My problem is a square tube, the complete inside wall is coated with active material (just like a normal catalyst), and mixture goes in from one side, goes out from the other.


Therefore I set Surface Area Washcoat Factor to 1, then "coverage dependent reaction" I used the simplest one, set ‘Mu’ to 1.


But I have further question about the conservation of site species number, as you say: NH3 + Pt(s) => NH3(s)


Is Pt(s) site species? Can I take it like the surface coating? So, I only need to entry a value to absorption capacity (site density) on this site.
Then NH3(s) is the second site species, I guess I don’t need to set site density for this one, as it is a product on the surface. Is my thought correct so far?


The problem here is the molar mass, if NH3 is 1, Pt(s) is 1, and NH3 should be 2.

Will this affect the later reactions, (e.g. NH3(s)+A; NH3(s)+B …)?? because the NH3(s) is 2 molar now, then more reactions happens.
When I set NH3 to 1, Pt(s) to 1, NH3(s) to 1, FLUENT says mass imbalance.
I’m confused here.

At the moment, it is steady state, when this works, I’ll move to transient to see surface coverage % vs time (also with temperature variation). But right now, the surface coverage is either 1 or 0. somewhere is wrong but I havn't figure out yet : (
Juun is offline   Reply With Quote

Old   April 1, 2014, 15:20
Post
  #4
Member
 
vlg
Join Date: Jul 2011
Location: My home :)
Posts: 81
Rep Power: 18
villager is on a distinguished road
Catalytic coating is a porous structure, therefore, the surface usually evolves (increase) after coating compared with the wall area. Chemisorption data for active surface determination may be taken into account.

Yes, Pt(s) is site species - imagine e.g. in initial state all surface is free of species - then the surface coverage of Pt(s) (free sites! - but it's indifferently for the model) is equal to 1.
Site density is set for the whole site! Don't mix site and species concepts. This model understands sites as some places (points) on catalyst. Every site can be filled with some species or it can be void.
Different sites = Pt and Ni e.g., if you have both (mixed catalyst). Then different species on different sites (like NH3(s) and NH3(s2) ). I think, it's not your case.

NH3(s) is NH3 linked with platinum, so its "molar mass" is equal to NH3 mass + Pt mass. I checked different m.m. for Pt (both consistent with NH3(s) and not), it seems like having no effect on solution. M.m. are important when simulating etching/deposition. But I prefer to maintain stoichometry balances - for Pt(s) I would input its real m.m. and so on.
but when you import kinetic mechanism from CHEMKIN files NH3 molar mass is important if you use sticking parametrization of adsorption (
see Theory Guide or "Surface Chemkin manual" for more info on that).

Try to move reaction parameters.
Juun likes this.
villager is offline   Reply With Quote

Old   April 3, 2014, 12:16
Default
  #5
New Member
 
Juun
Join Date: Aug 2013
Location: UK
Posts: 12
Rep Power: 13
Juun is on a distinguished road
Hi, Villager. I am so grateful that you have taken time to help me with the problems. The result is close to my expectation now.
Despite there is an imbalance in the reactions (the strange thing is that, when I use 1 molar mass for NH3 1, Pt(s) and NH3(s) on my laptop’s model, there is no error saying mass imbalance, unless I use 1, 1, 2 or something else. However, when it comes the desktop, the error appears in both case, 1,1,1 and 1,1,2, all other settings are the same. I just ignored it and run calculation, the result looks OK.) Not sure why, but adsorption works. Thank you.
Juun is offline   Reply With Quote

Old   November 21, 2014, 03:50
Default
  #6
New Member
 
Zaragoza
Join Date: Nov 2014
Posts: 4
Rep Power: 11
asfrallego is on a distinguished road
Hello Juun and villager.
I have a problem similar to described by Juun. I want to simulate CO2 adsorption with amine: CO2(bulk)+2DEA(site)-->product1(site)+product2(site). Im doing a transient simulation, the reaction rate looks working fine, DEA concentration reduces with time as products increases. Heat due the adsorption seems to work well.
The problem is that CO2 does not flow when transient+reaction is togheder running. I can stablish a CO2 concentration at fluid volume, but it does not change with time neither varies its concnetration (seems to be 'freezed').
I hope you both can help me with that issue.
asfrallego is offline   Reply With Quote

Old   November 21, 2014, 09:03
Default
  #7
Member
 
vlg
Join Date: Jul 2011
Location: My home :)
Posts: 81
Rep Power: 18
villager is on a distinguished road
Quote:
Originally Posted by asfrallego View Post
Hello Juun and villager.
I want to simulate CO2 adsorption with amine: CO2(bulk)+2DEA(site)-->product1(site)+product2(site).
Do you have an aBsorption or an aDsorption process?
MEA/DEA are traditionally used for aBsorption of CO2.
Adsorption of CO2 is also possible with amine-based adsorbents.
http://aaqr.org/vol12_no5_october201...32_745-769.pdf
So, you have, as I understand, you consider amine-based solid adsorbent with active surface for capturing CO2. Am I right?
villager is offline   Reply With Quote

Old   November 24, 2014, 02:57
Default
  #8
New Member
 
Zaragoza
Join Date: Nov 2014
Posts: 4
Rep Power: 11
asfrallego is on a distinguished road
Thanks for your reply villager.

The process is an adsorption (amine impregnated adsorbent). My support is a gamma-alumina with DEA. The paper you send me I have read before, but not much newest things it says.
I try to simulate the alumina support as a porous media and the DEA as a surface coverage specie. i have all the input data to feed Fluent (still waiting for real Arrhenius paramethers from the lab) but this seems to work all properly.

But my problem is that I cannot simulate at the same time the surface reaction and the flow for CO2 (and other non-reacting species). I have tried to change the products form site species to solid species (to observ the rate of deposition) but this does not change the problem.
If you have any new suggestion, I am waiting to read it!
asfrallego is offline   Reply With Quote

Old   June 7, 2015, 07:11
Default
  #9
New Member
 
Deependra
Join Date: Jun 2015
Location: Munich
Posts: 10
Rep Power: 11
Deependra is on a distinguished road
Hello Villager,

I am working on a same kind of problem that is as follows:

There is a cylindrical chamber in which some volatile compound being injected from inlet and there is wall inside of this chamber that is made of some material on which I need to simulate Adsorption and Desorption of this volatile compund.

Till this material continues adsorbing VOC compound or act as a sink, injecting VOC continues and once it reaches to Equilibrium when there is no Desorption and Adsorption or it does not act as source or sink,we stop injecting VOC and inject fresh air.So Desorption starts from wall.

And this condition I have to defined as user defined function on this wall material surface as boundary condition.

This whole procedure I need to do in Ansys Fluent.So if you have any idea about how to initiate and proceed for this Problem,Please share your ideas.

Thanks for Suggestions.
Deependra is offline   Reply With Quote

Old   April 26, 2016, 17:35
Default
  #10
New Member
 
Join Date: May 2015
Posts: 3
Rep Power: 11
Ebimo is on a distinguished road
Hello guys

Im trying to do species transport to my problem. I chose mixture-template for my reaction and also I defined my needed material. But in the mixture-template dialogue box when I am trying to add my material and remove N2, O2 and H2O (the default ones), I cannot remove them (I can add) and it shows me this error:

"Species n2 is used in one or more reactions." for example when Im trying to remove n2 from my mixture.

If someone have any experience plz help me. I really stuck there
Ebimo is offline   Reply With Quote

Old   June 20, 2016, 20:56
Post
  #11
Member
 
vlg
Join Date: Jul 2011
Location: My home :)
Posts: 81
Rep Power: 18
villager is on a distinguished road
Quote:
Originally Posted by Ebimo View Post
But in the mixture-template dialogue box when I am trying to add my material and remove N2, O2 and H2O (the default ones), I cannot remove them (I can add) and it shows me this error:

"Species n2 is used in one or more reactions." for example when Im trying to remove n2 from my mixture.
I confirm it. it seems at least strange!
I have found a workaround:
File -> import -> Chemkin Mechanism -> specify the Gas-phase CHEMKIN Reaction Mechanism File with your mixture. The file example is listed below:
Code:
ELEMENTS
O H C
END
SPECIES
CO CO2 O2 CH3OH H2O
END
The code defines 5-component gas mixture.
villager is offline   Reply With Quote

Old   August 31, 2016, 15:54
Default
  #12
Senior Member
 
Joe
Join Date: Feb 2012
Location: Canada
Posts: 112
Rep Power: 14
enayath is on a distinguished road
Hello Guys,

I am new to Fluent, so I'm sorry if my questions are basic:
I have a reaction (both forward and backward) with liquid and solid species in a domain that has a porous medium and Ammonia as a solvent. Basically, a 3-D geometry that I want to use LES model with species reactions. There is a porous medium that contain my nutrient basket and some solid sections in the domain that deposition of these nutrients would take place on them.(although both etching and deposition happen in nutrient basket and solid sections). I read the Fluent user guide and theory guide. Now have some questions:

1- What is the different between volumetric reaction and surface wall reaction?
2- In the "Species Model" box, in the "Reactions" section, can I just select wall surface and not volumetric? if so, what does that mean?
3- I think for my case, I should deal with species and solid species. no site species. right?
4- what does "abundant species by mass" mean in the user guide?...It says the last species should be the abundant species... How can I find/define my abundant species?

Thank you in advance!
enayath is offline   Reply With Quote

Old   February 8, 2017, 11:03
Default
  #13
New Member
 
Felipe Costa
Join Date: Jan 2017
Posts: 4
Rep Power: 9
grapheeene is on a distinguished road
Hi. I am trying to model a simple and single surface reaction

CH4 (g) -> C(s) + 2H2 (g)

on a "CVD reactor" to grow C(s) (graphene, or graphite). The inlets are gases inyections of Ar and CH4 at the top of the CVD, and two thin parallel Cu surfaces are in the center (the catalyst), where C would be deposited. I am trying to model the growth of C(s) and the amount of H2 on the entire region, with and without thermal difussion (Soret).

The question is that I have some trouble with setting up the mixture and the reaction. I have a confusion with 'site' and 'solid' species in this case, and with reactions and mechanisms (since it's only superficial). I need the system to stop generating H2 when the surface is full of C(s). And I think the reaction that needs to be computed as a 'wall surface' reaction is

CH4 (g) + C(site) -> C(site) + C(solid or bulk) + 2H2 (g)

I know the pyrolisis of CH4 is more complex but I want to start with something. Am I right?

My model with and without thermodifussion result in almost the same profile for H2 mole fraction for steady state, there is a tiny difference (less than 30%).. I think this is wrong, because I should find more H2 moles near heated surfaces with thermodif, it should be near one order of magnitude. Why does this happen?

I hope anyone can help me with anything.
dipoo likes this.
grapheeene is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 11:57.