CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

"reaction" threshold

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

LinkBack Thread Tools Display Modes
Old   October 3, 2009, 10:59
Default "reaction" threshold
New Member
Join Date: Oct 2009
Posts: 1
Rep Power: 0
yansongdi is on a distinguished road
Please, can you suggest which solvers I can combine to model the transport in pipes of dissolutions close to the solubility limit?

My plan is to have passive transport of a component X as long as it is below the solubility Xmax. Where X>Xmax, it should precipitate as a solid Xs, settle to the bottom of the pipe, and remain there until some fluid parcel with X<Xs dissolves the precipitate again.

I would define two scalar fields X and Xs which describe the dissolved concentration and the precipitated amount (VOF-method).

The question is, how to describe the precipitation process. I would define
where X>Xmax or Xsol>0 and dX/dt=dXs/dt=0 elsewhere. The rate k would be adjusted to about 1e3/runTime.deltaT() to force near-equilibirum conditions.

Can I use the chemistrySolver and the ODEchemistryModel to model "reactions" which are only activated above a solubility limit Xmax? If yes, please provide a few lines of code to show the essential steps or point me to some useful documentation.
Another question is, whether the combination of transport and reaction as it is implemented in reactingFoam is still valid, if the reaction time scale is much faster than the transport time scales?
As I understand reactingFoam, it solves first the chemical ODE and then the transport step. The disadvantage is that fast kinetics are frozen during the transport step. Does some solver exist which solves the chemical reaction implicitly during the transport step?

Additional processes, like temperature dependent solubility, heat release or volume change during dissolution or precipitation, will also be important at a later time. But for the moment I only need to find a way how to model the precipitation.

yansongdi is offline   Reply With Quote


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
Topo change doum OpenFOAM Running, Solving & CFD 0 June 22, 2007 16:26
How much is the refine threshold for grid adaption Jie FLUENT 0 February 21, 2002 01:43

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