@ Hannah and Perry:
I need some time to write equations to show you what seems wrong to me when we use energy equation, formulated with mass fluxes |
1 Attachment(s)
First of all, i would like to post derivation of mixture equation for temperature.
As you can see, it is possible to derive equations without using mass fluxes and heat capacities in fluxes expressions. Why i don't like mass fluxes? Because interFoam's pressure equation solves for volumetric fluxes. Until your densities are constant (doesn't depends on temperature), you can easilly construct mass flux which satisfies mass continuity. But for other cases you will get in "trouble" |
fireFoam
Hi,
I want to model small pool fire 2d case with methanol. for steady state flame propagation I want to implement interface boundary condition with on FireFoam. 1. conduction by flame = latent heat* velocity 2.stefan boundary condition with velovity == diffusion*mass fraction grad/1-mass fraction. 3.Clausius-Clapeyron equation. Can any body help how should I proceed ????????????????????? |
Hello Atul,
this thread's topic has nearly nothing in common with your post. Please only make posts where it fits to the topic and/or start a new thread. Also have a look whether you find tutorials on the wiki or somewhere else. |
3 Attachment(s)
@ Hannah
I had another look on my own and also your case setup of the 2D damTempBreak case. Apart from that we use our own solvers, which should be also the same except the OpenFOAM version, I adjusted the last details to your setup. And still when I view the results in paraview e.g. after 0.5 sec AND RESCALE the caption I can see temperatures below minimum temperature and above maximum temperature of the starting condition (which is physically impossible). Hannah, could you please review your damTempBreak case and tell if you make the same observation? I am attaching pictures of temperatur field for 0 and 0.5 sceonds and my case file. Attachment 44310 Attachment 44311 Attachment 44312 |
2 Attachment(s)
I just checked, and I do see that effect, although not nearly as extreme as that. With initial fluid temperatures of 270 and 300 K, I end up with fluids with maximum and minimum temperatures of ~267 and 301 K (see attached plot). I didn't notice this before, since ParaView didn't change the scale by default.
However, for my actual viscous fluid extrusion simulations the effect seems to be negligible (min 293.15 K ambient air, max 315.15 K input fluid). There's a brief excursion about halfway through the simulation, but the model recovers and the minimum and maximum temperatures are otherwise very steady and accurate (see attached plot). Perhaps the wandering temperatures are due to small numerical errors that add up from poor mesh or temporal resolution in the damBreakTemp case? For my purposes it's definitely given very good results. |
Well I cannot explain why the whole aproach seems to have worked in previous versions of OF an still works ok in 2.3 while it doesn't in 2.4.0.
Quote:
|
Quote:
Hi Perry, I know this is a really old post however i wonder if you managed to solve this problem as i have tried so many ways but i keep getting your error. Kind regards Rimsha |
1 Attachment(s)
Hello Rimsha,
Disclaimer: I did all this in OF 2.4.0 and got unphysical results. Seems like Hannah did not in version 2.3. Nevertheless, I did succesfully compile the solver and libraries. When I remember correctly the missing change was done in the TEqn.H file: "surfaceScalarField kappaf = incompressibleTwoPhaseMixture.kappaf();" was changed to "surfaceScalarField kappaf = mixture.kappaf();" I can't find my documentation with english commentary, but in spite of that the attached file should help you (there is also described how to setup a test case). As closing remark: I went over to using compressibleInterFoam, which solves the energy equation. |
Quote:
Thank you for your reply, I managed to complied the code yesterday, i am also getting strange results when the difference in cp is large in both fluids so i am looking into this right now. Thank you very much for the document it is highly appreciated. Kind Regards Rimsha |
Hi,Hannah and Perry
I'm working on the similar object. I've used the wolfdynamics page to make the same modifications as you to incorporate temperature into OpenFOAM (3.0).The damBreak case works well and the temperature field is physically.When I try to simulate the object which is a simulation of slug flow temperature field about 2.6-meters long pipeline with radius of 0.013 meter.Temperature of air and water are both 313 K.And the temperature of tube wall is 277 K.But some of final value of temperature is much higher than anticipated results(which is physically impossible).Some values are physical. The damBreak works well,but the pipeline can't. I really don't know what the problem is. Please help!:( Julia |
Hello Julia,
initially I don't know what the problem is, too. I would recommend you to start a new thread and supply more detailed information in there, like: turbulence model?, boundary, starting conditions, zip-file of your case folder, maybe a picture of physically impossible temperature field. Greetings Perry |
Thank you for your reply,I will strat a new thread:)
|
working solver
1 Attachment(s)
hey,
i just worked through the solver and made some little lib and class changes and it works fine with OF2.3.1. |
Quote:
Kind Regards Rimsha |
Can you let us all know what the problem was and how you fixed it?
|
1 Attachment(s)
Hi Smog,
sure thing, there were a number of things that was wrong with my libraries, 1) Attachment 48420 All the Foam:: incompresibleTwoPhaseMixture: : incompressibleTwoPhaseMixture needed to be changed to the picture attached above so to the name of the new library name i have selected. 2) error of Kaapaf in the Temp equation was also a concern which was solved using the suggestion by Perry above. Apart from that i cant remember what else however if i do remember i will post back. Kind Regards Rimsha |
Help please!!!
Hello Guys,
Can someone help me please, I need to model melting of a powder particle (solid to liquid, - MODELLING OF BLOWN-POWDER REPAIR PROCESS) and I have managed to modify, compile and use the new SOlver (Inter_Thermal_Foam) as done Perry, hannah and the attached PDFs (Thanks Perry), and it seems to work fine with the "damBreak" tutorial, but now need to adpot it to mine, but have no idea where to start. I was told that I would need to add the "SOURCE" term to the Energy euation TEqn.H, and a few other modifications to files. I also came across this link, but have no idea what he (Fabian) did; https://www.cfd-online.com/Forums/op...se-change.html Any help would be greatly appreciated. You can also contact me via email if that helps. Many thanks in advance. Ehi |
Help please!!!
Hello Guys,
Can someone help me please, I'm trying to model melting of a powder particle in space or unit volume (solid to liquid), and I have so far managed to modify, compile and use the new Solver (Inter_Thermal_Foam) as carried out and beautifully explained/detailed by Perry, Hannah and the attached PDFs (Thanks Perry!), and it seems to work fine with the "damBreak" tutorial, but now need to adopt it to mine, but have no idea where to start. I was told that I would now need to add the "SOURCE" term t(from "STEFAN"'s condition) to the Energy equation TEqn.H, and a few other modifications to the appropriate files. HOW DO I DO THIS PLEASE? I also came across this link, but have no idea what/how he (Fabian) did (it); https://www.cfd-online.com/Forums/op...se-change.html Any help would be greatly appreciated. |
adding "SOURCE TERM" to the modified interFoam solver and Shared Library
Hello guys,
I'm currently working on additive manufacturing (Blown Powder repair) for Aerospace Turbine components PhD, and successfully modified and ran a modified version of following the instructions. But now I need to incorporate the "SOURCE TERM" as you've described on here; https://www.cfd-online.com/Forums/op...se-change.html, my question is, HOW DO I GO ABOUT MODIFYING THE FILES PLEASE, and do I need to modify files in the src/transportModels/incompressible, below doing the same for the "SOLVER" folder? I JUST NEED TO MODEL A SOLID MELTING, and will hopefully take it from there. solve ( hEqn == - hs*fvc::ddt(rho, gamma) //Source term from Stefan condition ); Thanks in advance, and keep up the great work you guys have been doing. Kind regards, Ekoh E E |
All times are GMT -4. The time now is 09:57. |