Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

negative concentrations at flux boundary results in singularity error

Please login with a confirmed email address before reporting spam

greetings

I am working on a heterogenous catalysis model with fast adsorption reactions at one of the boundaries. Right near the boundary at time=0, the concentration of the species being adsorbed is very small or even negative, which causes a singularity error.

I have set my initial conditions for concentrations in transport of diluted species, so i am not sure why this is occuring at the boundary. i can avoid the singularity by setting the following initial condition in the rate expression:

( (t>0)* cCO2 + (t<=0)* 0.03 [mol/L])

(cCO2 is the concentration of the species being adsorbed, and 0.03 is the initial concentration. The rate of adsorption is modeled using Hertz-Knudsen kinetics.)

however, I was wondering if there is a better, more physically reasonable way to solve this issue?

thank you


0 Replies Last Post 12.06.2020, 00:13 GMT-4
COMSOL Moderator

Hello Kelsey Levine

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.