pith. sign in

arxiv: 2606.01654 · v1 · pith:XMNJIJUCnew · submitted 2026-06-01 · ⚛️ physics.flu-dyn

Interaction between vapor bubbles during flow boiling heat transfer in microchannels

Pith reviewed 2026-06-28 13:12 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords vapor bubblesflow boilingmicrochannelsbubble interactionheat transferphase changenumerical simulationReynolds number
0
0 comments X

The pith

In microchannel flow boiling, trailing vapor bubbles reduce leading bubble growth because they absorb heat through vaporization.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Numerical simulations examine multiple vapor bubbles in microchannel flow boiling and show that the leading bubble ends up smaller when rear bubbles are present than in an otherwise identical single-bubble case. The reduction occurs because vaporization at the rear bubbles takes up heat that would otherwise reach the leading bubble. The same simulations track how this effect changes with the initial volume ratio between bubbles, with flow Reynolds number, and with the thickness of the heated wall. The results indicate that bubble growth and coalescence timing depend on competition for the superheated thermal boundary layer.

Core claim

For the same initial size and position of the leading bubble, its final size in a single-bubble microchannel exceeds the size reached by the leading bubble when multiple bubbles are present, because vaporization at the rear bubbles absorbs heat. When the initial volume ratio of the leading bubble to the rear bubble is reduced, the leading bubble grows even smaller downstream because it contacts less of the superheated layer. Raising the Reynolds number produces modestly larger bubbles upstream since the bubbles reach the superheated fluid sooner, while increasing bottom-wall thickness raises upstream wall temperature by conduction and thereby speeds bubble growth and advances coalescence.

What carries the argument

Numerical simulation of the coupled heat transfer, phase change, and fluid flow around multiple vapor bubbles, used to isolate the effect of heat absorption by rear bubbles on leading-bubble growth.

If this is right

  • The leading bubble grows smaller as the number of trailing bubbles increases because of heat absorption at the rear sites.
  • A smaller initial volume ratio between leading and rear bubbles produces a smaller leading bubble downstream due to reduced contact with the superheated layer.
  • Higher Reynolds number yields slightly larger bubbles upstream because the bubbles encounter superheated fluid earlier.
  • Thicker bottom walls raise upstream wall temperature by conduction, increasing bubble growth rates and advancing coalescence.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The heat-competition mechanism could reduce net heat-removal rate in densely nucleating regimes compared with sparse nucleation.
  • Spacing or nucleation-site control might be used to tune the balance between bubble growth and coalescence in cooling designs.
  • The same interaction is likely to appear in other confined phase-change flows where multiple interfaces draw from a shared thermal boundary layer.

Load-bearing premise

The numerical model accurately captures the coupled heat transfer, phase change, and fluid flow physics without discretization or modeling errors large enough to change the reported bubble-interaction effects.

What would settle it

Side-by-side experimental images or volume measurements of the leading bubble in single-bubble versus multi-bubble flow boiling runs that match the simulated initial sizes, positions, Reynolds numbers, and wall thicknesses.

Figures

Figures reproduced from arXiv: 2606.01654 by Mengqi Ye, Odumuyiwa A. Odumosu, Tianyou Wang, Zhizhao Che.

Figure 1
Figure 1. Figure 1: Simulation setup for multiple bubble flow boiling in microchannels. (a) Sectional view of the channel wall to show two bubbles and the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: At the onset of the flow boiling process, the bubbles are transported in the flow direction by the saturated [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Time variation of the total bubble volume [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Snapshots of bubble growth and temperature fields at di [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Flow boiling in microchannels with different numbers of bubbles: (a, b) velocity fields at the middle cross-section in the z-direction; (c, d) temperature fields at the solid-fluid interface of the bottom wall. Panels (a, c) show results when the bubbles are in the adiabatic region (t = 1.2 ms), while panels (b, d) show results when the bubbles are in the heated region (t = 5 ms). (e) 4.8 ms (f) 5 ms y x O… view at source ↗
Figure 5
Figure 5. Figure 5: Temperature fields for microchannels with di [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Nusselt number Nubottom plotted versus the dimensionless position of the leading bubble nose xn/Dh for different numbers of bubbles. 0 1 2 3 4 5 0 5 10 15 20 25 30 35 Time (ms) One-bubble Two-bubble Three-bubble VL / V0 (-) (a) 0 1 2 3 4 5 0 5 10 15 20 25 30 Time (ms) One-bubble Two-bubble Three-bubble Vtotal / V0 (-) (b) 0 1 2 3 4 5 0 5 10 15 20 25 30 Time (ms) One-bubble Two-bubble Three-bubble xn / D h … view at source ↗
Figure 7
Figure 7. Figure 7: (a) Dimensionless leading bubble volume VL/V0, (b) dimensionless total bubble volume Vtotal/V0, and (c) dimensionless position of the leading bubble nose xn/Dh, plotted versus time for different numbers of bubbles. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Pressure distribution contour at the middle cross-section in z-direction for microchannels with di [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Snapshots of bubble growth and temperature fields at di [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Velocity field at middle cross-section in the [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Nusselt number Nubottom plotted versus the dimensionless position of the leading bubble nose xn/Dh for different bubble volume ratios. 11 [PITH_FULL_IMAGE:figures/full_fig_p011_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: (a) Dimensionless total bubble volume Vtotal/V0, (b) dimensionless position of the leading bubble nose xn/Dh, and (c) dimensionless trailing bubble volume VT /VT,0, plotted versus time for different bubble volume ratios. The sudden increase of VT /VT,0 indicates the occurrence of the bubble coalescence and is marked with a square symbol. 3.3. Effect of inlet Reynolds number The influence of the inlet Reyn… view at source ↗
Figure 13
Figure 13. Figure 13: Snapshots of bubble growth and temperature fields at di [PITH_FULL_IMAGE:figures/full_fig_p013_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Flow boiling in microchannels with different Re numbers for two-bubble cases: (a, b) temperature fields at the solid-fluid interface of the bottom wall; (c, d) velocity field at the middle cross-section in the z-direction. Panels (a, c) show results when the bubbles are in the adiabatic region (t = 1.5 ms), while panels (b, d) show results when the bubbles are in the heated region (t = 5.6 ms). Because of… view at source ↗
Figure 15
Figure 15. Figure 15: Nusselt number Nubottom plotted versus the dimensionless position of the leading bubble nose xn/Dh for different Re numbers. 0 1 2 3 4 5 6 0 2 4 6 8 10 12 14 16 18 Time (ms) (a) Re = 360 Re = 400 Re = 440 Re = 480 Vtotal / V0 (-) 0 1 2 3 4 5 6 0 5 10 15 20 25 30 xn / Dh (-) Time (ms) (b) Re = 360 Re = 400 Re = 440 Re = 480 [PITH_FULL_IMAGE:figures/full_fig_p014_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: (a) Dimensionless total bubble volume Vtotal/V0, (b) dimensionless position of the leading bubble nose xn/Dh, plotted versus time for different Re numbers. they can quickly catch up in the later stage because of the bubble expansion. 3.4. Effect of bottom wall thickness In practical applications, varying the bottom wall thickness determines the microchannel compactness, the ease of manufacturing process, … view at source ↗
Figure 17
Figure 17. Figure 17: Snapshots of bubble growth and temperature fields at di [PITH_FULL_IMAGE:figures/full_fig_p015_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Temperature distribution for different bottom wall thicknesses along the centre line of the solid-fluid interface of the heated bottom wall (z = 0). (a) Temperature distribution when the bubbles are in the adiabatic region (t = 1.8 ms), (b) Temperature distribution when the bubbles are in the heated region (t = 5.4 ms). thicknesses, the instants of the bubble coalescence are not the same. For microchannel… view at source ↗
Figure 19
Figure 19. Figure 19: Flow boiling in microchannels with different bottom wall thicknesses for two-bubble case: (a, b) temperature fields at the solid-fluid interface of the bottom wall; (c, d) velocity field at the middle cross-section in the z-direction. Panels (a, c) show results when the bubbles are in the adiabatic region (t = 1.8 ms), while panels (b, d) show results when the bubbles are in the heated region (t = 5.4 ms)… view at source ↗
Figure 20
Figure 20. Figure 20: (a) Dimensionless total bubble volume Vtotal/V0, (b) dimensionless position of the leading bubble nose xn/Dh, plotted versus time for different bottom wall thicknesses. to those with thinner bottom walls. The difference in the magnitude of the Nusselt number is small because the microchannel material adopted (copper) has both higher thermal conductivity and diffusivity, which make it to transfer heat effe… view at source ↗
Figure 21
Figure 21. Figure 21: Nusselt number Nubottom versus the dimensionless position of the leading bubble nose xn/Dh for different bottom wall thicknesses. higher Reynolds number results in a higher bubble moving speed, a stronger effect of convection, and a higher Nusselt number, in contrast at the downstream of the heated region, a higher Reynolds number results in a lower superheat degree, a slower expansion of the bubble size,… view at source ↗
read the original abstract

Microchannel flow boiling is an efficient cooling solution for high-power-density miniaturized systems. Many studies on microchannel flow boiling focused on the dynamics of single vapor bubbles, while neglecting the interaction between bubbles, which is important in relevant applications. Here, numerical simulations are carried out to study the interaction between multiple vapor bubbles in microchannel flow boiling. The results show that for different numbers of bubbles in the microchannels with the same initial size and position of leading bubbles, the bubble size in a single-bubble microchannel is larger compared to the leading bubble of multiple-bubble cases because of heat absorption by the vaporization at the rear bubbles. As the initial volume ratio between the leading bubble and the rear bubble decreases, the leading bubble size in the downstream becomes smaller because of the reduced contact with the superheated thermal boundary layer. With increasing the Reynolds number, both the leading and the trailing bubbles increase slightly in size in the upstream of the heated region, because the bubbles at higher Reynolds number move faster and firstly get in contact with the superheated fluid. The increase in the bottom wall thickness increases the growth rate of the multiple bubble sizes with earlier bubble coalescence because of the higher upstream wall temperature by heat conduction in the solid wall.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper reports numerical simulations of vapor bubble interactions during flow boiling in microchannels. It claims that, for fixed initial leading-bubble size and position, the leading bubble grows larger in single-bubble cases than in multi-bubble cases because rear bubbles absorb heat via vaporization; additional parametric trends are reported for initial volume ratio, Reynolds number, and bottom-wall thickness, including earlier coalescence at larger wall thickness.

Significance. If the numerical model is shown to be free of significant discretization or constitutive errors, the work would supply concrete mechanistic insight into multi-bubble heat-transfer coupling that is currently missing from the predominantly single-bubble literature, with direct relevance to microchannel cooling design.

major comments (2)
  1. [Abstract/Results] Abstract and Results: the central claim that rear-bubble vaporization reduces leading-bubble growth rests entirely on the numerical data, yet the manuscript supplies no experimental validation, mesh-independence study, or boundary-condition verification; without these the reported interaction mechanism cannot be assessed for numerical artifact.
  2. [Numerical Methods] Numerical Methods (assumed section describing the solver): no information is given on the phase-change model (e.g., sharp-interface vs. diffuse), the treatment of the solid-fluid conjugate heat transfer, or the time-step and spatial discretization criteria; these choices directly affect whether the reported Reynolds-number and wall-thickness trends are physically reliable.
minor comments (2)
  1. [Abstract] Abstract: the long sentence describing Reynolds-number effects is difficult to parse; splitting it would improve clarity.
  2. [Results] The manuscript does not state the initial bubble shapes or the precise definition of 'bubble size' used for the reported comparisons.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments on our numerical study of vapor bubble interactions in microchannel flow boiling. We address each major comment below.

read point-by-point responses
  1. Referee: [Abstract/Results] Abstract and Results: the central claim that rear-bubble vaporization reduces leading-bubble growth rests entirely on the numerical data, yet the manuscript supplies no experimental validation, mesh-independence study, or boundary-condition verification; without these the reported interaction mechanism cannot be assessed for numerical artifact.

    Authors: Our work is a purely numerical investigation intended to provide mechanistic insight into multi-bubble coupling that is difficult to resolve experimentally. We agree that mesh-independence and boundary-condition verification are essential for credibility and will add these studies (including grid convergence plots and domain-size checks) to the revised manuscript. Experimental validation lies outside the current scope. revision: partial

  2. Referee: [Numerical Methods] Numerical Methods (assumed section describing the solver): no information is given on the phase-change model (e.g., sharp-interface vs. diffuse), the treatment of the solid-fluid conjugate heat transfer, or the time-step and spatial discretization criteria; these choices directly affect whether the reported Reynolds-number and wall-thickness trends are physically reliable.

    Authors: We apologize for the insufficient detail. The revised manuscript will expand the Numerical Methods section to specify the sharp-interface phase-change model (volume-of-fluid with explicit mass transfer), the conjugate heat-transfer treatment at the solid-fluid interface, and the chosen time-step and spatial discretization criteria together with convergence checks. revision: yes

standing simulated objections not resolved
  • Experimental validation of the reported bubble-interaction mechanism

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper reports outcomes from direct numerical simulations of multi-bubble flow boiling without any analytical derivation chain, fitted parameters, or self-citations that reduce claims to inputs by construction. All reported trends (bubble-size differences due to heat absorption, Reynolds-number effects, wall-thickness effects) are stated as simulation results rather than predictions derived from prior equations or ansatzes within the work. No load-bearing step matches any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the simulation framework itself is not described.

pith-pipeline@v0.9.1-grok · 5761 in / 1008 out tokens · 23533 ms · 2026-06-28T13:12:43.903307+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 31 canonical work pages

  1. [1]

    O. A. Odumosu, H. Xu, T. Wang, Z. Che, Growth of elongated vapor bubbles during flow boiling heat transfer in wavy microchannels, Applied Thermal Engineering 223 (2023) 119987.doi:10.1016/j.applthermaleng. 2023.119987

  2. [3]

    Zhang, Y

    X. Zhang, Y . Zhang, C. Yang, J. Liu, B. Liu, Z. Zhou, Study of flow and heat transfer performance of nanofluids in curved microchannels with different curvatures, Numerical Heat Transfer, Part A: Applications (2024).doi: 10.1080/10407782.2024.2321520. 18

  3. [4]

    Baldassari, M

    C. Baldassari, M. Marengo, Flow boiling in microchannels and microgravity, Progress in Energy and Combustion Science 39 (1) (2013) 1–36.doi:10.1016/j.pecs.2012.10.001

  4. [5]

    S. S. Bertsch, E. A. Groll, S. V . Garimella, Review and comparative analysis of studies on saturated flow boiling in small channels, Nanoscale and Microscale Thermophysical Engineering 12 (3) (2008) 187–227.doi:10. 1080/15567260802317357

  5. [6]

    S. V . Garimella, C. B. Sobhan, Transport in microchannels - a critical review, Annual Review of Heat Transfer 13 (2003) 1–50.doi:10.1615/annualrevheattransfer.v13.30

  6. [7]

    Z. Guo, D. Fletcher, B. Haynes, A review of computational modelling of flow boiling in microchannels, The Journal of Computational Multiphase Flows 6 (2) (2014) 79–110.doi:10.1260/1757-482x.6.2.79

  7. [8]

    Harirchian, S

    T. Harirchian, S. V . Garimella, Boiling heat transfer and flow regimes in microchannels—a comprehensive understanding, Journal of Electronic Packaging 133 (1) (2011) 011001.doi:10.1115/1.4002721

  8. [9]

    S. T. Kadam, I. Hassan, R. Kumar, M. A. Rahman, Bubble dynamics in microchannel: An overview of the state-of-the-art, Meccanica 56 (3) (2021) 481–513.doi:10.1007/s11012-020-01300-4

  9. [10]

    M. M. Nahar, B. Ma, K. Guye, Q. H. Chau, J. Padilla, M. Iyengar, D. Agonafer, Review article: Microscale evaporative cooling technologies for high heat flux microelectronics devices: Background and recent advances, Applied Thermal Engineering 194 (2021) 117109.doi:10.1016/j.applthermaleng.2021.117109

  10. [11]

    J. R. Thome, Boiling in microchannels: a review of experiment and theory, International Journal of Heat and Fluid Flow 25 (2) (2004) 128–139.doi:10.1016/j.ijheatfluidflow.2003.11.005

  11. [13]

    Gedupudi, Y

    S. Gedupudi, Y . Q. Zu, T. G. Karayiannis, D. B. R. Kenning, Y . Y . Yan, Confined bubble growth during flow boiling in a mini/micro-channel of rectangular cross-section part i: Experiments and 1-d modelling, International Journal of Thermal Sciences 50 (3) (2011) 250–266.doi:10.1016/j.ijthermalsci.2010.09.001

  12. [14]

    Halon, Z

    S. Halon, Z. Krolicki, R. Revellin, B. Zajaczkowski, Heat transfer characteristics of flow boiling in a micro channel array with various inlet geometries, International Journal of Heat and Mass Transfer 187 (2022) 122549. doi:10.1016/j.ijheatmasstransfer.2022.122549

  13. [15]

    Kumar, G

    R. Kumar, G. Singh, D. Mikielewicz, Effect of asymmetric fluid flow distribution on flow boiling in a microchannel heat sink – an experimental investigation, Applied Thermal Engineering 213 (2022) 118710. doi:10.1016/j.applthermaleng.2022.118710

  14. [18]

    Zhang, J

    Y . Zhang, J. Wang, W. Liu, Z. Liu, Heat transfer and pressure drop characteristics of r134a flow boiling in the parallel/tandem microchannel heat sinks, Energy Conversion and Management 148 (2017) 1082–1095.doi: 10.1016/j.enconman.2017.06.067

  15. [19]

    S. Zhou, B. Shu, Z. Yu, Y . Huang, Y . Zhang, Experimental study and mechanism analysis of the flow boiling and heat transfer characteristics in microchannels with different surface wettability, Micromachines 12 (8) (2021) 881.doi:10.3390/mi12080881. 19

  16. [20]

    Y . Luo, W. Li, K. Zhou, K. Sheng, S. Shao, Z. Zhang, J. Du, W. J. Minkowycz, Three-dimensional numerical simulation of saturated annular flow boiling in a narrow rectangular microchannel, International Journal of Heat and Mass Transfer 149 (2020) 119246.doi:10.1016/j.ijheatmasstransfer.2019.119246

  17. [21]

    Mukherjee, S

    A. Mukherjee, S. G. Kandlikar, Z. J. Edel, Numerical study of bubble growth and wall heat transfer during flow boiling in a microchannel, International Journal of Heat and Mass Transfer 54 (15-16) (2011) 3702–3718. doi:10.1016/j.ijheatmasstransfer.2011.01.030

  18. [22]

    Zhuan, W

    R. Zhuan, W. Wang, Flow pattern of boiling in micro-channel by numerical simulation, International Journal of Heat and Mass Transfer 55 (5) (2012) 1741–1753.doi:10.1016/j.ijheatmasstransfer.2011.11.029

  19. [23]

    Ferrari, M

    A. Ferrari, M. Magnini, J. R. Thome, Numerical analysis of slug flow boiling in square microchannels, International Journal of Heat and Mass Transfer 123 (2018) 928–944.doi: 10.1016/j.ijheatmasstransfer.2018.03.012

  20. [24]

    Kaushik, N

    M. Magnini, O. K. Matar, Numerical study of the impact of the channel shape on microchannel boiling heat transfer, International Journal of Heat and Mass Transfer 150 (2020) 119322.doi:10.1016/j. ijheatmasstransfer.2020.119322

  21. [25]

    Magnini, B

    M. Magnini, B. Pulvirenti, J. R. Thome, Numerical investigation of hydrodynamics and heat transfer of elongated bubbles during flow boiling in a microchannel, International Journal of Heat and Mass Transfer 59 (2013) 451– 471.doi:10.1016/j.ijheatmasstransfer.2012.12.010

  22. [26]

    Z. Guo, B. S. Haynes, D. F. Fletcher, Numerical simulation of annular flow boiling in microchannels, The Journal of Computational Multiphase Flows 8 (1) (2016) 61–82.doi:10.1177/1757482X16634205

  23. [27]

    A. Priy, I. Ahmad, M. K. Khan, M. Pathak, Bubble interaction and heat transfer characteristics of microchannel flow boiling with single and multiple cavities, Journal of Thermal Science and Engineering Applications 16 (6) (2024) 061010.doi:10.1115/1.4065187

  24. [28]

    Widyatama, M

    A. Widyatama, M. Venter, D. Orejon, K. Sefiane, Experimental investigation of bubble dynamics and flow patterns during flow boiling in high aspect ratio microchannels with the effect of flow orientation, International Journal of Thermal Sciences 189 (2023) 108238.doi:10.1016/j.ijthermalsci.2023.108238

  25. [29]

    Q. Liu, B. Palm, Numerical study of bubbles rising and merging during convective boiling in micro-channels, Applied Thermal Engineering 99 (2016) 1141–1151.doi:10.1016/j.applthermaleng.2016.01.116

  26. [30]

    Consolini, J

    L. Consolini, J. R. Thome, A heat transfer model for evaporation of coalescing bubbles in micro-channel flow, International Journal of Heat and Fluid Flow 31 (1) (2010) 115–125.doi:10.1016/j.ijheatfluidflow. 2009.10.005

  27. [31]

    Magnini, B

    M. Magnini, B. Pulvirenti, J. R. Thome, Numerical investigation of the influence of leading and sequential bubbles on slug flow boiling within a microchannel, International Journal of Thermal Sciences 71 (2013) 36–52. doi:10.1016/j.ijthermalsci.2013.04.018

  28. [32]

    Q. Liu, W. Wang, B. Palm, C. Wang, X. Jiang, On the dynamics and heat transfer of bubble train in micro- channel flow boiling, International Communications in Heat and Mass Transfer 87 (2017) 198–203.doi: 10.1016/j.icheatmasstransfer.2017.07.002

  29. [33]

    Lombaard, M

    L. Lombaard, M. A. Moghimi, P. Valluri, J. P. Meyer, Interaction between multiple bubbles in microchannel flow boiling and the effects on heat transfer, International Communications in Heat and Mass Transfer 129 (2021) 105703.doi:10.1016/j.icheatmasstransfer.2021.105703

  30. [34]

    Scheufler, J

    H. Scheufler, J. Roenby, Twophaseflow: A framework for developing two phase flow solvers in openfoam, OpenFOAM Journal 3 (2023) 200–224.doi:10.51560/ofj.v3.80. 20

  31. [35]

    J. U. Brackbill, D. B. Kothe, C. Zemach, A continuum method for modeling surface tension, Journal of Computational Physics 100 (2) (1992) 335–354.doi:10.1016/0021-9991(92)90240-y

  32. [36]

    Hardt, F

    S. Hardt, F. Wondra, Evaporation model for interfacial flows based on a continuum-field representation of the source terms, Journal of Computational Physics 227 (11) (2008) 5871–5895.doi:10.1016/j.jcp.2008.02. 020

  33. [37]

    Scriven, On the dynamics of phase growth, Chemical Engineering Science 10 (1) (1959) 1–13.doi:10

    L. Scriven, On the dynamics of phase growth, Chemical Engineering Science 10 (1) (1959) 1–13.doi:10. 1016/0009-2509(59)80019-1. 21