CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

boundedness of the result (negative solute concenteration)

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 1, 2010, 18:54
Default boundedness of the result (negative solute concenteration)
  #1
New Member
 
Muhammad reza hassani
Join Date: Apr 2010
Posts: 29
Rep Power: 16
mhassani is on a distinguished road
Dear foamers
I have added a scalar transport equation to nonNewtonianIcoFoam solver to solve mass transport equation. I expected the results of concentration to be positive all over the domain, however I found out that there are some negative value for the concentration at some points.
anyone have any idea of how I can change the setting of the solver to get more realistic results? may be changes in the interpolation scheme or time marching (implicit or explicit) can be effective?
by the way, the case is the release of a solute at the inlet of an 2D expansion of a tube.
muhammad
mhassani is offline   Reply With Quote

Old   September 2, 2010, 19:18
Default
  #2
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 452
Rep Power: 24
santiagomarquezd will become famous soon enough
Hi Muhammad. What div and temporal schemes are you using? Try gauss upwind for div and euler for time, both are first order accurate but ever bounded. If you have problems again you have to look at in another place.

Good luck.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   November 30, 2010, 05:30
Default hi
  #3
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 16
siavash_abadeh is on a distinguished road
Quote:
Originally Posted by mhassani View Post
Dear foamers
I have added a scalar transport equation to nonNewtonianIcoFoam solver to solve mass transport equation. I expected the results of concentration to be positive all over the domain, however I found out that there are some negative value for the concentration at some points.
anyone have any idea of how I can change the setting of the solver to get more realistic results? may be changes in the interpolation scheme or time marching (implicit or explicit) can be effective?
by the way, the case is the release of a solute at the inlet of an 2D expansion of a tube.
muhammad
hi mhassani
i have a problem and i think you can help me.
my problem is mass transfer from solid to liquid with electrochemical method and i want simulate this problem with openfoam and i cant find solver for it.
how resolvent your problem?
use a solver or write a code for it.
if you write a code : can i have your code?
very very thank in advance
siavash_abadeh is offline   Reply With Quote

Old   December 2, 2010, 14:29
Default
  #4
New Member
 
Muhammad reza hassani
Join Date: Apr 2010
Posts: 29
Rep Power: 16
mhassani is on a distinguished road
Dear Siavash
sorry for late response of mine. I don't know about electrochemical method. the problem of mine involved with a simple convection-diffusion equation. the method of adding the equation was described through a powerpoint tutorials available on internet. it is no difficult. you can find it so easy.
you just create a file with "*.h" extension and include it in the solver code. you have to compile the new solver after modification of the make files.
the boundedness of the concentration still remains unsolved.
regards
hassani
mhassani is offline   Reply With Quote

Old   December 4, 2010, 04:12
Default
  #5
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 16
siavash_abadeh is on a distinguished road
Quote:
Originally Posted by mhassani View Post
Dear Siavash
sorry for late response of mine. I don't know about electrochemical method. the problem of mine involved with a simple convection-diffusion equation. the method of adding the equation was described through a powerpoint tutorials available on internet. it is no difficult. you can find it so easy.
you just create a file with "*.h" extension and include it in the solver code. you have to compile the new solver after modification of the make files.
the boundedness of the concentration still remains unsolved.
regards
hassani
dear hassani
thanks for your attention
i have prolem same you.
i have mass transport eq in a domain and i want to know:
do you have code for mass transport eq with openfoam?
can you give me your code? maybe i can change it and use it in my case.
very thanks in advance.
siavash_abadeh is offline   Reply With Quote

Reply


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
negative element volume (CFX-10.0) CFDworker CFX 8 September 27, 2011 19:16
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11
[ICEM] Hex refinement negative quality afterwards Anorky ANSYS Meshing & Geometry 8 March 22, 2010 04:25
Contour plot with negative velocity J CFX 11 November 3, 2008 17:46
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00


All times are GMT -4. The time now is 12:59.