CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

ScalartranportFoam solving with chemical reaction

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

LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2017, 03:43
Default ScalartranportFoam solving with chemical reaction
New Member
Join Date: Apr 2014
Posts: 8
Rep Power: 11
Christian86 is on a distinguished road
Dear Foamers,

i have some trouble with solving the scalar transport equation in conjunction with a chemical reaction in the concentration field.

The simulation task sketches as follows:
Two Species are transported along a computational domain. When some amount of species A touches Species B a chemical reaction (sink in the concentration field) takes place. The simulation works for moderate Diffusion and chemical reaction strength. When the reaction strength is rised, the simulation is still stable, but there are heavy unphysical oscillations inside the domain. The oscillations come along the conncection lines with the underlying blocks of the block structured mesh.

The code is actually working, so in my opinion is seems like i took the wrong options for the simulation.

A picture of the error and of the used descritization schemes is included.

Has anyone a hint, how what could be the problem?


Attached Images
File Type: jpg mistake.jpg (175.8 KB, 30 views)
File Type: jpg schemes.jpg (82.4 KB, 24 views)
Christian86 is offline   Reply With Quote

Old   March 17, 2017, 04:19
New Member
Join Date: Feb 2017
Posts: 22
Rep Power: 9
sepp.zell is on a distinguished road
Hi Christian,

in which way you modeled the reaction? Is the sink term modeled as some reaction rate times the concentration (e.g. r_c*C, with C the concentration)?
Am I correct that the reaction rate is only non-zero at the interface?
If this is the case, you might have large jumps in the coefficients at the interface (when you switch on reaction) and this jumping coefficients can cause oscillations.

I started to model flow through a porous medium.
Here the resistance is added as source term in part of the domain. porousSimpleFoam produces velocity oscillations around the interface (see e.g. ).

Therefore I guess that you might get the same problems with oscillations at your interface. I don't know if those oscillations at interface are a general problem of openFOAM. Maybe some more experienced user can give here more insight.


sepp.zell is offline   Reply With Quote


chemical reactions, mass transfer, numerical error

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 03:20
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 08:39.