Pith. sign in

REVIEW 4 major objections 6 minor 1 cited by

Thread-safe multiphase lattice Boltzmann model for droplet and bubble dynamics at high density and viscosity contrasts

T0 review · 4 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A thread-safe, third-order lattice Boltzmann method reproduces droplet and bubble dynamics at density ratios up to 1000 and viscosity ratios up to 100.

desk verdict A genuinely useful extension of thread-safe LBM to high density contrasts, but the printed Laplacian stencil in Eq. 19 does not annihilate constants, so the method as written is not reproducible until that is fixed. read the letter →

arxiv 2501.00846 v1 pith:SYBPPJCP submitted 2025-01-01 physics.flu-dyn physics.comp-ph

Marco Lauricella , Adriano Tiribocchi , Sauro Succi , Luca Brandt , Aritra Mukherjee , Michele La Rocca , Andrea Montessori This is my paper · ORCID
classification physics.flu-dynphysics.comp-ph MSC 76M2876T10 PACS 47.11.-j47.55.D
keywords latticeBoltzmannmethodthread-safeGPUcomputingAllen-Cahnphasefieldhighdensitycontrastmultiphaseflowrisingbubbledropletcollisionregularizedair-watersimulation
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 claims that a thread-safe lattice Boltzmann implementation, which reconstructs distribution populations from macroscopic fields to avoid race conditions, can be extended to two-fluid flows with density ratios up to 1000 and viscosity ratios up to 100 by adding an Allen-Cahn interface tracker and three extra force terms. The resulting solver reproduces the rising-bubble deformation regimes on the Eötvös–Galilei map for an air-water-like contrast and matches experimental sequences of head-on and off-axis water-droplet collisions, including satellite formation. It also reports weak-scaling efficiency near 0.95 up to 32 GPUs, so high-contrast multiphase simulations become practical on shared-memory and GPU architectures. A sympathetic reader would care because this pushes diffuse-interface lattice Boltzmann methods into the parameter range of realistic environmental, microfluidic, and biological flows.

What carries the argument

The central object is the fused streaming-collision update $f_i(\mathbf{x}+\mathbf{c}_i\Delta t,t+\Delta t)=f_i^{eq}(\mathbf{x},t)+(1-\omega)f_i^{neq}(\mathbf{x},t)+S_i(\mathbf{x},t)$, in which both the equilibrium and non-equilibrium parts are reconstructed from macroscopic fields at each node rather than read from memory, using the recursive properties of Hermite polynomials up to third order on a D3Q27 (27-velocity cubic) stencil. This reconstruction is what removes race conditions, letting streaming and collision be fused into one kernel. The density and viscosity contrast is then injected through the three force terms, while the conservative Allen-Cahn equation $\partial_t\phi+u_\alpha\partial_\alpha\phi=D\partial_\alpha\partial_\alpha\phi-\kappa\partial_\alpha(\phi(1-\phi)n_\alpha)$ is advanced by a forward-time centred-space finite-difference scheme to track the interface.

What would settle it

A static planar interface with $\rho_L/\rho_G=1000$ and no body force should remain perfectly quiescent: if the spurious velocity produced by the unbalanced discrete pressure and viscous force terms does not vanish, or does not decrease as the interface width is refined, the high-density-ratio capability claim is false. The same test with $\nu_L/\nu_G=100$ isolates the viscous correction.

Watch

Extended reading notes

Core claim

The paper's central claim is that the thread-safe lattice Boltzmann update—rebuilding post-collision populations from macroscopic moments instead of reading neighbouring populations—can carry a high-contrast multiphase model if the non-equilibrium part is reconstructed to third order in Hermite polynomials and the density/viscosity contrast is injected through three explicit forces: surface tension $F_s=\mu_\phi\partial_\alpha\phi$, pressure correction $F_p=-p^*c_s^2\partial_\alpha\rho$, and viscous correction $F_\nu=-\frac{\nu\omega}{c_s^2\Delta t}[\sum_i(f_i-f_i^{eq})c_{i\alpha}c_{i\beta}]\partial_\alpha\rho$. With this construction, the paper demonstrates that $\rho_L/\rho_G=1000$ and $\nu_L/\nu_G=100$ are reachable on a 27-velocity cubic lattice, recovering the rising-bubble regime map and reproducing experimental water-droplet collision outcomes at $\mathit{Oh}=0.0044$.

Load-bearing premise

The method's high-density-ratio reach rests on the assumption that the pressure-gradient and viscous force corrections (Eqs. 16–17) exactly recover the target momentum equation at the discrete interface, using the stencils of Eqs. 18–19; if those forces are not balanced at the sharpest interface cells, the claim of density ratio 1000 collapses even though the thread-safe streaming itself is sound.

Editorial extensions

If this is right

  • At density ratio 1000 and viscosity ratio 100, the method reproduces the air-water rising-bubble deformation regimes on the Eötvös–Galilei map, including axisymmetric, hat-shaped, peripheral-breakup, central-breakup, and oscillatory regimes, matching the benchmark reference when the interface is resolved.
  • Head-on and off-axis water-droplet collisions with Weber numbers from 10 to 96 and Reynolds numbers from roughly 700 to 2200 reproduce the experimental coalescence and separation boundaries and the formation of satellite droplets at Ohnesorge number 0.0044.
  • Because equilibrium and non-equilibrium populations are reconstructed from macroscopic fields, the fused stream-and-collision kernel is free of read-after-write race conditions, which is what permits weak-scaling efficiency near 0.95 up to 32 GPUs.
  • A train of five raindrops impacting a solid substrate reproduces the observed lamella spreading, recoil, and growth of the liquid film with successive impacts, giving a direct numerical route to studying microplastic and pollutant dispersion by rain.
  • Third-order recursive regularization raises the stability ceiling of the thread-safe scheme, so a droplet diameter of 100 lattice units with a Cahn number of 0.04 is enough to capture reflexive separation and satellite droplets.

Reading between the lines

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

  • The high-density-ratio capability ultimately depends on the balance of the pressure and viscous force corrections at the sharpest interface cells; a direct test the paper does not report is to compute the discrete momentum balance across a static interface at density ratio 1000 before running any dynamic simulation.
  • The paper's assertion that machine learning is non-essential for complex multiphase flow is a broader opinion than the benchmarks prove; the measurements establish scalability and fidelity, not a comparison with learned surrogate models.
  • The raindrop-train simulation points toward a testable extension: quantifying how impact spacing, surface wettability, and roughness control the breakup statistics and the ejection of microplastics from the lamella rim.
  • The paper flags dual-grid resolution of the phase field as future work; if implemented, the interface width could be controlled independently of flow resolution, which would further relax the resolution constraint seen in the thin-film breakup case.
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

4 major / 6 minor

Summary. The manuscript presents a hybrid lattice Boltzmann method for two-component flows with high density and viscosity contrasts, combining a thread-safe, high-order regularized collision-streaming scheme with a conservative Allen-Cahn interface-tracking equation solved by finite differences. The surface tension, pressure-gradient, and viscous forces are added as external forcings, with the pressure and viscous contributions intended to extend the method to density and viscosity ratios up to 1000 and 100. Validation is carried out against reference simulations for rising bubbles and against experiments for head-on and off-axis droplet collisions, and GPU strong/weak scaling is reported on the Leonardo supercomputer. An application to successive raindrop impacts on a solid surface is also presented.

Significance. If the method performs as claimed, it would be a useful contribution to GPU-oriented multiphase LBM, since the thread-safe formulation removes race conditions while retaining third-order Hermite reconstruction. The paper has concrete strengths: it benchmarks against external experiments and independent simulations rather than only self-comparisons, it demonstrates the challenging air-water density ratio, and it reports scaling data on a modern GPU cluster. However, the manuscript as written contains an inconsistency in the discrete Laplacian operator used in several load-bearing equations, and the validation is largely qualitative, with no grid-convergence study and at least one visible discrepancy with the cited experiment. Because the central claim depends on these points, the paper cannot be accepted in its present form.

major comments (4)
  1. [II.C, Eq. (19)] The discrete second-derivative operator printed in Eq. (19) does not annihilate constant fields. For a constant field Ψ=C, the right-hand side evaluates to (1-2w0)C/c_s^2, which is (11/27)C/c_s^2 for the D3Q27 weights (w0=8/27), rather than zero. Since this operator is used for the chemical potential in Eq. (3) and for the Allen-Cahn diffusion term in Eq. (20), the method as written contains a spurious constant-field source in every time step, independent of density ratio. The standard isotropic Laplacian has the form (2/c_s^2)(Σ_i w_i Ψ(x+c_i)-Ψ(x)) or an equivalent expression that annihilates constants. The authors should either correct Eq. (19) or explicitly define the operator actually implemented, and they should show that this operator passes the constant-annihilation test. This is load-bearing because the high-density-contrast claim relies on the discrete forces being balanced at the interface.
  2. [III.A, Fig. 4 and surrounding text] The validation of the rising-bubble case is qualitative: the center-of-mass velocity comparison with Adelsberger et al. is presented graphically with no error metrics, and the statement that 'good agreement' holds 'provided that the interface width is adequately resolved' is not supported by a grid-convergence or interface-width-convergence study. Given that the main novelty is the ability to handle density ratio 1000 and viscosity ratio 100, the absence of a resolution study is a significant gap: it leaves open whether the observed agreement arises from correct force balance or from sufficiently diffuse interfaces masking errors.
  3. [III.B, We=25 case] The paper acknowledges that the simulation produces a satellite droplet after thread breakup at We=25, while the cited experiments of Ashgriz and Poo do not report one. The explanation invoking Huang et al. is plausible, but it is not accompanied by a quantitative comparison (satellite size, pinch-off time, or a resolution study showing that this satellite is converged). In addition, Fig. 7 attributes a thin-film breakup at We=40 to insufficient resolution. Together these points undermine the abstract's claim of results 'in agreement with experiments' and need to be addressed with either a systematic resolution/parameter study or a more nuanced statement of agreement.
  4. [II.C, Eqs. (15)-(17)] The force-splitting construction assumes that the discrete pressure-gradient force F_p and viscous force F_ν, together with the surface tension force, reproduce the target momentum equation at density ratio 1000. No test is reported for spurious currents at a stationary droplet or for the discrete balance of these forces at a sharp interface. Such a test is directly relevant because the discrete stencils of Eqs. (18)-(19) determine whether the pressure-gradient and viscous forces are accurately represented. The authors should add a quantitative assessment, e.g., maximum spurious velocity as a function of grid resolution and interface width for a stationary droplet at high density ratio.
minor comments (6)
  1. [II.C, Eq. (19)] The left-hand side of Eq. (19) has two indices (∂α∂βΨ) while the right-hand side is a scalar; please clarify whether this equation is intended as the Laplacian or as the full Hessian, and supply the correct lattice weights and central term.
  2. [III.B, Fig. 5 caption] The caption of Fig. 5 states 'We=23' whereas the text discussing this case refers to We=25; the inconsistency should be reconciled.
  3. [II.C, Eq. (3)] The text spells the name 'Jaqmin' but the reference [22] is by Jacqmin; please correct the spelling.
  4. [II.E, Performance section] The paragraph claiming that machine learning is 'non-essential' for studying complex flows is an opinion not supported by the scaling benchmarks or by any comparison with ML-based methods; it should be removed or supported with quantitative evidence.
  5. [III.B, Fig. 13] The collision map reports only the numerical simulation outcomes as triangles; overlaying the experimental regime boundaries from [34] would make the claimed agreement quantitatively assessable.
  6. [II.A, Eq. (5)] The equilibrium profile in Eq. (5) uses absolute values of coordinate differences, which is not the standard planar interface profile and is introduced without derivation; please clarify this choice and its relation to the coordinate axes.

Circularity Check

0 steps flagged · score 2.0 of 10

No material circularity: the thread-safe and regularized LBM components are prior-work lineage, while the high-density-contrast capability is benchmarked against external experiments and independent simulations rather than derived from fitted inputs.

full rationale

The paper's central claim is that a hybrid Allen-Cahn / high-order thread-safe LBM reproduces droplet and bubble dynamics at density and viscosity contrasts up to 1000 and 100. That claim is tested against two external references: the rising-bubble center-of-mass velocities of Adelsberger et al. (ref. 33) and the experimental droplet-collision morphologies of Ashgriz and Poo (ref. 34). The dimensionless parameters (Eo, Ga, We, Re, Oh) are set from the reference cases, and the agreement is an outcome of the simulation rather than a fitted quantity. The force terms in Sec. II.C are derived from the target momentum equation: F_p removes the p* c_s^2 grad rho contribution of the pressure gradient, and F_nu is obtained from the existing relation between the second-order moment of the discrete distribution and the deviatoric stress tensor (ref. 30), so the momentum balance is not imported as an ansatz. The self-citations (refs. 18-21 and the accLB code) supply the thread-safe collision, Hermite regularization, and code lineage; they do not by themselves establish the high-density-contrast result, which is why the self-citation density does not make the derivation circular. Interface width and lattice resolution are numerical choices, not calibration constants. One independent concern: the printed discrete Laplacian in Eq. 19 does not annihilate constant fields on D3Q27, since for constant Psi it returns (1-2w0)Psi/c_s^2 = (11/27)Psi/c_s^2 with w0=8/27, which would affect the chemical potential and Allen-Cahn diffusion terms as written. This is a correctness and reproducibility issue rather than a circularity issue, and it does not raise the circularity score.

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

The central claim rests on standard LBM multiscale analysis, a third-order Hermite closure, a phase-field model of interfacial forces, and discrete stencils for gradients. The only hand-chosen numerical parameter that affects all benchmarks is the interface width delta, with the Allen-Cahn mobility left unspecified. No new physical entities are introduced.

free parameters (2)
  • Interface width delta = 4 lattice units in droplet benchmarks; not reported for bubble runs
    A numerical parameter chosen by hand that controls the smearing of density and viscosity across the interface; it directly affects all benchmark results.
  • Allen-Cahn diffusivity/mobility D
    Appears in Eqs.4 and 20 and controls interface relaxation, but no value or calibration procedure is reported in the paper.
assumptions (5)
  • standard math The discrete LBM with the chosen equilibrium and forcing recovers the incompressible Navier-Stokes equations via Chapman-Enskog analysis.
    Invoked implicitly in Section II.B; no multiscale expansion is shown in this paper.
  • standard math Third-order non-equilibrium moments are reconstructed from the second-order moment and velocity via the recursive Hermite relation a3_neq = u a2_neq + ...
    Taken from Malaspinas and used in Eq.10; this closure is load-bearing for the high-order TSLB claim.
  • domain assumption The conservative Allen-Cahn equation with the tanh equilibrium profile captures two-phase interface dynamics, and the chosen chemical potential yields the surface tension force.
    The phase-field model is assumed to reproduce sharp-interface capillary physics at density ratio 1000; no sharp-interface limit is derived.
  • standard math The discrete gradient and Laplacian stencils of Eqs.18-19 are sufficiently isotropic for the surface tension, pressure, and viscous forces at high density contrast.
    Borrowed from Thampi et al.; accuracy on coarse interface resolutions is not analyzed.
  • standard math The viscous force in Eq.17 can be expressed through the non-equilibrium second-order moment as in Kruger et al.
    Used without derivation in Section II.C.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thread-safe multiphase lattice Boltzmann model for droplet and bubble dynamics at high density and viscosity contrasts." pith.science (2026). https://pith.science/paper/SYBPPJCP

@misc{pith2026250100846,
  author       = {Pith},
  title        = {Pith review of: Thread-safe multiphase lattice Boltzmann model for droplet and bubble dynamics at high density and viscosity contrasts},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SYBPPJCP}},
  note         = {Machine review of arXiv:2501.00846}
}
read the original abstract

This study presents a high-order, thread-safe version of the lattice Boltzmann (LBM) method, incorporating an interface-capturing equation, based on the conservative Allen-Cahn equation, to simulate incompressible two-component systems with high-density and viscosity contrasts. The method utilizes a recently proposed thread-safe implementation optimized for shared memory architectures and it is employed to reproduce the dynamics of droplets and bubbles on several test cases with results in agreement with experiments and other numerical simulations from the literature. The proposed approach offers promising opportunities for high-performance computing simulations of realistic fluid systems with high-density and viscosity contrasts for advanced applications in environmental, atmospheric and meteorological flows, all the way down to microfluidic and biological systems, particularly on GPU-based architectures.

Figures

Figures reproduced from arXiv: 2501.00846 by the authors.

Figure 1
Figure 1. FIG. 1. Plot of the strong scalings for cubic boxes of sides 512(panels (a) to (c) blue lines with plusses) and 1024 (panels (a) to [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Sketch of the computational domain for the rising bubble test. [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Deformation regimes of a rising bubble as a function of [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (11 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Center-of-mass velocity ( [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Reflexive head-on separation at [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Reflexive head-on separation with satellite droplet formation at [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Top view of head-on collision at [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Reflexive head-on separation with formation of multiple satellite droplets at [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]
Figure 9
Figure 9. Figure 9: and 10 illustrate two representative cases at relatively low Weber numbers, corresponding to x = 0.5 and W e = 10, and x = 0.28 and W e = 53, respectively. x is impact parameter, i.e. the ratio of the distance between the centers of the droplets and their diameter. In …
Figure 10
Figure 10. Figure 10: FIG. 10. Coalescence collision at [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11. Collision with satellite droplet formation of at [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12. Collision with satellite droplet formation of at [PITH_FULL_IMAGE:figures/full_fig_p013_12.png]
Figure 13
Figure 13. Figure 13: FIG. 13 [PITH_FULL_IMAGE:figures/full_fig_p014_13.png]
Figure 14
Figure 14. Figure 14: FIG. 14. (a-b) sequence of multiple impacts between falling droplets and a neutral solid substrate at [PITH_FULL_IMAGE:figures/full_fig_p015_14.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. accLB: A High-Performance Lattice Boltzmann Code for Multiphase Turbulence on Multi-Gpu Architectures

    physics.flu-dyn 2025-05 conditional novelty 5.0 of 10

    accLB is a multi-GPU lattice Boltzmann code that reaches more than 150 GLUPS and reproduces single-phase and bubble-laden HIT energy spectra.

Reference graph

Works this paper leans on

39 extracted references · 33 canonical work pages · cited by 1 Pith paper

  1. [1]

    Allen-Cahn equation for interface tracking The interface tracking is performed by capturing the dynamic evolution of a phase field , i.e. ϕ, via the following conservative Allen-Cahn equation: ∂tϕ + uα∂αϕ = D∂α∂αϕ − κ∂α(ϕ(1 − ϕ)nα) (4) 3 where D is the diffusivity of the interface and κ = 4D/δ (δ being the interface width). In this model, ϕ can assume any...

    work page

  2. [2]

    M. W¨ orner, Numerical modeling of multiphase flows in microfluidics and micro process engineering: a review of methods and applications, Microfluidics and nanofluidics 12 (6) (2012) 841–886

    work page 2012

  3. [3]

    Q. T. Ho, J. Carmeliet, A. K. Datta, T. Defraeye, M. A. Delele, E. Herremans, L. Opara, H. Ramon, E. Tijskens, R. Van Der Sman, et al., Multiscale modeling in food engineering, Journal of food Engineering 114 (3) (2013) 279–291

    work page 2013

  4. [4]

    Parker, Multiphase flow and transport in porous media, Reviews of Geophysics 27 (3) (1989) 311–328

    J. Parker, Multiphase flow and transport in porous media, Reviews of Geophysics 27 (3) (1989) 311–328

    work page 1989

  5. [5]

    W. W. Grabowski, L.-P. Wang, Growth of cloud droplets in a turbulent environment, Annual review of fluid mechanics 45 (1) (2013) 293–324

    work page 2013

  6. [6]

    Bogdan, A

    M. Bogdan, A. Montessori, A. Tiribocchi, F. Bonaccorso, M. Lauricella, L. Jurkiewicz, S. Succi, J. Guzowski, Stochastic jetting and dripping in confined soft granular flows, Physical Review Letters 128 (12) (2022) 128001

    work page 2022

  7. [7]

    M. T. Guo, A. Rotem, J. A. Heyman, D. A. Weitz, Droplet microfluidics for high-throughput biological assays, Lab on a Chip 12 (12) (2012) 2146–2155

    work page 2012

  8. [8]

    Montessori, M

    A. Montessori, M. Lauricella, N. Tirelli, S. Succi, Mesoscale modelling of near-contact interactions for complex flowing interfaces, Journal of Fluid Mechanics 872 (2019) 327–347

    work page 2019

Show all 39 references
  1. [9]

    E. E. Michaelides, M. Sommerfeld, B. van Wachem, Multiphase flows with droplets and particles, CRC Press, 2022

    2022

  2. [10]

    C. R. Constante-Amores, L. Kahouadji, A. Batchvarov, S. Shin, J. Chergui, D. Juric, O. K. Matar, Direct numerical simulations of transient turbulent jets: vortex-interface interactions, Journal of Fluid Mechanics 922 (2021) A6. doi: https://doi.org/10.1017/jfm.2021.519

    2021 doi

  3. [11]

    Montessori, M

    A. Montessori, M. Lauricella, E. Stolovicki, D. A. Weitz, S. Succi, Jetting to dripping transition: Critical aspect ratio in step emulsifiers, Physics of Fluids 31 (2) (2019). doi:https://doi.org/10.1063/1.5084797

    2019 doi

  4. [12]

    Montessori, M

    A. Montessori, M. Lauricella, S. Succi, E. Stolovicki, D. Weitz, Elucidating the mechanism of step emulsification, Physical Review Fluids 3 (7) (2018) 072202. doi:https://doi.org/10.1103/PhysRevFluids.3.072202

    2018 doi

  5. [13]

    L. Yi, C. Wang, S. G. Huisman, C. Sun, Recent developments of turbulent emulsions in taylor–couette flow, Philosophical Transactions of the Royal Society A 381 (2243) (2023) 20220129. doi:https://doi.org/10.1098/rsta.2022.0129

    2023

  6. [14]

    L. G. Martinez, B. Duret, J. Reveillon, F. X. Demoulin, A new dns formalism dedicated to turbulent two-phase flows with phase change, International Journal of Multiphase Flow 143 (2021) 103762. doi:https://doi.org/10.1016/j. ijmultiphaseflow.2021.103762

    2021

  7. [15]

    Succi, G

    S. Succi, G. Amati, M. Bernaschi, G. Falcucci, M. Lauricella, A. Montessori, Towards exascale lattice boltzmann computing, Computer & Fluids 181 (2019) 107–115. doi:https://doi.org/10.1016/j.compfluid.2019.01.005

    2019 doi

  8. [16]

    Succi, G

    S. Succi, G. Amati, F. Bonaccorso, M. Lauricella, M. Bernaschi, A. Montessori, A. Tiribocchi, Toward exascale design of soft mesoscale materials, J. Comput. Sci. 46 (2020) 101175. doi:https://doi.org/10.1016/j.jocs.2020.101175

    2020

  9. [17]

    Kr¨ uger, H

    T. Kr¨ uger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E. M. Viggen, The lattice boltzmann method, Springer International Publishing 10 (978-3) (2017) 4–15

    2017

  10. [18]

    Montessori, G

    A. Montessori, G. Falcucci, Lattice Boltzmann modeling of complex flows for engineering applications, Morgan & Claypool Publishers, 2018

    2018

  11. [19]

    Montessori, M

    A. Montessori, M. Lauricella, A. Tiribocchi, M. Durve, M. La Rocca, G. Amati, F. Bonaccorso, S. Succi, Thread-safe lattice boltzmann for high-performance computing on gpus, Journal of Computational Science 74 (2023) 102165

    2023

  12. [20]

    Montessori, M

    A. Montessori, M. La Rocca, G. Amati, M. Lauricella, A. Tiribocchi, S. Succi, High-order thread-safe lattice Boltzmann model for high performance computing turbulent flow simulations, Physics of Fluids 36 (3) (2024) 035171. arXiv:https: //pubs.aip.org/aip/pof/article-pdf/doi/1...

    2024 doi

  13. [21]

    Montessori, L

    A. Montessori, L. A. Hegele Jr, M. Lauricella, A high-performance lattice boltzmann model for multicomponent turbulent jet simulations, arXiv preprint arXiv:2403.15773 (2024)

    2024 arXiv

  14. [22]

    Malaspinas, Increasing stability and accuracy of the lattice boltzmann scheme: recursivity and regularization, arXiv preprint arXiv:1505.06900 (2015)

    O. Malaspinas, Increasing stability and accuracy of the lattice boltzmann scheme: recursivity and regularization, arXiv preprint arXiv:1505.06900 (2015)

    2015 arXiv

  15. [23]

    Jacqmin, Calculation of two-phase navier–stokes flows using phase-field modeling, Journal of computational physics 155 (1) (1999) 96–127

    D. Jacqmin, Calculation of two-phase navier–stokes flows using phase-field modeling, Journal of computational physics 155 (1) (1999) 96–127

    1999

  16. [24]

    Kim, A continuous surface tension force formulation for diffuse-interface models, Journal of computational physics 204 (2) (2005) 784–804

    J. Kim, A continuous surface tension force formulation for diffuse-interface models, Journal of computational physics 204 (2) (2005) 784–804

    2005

  17. [25]

    Schenk, G

    M. Schenk, G. Giamagas, A. Roccon, A. Soldati, F. Zonta, Computationally efficient and interface accurate dual-grid phase-field simulation of turbulent drop-laden flows, Journal of Fluids Engineering 146 (12) (2024)

    2024

  18. [26]

    Ding, C.-j

    H. Ding, C.-j. Yuan, On the diffuse interface method using a dual-resolution cartesian grid, Journal of Computational Physics 273 (2014) 243–254

    2014

  19. [27]

    Grad, On the kinetic theory of rarefied gases, Communications on pure and applied mathematics 2 (4) (1949) 331–407

    H. Grad, On the kinetic theory of rarefied gases, Communications on pure and applied mathematics 2 (4) (1949) 331–407

    1949

  20. [28]

    Coreixas, G

    C. Coreixas, G. Wissocq, G. Puigt, J.-F. Boussuge, P. Sagaut, Recursive regularization step for high-order lattice boltzmann methods, Physical Review E 96 (3) (2017) 033306. 17

    2017

  21. [29]

    Grad, Note on n-dimensional hermite polynomials, Communications on Pure and Applied Mathematics 2 (4) (1949) 325–330

    H. Grad, Note on n-dimensional hermite polynomials, Communications on Pure and Applied Mathematics 2 (4) (1949) 325–330

    1949

  22. [30]

    Z. Guo, C. Zheng, B. Shi, Discrete lattice effects on the forcing term in the lattice boltzmann method, Physical review E 65 (4) (2002) 046308

    2002

  23. [31]

    Kr¨ uger, F

    T. Kr¨ uger, F. Varnik, D. Raabe, Shear stress in lattice boltzmann simulations, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics 79 (4) (2009) 046704

    2009

  24. [32]

    S. P. Thampi, S. Ansumali, R. Adhikari, S. Succi, Isotropic discrete laplacian operators from lattice hydrodynamics, Journal of Computational Physics 234 (2013) 1–7. doi:https://doi.org/10.1016/j.jcp.2012.07.037. URL https://www.sciencedirect.com/science/article/pii/S0021999112004226

    2013 doi

  25. [33]

    M. K. Tripathi, K. C. Sahu, R. Govindarajan, Dynamics of an initially spherical bubble rising in quiescent liquid, Nature communications 6 (1) (2015) 6268

    2015

  26. [34]

    Adelsberger, P

    J. Adelsberger, P. Esser, M. Griebel, S. Groß, M. Klitz, A. R¨ uttgers, 3d incompressible two-phase flow benchmark com- putations for rising droplets (2014)

    2014

  27. [35]

    Ashgriz, J

    N. Ashgriz, J. Poo, Coalescence and separation in binary collisions of liquid drops, Journal of Fluid Mechanics 221 (1990) 183–204

    1990

  28. [36]

    The sharp-interface limit of the cahn–hilliard/navier–stokes model for binary fluids, Journal of Fluid Mechanics 714 (2013) 95–126

    2013

  29. [37]

    Huang, K.-L

    K.-L. Huang, K.-L. Pan, C. Josserand, Pinching dynamics and satellite droplet formation in symmetrical droplet collisions, Physical review letters 123 (23) (2019) 234502

    2019

  30. [38]

    Amani, N

    A. Amani, N. Balcazar, E. Guti´ errez, A. Oliva, Numerical study of binary droplets collision in the main collision regimes, Chemical Engineering Journal 370 (2019) 477–498

    2019

  31. [39]

    J. Li, H. Zhang, Q. Liu, Dynamics of a successive train of monodispersed millimetric-sized droplets impact on solid surfaces at low weber number, Experimental Thermal and Fluid Science 102 (2019) 81–93

    2019

Pith tools

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