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   September 2, 2019, 18:34
Default Bubble is not increasing in size in Nucleate boiling
  #1
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
Hello everyone,
I am running my simulation on Nucleate Boiling. My bubble has been created successfully like a semi-sphere on a heated surface but for some reason it is not increasing its size. What can be the reason behind this?

Any idea would be appreciated. Thank you.
Fahmida is offline   Reply With Quote

Old   September 12, 2019, 18:43
Default
  #2
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10
saddy is on a distinguished road
what solver are you using ?
saddy is offline   Reply With Quote

Old   September 17, 2019, 02:37
Smile
  #3
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 8
Bali is on a distinguished road
Quote:
Originally Posted by Fahmida View Post
Hello everyone,
I am running my simulation on Nucleate Boiling. My bubble has been created successfully like a semi-sphere on a heated surface but for some reason it is not increasing its size. What can be the reason behind this?

Any idea would be appreciated. Thank you.
hi FaRa,
could you upload your case files so that we can see where the problem is
Bali is offline   Reply With Quote

Old   September 17, 2019, 16:12
Default
  #4
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
hi Saddy,

Thank you for your response. I am using interThermalPhaseChangeFoam solver which has been extended from interFoam.

Here is the paper link:

https://www.sciencedirect.com/scienc...52711016300309
Fahmida is offline   Reply With Quote

Old   September 17, 2019, 16:48
Default Case files
  #5
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
Hi Bali, thank you for your response.

Here is my blockMesh file:

Code:
convertToMeters 1;

vertices
(
    (0 0 0)
    (0.031 0 0)
    (0.031 0.031 0)
    (0 0.031 0)
    (0 0 0.031)
    (0.031 0 0.031)
    (0.031 0.031 0.031)
    (0 0.031 0.031)
);

blocks
(
     hex (0 1 2 3 4 5 6 7) (120 120 120) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    bottom
    {
        type wall;
        faces
        (
            (1 5 4 0)
        );
    }
    outlet
    {
        type patch;
        faces
        (
            (3 7 6 2)
        );
    }
    walls
    {
        type wall;
        faces
        (
            (0 4 7 3)
            (2 6 5 1)
            (0 1 2 3)
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);
Here is the toposetDict:

Code:
actions
(
    {
        name c0;
        type cellSet;
        action new;
        source sphereToCell;
        sourceInfo
        {
            centre (0.0155 0 0.0155);
            radius 0.001;
        }
    }

    {
        name c0;
        type cellSet;
        action delete;
        source boxToCell;
        sourceInfo
        {
            
                  box (0.0145 -0.001 0.0145) (0.0165 0 0.0165);
        }
    }

);
Here is the setFieldDict:

Code:
defaultFieldValues
(
        volScalarFieldValue alpha1 0
        volScalarFieldValue T 559.95
        volVectorFieldValue U (0 0 0)
);

regions
(


   boxToCell
    {
        box (0 0 0) (0.031 0.026 0.031);
        fieldValues
        (
            volScalarFieldValue alpha1 1
            volScalarFieldValue T 559.95
            volVectorFieldValue U (0 0 0)
        );
    }
   
    cellToCell
    {
        set c0;

        fieldValues
        (
            volScalarFieldValue alpha1 0
            volScalarFieldValue T 559.95
            volVectorFieldValue U (0 0 0)
        );
    }
);

ParaView file is also attached.

Here, I would like to mention that I was able to grow the bubble size when performing in 2D with meshes of 150x150. For 3D, I have tried several meshes upto 120x120x120 with different values of Evaporation threshold number but the bubble size did not increase. I think the mesh of 150x150x150 may solve the issue. But since I am performing my simulation on my laptop with only 6 cores parallel processing, it is not possible for me to turn my pc on for a week.

So any alternate suggestion is needed.
Attached Files
File Type: zip anim.zip (51.9 KB, 29 views)
Fahmida is offline   Reply With Quote

Old   November 17, 2019, 05:26
Default Bubble detachment diameter
  #6
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
Hello everyone, I came out successfully in growing my bubble .

Now I have got the bubble detachment time by watching it in paraview. But I also need to calculate the detachment diameter in controlDict. Does anyone have any idea how to do this?
Attached Images
File Type: png detachTime2_.24s.png (9.0 KB, 54 views)
micheall likes this.
Fahmida is offline   Reply With Quote

Old   November 17, 2019, 13:07
Default
  #7
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
I have been able to calculate the vapor surface area for each timestep using swak4foam. But this is including the vapor over the liquid surface . Is there any way to remove them?
Fahmida is offline   Reply With Quote

Old   November 17, 2019, 13:33
Default
  #8
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hi Fahmida,

I am also using the interThermalPhaseChangeFoam framework for modelling bubble sliding. I was curious what boundary conditions for the p-rgh field you have been using for your simulations.

Kind regards,

Sam
saj216 is offline   Reply With Quote

Old   November 17, 2019, 14:07
Default
  #9
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
Hello,

In my computational domain consisting air and liquid has a fixedFluxPressure boundary condition for the inlet and side walls and a total pressure with p0 0 for the outflow.
Fahmida is offline   Reply With Quote

Old   November 19, 2019, 20:07
Default
  #10
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Thanks,

Did you find out why the bubble size did not increase in your simulation?
saj216 is offline   Reply With Quote

Old   November 20, 2019, 06:01
Default
  #11
Member
 
FaRa
Join Date: Jul 2019
Posts: 33
Rep Power: 7
Fahmida is on a distinguished road
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.
Fahmida is offline   Reply With Quote

Old   December 3, 2019, 23:41
Default
  #12
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Thank you for the help.

How did you initialise the perfect bubble shape especially on a perfectly smooth surface.

Sorry for asking, I'm not very good at OpenFoam
saj216 is offline   Reply With Quote

Old   December 4, 2019, 00:10
Default
  #13
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 8
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Thank you for the help.

How did you initialise the perfect bubble shape especially on a perfectly smooth surface.

Sorry for asking, I'm not very good at OpenFoam
I think there is a tutorial case in interThermalPhaseChangeFoam for single bubble growth.
BTW, the solver evapVOFHardt based on of231 is also useful for bubble growth simulation.
It should be noted that spurious current might be an issue in this kind of simulation
Bali is offline   Reply With Quote

Old   December 4, 2019, 01:15
Default
  #14
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hi,

Thanks, yeah there is but the simulation relies on a cavity in the mesh rather than a smooth surface. Btw have you had any issues changing boundary conditions of the thermalphasechangefoam solver. When I have the simulation has quickly collapsed.

For example for a vapour surrrounded in a superheated liquid scenario. I can provide for detail if required. Thanks again.
saj216 is offline   Reply With Quote

Old   December 4, 2019, 01:22
Default
  #15
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 8
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hi,

Thanks, yeah there is but the simulation relies on a cavity in the mesh rather than a smooth surface. Btw have you had any issues changing boundary conditions of the thermalphasechangefoam solver. When I have the simulation has quickly collapsed.

For example for a vapour surrrounded in a superheated liquid scenario. I can provide for detail if required. Thanks again.
Hi, Sam

it will be good to check your case if can upload it.

Best
Bali is offline   Reply With Quote

Old   December 8, 2019, 00:03
Default
  #16
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Hi Bali,

Sorry, for the late reply. This is the common error, I repeatedly get. I assume it is an issue of my pressure boundaries. If you understand why, could you explain to me as simply as possible .

--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/joseph/OpenFOAM/joseph-2.4.0/run/CFD-PC-master/CFD-PC/interThermalPhaseFoam/tutorials/VapourGrowthWithoutLayer/system/data.solverPerformance.p_rgh at line 0.

From function operator>>(Istream&, Scalar&)
in file lnInclude/Scalar.C at line 93.

Thank you for the help.

Sam
Attached Files
File Type: gz blockMeshDict.m4.tar.gz (893 Bytes, 7 views)
File Type: gz initialboundaries.tar.gz (922 Bytes, 9 views)
File Type: gz system.tar.gz (3.8 KB, 8 views)
saj216 is offline   Reply With Quote

Old   December 8, 2019, 00:33
Default
  #17
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 8
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Hi Bali,

Sorry, for the late reply. This is the common error, I repeatedly get. I assume it is an issue of my pressure boundaries. If you understand why, could you explain to me as simply as possible .

--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/joseph/OpenFOAM/joseph-2.4.0/run/CFD-PC-master/CFD-PC/interThermalPhaseFoam/tutorials/VapourGrowthWithoutLayer/system/data.solverPerformance.p_rgh at line 0.

From function operator>>(Istream&, Scalar&)
in file lnInclude/Scalar.C at line 93.

Thank you for the help.

Sam

hi Sam,

I have checked your files. Here some points:

1. why the "base" boundary is applied with "symmetry"? if it is a base wall boundary, "symmetry" is not appropriate.

2. I guess you can not use "fixedvalue 0" for "top" and "right" boundary at the same time. one of them could be changed into "zeroGradient"
Bali is offline   Reply With Quote

Old   December 8, 2019, 00:39
Default
  #18
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
Sorry for the confusion, the base in this case is symmetry because the bubble is initialised in the simulation so that it is a 1/4 of a bubble. So that the computing requirement can be reduced.

Okay I will try that and come back to you if it works.

Thanks again.
saj216 is offline   Reply With Quote

Old   December 8, 2019, 00:40
Default
  #19
New Member
 
Bali
Join Date: Oct 2018
Posts: 14
Rep Power: 8
Bali is on a distinguished road
Quote:
Originally Posted by saj216 View Post
Sorry for the confusion, the base in this case is symmetry because the bubble is initialised in the simulation so that it is a 1/4 of a bubble. So that the computing requirement can be reduced.

Okay I will try that and come back to you if it works.

Thanks again.
OK. but I guess "symmetryPlane" might be better
Bali is offline   Reply With Quote

Old   December 8, 2019, 22:47
Default
  #20
Member
 
Sam
Join Date: May 2019
Posts: 64
Rep Power: 7
saj216 is on a distinguished road
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 ?
saj216 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 07:21
TwoPhaseEulerFoam - Bubble diameter size Problem BlnPhoenix OpenFOAM Running, Solving & CFD 8 September 17, 2019 16:59
How to make UDFs for Nucleate Pool Boiling Simulation? SIKJAE Fluent UDF and Scheme Programming 1 August 4, 2018 08:07
OpenFOAM UpdateCoeffs Nucleate Wall Boiling suneth.warna OpenFOAM Programming & Development 0 July 31, 2016 18:42
Flow boiling & bubble formation Venkat CFX 1 July 23, 2009 10:20


All times are GMT -4. The time now is 11:01.