|
[Sponsors] |
![]() |
![]() |
#281 |
Senior Member
|
Hello Foamers,
# Post 19 --> Mr. Fabian Roesler had constructed the solver for melting problem in OpenFOAM 2.1.1 based on Stefan Task . I have modified this solver in simplest way using SIMPLE algorithm based on only Temperature (TEqn.H) as shown in the attachment. Here, I have neglected pressure (pEqn.H) and velocity (UEqn.H) terms to fit the conditions based on my solver as explained in next comment section. In this case, by default 'alpha' is a function of temperature, 'lamda' and 'cp' is the function of 'alpha' This solver explains the movement of interface from right to left for non-isothermal phase change where melting occurs. [Stefan Problem - moving boundary is met] |
|
![]() |
![]() |
![]() |
![]() |
#282 | ||
Senior Member
|
Hello Foamers,
Following my above comment, I would like to explain my problem conditions to fit into the phase change solver. After solving the simple energy equation, I want to include the phase transition phenomenon - Boiling, Condensation and Convection [ATTACHMENT 1]. ATTACHMENT 2 explains the schematic 1D model with phenomenon. Here the main equation is only ENERGY equ. The problem is basically the coal pyrolysis. Initially, the coal is considered as wet and 2 phases of moisture and steam need to be implemented. After evaporation of moisture, once if wet coal is converted into dry ~ pyrolysis will take place, which is a different topic of discussion. Position of interface (boiling plane) is explicitly defined based on mass and heat balance ~ as not defined in previous Fabian Roesler problem. Quote:
(1) Boiling: When the heating wall (left side) reaches the temperature of 100deg, the moisture evaporates [the condition (alpha) w = 0] instantly. [ATTACHMENT 3] (2) Condensation: The steam condenses (right side of chamber) until the temperature reaches 100deg on the right side. [ATTACHMENT 4] (3) Convection: When T = 100deg, the steam flow after moisture evaporation. [ATTACHMENT 5] How to fix such condition at interface quoted below ==> Quote:
Although the equations to solve looks simple, I couldn't able to figure out the exact path to find my solution (considering my queries mainly marked in "RED mark") Kindly someone share their ideas please. Thank you !!! |
|||
![]() |
![]() |
![]() |
![]() |
#283 | |||
Senior Member
|
Hello Foamers,
To make it simple, I will explain the main core of the problem about phase change. alpha as a function of T in melting problem ~ as mentioned in above explanation Quote:
Quote:
Heat is transferred by wall conduction only. When the temperature reaches T=100deg (heating wall), the surface reaches above the saturation temperature, and so liquid (moisture) evaporates leading to the motion of vapor-liquid flat interface. If the vapor-liquid interface moves to the complete right, then the chamber is saturated with only evaporated vapor. Interface stays flat for 1D PROBLEM. //************************************************** *****************// Following the above condition, mass and heat balance are calculated as, Quote:
Kindly someone help me by sharing ideas.. Thank you Last edited by Kummi; July 17, 2019 at 00:11. |
||||
![]() |
![]() |
![]() |
![]() |
#284 |
Member
Hasan Celik
Join Date: Sep 2016
Posts: 64
Rep Power: 9 ![]() |
Is there any CHT melting solver that you know?
Thanks! |
|
![]() |
![]() |
![]() |
![]() |
#285 |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Hello Everyone,
First of all, I thank everyone who has contributed to this wonderful thread. I downloaded the solver posted by Typhian in post #263. I am able to compile it successfully in OF V6 and perform the test case too. But, when I changed the mesh resolution to 132X100 in the blockMesh dictionary. The case gives an abnormal result. After changing the entries int he blockMeshdict, I ran blockMesh and meltFoam. I believe it has something to do with the cellToRegion file. Someone please explain me the use of cellToRegion file. I mean what it does exactly. Sorry if my question is too silly. Thank You. |
|
![]() |
![]() |
![]() |
![]() |
#286 | ||
Senior Member
|
Hello Pavithra,
I downloaded the solver and compiled in version 5.0 (post #263). U can find the same case in the attachment. I don't find any trouble while compiling and solving this case. Quote:
You have compiled in OF V6 - thats where the problem is, I hope so. Quote:
cellToRegion file is like any other field file but not associated to a specific field. The BC type is either zeroGradient or calculated. calculated (if defined) --> in case of coupling patches to other regions zeroGradient (if defined) -->boundary patches (outside of domain) Hope it helps. Thank you ^^ |
|||
![]() |
![]() |
![]() |
![]() |
#287 |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Respected Sir,
Thank you so much for your kind reply. I have attached the liquid fraction and temperature distribution at 2 seconds. I find the liquid fraction and temperature distribution at 2 secs as abnormal. Kindly please give your comments on this. So, in normal case, cellToRegion file does not have any use. The default case should produce the results even without the cellToRegion file. Sir, please correct me, if I am wrong. Thank You. |
|
![]() |
![]() |
![]() |
![]() |
#288 |
Senior Member
|
In 1D case of mesh size (100 1 1) - the results looks reasonable.
However, for 2D case, the behavior is weird. Have you tried increasing the width size or decreasing the default time size less than 0.005 ? |
|
![]() |
![]() |
![]() |
![]() |
#289 | |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Quote:
But, Sir please help me understand that how increasing the width can resolve this issue. Thank You |
||
![]() |
![]() |
![]() |
![]() |
#290 | ||
Senior Member
|
Hello Pavithra,
Quote:
Quote:
I am not proficient in multiphase problems. My field of work is different. However, I am sharing ideas to the best of my knowledge. In the attachment given below, both images are captured at 170s. When the temperature range is between 300-310K, the contour plot looks apparent. But when the range is rescaled btw 300-340K, there found strange behavior - bottom left corner temperature rises. Time reduction doesn't helps. I am not sure where could be the problem. ~Have you gone through this post completely from the beginning. I am sure it will give you the spark. ~Have you tried the extended version coding by Fabian Roesler [post #81] ~Posts #130 #151 #154 draw contour results - might be helpful for you to validate ~Fabian Roesler - Ph.D thesis on post #180 ~Ole Richter posts from #201-205 - about incompressible two phase mixture might help too.. Keep updating your progress. You will find a way one day!! Thank you |
|||
![]() |
![]() |
![]() |
![]() |
#291 | |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Quote:
Sir, Thank you so much for your detailed reply. I will go through all the posts carefully and will update my progress here. Thank you so much Sir. |
||
![]() |
![]() |
![]() |
![]() |
#292 |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Hello Everyone,
I am trying to use the solver by Fabian posted in this thread. Could someone please clarify me the following points? 1) What is Tdim in constant/transportProperties? 2) From his paper, I can see that Tdim = (T_liquidus) - (T_solidus). If this is true, for isothermal melting, Tdim = 0. If I set Tdim = 0, the solver crashes due to division by 0. 3) Moreover the value of Tdim affects the results, drastically. Smaller values of Tdim leads to accelerated melting and viceversa. Please help me in choosing a right value of Tdim for isothermal melting. Fabian's solver and link to his paper are attached. https://link.springer.com/article/10...231-011-0866-9 Thank You. |
|
![]() |
![]() |
![]() |
![]() |
#293 |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Hello Everyone,
Finally, I have managed to modify buoyantBoussinesqPimpleFoam to simulate melting problems. I have attached the solver and Gallium melting testcase. The solver is tested in OF v6. Any feedback or suggestions are welcome. Also, please give me a direction or any guidance to parallelize this solver. The solver works well when running on a single core. But, blows up while running in parallel. Thank You. - Pavithra. |
|
![]() |
![]() |
![]() |
![]() |
#294 |
Member
Join Date: Apr 2019
Location: India
Posts: 81
Rep Power: 7 ![]() |
Hi Everyone,
I am attaching the results of some validation cases, herewith. I have attached the solver, cases and results. Any suggestions for improvements are welcome. Thank You. - Pavithra |
|
![]() |
![]() |
![]() |
![]() |
#295 |
New Member
Luca
Join Date: Jun 2019
Location: Colorado, USA
Posts: 6
Rep Power: 7 ![]() |
Kummi,
For a finite volume solver shouldn't the energy balance always be written in W/m^3? In this case there could be some confusion due to the length attributed to dx/dt (the velocity of the interface, correct me if I am mistaken). I think you want to divide dx/dt by an appropriate length scale for your vapor-liquid interface transport phenomena, then the source term can be added to the TEqn with correct units. It seems simple but in my experience careful dimensional analysis can solve/avoid many issues when writing equations for a finite volume solver. Cheers |
|
![]() |
![]() |
![]() |
![]() |
#296 | |
Senior Member
|
Hello Luca,
Thank you for your comments. Energy balance will be written in terms of W/m3 in finite volume approach basically. My problem is a kind of 1D pyrolysis (Coal) modelling, which is focussed mainly on heat loss - thus only Energy equation is included in my work. ENERGY EQU: rho*Cp*(dT/dt) = del/delx[K*dT/dx ] + rho*dQ/dT + r*cp*dT/dx (W/m3) Unsteady conduction = Diffusion term + SOURCE TERM 1 + SOURCE TERM 2 SOURCE TERM 1 = heat loss due to pyrolysis SOURCE TERM 2 = phase change of moisture (moisture embedded in the wet coal) My query is all about the SOURCE TERM 2. Concerning it, I posted my query here - because this post discusses phase change problem - thought of gaining some assistance here for my work. Pyrolysis - Arrhenius-like degradation chemistry Drying - surface modelling at boiling plane (drying when moisture boils at 100deg) Since, drying is based on surface modelling, the mass and heat balances unit is calculated at the surface. The calculated mass balance (r) is multiplied with cp and dT/dx as SOURCE TERM 2 gives the unit of W/m3. Quote:
How the source terms can be included in such problems? I have contacted certain people, few replied as they are not familier with pyrolysis, others asked me to look into interFoam solvers. But the approach of interFoam seems to different with separate transport equation for volume fraction, which is not the same in my case. I'm trying to figure out and learn as how the problem can be attacked based on algorithm and importantly how such problems can be approached step by step in OpenFOAM// Please share your ideas, it will be highly helpful. Thank you. |
||
![]() |
![]() |
![]() |
![]() |
#297 |
New Member
Luca
Join Date: Jun 2019
Location: Colorado, USA
Posts: 6
Rep Power: 7 ![]() |
I should preface by saying I am relatively new to OpenFOAM and am currently working on a cht problem with solidification and melting, hence my monitoring of this thread. I do have more extensive experience with CFD modeling in general.
Based on your last post I don't see a problem. Perhaps I am misunderstanding, but it looks as if you have a valid energy equation in W/m^3. The 2d surface model is an additional function used to define your source term, the units of this don't necessarily matter so long as they are correct within the function definition. Do you know dx/dt based on results from the previous timestep? The way I see it, you should be able to calculate (r) with a function object, but I am not sure if you could then use that term to your TEqn. |
|
![]() |
![]() |
![]() |
![]() |
#298 | |||
Senior Member
|
Hello Luca,
Sorry for my late response. Quote:
Quote:
Quote:
Thank you |
||||
![]() |
![]() |
![]() |
![]() |
#299 | ||
Senior Member
|
Hello Luca,
In previous message, I have given a condition for moisture evaporation (phase change) quoted in a box. Quote:
I'm hereby attaching the solver (OF211). In this solver, under TEqn.H, I solved a if-else loop based on the above condition. Quote:
If you have any ideas about it, kindly do share. It will be helpful. Thank you |
|||
![]() |
![]() |
![]() |
![]() |
#300 |
New Member
d durga prasad
Join Date: Nov 2019
Posts: 5
Rep Power: 6 ![]() |
I was trying to make solidification model of pure metal.
can I use erfConvectiveMeltingPimpleFoam solver for solidification just by changing initial and boundary conditions or do I need to make any changes to the solver part? |
|
![]() |
![]() |
![]() |
Tags |
melting openfoam |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Melting and solidification with free surface problem? | cqlwj123 | CFX | 6 | July 25, 2013 02:46 |
Can I solve this problem by Fluent? | Kai_kc | FLUENT | 1 | October 27, 2010 05:29 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 06:29 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 19:13 |
Melting Problem | M | FLUENT | 0 | April 29, 2007 16:07 |