|
[Sponsors] |
May 28, 2009, 12:18 |
Solving Multispecies Reactive Flow Model
|
#1 |
New Member
Sanjib Das Sharma
Join Date: May 2009
Posts: 22
Rep Power: 17 |
Hi,
I am new to OpenFOAM. However, I am undertaking a difficult project which involves the following: 1) Steady state, 3D Model, 2) Single Phase Turbulent Flow 3) Multiple species, Multiple Reaction 4) Density change due to liquid-to-vapor flashing 5) Local Grid refinement The mesh is already created using GAMBIT. So need to transfer that in OpenFOAM format. Any suggestions will help me immensely. If there are similar cases solved already, please send me a link or a copy of that. My email address is: sanjibdsharma@gmail.com Best regards, Sanjib |
|
June 2, 2009, 14:37 |
|
#2 |
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 17 |
Sanjib
I would suggest to read more about two solvers: reactingFoam and dieselFoam. I guess they are not steady state solvers. Refer to Chalmers OpenFoam tutorials for some details. Local grid refiniment is possible only with static mesh. Guess adaptive local mesh refinement is not implemented in them. Kindly keep us updated with this thread about the progress, as it seems like an interesting topic. Rachel |
|
June 2, 2009, 14:43 |
|
#3 |
New Member
Sanjib Das Sharma
Join Date: May 2009
Posts: 22
Rep Power: 17 |
Thanks Rachel. Where do I find the Chalmers OpenFoam Tutorials ? Please send me the link.
Regards, Sanjib |
|
June 2, 2009, 14:46 |
|
#4 |
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 17 |
Hi Sanjib,
A little bit of googling would have revealed the result Anyways, here are the link. http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/ http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2008/ |
|
August 12, 2009, 08:26 |
|
#5 |
New Member
Sanjib Das Sharma
Join Date: May 2009
Posts: 22
Rep Power: 17 |
Hi Rachel,
I tried to use the reactingfoam solver for a very simple case. Since I wanted to include only species transport without any reaction, so what made the changes in the code accordingly (turning off the chemistry header files and chemistryProperties input file). However, since the thermodynamic data are required, so I kept the thermophysicalProperties input file. But this file access data from chemkin database with chemkin.inp and therm.dat file. Although I could compile the code, but while running it the therm.dat file could not be accessed. 1) Do I have to compile the chemkin first ? 2) How do I use thermodynamic data without the chemkin format in the thermophysicalProperties input file ? Thanks for your help, Sanjib |
|
August 14, 2009, 04:51 |
|
#6 |
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 17 |
Can you post the output of your log file?
I do not think that chemkin has to be compiled. IT is just the Chemkin format. OF-1.6 does not use Chemkin format for some solvers. You might want to check new solvers in OF, e.g.- rhoReactingFoam I have some problems with this new solver http://www.cfd-online.com/Forums/ope...ctingfoam.html Let me know if you succeed with this solver. |
|
August 14, 2009, 09:57 |
Some issues with OpenFoam
|
#7 |
New Member
Sanjib Das Sharma
Join Date: May 2009
Posts: 22
Rep Power: 17 |
Hi Rachel,
Thanks for your reply. I understood what changes are necessary. Finally I could compile and run the reactingFoam module for a non-reactive species transport case. However, there are more serious issues which I am not able to understand. I am copying an earlier mail below: ************************************************** ***** I am new to OpenFoam, however have been using FLUENT for a long time. To start a new project in OpenFoam, I thought it might be good to compare a benchmark case between OpenFoam and FLUENT. So, I took turbulent flow over a BACKWARD FACING STEP. Following are my conclusions: 1) Substantial difference in the velocity, turbulent-kinetic-energy and turbulent-dissipation fields. 2) The reattachment length is shorter in OPENFOAM prediction; While there are two circulation zones, top and bottom of the jet, predicted by OPENFOAM, there is only one predicted by FLUENT, the bottom one. 3) While I can use higher order discretization schemes in FLUENT, in OpenFoam I can use only "upwind" which is a first-order, bounded scheme. When I used 2nd order schemes, the results showed instability and no convergence was achieved. Does this mean we have to use only "upwind" and cannot use the higher order schemes ? 4) I also wanted to test the linear solvers by replacing default PCG/PBiCG with GAMG. However, this did not work and the simulation stopped abruptly. Does this mean we can use only PCG/PBiCG ? These are critical to my project and would help immensely to understand how OpenFoam works. ************************************************** ****** Regards, Sanjib |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 07:20 |
BuoyantSimpleFoam | amitshah | OpenFOAM Running, Solving & CFD | 4 | June 3, 2009 07:04 |
What is the difference between liquid reactive flow and gas reactive flow? | James | Main CFD Forum | 6 | May 15, 2009 13:14 |
Problems with sonicTurbFoam | tangd | OpenFOAM Running, Solving & CFD | 0 | May 29, 2006 10:07 |
fluid flow fundas | ram | Main CFD Forum | 5 | June 17, 2000 22:31 |