Pith. sign in

REVIEW 3 major objections 5 minor 81 references

A Population Balance Model for Large Eddy Simulation of Polydisperse Droplet Evolution

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read A breakup model extends droplet-size prediction into the viscous range and matches jet-in-crossflow data.

desk verdict The viscous-range kernel extension and LES-PBE application are genuinely new and useful; the validation is shape-only and the SGS closure is untested, but the paper deserves review. read the letter →

arxiv 1908.02397 v1 pith:Y77CJH2P submitted 2019-08-06 physics.flu-dyn physics.comp-ph

classification physics.flu-dynphysics.comp-ph
keywords dropletbreakuppopulationbalanceequationlargeeddysimulationpolydispersedropletsviscousrangeBatchelorstructurefunctionjetincrossflowoilspillmodeling
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper develops a population balance model for polydisperse droplet evolution that is built to run inside large eddy simulations, and tests it against lab experiments of oil droplets in turbulence. The central claim is that a breakup kernel based on droplet-eddy collisions, extended from the inertial to the viscous range with a Batchelor structure function and one fitted constant, can predict the relative droplet size distribution of a turbulent jet in crossflow. If true, the model gives engineers and oceanographers a practical way to compute how turbulent transport and breakup reshape oil droplet distributions, including the intermittent variability that mean-flow models miss.

What carries the argument

The central object is the breakup frequency integral g(di) = K ∫ π/4(di+de)^2 u_e(de) Ω(di,de) dn_e(de), evaluated with droplet-eddy collisions treated kinetically. The key extension is replacing the inertial-range eddy velocity with the Batchelor structure function S2(r) = C2 $ε^{{2/3}}$ $r^{{2/3}}$ [1 + (r/(γ2 η))^{-2}]^{-2/3}, which transitions smoothly to the viscous range; the parameter γ2 = (15 C2)^{3/4} ≈ 13 sets the crossover scale. A fast analytic fit for the integral makes the breakup rate cheap enough to evaluate at every LES grid point and timestep.

What would settle it

Run the same LES breakup model at a much finer grid resolution (or in a DNS-resolved setting) for the identical jet-in-crossflow case and compare the predicted size distributions; if the coarse-grid and fine-grid results diverge by more than the reported experimental scatter, the subgrid closure and the fitted K* are not transferable across resolutions.

Watch

Extended reading notes

Core claim

The paper claims that turbulent breakup of droplets in the viscous range, where the Kolmogorov scale is larger than the droplet, can be captured with a single breakup-frequency model that smoothly blends inertial and viscous scalings. Using a Batchelor blending function for the eddy velocity, the breakup frequency becomes a function of a local Reynolds number, Ohnesorge number, and density-viscosity ratio, with one fitted constant K* = 0.2 tuned on a breaking-wave experiment. Applied to a three-dimensional LES of a jet in crossflow with 15 droplet size bins, the model reproduces the measured relative size distribution and reveals that the Sauter mean diameter and total surface area fluctuate strongly and non-Gaussianly in time and space.

Load-bearing premise

The paper assumes that subgrid fluctuations of dissipation rate and droplet concentration are uncorrelated, so the filtered breakup source can be computed from filtered quantities alone; if these correlations matter at high Reynolds numbers, the predicted breakup rates will be biased.

Editorial extensions

If this is right

  • The LES model can predict polydisperse droplet transport and breakup in flows where the Kolmogorov scale exceeds the droplet size, which earlier inertial-range-only models overestimate.
  • The fitted constant K* = 0.2 transfers from a breaking-wave experiment to a round jet and then to a jet in crossflow, suggesting the breakup kernel is flow-independent to first order.
  • LES outputs yield temporal and spatial statistics of droplet surface area and Sauter mean diameter, quantities that control biodegradation rates and optical properties of oil plumes.
  • The model's insensitivity to the assumed initial droplet injection size (monodisperse vs bidisperse) suggests far-field size distributions are set by breakup dynamics rather than nozzle details.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the subgrid closure holds, the same kernel could be applied to bubble plumes and spray atomization downstream of primary breakup, where coalescence remains negligible.
  • The fitted K* may absorb the uncertainty in the upper integration cutoff and daughter-size probability; a systematic study varying both could separate physical breakup physics from closure tuning.
  • The strong non-Gaussian variability of d32 and total surface area implies that single-point mean predictions of oil fate may understate the risk of locally high surface area available for biodegradation.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper develops a population balance model for polydisperse droplet evolution in an LES framework. The breakup frequency is modeled as a droplet–eddy collision integral with a Batchelor structure function that smoothly bridges inertial and viscous ranges, and is parameterized by algebraic fits. The one free constant K* is calibrated against the Li et al. (2017) breaking-wave experiment, tested against Eastwood et al. (2004) jet data, and used in a 3D LES of a jet in crossflow. The LES predictions are compared with Murphy et al. (2016) data using volume-normalized relative size/volume distributions. The paper also reports intermittency statistics of Sauter mean diameter, total surface area, and breakup source terms.

Significance. The viscous-range extension via a structure function is a useful physical improvement over standard inertial-range kernels, and the algebraic fits make LES application affordable. The paper is honest about limitations: it notes the subgrid-closure assumption, the limited statistical convergence of Murphy et al., and tests sensitivity to the injection distribution and daughter-probability model. If the remaining validation issues are resolved, the model would be a practical tool for oil-spill and spray applications, where polydisperse droplet evolution matters. The current evidence supports the relative-shape predictions in the specific Murphy et al. configuration, but not yet a generally transferable quantitative prediction.

major comments (3)
  1. [§4.1, Eq. (4.3)] The closure \widetilde{g n}\approx\tilde{g}\tilde{n} and g(\epsilon)\approx g(\tilde{\epsilon}) is load-bearing, since the grid scale is 5.2 mm while \eta\approx 13–60 \mu m (§4.3), so breakup rates are set by unresolved eddies. The breakup frequency depends nonlinearly on \epsilon via Eqs. (2.17)–(2.18), and Fig. 17 shows grid-scale \epsilon varying over orders of magnitude, so Jensen-type bias is expected. The 1D calibrations (§3.1, §3.2) use mean or parameterized dissipation and do not exercise this closure. Because no grid-resolution or SGS-model sensitivity is reported, the LES comparisons in Figs. 11–13 cannot yet be regarded as a clean test of the model; the authors should quantify the closure error or at least demonstrate resolution insensitivity.
  2. [§3.2, Fig. 6] The claim that K*=0.2 is robust is weakened by the Eastwood comparison: for heptane and 10 cSt silicone oil the model with K*=0.2 overpredicts decay, and K*=0.1 and 0.15 are needed. Since no parameter dependence is identified, the fitted constant is not shown to be transferable across fluids. Please either explain why these discrepancies are within model uncertainty or restrict the transferability claim and assess the impact on the Murphy comparison in §4.
  3. [§4.3, Figs. 11–13] The total oil concentrations differ between simulation and experiment by factors of 1.4 and 3.7, and the comparisons are normalized by total volume concentration. Therefore the abstract's claim of 'good agreement for the relative size distribution' is accurate only for the shape, not the absolute concentration. The model also predicts absolute concentration, so this discrepancy is a substantive part of the validation and should be discussed quantitatively (e.g., statistical convergence of Murphy et al., measurement volume, effective injection in Appendix B) rather than only used to justify a shape-only comparison.
minor comments (5)
  1. [Figure 13 caption] The caption states the comparison is 'at x = 1.3 m', but the text describes the total size distribution over the window x = 0.76–1.66 m; please align the caption with the text.
  2. [Eqs. (2.20)–(2.21)] Using the symbol e both as Euler's number and as the coefficient e(y) in the fit is confusing; consider renaming the coefficient (for example, f or h).
  3. [Eq. (3.4)] The error measure E uses the squared logarithmic difference but does not state whether absolute values or signed differences are intended; please clarify.
  4. [Appendix A] The notation log10[-c(y)] = ... is unusual because c(y) appears in the fit G = a x^b + c x^d - e; please clarify the signs and the domain of validity of the fit.
  5. [Acknowledgements heading] The word 'Acknowledegements' is misspelled; it should be 'Acknowledgements'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: K* is fitted to an external breaking-wave experiment and the LES size-distribution comparison is an out-of-sample relative-shape prediction.

full rationale

The paper's central claim is that a population balance breakup model with one fitted constant K*=0.2 predicts relative droplet size distributions in an LES of a jet in crossflow. The constant is calibrated in §3.1 against the Li et al. (2017) breaking-wave experiment by minimizing the error measure E in Eq. (3.4) within a 1D vertical-diffusion model, and then used without refitting in the Eastwood et al. comparison (§3.2) and the Murphy et al. LES (§4). The Murphy et al. comparison therefore tests the model's shape prediction outside the calibration dataset; nothing in the fit fixes the shape of the normalized size distribution N_i* at x=0.76–1.66 m. The normalization of both experiment and simulation by total volume in Eq. (4.8) weakens the comparison but does not make it an identity, since the relative bin-to-bin shape is still produced by the coupled transport and breakup dynamics. The §4.1 closure ∼gn≈g̃ñ and g(ε)≈g(ε̃) is an acknowledged approximation explicitly flagged by the authors as needing further subgrid modeling at very high Reynolds numbers; that is a robustness/correctness limitation, not a circular step. The Appendix A fits are algebraic surrogates for numerical evaluation of the authors' own integral (A 1), used for computational speed; fitting a fast formula to one's own model is not importing the target result. Eσ in Eq. (2.13) is computed from the same formation-energy and daughter-probability model, but this is internal consistency among model components, not a derivation of a prediction from itself. Self-citations to Yang et al. (2014, 2015, 2016) and Bou-Zeid et al. (2005) concern the LES code, SGS model, and droplet velocity expansion; they are infrastructure citations and are not load-bearing for the breakup-model claim. The Eastwood et al. comparison showing that K*=0.2 is suboptimal for heptane and 10cSt silicone oil is reported by the authors themselves (§3.2) and indicates transferability risk, but that is a model-adequacy issue, not circularity. On the supplied evidence, there is no step in which a 'prediction' or 'first-principles result' is equivalent by construction to its input data or to a self-citation chain.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The model's predictive content reduces to one fitted parameter K*, plus a set of numerical curve-fit coefficients, and relies on standard turbulence closures (K41, Batchelor blending, Smagorinsky/LASD SGS) and phenomenological breakup assumptions from the literature. No new physical entities are introduced.

free parameters (5)
  • K* = 0.2
    Dimensionless breakup prefactor fitted by minimizing the log-density error E in eq (3.4) against the Li et al. (2017) breaking-wave data. The same paper's Eastwood jet test shows 0.1 and 0.15 fit better for heptane and 10 cSt silicone oil.
  • kD = 0.3
    Coefficient in the turbulent diffusion model D(t)=kD u'(t)L(t) used in the 1D wave-breakup fit; chosen in accordance with Li et al. (2017), not independently determined here.
  • dmin = 1 µm
    Minimum daughter droplet diameter introduced in eq (2.6) to bound the breakup-probability model; chosen by hand.
  • alpha = approx 2
    Prefactor in the viscous resistive energy eq (2.14), taken from Calabrese et al. (1986) and Skartlien et al. (2013); approximate constant.
  • breakup-integral fit coefficients (A3) = Tables 4-5
    Coefficients a1-a4, c1-c4, d1-d4, e1-e4 fitted with MATLAB curve fitting to the numerically evaluated integral gf for fixed Gamma; they are numerical surrogates rather than physical parameters.
assumptions (6)
  • domain assumption K41 inertial-range scaling and the second-order structure function S2(r) with C2=2.1 and gamma2=13
    Standard turbulence theory used for eddy fluctuation velocities; invoked in eq (2.10) and throughout the breakup kernel.
  • ad hoc to paper Droplets and eddies interact as in kinetic theory of gases, with a sharp cutoff at eddy size equal to droplet size
    Stated in §2.2 as a crucial assumption; the cutoff dependence is lumped into K*.
  • domain assumption Daughter size distribution follows the Tsouris and Tavlarides energy-based U-shaped model, with volume conservation by complementary bins
    Used in eq (2.7) to define the breakup probability; the paper shows a Gaussian alternative does not fit the wave data.
  • domain assumption Dispersed phase volume fraction is low enough to neglect coalescence and treat droplets as spherical active scalars
    Sets Sc=0 in eq (1.1) and justifies the Eulerian scalar treatment in §1 and §4.1.
  • ad hoc to paper Subgrid closure \widetilde{g n}≈\tilde{g}\tilde{n} and g(\tilde{ε}) for the filtered breakup source
    Explicitly stated in §4.1 as an approximation that may need refinement at high Reynolds numbers.
  • domain assumption Ferry and Balachandar expansion valid for grid Stokes number St_delta << 1
    Used in eq (4.4) for droplet transport velocity; validity requires small droplet response time relative to resolved eddies.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Population Balance Model for Large Eddy Simulation of Polydisperse Droplet Evolution." pith.science (2026). https://pith.science/paper/Y77CJH2P

@misc{pith2026190802397,
  author       = {Pith},
  title        = {Pith review of: A Population Balance Model for Large Eddy Simulation of Polydisperse Droplet Evolution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y77CJH2P}},
  note         = {Machine review of arXiv:1908.02397}
}
abstract

In the context of many applications of turbulent multi-phase flows, knowledge of the dispersed phase size distribution and its evolution is critical to predicting important macroscopic features. We use a population dynamics model for polydisperse droplet distributions specifically adapted to a LES framework including a model for droplet breakup due to turbulence, neglecting coalescence. Following earlier methods used in the Reynolds averaged Navier--Stokes framework, the droplet breakup due to turbulent fluctuations is modelled by treating droplet-eddy collisions as in kinetic theory of gases. In order to also model smaller droplets comparable to or smaller than the Kolmogorov scale we extend the breakup kernels using a structure function model that smoothly transitions from the inertial to the viscous range. The model includes a dimensionless coefficient that is fitted by comparing predictions in a one-dimensional version of the model with a laboratory experiment of oil droplet breakup below breaking waves. After initial comparisons of the one-dimensional model to measurements of oil droplets in an axisymmetric jet, it is then applied in a three-dimensional LES of a jet in crossflow with large oil droplets of a single size being released at the source of the jet. We model the concentration fields using $N_d =15$ bins of discrete droplet sizes and solve scalar transport equations for each bin. The resulting droplet size distributions are compared with published experimental data, and good agreement for the relative size distribution is obtained. The LES results also enable us to quantify size distribution variability. We find that the probability distribution functions of key quantities such as the total surface area and the Sauter mean diameter of oil droplets are highly variable, some displaying strong non-Gaussian intermittent behavior.

Figures

Figures reproduced from arXiv: 1908.02397 by the authors.

Figure 1
Figure 1. Left: Sketch of the wave breaking experiment of Li [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Time evolution of the droplet size distribution for two values [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Average square error between predicted and measured [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figures from the paper (22 more)
Figure 4
Figure 4. Figure 4: Time evolution of the droplet size distribution for [PITH_FULL_IMAGE:figures/full_fig_p014_4.png]
Figure 5
Figure 5. Figure 5: We depict the effect of changing the maximum diameter size o [PITH_FULL_IMAGE:figures/full_fig_p015_5.png]
Figure 6
Figure 6. Figure 6: Evolution of the N(dmax) with downstream distance for the four different oils in 1. The symbols correspond to the experimental data of the different oils, 50 cSt Silicone oil ( ), Olive oil ( ), 10 cSt Silicone oil ( ) and Heptane (∗). The different lines correspond to…
Figure 7
Figure 7. Figure 7: Sketch of the simulation domain and dimensions. [PITH_FULL_IMAGE:figures/full_fig_p020_7.png]
Figure 8
Figure 8. Figure 8: Contour plots of instantaneous concentration fields at t [PITH_FULL_IMAGE:figures/full_fig_p021_8.png]
Figure 9
Figure 9. Figure 9: Time averaged concentration fields at the midplane of the j [PITH_FULL_IMAGE:figures/full_fig_p022_9.png]
Figure 10
Figure 10. Figure 10: Contour plot of the logarithm of the averaged total volu [PITH_FULL_IMAGE:figures/full_fig_p023_10.png]
Figure 11
Figure 11. Figure 11: Comparison of LES model at x = 0.76 m for mono-dispersed injection ( ) and bi-dispersed injection ( ), and experimental data from (Murphy et al. 2016) measured at the corresponding time (∗). Left panel: Relative size distribution from LES, Right panel: Relative volume…
Figure 12
Figure 12. Figure 12: Comparison of LES model at x = 1.3 m ( ) , and experimental data from (Murphy et al. 2016) measured at the corresponding time (∗). Left panel: Relative size distribution from LES, Right panel: Relative volume distribution. to a higher concentration at the top end. We …
Figure 13
Figure 13. Figure 13: Comparison of LES model at x = 1.3 m ( ) , and experimental data from (Murphy et al. 2016) measured at the corresponding time (∗). Left panel: Relative size distribution from LES, Right panel: Relative volume distribution. 0 0.5 1 1.5 2 x (m) 0 0.2 0.4 0.6 0.8 1 z c 5…
Figure 14
Figure 14. Figure 14: Left panel: evolution of centroid of various droplet plume [PITH_FULL_IMAGE:figures/full_fig_p025_14.png]
Figure 15
Figure 15. Figure 15: Average d32 diameter as function of downstream distance x at various plume heights. The lines correspond to z = 0.55 m ( ), z = 0.50 m( ) and z = 0.45 m( ) breakup, though the rate of change diminishes and appears to reach a stationary scale of about d32 ≈ 300 µm at l…
Figure 16
Figure 16. Figure 16: Representative time signals (left panels) and histograms [PITH_FULL_IMAGE:figures/full_fig_p027_16.png]
Figure 17
Figure 17. Figure 17: Representative time signals (left panels) and histograms [PITH_FULL_IMAGE:figures/full_fig_p028_17.png]
Figure 18
Figure 18. Figure 18: Representative time signals (left panels) and histograms [PITH_FULL_IMAGE:figures/full_fig_p028_18.png]
Figure 19
Figure 19. Figure 19: Representative time signals (left panels) and histograms [PITH_FULL_IMAGE:figures/full_fig_p029_19.png]
Figure 20
Figure 20. Figure 20: Time history of S˜ b,i normalized by concentration for different droplet sizes at z = 0.56 m on the centerline. (a) and (b) represent the droplet of size 1000 µm at two different x locations, (c) is the time history for d = 432 µm and (d) is for d = 20 µm. Dotted line…
Figure 21
Figure 21. Figure 21: The fit is represented by the dashed lines whereas the nu [PITH_FULL_IMAGE:figures/full_fig_p032_21.png]
Figure 22
Figure 22. Figure 22: The fit is represented by the dashed lines whereas the nu [PITH_FULL_IMAGE:figures/full_fig_p032_22.png]
Figure 23
Figure 23. Figure 23: The fit is represented by the dashed lines computed by int [PITH_FULL_IMAGE:figures/full_fig_p034_23.png]
Figure 24
Figure 24. Figure 24: Left panel: Variation of the mean centerline velocity with a [PITH_FULL_IMAGE:figures/full_fig_p036_24.png]
Figure 25
Figure 25. Figure 25: Sketch depicting the nozzle placement in LES (at vertical [PITH_FULL_IMAGE:figures/full_fig_p036_25.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

81 extracted references · 77 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    Albertson, J. D. & Parlange, M. B. 1999 Surface length scales and shear stress: Implications for land-atmosphere interaction over complex terrain . Water Resour. Res. 35 (7), 2121--2132

  4. [4]

    1981 Two Phase Flows in Chemical Engineering

    Azbel , D. 1981 Two Phase Flows in Chemical Engineering . Cambridge University Press

  5. [5]

    Bandara, U. C. & Yapa, P. D. 2011 Bubble sizes, breakup, and coalescence in deepwater gas/oil plumes . J. Hydraul. Eng. 137 (7), 729--738

  6. [6]

    Batchelor, G. K. 1951 Pressure fluctuations in isotropic turbulence . Proc. Cambr. Phil. Soc. 47 (July 1950), 359--374

  7. [7]

    & Peck, R.E

    Bossard, J.A. & Peck, R.E. 1996 Droplet sie distribution effects in spray combustion . Symposium (International) on Combustion 26 (1), 1671 -- 1677

  8. [8]

    , Meneveau, C

    Bou-Zeid, E. , Meneveau, C. & Parlange, M. 2005 A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows . Phys. Fluids 17 (2), 025105

Show all 81 references
  1. [9]

    Brandvik, P. J. , Johansen, . , Leirvik, F. , Farooq, U. & Daling, P. S. 2013 Droplet breakup in subsurface oil releases - Part 1: Experimental study of droplet breakup and effectiveness of dispersant injection . Mar. Pollut. Bull. 73 (1), 319--326

  2. [10]

    Calabrese, R. V. , Wang, C. Y. & Bryner, N. P. 1986 Drop breakup in turbulent stirred tank contactors. Part I: Effect of Dispersed-Phase Viscosity . AlChE J. 32 (4), 677--681

  3. [11]

    , Meneveau, C

    Chamecki, M. , Meneveau, C. & Parlange, M. B. 2008 A hybrid spectral/finite-volume algorithm for large-eddy simulation of scalars in the atmospheric boundary layer . Boundary Layer Meteorol. 128 (3), 473--484

  4. [12]

    & James, M

    Chatzi, E. & James, M. L. 1987 Analysis of Interactions for Liquid—Liquid Dispersions in Agitated Vessels . Ind. Eng. Chem. Res. 26 (11), 2263--2267

  5. [13]

    , Yang, D

    Chen, B. , Yang, D. , Meneveau, C. & Chamecki, M. 2018 Numerical study of the effects of chemical dispersant on oil transport from an idealized underwater blowout . Phys. Rev. Fluids 3 , 083801

  6. [14]

    & Karagozian, A

    Cortelezzi, L. & Karagozian, A. R. 2001 On the formation of the counter-rotating vortex pair in transverse jets . J. Fluid Mech. 446 , 347--373

  7. [15]

    Coulaloglou, C. A. & Tavlarides, L. L. 1977 Description of interaction processes in agitated liquid-liquid dispersions . Chem. Eng. Sci. 32 (11), 1289--1297

  8. [16]

    , Moureau, V

    Desjardins, O. , Moureau, V. & Pitsch, H. 2008 An accurate conservative level set/ghost fluid method for simulating turbulent atomization . J. Comput. Phys. 227 (18), 8395 -- 8416

  9. [17]

    , Luret, G

    Duret, B. , Luret, G. , Reveillon, J. , Menard, T. , Berlemont, A. & Demoulin, F.X. 2012 Dns analysis of turbulent mixing in two-phase flows . Int. J. Multiphase Flow 40 , 93 -- 105

  10. [18]

    Eastwood, C. D. , Armi, L. & Lasheras, J. C. 2004 The breakup of immiscible fluids in turbulent flows . J. Fluid Mech. 502 , 309–333

  11. [19]

    & Balachandar, S

    Ferry, J. & Balachandar, S. 2001 A fast eulerian method for disperse two-phase flow . Int. J. Multiphase Flow 27 (7), 1199 -- 1226

  12. [20]

    Gaskell, P. H. & Lau, A. K C 1988 Curvature compensated convective transport: SMART, A new boundedness preserving transport algorithm . Int. J. Numer. Methods Fluids 8 (6), 617--641

  13. [21]

    , Malkiel, E

    Gopalan , B. , Malkiel, E. & Katz, J. 2008 Experimental investigation of turbulent diffusion of slightly buoyant droplets in locally isotropic turbulence . Phys. Fluids 20 (9), 095102

  14. [22]

    & Herrmann, M

    Gorokhovski, M. & Herrmann, M. 2008 Modeling primary atomization . Annu. Rev. Fluid Mech. 40 (1), 343--366

  15. [23]

    2010 Detailed numerical simulations of the primary atomization of a turbulent liquid jet in crossflow

    Herrmann, M. 2010 Detailed numerical simulations of the primary atomization of a turbulent liquid jet in crossflow . J. Eng. Gas Turbines Power 132 (6)

  16. [24]

    2013 A sub-grid surface dynamics model for sub-filter surface tension induced interface dynamics

    Herrmann, M. 2013 A sub-grid surface dynamics model for sub-filter surface tension induced interface dynamics . Computers and Fluids 87 , 92--101

  17. [25]

    Hinze, J. O. 1955 Fundamentals of the hydrodynamic mechanism of splitting in dispersion processes . AlChE J. 1 (3), 289--295

  18. [26]

    & Katz, S

    Hulburt, H.M. & Katz, S. 1964 Some problems in particle technology . Chem. Eng. Sci. 19 (8), 555--574

  19. [27]

    Hussein, H. J. , Capp, S. P. & George, W. K. 1994 Velocity measurements in a high-Reynolds number, momentum-conserving, axisymmetric, turbulent jet . J. Fluid Mech. 258 (1994), 31--75

  20. [28]

    Igel, A. L. & van den Heever, S. C. 2016 The Importance of the Shape of Cloud Droplet Size Distributions in Shallow Cumulus Clouds. Part I: Bin Microphysics Simulations . J. Atmos. Sci. pp. JAS--D--15--0382.1

  21. [29]

    Jakobsen, H. A. 2014 The Population Balance Equation\/ , pp. 937--1003 . Cham: Springer International Publishing

  22. [30]

    , Brandvik, P

    Johansen, . , Brandvik, P. J. & Farooq, U. 2013 Droplet breakup in subsea oil releases - Part 2: Predictions of droplet size distributions with and without injection of chemical dispersants . Mar. Pollut. Bull. 73 (1), 327--335

  23. [31]

    & Meneveau, C

    Johnson, P.L. & Meneveau, C. 2018 Predicting viscous-range velocity gradient dynamics in large-eddy simulations of turbulence . Journal of Fluid Mechanics 837 , 80--114

  24. [32]

    Kolmogorov, A N 1949 On the breakage of drops in a turbulent flow. Dokl. Akad. Nauk SSSR 66 , 825--828

  25. [33]

    , Aoki, M

    Konno, M. , Aoki, M. & Saito, S. 1982 Scale effect on breakup process in liquid-liquid agitated tanks . J. Chem. Eng. Jpn. 16 , 312--319

  26. [34]

    , Peglow, M

    Kumar, J. , Peglow, M. , Warnecke, G. , Heinrich, S. & M \" o rl, L. 2006 Improved accuracy and convergence of discretized population balance for aggregation: The cell average technique . Chem. Eng. Sci. 61 (10), 3327--3342

  27. [35]

    & Ramakrishna, D

    Kumar, S. & Ramakrishna, D. 1996 On the solution of population balance equations by discretization -I. A fixed pivot technique . Chem. Eng. Sci. 51 (8), 1311--1332

  28. [36]

    Lasheras, J. C. , Eastwood, C. , Mart \' i nez-Baz \' a n, C. & Montaes, J. L. 2002 A review of statistical models for the break-up an immiscible fluid immersed into a fully developed turbulent flow . Int. J. Multiphase Flow 28 (2), 247--278

  29. [37]

    & Mewes, D

    Lehr, F. & Mewes, D. 1999 A transport equation for the interfacial area density applied to bubble columns . Chem. Eng. Sci. 56 , 1159--1166

  30. [38]

    , Miller, J

    Li, C. , Miller, J. , Wang, J. , Koley, S. S. & Katz, J. 2017 Size distribution and dispersion of droplets generated by impingement of breaking waves on oil slicks . J. Geophys. Res. Oceans 122 (10), 7938--7957

  31. [39]

    & Lucas, D

    Liao, Y. & Lucas, D. 2009 A literature review of theoretical models for drop and bubble breakup in turbulent dispersions . Chem. Eng. Sci. 64 (15), 3389--3406

  32. [40]

    Liu, Y. G. , You, L. G. , Yang, W. N. & Liu, F. 1995 On the Size Distribution of Cloud Droplets . Atmos. Res. 35 (2-4), 201--216

  33. [41]

    & Svendsen, H

    Luo, H. & Svendsen, H. F. 1996 Theoretical model for drop and bubble breakup in turbulent dispersions . AICHE J 42 (5), 1225--1233

  34. [42]

    , David, R

    Marchal, P. , David, R. , Klein, J. P. & Villermaux, J. 1988 Crystallization and precipitation engineering-I. An efficient method for solving population balance in crystallization with agglomeration . Chem. Eng. Sci. 43 (1), 59--67

  35. [43]

    , Deane, G

    Mart \' i nez-Baz \' a n , Rodr \' i guez-rodr \' i guez, J. , Deane, G. B. , Monta \ n \' e s, J. L. & Lasheras, J. C. 2010 Considerations on bubble fragmentation models . J. Fluid Mech. 661 , 159–177

  36. [44]

    , Monta \ n \' e s, J

    Mart \' i nez-Baz \' a n, C. , Monta \ n \' e s, J. L. & Lasheras, J. C. 1999 a\/ On the breakup of an air bubble injected into a fully developed turbulent flow. part 1. breakup frequency . J. Fluid Mech. 401 , 157--182

  37. [45]

    , Monta \ n \' e s, J

    Mart \' i nez-Baz \' a n, C. , Monta \ n \' e s, J. L. & Lasheras, J. C. 1999 b\/ On the breakup of an air bubble injected into a fully developed turbulent flow. part 2. size pdf of the resulting daughter bubbles . J. Fluid Mech. 401 , 183–207

  38. [46]

    , Xue, X

    Murphy, D W. , Xue, X. , Sampath, K. & Katz, J. 2016 Crude oil jets in crossflow: Effects of dispersant concentration on plume behavior . J. Geophys. Res. Oceans 121 (6), 4264--4281

  39. [47]

    , Gupta, J

    Narsimhan, G. , Gupta, J. P. & Ramkrishna, D. 1979 A model for transitional breakage probability of droplets in agitated lean liquid-liquid dispersions . Chem. Eng. Sci. 34 (2), 257--265

  40. [48]

    , Ramkrishna, D

    Narsimhan, G. , Ramkrishna, D. & Gupta, J. P. 1980 Analysis of drop size distributions in lean liquid-liquid dispersions . AlChE J. 26 (6), 991--1000

  41. [49]

    , Kronenburg, A

    Neuber, G. , Kronenburg, A. , Stein, O. T. & Cleary, M. J. 2017 MMC-LES modelling of droplet nucleation and growth in turbulent jets . Chem. Eng. Sci. 167 , 204--218

  42. [50]

    Nissanka, I. D. & Yapa, P. D. 2016 Calculation of oil droplet size distribution in an underwater oil well blowout . Journal of Hydraulic Research 54 (3), 307--320

  43. [51]

    Nissanka, I. D. & Yapa, P. D. 2018 Calculation of oil droplet size distribution in ocean oil spills: A review . Mar. Pollut. Bull. 135 , 723 -- 734

  44. [52]

    North, E. W. , Adams, E. E. , Thessen, A. E. , Schlag, Z. , He, R. , Socolofsky, S. A. , Masutani, S. M. & Peckham, S. D. 2015 The influence of droplet size and biodegradation on the transport of subsurface oil droplets during the Deepwater Horizon spill: A model sensitivity s...

  45. [53]

    , Thornock, J

    Pedel, J. , Thornock, J. N. , Smith, S. T. & Smith, P. J. 2014 Large eddy simulation of polydisperse particles in turbulent coaxial jets using the direct quadrature method of moments . Int. J. Multiphase Flow 63 , 23--38

  46. [54]

    Cambridge: Cambridge Univ

    Pope, S B 2011 Turbulent flows \/ . Cambridge: Cambridge Univ. Press

  47. [55]

    , Meneveau, C

    Port \' e -Agel, F. , Meneveau, C. & Parlange, M. B. 2000 A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer . J. Fluid Mech. 415 , 261--284

  48. [56]

    Prince, M. J. & Blanch, H. W. 1990 Bubble Coalescence and Break-up in Air-Sparged Bubble Columns . AlChE J. 36 (10), 1485--1499

  49. [57]

    1985 The status of population balances

    Ramkrishna, D. 1985 The status of population balances . Rev. Chem. Eng. 3 (1), 49--95

  50. [58]

    Randolph, A. D. 1964 A population balance for countable entities . Can. J. Chem. Eng. 42 (6), 280--281

  51. [59]

    , Cleary, M

    Salehi, F. , Cleary, M. J. & Masri, A. R. 2017 Population balance equation for turbulent polydispersed inertial droplets and particles . J. Fluid Mech. 831 , 719--742

  52. [60]

    , Ramkrishna, D

    Sathyagal, A.N. , Ramkrishna, D. & Narsimhan, G. 1996 Droplet breakage in stirred dispersions. breakage functions from experimental drop-size distributions . Chem. Eng. Sci. 51 (9), 1377 -- 1391

  53. [61]

    & Yamamoto, K

    Sato, Y. & Yamamoto, K. 1987 Lagrangian measurement of fluid-particle motion in an isotropic turbulent field . J. Fluid Mech. 175 , 183–199

  54. [62]

    ICLASS 2012, 12th Triennial Int

    Seubert, N , Kronenburg, A , Stein, O T , Ge, Y & Cleary, M J 2012 Large Eddy Simulation-Probability Density Function modelling of nucleation and condensation of DBP droplets in a turbulent jet . ICLASS 2012, 12th Triennial Int. Conference on Liquid Atomization and Spray Syste...

  55. [63]

    & Rigopoulos, S

    Sewerin, F. & Rigopoulos, S. 2017 An LES-PBE-PDF approach for modeling particle formation in turbulent reacting flows . Phys. Fluids 29 (10)

  56. [64]

    , Sollum, E

    Skartlien, R. , Sollum, E. & Schumann, H. 2013 Droplet size distributions in turbulent emulsions: Breakup criteria and surfactant effects from direct numerical simulations . J. Chem. Phys. 139 (17), 174901

  57. [65]

    1963 General Circulation Experiments With the Primitive Equations

    Smagorinsky, J. 1963 General Circulation Experiments With the Primitive Equations . Monthly Weather Review 91 (3), 99--164

  58. [66]

    Smoluchowski , M. V. 1916 Drei Vortrage uber Diffusion, Brownsche Bewegung und Koagulation von Kolloidteilchen . Zeitschrift fur Physik 17 , 557--585

  59. [67]

    & Jakobsen, H

    Solsvik, J. & Jakobsen, H. A. 2016 Development of Fluid Particle Breakup and Coalescence Closure Models for the Complete Energy Spectrum of Isotropic Turbulence . Ind. Eng. Chem. Res. 55 (5), 1449--1460

  60. [68]

    , Maa , S

    Solsvik, J. , Maa , S. & Jakobsen, H. A. 2016 Definition of the Single Drop Breakup Event . Ind. Eng. Chem. Res. 55 (10), 2872--2882

  61. [69]

    Tseng, Y. H. , Meneveau, C. & Parlange, M. B. 2006 Modeling flow around bluff bodies and predicting urban dispersion using large eddy simulation . Environ. Sci. Technol. 40 (8), 2653--2662

  62. [70]

    & Tavlarides, L

    Tsouris, C. & Tavlarides, L. L. 1994 Breakage and coalescence models for drops in turbulent dispersions . AlChE J. 40 (3), 395--406

  63. [71]

    & Wang, J

    Wang, T. & Wang, J. 2007 Numerical simulations of gas-liquid mass transfer in bubble columns with a CFD-PBM coupled model . Chem. Eng. Sci. 62 (24), 7107--7118

  64. [72]

    & Antonia, R

    Xu, G. & Antonia, R. 2002 Effect of different initial conditions on a turbulent round free jet . Exp. Fluids 33 (5), 677--683

  65. [73]

    , Chen, B

    Yang, D. , Chen, B. , Chamecki, M. & Meneveau, C. 2015 Oil plumes and dispersion in langmuir, upper-ocean turbulence: Large-eddy simulations and k-profile parameterization . J. Geophys. Res. Oceans 120 (7), 4729--4759

  66. [74]

    , Chen, B

    Yang, D. , Chen, B. , Socolofsky, S. A. , Chamecki, M. & Meneveau, C. 2016 Large-eddy simulation and parameterization of buoyant plume dynamics in stratified flow . J. Fluid Mech. 794 , 798--833

  67. [75]

    , Meneveau, C

    Yang, D. , Meneveau, C. & Shen, L. 2014 Effect of downwind swells on offshore wind energy harvesting – A large-eddy simulation study . Renewable Energy 70 , 11--23

  68. [76]

    , Boufadel, M C

    Zhao, L. , Boufadel, M C. , Socolofsky, S A. , Adams, E. , King, T. & Lee, K. 2014 a\/ Evolution of droplets in subsea oil and gas blowouts: Development and validation of the numerical model VDROP-J . Mar. Pollut. Bull. 83 (1), 58 -- 69

  69. [77]

    , Shaffer, F

    Zhao, L. , Shaffer, F. , Robinson, B. , King, T. , D’Ambrose, C. , Pan, Z. , Gao, F. , Miller, R. S. , Conmy, R. N. & Boufadel, M. C. 2016 Underwater oil jet: Hydrodynamics and droplet size distribution . Chem. Eng. J. 299 , 292 -- 303

  70. [78]

    , Torlapati, J

    Zhao, L. , Torlapati, J. , Boufadel, M C. , King, T. , Robinson, B. & Lee, K. 2014 b\/ VDROP : A comprehensive model for droplet formation of oils and gases in liquids - incorporation of the interfacial tension and droplet viscosity . Chem. Eng. J. 253 , 93 -- 106

  71. [79]

    @esa (Ref

    \@ifclassloaded aguplus natbib The aguplus class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later remove the command natbib from the document \@ifclassloaded nlinproc natbib The nlinproc class already includes natbib cod...

  72. [80]

    @stdbsttrue NAT@ctr \@lbibitem[ NAT@ctr ] \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 [ @natanchorstart #2\@extra@b@citeb \@biblabel @num @natanchorend] @ifc...

  73. [81]

    A fast Eulerian method for disperse two-phase flow

    @open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifundefined NAT@sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifundefined bib@heading @heading NAT@ctr thebibliography [1] @ \@biblabel NAT@ctr \@bib...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.