natural convection in a melting furnace

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

 April 24, 2016, 06:55 natural convection in a melting furnace #1 New Member   donny Join Date: Feb 2016 Posts: 3 Rep Power: 2 hey guys, im new in openfoam. i have to simulate two gases in a melting furnance. i know i have to use the buoyantBoussinesqSimpleFoam, but i dont know how to integrate two different gases in the hotRoom tutorial. the gases are not reacting. The temperature difference accure natural convection. The Temperature of the floor is 1023 K an the ceiling 623 K. Please i need your help. regards Donny

 April 25, 2016, 03:09 #2 Senior Member     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,548 Blog Entries: 6 Rep Power: 27 Dear donny, first hint, buossinesq approximation is in my opinion the wrong choice. You should take the buoyantSimpleFoam but to clear my mind just one question, why you think boussinesq is the right approach? __________________ Best regards, Tobias Holzmann Some interesting OpenFOAM tutorials, publications and videos on www.Holzmann-cfd.de OpenFOAM Beginners should check out the new wiki on wiki.openfoam.com A list of some active OpenFOAM contributers can be found »here« A book about the basics of »Mathematics, Numerics, Derivations and OpenFOAM« can be found on www.Holzmann-cfd.de

 April 25, 2016, 06:02 #3 New Member   donny Join Date: Feb 2016 Posts: 3 Rep Power: 2 Hey Tobi, thank you for the answer. i read in another forum that this solver is the right choice for natural convection. Therefore i tried it just for one gas and it worked, but i do not know how to do it with two gases. By the way, you have a privat message from me. best regards Donny

 April 25, 2016, 14:38 #4 Senior Member     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,548 Blog Entries: 6 Rep Power: 27 Hey guys, I was writing via messages with donny and now I will summarize. buoyantBoussinesq...Foam are using the Boussinesq-Approximation , you only should use them if the density change in your fluid is not high and the temperature range (dT) is low. Check out some references... e.g. Ferziger & Perić Donny want to implement a second gas, therefore the very simple approach would be to implement a passive scalar "S" that is limited between 0 and 1. (like alpha in VOF). This scalar represent the two gases (S = 0 -> gas1; S = 1 -> gas2; 0 < S < 1 -> mixture of gas1 and gas2) After that, you have to implement all thermodynamic data for gas1 and gas2 into the solver. Therefore checkout the createField.H and create new fields for rho1, rho2, nu1, nu2 (if incompressible assumption) or mu1, mu2 etc. The next step is to create the actual thermodynamic variables based on gas1 and gas2 and the passive scalar S: Code: `rho = rho1 * S + rho2 * (1-S);` You have to choose if you use openFOAM thermodynamics (then you normally should have two) or you implement your own polynoms (I think you also could use icoPoly.... for that). That means you have either to create a second thermodynamic object or you decouple all thermodynamic calculation with your polynoms. This model would be nice to handle and very easy. Another idea at the moment would be, to use some combustion solver and unset the reaction (kinetics). Once I did something like that with some old collegue who was exactly doing the same. The result between the combustion solver and the solver mentioned above was very close and a lot of faster. The easiest approach would be: implement only one passive scalar and use one fluid inside. Doing this leads to a first guess but hence we only use thermodynamics out of one gas, you will be not as accurate as in the method I mentioned above. If the thermodynamic data are very similar (cp, rho, nu, mu ...) then the passive scalar should work fine. This would be my suggestions. donny1991 likes this. __________________ Best regards, Tobias Holzmann Some interesting OpenFOAM tutorials, publications and videos on www.Holzmann-cfd.de OpenFOAM Beginners should check out the new wiki on wiki.openfoam.com A list of some active OpenFOAM contributers can be found »here« A book about the basics of »Mathematics, Numerics, Derivations and OpenFOAM« can be found on www.Holzmann-cfd.de Last edited by Tobi; April 26, 2016 at 03:06.

 April 26, 2016, 07:00 #5 New Member   donny Join Date: Feb 2016 Posts: 3 Rep Power: 2 Hey Tobi, thank you for this helpful answer. Im now using the reactionFoam solver without reaction and only with O2 and N2. I changed all BC's and the properties for my case. I tried to run the simulation but i got this. Reading g Creating reaction model Selecting combustion model PaSR Selecting chemistry type { chemistrySolver ode; chemistryThermo psi; } Selecting thermodynamics package { type hePsiThermo; mixture reactingMixture; transport sutherland; thermo janaf; energy sensibleEnthalpy; equationOfState perfectGas; specie specie; } Selecting chemistryReader foamChemistryReader chemistryModel: Number of species = 2 and reactions = 1 Selecting ODE solver KRR4 Reading field U Reading/calculating face flux field phi Creating turbulence model. Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon Floating point exception: 8 i dont know how to fix this problem. Under the following link you can download my case. thank you!! meltingfurnance.zip

 April 26, 2016, 08:01 #6 Senior Member     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,548 Blog Entries: 6 Rep Power: 27 Hi, maybe I will check it but as I told above, the bad thing using a reaction solver is, that it is soooo slow. Finally I am not sure if you have buoyancy included. So you have to check it yourself. My suggestion to you, build your own solver. Its so much faster, clean (you know whats going on) and really not a big deal to implement sth. like that into OpenFOAM. __________________ Best regards, Tobias Holzmann Some interesting OpenFOAM tutorials, publications and videos on www.Holzmann-cfd.de OpenFOAM Beginners should check out the new wiki on wiki.openfoam.com A list of some active OpenFOAM contributers can be found »here« A book about the basics of »Mathematics, Numerics, Derivations and OpenFOAM« can be found on www.Holzmann-cfd.de

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post flex00 FLUENT 2 January 16, 2014 12:50 Ciefdi OpenFOAM Running, Solving & CFD 0 November 7, 2013 12:44 jorien CFX 0 October 14, 2011 09:26 Alex CD-adapco 5 December 12, 2007 05:58 mauricio FLUENT 2 February 23, 2005 20:43

All times are GMT -4. The time now is 21:34.