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

Bubble is not increasing in size in Nucleate boiling

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 8, 2019, 21:56
Default
  #21
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Sad news, your suggestions didn't work. The symmetryPlane BC on the base affects the axisymmetric setup and prevents it running. Do you have any other suggestions ?
Hi Sam,

sorry I am a little confused about your case. 1/4 of the bubble will be simulated. So why using axis-symmetry? when using the axis boundary, I think the wedge boundary condition might be more suitable. And I am still not sure your case is in 2d or 3d.

If you want to setup the single bubble growth case just like the default case in the solver, using the makeaxis tool, I guess you have to reconsider your geometry. Maybe You have to draw the bottom half of the structure as well if the symmetry boundary is conflict on the base.

sorry for my poor English.

Best
Bali is offline   Reply With Quote

Old   December 8, 2019, 22:32
Default
  #22
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Your English is great, don't worry about it. It is a 2D axisymmetric simulation, i'm trying to replicate this scenario based on this paper (https://www.mdpi.com/1996-1073/10/3/272/htm).


initialising a thermal boundary layer at interface using funkySetFields

However, the BCs used for the solver used in this simulation does not work for mine.

Thanks for the help again.
saj216 is offline   Reply With Quote

Old   December 8, 2019, 23:35
Default
  #23
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Your English is great, don't worry about it. It is a 2D axisymmetric simulation, i'm trying to replicate this scenario based on this paper (https://www.mdpi.com/1996-1073/10/3/272/htm).


initialising a thermal boundary layer at interface using funkySetFields

However, the BCs used for the solver used in this simulation does not work for mine.

Thanks for the help again.
oh, I've seen this paper before, the solver they used is the evapVOFHardt I think. But you can use the validation case for the interThermalPhaseFoam solver. In Kunkelmann's PhD dissertation this case could be found.

I am really busy these days, perhaps you should check Kunkelmann's dissertation, or just try the simulation in a 2d plane without wedge boundary first, determining the boundary conditions. The files you uploaded were not complete for running in OpenFOAM so I cannot give targeted advice.

funkysetfields is good for the temperature fields initialization but, I recommend you simplify the case by using the default setFields tool first. When everything is settled, use funkysetfields tool to obatin a better initial temperature field.

Besides, in this case spurious currents should be avoided, SST model in interThemalPhaseFoam or surface tension model in evapVOFHardt solver is better than the default CSF model in interFoam.

good luck
Bali is offline   Reply With Quote

Old   December 8, 2019, 23:44
Default
  #24
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
No problem, thanks for the advice . I did try to upload the complete file but it was too big.
saj216 is offline   Reply With Quote

Old   December 8, 2019, 23:53
Default
  #25
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
No problem, thanks for the advice . I did try to upload the complete file but it was too big.
Hi Sam, I find a case regarding single bubble growth in superheated liquid,which can be run with interThermalPhaseChangeFoam. 1/8 of the bubble is simulated without the wedge boundary. check it maybe it will inspire you.

good luck
Attached Files
File Type: zip singleBubble.zip (8.9 KB, 23 views)
Bali is offline   Reply With Quote

Old   December 8, 2019, 23:55
Default
  #26
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hello again,

Do you have the link for this?
saj216 is offline   Reply With Quote

Old   December 8, 2019, 23:57
Default
  #27
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hello again,

Do you have the link for this?
ha? I have uploaded it in the last reply. maybe pm me your email address so I can send it to you
Bali is offline   Reply With Quote

Old   December 8, 2019, 23:59
Default
  #28
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Sorry I didn't notice. I'm a little tired . Thanks for all your help. If I find out why it didn't work i'll let you know.
saj216 is offline   Reply With Quote

Old   January 3, 2020, 14:37
Default
  #29
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
I believe I have resolved this issue. It was the way in which i was treating the temperature in terms of my initialisation.

I do have another question regarding the setting of the reference temperature in the transport properties. I know it is part of the thermal energy transport equation in the governing equations. However, how exactly should this be treated in reference to a bubble?
saj216 is offline   Reply With Quote

Old   June 2, 2020, 15:14
Default
  #30
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hey guys,

Does anyone know how to code for the bubble radius in postprocessing?
saj216 is offline   Reply With Quote

Old   June 3, 2020, 04:49
Default
  #31
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hey guys,

Does anyone know how to code for the bubble radius in postprocessing?
Hi Sam, if there is only one bubble in your case, you can use paraview to get the bubble diameter with respect to time easily. of course, write some codes in the controlDict file can plot the diameter when running the simulation.


functions
(
DataSummary
{
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
redirectType DataSummary;
outputControl timeStep;
outputInterval 10;
code
#{

//************************************************** *********************
//First get t and dt
scalar t = mesh().time().value();
// scalar dt = mesh().time().deltaTValue();
//************************************************** *********************

//************************************************** *********************
//Now get void fraction
const volScalarField& alpha1 = mesh().lookupObject<volScalarField>("alpha.phase1" );

const scalar V_Vapor = gSum( (1.0-alpha1.internalField()) * mesh().V() );
//************************************************** *********************
// Calculate the equilibrium diameter of the bubble
const scalar D_bubble = 2*cbrt(4*V_Vapor*3/4/3.141592654);
//************************************************** *********************
//Finally print out results:
//Now write out data:
if( Pstream::master() == true )
{
std:fstream fs;
fs.open ("Diameter.dat", std::fstream::app);
fs.precision(6);
fs << t << "\t" << D_bubble << "\t" << "\n" ;
fs.close();
}

#};

codeInclude
#{
#include <fstream>
#};
}


);
Bali is offline   Reply With Quote

Old   June 3, 2020, 05:12
Default
  #32
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hi Bali,

Thanks for the reply. Wow I should have realised that. My bad. The only problem is that I had to upgrade the solver to openfoam version 7 and I have not been able to run Data Summary in it . Sorry to ask but do you have any idea?
saj216 is offline   Reply With Quote

Old   June 3, 2020, 05:17
Default
  #33
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hi Bali,

Thanks for the reply. Wow I should have realised that. My bad. The only problem is that I had to upgrade the solver to openfoam version 7 and I have not been able to run Data Summary in it . Sorry to ask but do you have any idea?
I havent used of-7 but maybe add this line at the end of the controldict file is helpful

libs ("libutilityFunctionObjects.so");
Bali is offline   Reply With Quote

Old   June 3, 2020, 05:21
Default
  #34
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
It is already there, thanks for the help though. I had to do it so it can be used on a cluster. I wish I didn't have to. Such a pain to do.

Thanks again,

Sam
saj216 is offline   Reply With Quote

Old   June 3, 2020, 06:52
Default
  #35
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hi again Bali, I reread your comments and you stated that you can get bubble radius as a function of time through paraview. Could you expand on how to achieve this. Step by step if possible . Thank
saj216 is offline   Reply With Quote

Old   June 3, 2020, 07:47
Default
  #36
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 7
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hi again Bali, I reread your comments and you stated that you can get bubble radius as a function of time through paraview. Could you expand on how to achieve this. Step by step if possible . Thank
Hi Sam, it is really simple, just do some volume integral with the volume fraction field. Assume that alpha = 0 is gas when alpha = 1 is liquid.

1. calculator. input 1-alpha.
2. open the spreadsheet view and ctrl+a all the cell data. use integrate variables to obtain the volume of the vapor gas.
3. select the integrated data in spreadsheet and use plot selection over time
Bali is offline   Reply With Quote

Old   November 10, 2020, 02:00
Default
  #37
New Member
 
M Naarendharan
Join Date: Aug 2020
Posts: 21
Rep Power: 5
Naaren is on a distinguished road
Quote:
Originally Posted by Fahmida View Post
I had to reduce the domain size and increase the evaporation threshold number although I am not clear enough how the evaporation number was effecting this issue.

Hi,
I am facing a similar issue of bubble size not growing when I go for a 3D geometry in this solver. I would be happy if you can share the evaporation threshold number which was used in your simulation.


Thank you,
M Naarendharan
Naaren is offline   Reply With Quote

Old   November 11, 2020, 13:25
Default Movement of BUbble
  #38
New Member
 
M Naarendharan
Join Date: Aug 2020
Posts: 21
Rep Power: 5
Naaren is on a distinguished road
Hi everyone,
I have one more doubt. In a 2D Axisymmetric simulation, my bubble growth is proper, but when I go for a 3D simulation in a smaller domain, the bubble placed at the center moves towards a side wall and eventually hits the wall.

I am not sure about what is causing the movement of the bubble. I would be grateful if someone could help me with this issue.


Thanking you,
M Naarendharan
Naaren is offline   Reply With Quote

Old   November 18, 2020, 03:22
Default
  #39
New Member
 
M Naarendharan
Join Date: Aug 2020
Posts: 21
Rep Power: 5
Naaren is on a distinguished road
I am still facing the issue. Any suggestions would be of great help.


Thank you,
M Naarendharan
Naaren is offline   Reply With Quote

Old   November 25, 2020, 22:29
Default
  #40
New Member
 
M Naarendharan
Join Date: Aug 2020
Posts: 21
Rep Power: 5
Naaren is on a distinguished road
Hi,
I am still facing the issue. Any suggestions would be very much appreciated.


Thank you,
M Naarendharan
Naaren is offline   Reply With Quote

Reply

Tags
bubble, nucleate boiling

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Why pressure is increasing when Degassing BC is used for a pool boiling case? soumitra2102 CFX 5 August 4, 2020 06:21
TwoPhaseEulerFoam - Bubble diameter size Problem BlnPhoenix OpenFOAM Running, Solving & CFD 8 September 17, 2019 15:59
How to make UDFs for Nucleate Pool Boiling Simulation? SIKJAE Fluent UDF and Scheme Programming 1 August 4, 2018 07:07
OpenFOAM UpdateCoeffs Nucleate Wall Boiling suneth.warna OpenFOAM Programming & Development 0 July 31, 2016 17:42
Flow boiling & bubble formation Venkat CFX 1 July 23, 2009 09:20


All times are GMT -4. The time now is 05:41.