Pith. sign in

REVIEW 3 major objections 6 minor 43 references

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

T0 review · 3 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

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

Pith's one-line read The accLB GPU lattice Boltzmann solver sustains over 150 GLUPS on 64 GPUs and reproduces the -3 energy spectrum of bubble-laden turbulence.

desk verdict A solid, incremental GPU lattice-Boltzmann code paper with real scaling data and plausible physics, but the headline -3 spectrum needs a grid-convergence check before I'd trust it. read the letter →

arxiv 2505.01126 v1 pith:7SEE2QUI submitted 2025-05-02 physics.flu-dyn physics.comp-ph

Marco Lauricella , Aritra Mukherjee , Luca Brandt , Sauro Succi , Daulet Izbassarov , Andrea Montessori This is my paper · ORCID
classification physics.flu-dynphysics.comp-ph
keywords turbulentflowslatticeBoltzmannmethodphase-fieldGPUcomputingOpenACChomogeneousisotropicturbulencebubble-ladenenergyspectrum
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

The paper presents accLB, a Fortran lattice Boltzmann solver built for multiphase turbulence on GPU clusters. It claims sustained throughput above 150 GLUPS (billions of lattice-update operations per second) on 64 GPUs, with efficient strong and weak scaling on two pre-exascale GPU systems. The physical core is validation: single-phase homogeneous isotropic turbulence reproduces the $-5/3$ energy and $-7/3$ pressure spectra, while bubble-laden runs at a 6% volume fraction produce a transition to a $-3$ energy spectrum, matching experiments and theory, along with stronger small-scale intermittency. If these results hold, accLB provides a portable platform for direct numerical simulation of bubble-induced turbulence modulation at scale.

What carries the argument

The load-bearing object is the thread-safe, regularized lattice Boltzmann update, in which post-collision populations are reconstructed from macroscopic fields and the second-order Hermite coefficient $a^{(2)}_{1,\alpha\beta}$ of the non-equilibrium distribution instead of from stored full distributions; this removes race conditions and permits coalesced memory access on GPUs. The interface is advanced by a conservative Allen-Cahn phase-field equation with a hyperbolic-tangent equilibrium profile, and the advective term is discretized with a fifth-order WENO scheme. MPI handles domain decomposition with a two-lattice halo, while OpenACC directives provide GPU offloading from a single Fortran source, which is what allows the same code to run on NVIDIA and AMD accelerators.

What would settle it

Run the bubble-laden case at $\mathrm{Re}_\lambda \approx 114$ on grids of $512^3$, $1024^3$, and $2048^3$, and with the interface thickness cut by factors of two while keeping the Weber and Froude numbers fixed; if the $-3$ scaling systematically weakens or moves as resolution improves, the attribution to physical bubble-induced dissipation collapses. A complementary check is to compare the measured energy dissipation with resolved viscous dissipation: an excess that grows with refinement would indicate a numerical origin.

Watch

Extended reading notes

Core claim

The central claim is that a thread-safe regularized lattice Boltzmann method, coupled to a conservative Allen-Cahn phase-field interface tracker, gives a scalable and accurate solver for multiphase turbulence on GPU architectures. On the performance side, the paper reports weak-scaling peaks of about 155 GLUPS for single-phase and 105 GLUPS for two-phase flow on 64 GPU devices, with strong-scaling efficiencies that decline in the expected way as communication dominates. On the physics side, single-phase simulations recover the classical $-5/3$ energy and $-7/3$ pressure spectra and the expected structure-function scalings, while bubble-laden homogeneous isotropic turbulence shows a transition to a $-3$ energy spectrum in the intermediate wavenumber range and a rise in velocity-gradient flatness from roughly 5 to 16, which the paper interprets as bubble-induced dissipation and enhanced intermittency.

Load-bearing premise

The central assumption is that the steepened -3 energy spectrum seen in the bubble-laden runs is caused by the bubbles themselves rather than by numerical dissipation; the paper supports this with one parameter set and no test of how the result changes with grid spacing or interface thickness.

Editorial extensions

If this is right

  • Multiphase direct numerical simulations at Taylor Reynolds numbers above 100 become practical on current GPU supercomputers, since the code sustains more than 150 GLUPS and keeps weak-scaling efficiency high up to 64 GPUs.
  • The recurrence of $-5/3$ and $-7/3$ spectra in single-phase tests supports the thread-safe regularized LB formulation as a valid DNS tool for turbulent flows, not only for interface-dominated regimes.
  • The bubble-laden transition to $-3$ scaling, if physical, gives simulations a handle on pseudo-turbulence problems such as bubbly wake dynamics and mass transfer across interfaces.
  • Because the code relies on OpenACC, the same implementation can be benchmarked on NVIDIA and AMD GPU clusters without rewriting kernels, making it a baseline for portable lattice Boltzmann solvers.

Reading between the lines

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

  • Editorial inference: the thread-safe reconstruction from Hermite coefficients could be combined with adaptive mesh refinement or subgrid-scale models; the paper shows the halo and stencil machinery, but does not explore these extensions.
  • Editorial inference: a systematic sweep of bubble volume fraction, Weber number, and Froude number would reveal whether the $-3$ spectrum is a generic feature of buoyant bubble-laden turbulence or a feature of the near-critical regime studied here (Weber number near 1, Froude number near 2).
  • Editorial inference: measuring the spectral energy flux in the same simulations would test whether bubbles act as a forward energy sink or merely re-route energy between scales; the paper does not report that flux.
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 / 6 minor

Summary. The paper presents accLB, a Fortran-OpenACC lattice Boltzmann solver for multiphase turbulence on multi-GPU systems. The code couples a conservative Allen-Cahn phase-field equation with a regularized, thread-safe LBM, and uses MPI+OpenACC for hybrid parallelism. Performance is reported on Leonardo (NVIDIA A100) and LUMI (AMD MI250X), with weak scaling up to 64 GPUs and sustained throughput above 150 GLUPS. Physical validation includes single-phase HIT at Reλ=150, 230, 370 (energy and pressure spectra, structure functions, gradient statistics) and bubble-laden HIT at Reλ=65, 86, 114 with 6% void fraction, where a steeper -3 energy spectrum and enhanced intermittency are reported. The central claim is that accLB is a scalable and accurate platform for multiphase turbulence, with the -3 spectrum attributed to bubble-induced dissipation as seen in experiments and Risso's theory.

Significance. If the stated results hold, the paper provides a useful data point for portable GPU programming of multiphase LB methods, since OpenACC portability across NVIDIA and AMD hardware is demonstrated with concrete scaling data on two pre-exascale systems. The single-phase validation reproduces standard HIT benchmarks (energy and pressure spectra, structure functions, skewness/flatness), and the bubble-laden simulations capture a -3 spectral range and increased flatness consistent with existing experiments and Risso's model. The physical claims are testable and the code's performance numbers are reproducible in principle from the described benchmarks. The main value is in the combination of portable performance and a physically relevant demonstration, rather than a new theoretical result.

major comments (3)
  1. [§5.2, Figure 5] The central physical validation—the -3 energy spectrum—is attributed to bubble-induced dissipation, but the manuscript reports only a single parameter set (6% void fraction, γ=0.01, We_T≈1, Fr_T≈2) and provides no grid-convergence or interface-width sensitivity study. Since the single-phase runs at Reλ=230 and 370 are explicitly under-resolved (Kolmogorov scale 0.4Δx and 0.15Δx), the numerical dissipation behavior of the regularized LB/WENO-5 scheme is not characterized at the resolutions relevant to the bubble-laden runs (Reλ=65–114, with no resolution statement). Without evidence that the -3 plateau persists under grid refinement and interface-width variation at fixed physical parameters, the claim that the steep spectrum is physical rather than numerical is not established. Please add such convergence tests, or explicitly temper the physical attribution in the conclusions.
  2. [Abstract, §5.1, §6] The manuscript repeatedly describes the single-phase simulations as direct numerical simulations (Abstract; Section 6), yet Section 5.1 states that only the Reλ=150 case has Kolmogorov scale approximately equal to Δx, while Reλ=230 and 370 have Kolmogorov scale 0.4Δx and 0.15Δx, i.e., they are beyond the strict DNS regime. The DNS label should be reserved for the Reλ=150 case, or the under-resolved cases should be explicitly labeled as under-resolved or implicit-LES-like; otherwise the accuracy claims in the abstract and conclusions are overstated.
  3. [§5.2, Figure 6] The comparison between single-phase and bubble-laden cases is not made at matched Reynolds numbers. In Figure 6 the single-phase reference is at Reλ=150 while the bubble-laden case is at Reλ=114; the text says the simulation parameters are 'otherwise identical', but the resulting Taylor Reynolds numbers differ. Because the bubble-laden spectra are reported at Reλ=65–114 and no single-phase spectra at the same Reλ and forcing parameters are shown, the possibility that the steeper spectrum is a low-Reynolds-number effect or a consequence of different effective forcing is not excluded. Please provide matched-Reλ single-phase reference runs with identical forcing amplitude and wavenumber.
minor comments (6)
  1. [§2.1.1] The sentence 'the convective term in Eq. (4) is discretized using a fifth-order WENO scheme' refers to Eq. (4), but in the manuscript Eq. (4) is the lattice Boltzmann equation; the phase-field advection term appears in Eq. (3). The cross-reference should point to Eq. (3).
  2. [Figure 3 caption] The caption for panel (c) states that the third-order structure function S3(r) is compared with r^2 scaling in the dissipation range and r^{2/3} scaling in the inertial range. These are the scalings for S2(r); the text and Kolmogorov theory give S3(r) = -4/5 ε r in the inertial range. The caption should be corrected.
  3. [§5.2] The bubble-laden runs do not report the bubble diameter, interface width δ, or the number of bubbles in the domain. These parameters are needed for reproducibility and to interpret the Weber and Froude numbers, especially since the interface width controls diffuse-interface dissipation.
  4. [Eq. (5)] The Guo forcing term in Eq. (5) is written with an ambiguous bracket structure: S_i = w_i ( (c_{iα} - u_α)/c_s^2 + (c_{iβ} u_β)/c_s^4 c_{iα} ) F_α. Please verify the indices and parentheses; the standard form is S_i = w_i [ (c_{iα} - u_α)/c_s^2 + (c_{iα} u_β c_{iβ})/(2 c_s^4) ] F_α.
  5. [References] Reference [23] is cited as an arXiv preprint (arXiv:2501.00846); if the companion paper has been published, the reference should be updated with the journal citation.
  6. [Figure 2(b)] The pressure spectrum E_p(k) is not defined explicitly; please state the normalization (e.g., E_p(k) from the angle-averaged pressure field) so that the comparison with the -7/3 scaling can be checked quantitatively.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation; central claims are validated against external benchmarks, theory, and experiments.

full rationale

The paper's central claims—GPU scaling and physical fidelity—are tested against independent external references. Single-phase HIT spectra are compared with Kolmogorov −5/3 and pressure −7/3 scalings, structure functions against r2/3 and 4/5 laws, and bubble-laden spectra against experimental observations [27] and Risso's theoretical model [36]. These are not outputs of the solver's own fit; they are external criteria. The thread-safe/regularized LB formulation is inherited from the authors' earlier papers [29,30,23], but this is methodological provenance, not a derivation whose conclusion equals its premise: the paper's contribution is the accLB implementation, scaling measurements, and validation runs, all of which stand or fall on the reported benchmarks and comparisons. The lack of grid-convergence/interface-width sensitivity noted for the −3 spectrum is an evidentiary limitation concerning numerical dissipation, not a circularity: the prediction is not constructed from the data it is compared with. No load-bearing self-citation chain or fitted-input-renamed-as-prediction is present.

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

The solver depends on a set of numerical and physical parameters that are either not reported or chosen to set target Reynolds and Weber numbers, plus several domain assumptions inherited from prior work. No new physical entities are introduced.

free parameters (5)
  • ABC forcing amplitude A = not reported (tuned to target Re_lambda)
    Set per case to reach Re_lambda=150/230/370 and 65/86/114; without A the turbulence state cannot be reproduced.
  • ABC forcing wavenumber k_f = not reported
    Defines the large scale energy injection; needed to fix the inertial range in spectra.
  • Surface tension coefficient gamma = 0.01 lattice units
    Fixed in bubble-laden runs to obtain turbulent Weber number near 1; no sensitivity study.
  • Phase-field interface width delta and diffusivity D = not reported
    Control Allen-Cahn interface width and numerical dissipation; no values or convergence test given for production runs.
  • Bubble/droplet density and viscosity ratios = not reported
    Section 4 benchmark uses density ratio 100, but the bubble-laden HIT parameters are not tabulated; needed to reproduce the -3 spectrum.
assumptions (5)
  • domain assumption Regularized thread-safe LBM recovers incompressible Navier-Stokes with the quoted forcing terms
    Basis of all simulations; inherited from Refs. [25,29,30], not re-derived here.
  • domain assumption Kolmogorov -5/3 and -7/3 scalings are appropriate validation targets for single-phase HIT
    Used as benchmarks in Section 5.1.
  • domain assumption The -3 spectrum in dispersed flows is a physical signature of bubble-induced dissipation
    The paper's main physical result is interpreted through Risso's model and Martinez Mercado's experiments [36,27].
  • domain assumption ABC forcing sustains statistically stationary HIT without significant anisotropy
    Assumed to maintain the turbulent state; no independent check reported.
  • domain assumption Conservative Allen-Cahn phase-field accurately tracks interfaces at the chosen numerical parameters
    Required for multiphase interface dynamics; no dedicated convergence or validation test in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of accLB: A High-Performance Lattice Boltzmann Code for Multiphase Turbulence on Multi-Gpu Architectures." pith.science (2026). https://pith.science/paper/7SEE2QUI

@misc{pith2026250501126,
  author       = {Pith},
  title        = {Pith review of: accLB: A High-Performance Lattice Boltzmann Code for Multiphase Turbulence on Multi-Gpu Architectures},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7SEE2QUI}},
  note         = {Machine review of arXiv:2505.01126}
}
abstract

In this work, we present accLB, a high-performance Fortran-based lattice Boltzmann (LB) solver tailored to multiphase turbulent flows on multi-GPU architectures. The code couples a conservative phase-field formulation of the Allen-Cahn equation with a thread-safe, regularized LB method to capture complex interface dynamics. Designed from the ground up for HPC environments, accLB employs MPI for distributed memory parallelism and OpenACC for GPU acceleration, achieving excellent portability and scalability on leading pre-exascale systems such as Leonardo and LUMI. Benchmark tests demonstrate strong and weak scaling efficiencies on multiple GPUs. Physical validation includes direct numerical simulations of homogeneous isotropic turbulence (HIT). Further, we examine bubble-laden HIT and observe a transition to a $-3$ energy scaling, as in experiments and theoretical predictions, due to bubble-induced dissipation, along with enhanced small-scale intermittency. These results highlight accLB as a robust and scalable platform for the simulation of multiphase turbulence in extreme computational regimes.

Figures

Figures reproduced from arXiv: 2505.01126 by the authors.

Figure 1
Figure 1. Panels (a) and (b) show the speed up S p and parallel efficiency Ep for the single component case measured in strong and weak scaling on the NVIDIA Ampere™ A100 and AMD Instinct™ MI250 GPU-based HPC infrastructures. Panels (c) and (d) report the corresponding data for the two-component benchmark case. takes the form: f(x) = A   sin(kfz) + cos(kf y) sin(kf x) + cos(kfz) sin(kf y) + cos(kf x)   , w… view at source ↗
Figure 2
Figure 2. Panel (a) shows the energy spectra for Reλ, 150, 230 and 370. The expected −5/3 power-law in the inertial range is well captured across a broad range of wavenumbers. Panel (b) displays the corresponding pressure spectra, which follow the −7/3 scaling in the inertial range. 100 101 102 r 10-7 10-6 10-5 10-4 10-3 10-2 S2(r) S2(r), Re =150 r 2 r 2/3 S2(r), Re =230 S2(r), , Re =370 100 101 102 r 10-11 10-10 10-9 10-8 10… view at source ↗
Figure 3
Figure 3. Panel (a): Normalized probability density function (PDF) of the longitudinal velocity gradients. Panel (b): Second-order structure [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Instantaneous snapshots of the bubble field panel (a) and velocity field (b) at [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Energy spectra for bubble-laden HIT at di [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: Comparison of velocity increment PDFs between single-phase and bubble-laden HIT. The enhanced tails in the latter case indicate [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 37 canonical work pages

  1. [1]

    Two-dimensional turbulence with the lattice boltzmann equation

    R Benzi and S Succi. Two-dimensional turbulence with the lattice boltzmann equation. Journal of Physics A: Mathematical and General, 23 (1):L1, 1990

    work page 1990

  2. [2]

    Streams: A high-fidelity accelerated solver for direct numerical simulation of compressible turbulent flows

    Matteo Bernardini, Davide Modesti, Francesco Salvadore, and Sergio Pirozzoli. Streams: A high-fidelity accelerated solver for direct numerical simulation of compressible turbulent flows. Computer Physics Communications, 263:107906, 2021. ISSN 0010-4655

    work page 2021

  3. [3]

    Lbcuda: A high-performance cuda port of lbsoft for simulation of colloidal systems.Computer Physics Communications, 277:108380,

    Fabio Bonaccorso, Marco Lauricella, Andrea Montessori, Giorgio Amati, Massimo Bernaschi, Filippo Spiga, Adriano Tiribocchi, and Sauro Succi. Lbcuda: A high-performance cuda port of lbsoft for simulation of colloidal systems.Computer Physics Communications, 277:108380,

    work page

  4. [4]

    Particle-laden turbulence: Progress and perspectives

    Luca Brandt and Filippo Coletti. Particle-laden turbulence: Progress and perspectives. Annual Review of Fluid Mechanics, 54(V olume 54, 2022):159–189, 2022. ISSN 1545-4479

    work page 2022

  5. [5]

    Extended self-similarity in the numerical simulation of three-dimensional homogeneous flows

    M Briscolini, P Santangelo, S Succi, and R Benzi. Extended self-similarity in the numerical simulation of three-dimensional homogeneous flows. Physical Review E, 50(3):R1745, 1994

    work page 1994

  6. [6]

    Gpu acceleration of cans for massively-parallel direct numerical simulations of canonical fluid flows

    Pedro Costa, Everett Phillips, Luca Brandt, and Massimiliano Fatica. Gpu acceleration of cans for massively-parallel direct numerical simulations of canonical fluid flows. Computers and Mathematics with Applications, 81:502–511, 2021. ISSN 0898-1221. Development and Application of Open-source Software for Problems with Numerical PDEs

    work page 2021

  7. [7]

    The interaction of droplet dynamics and turbulence cascade

    Marco Crialesi-Esposito, Sergio Chibbaro, and Luca Brandt. The interaction of droplet dynamics and turbulence cascade. Communications Physics, 6(1):5, 2023

    work page 2023

  8. [8]

    Demou, Marco Edoardo Rosti, Pedro Costa, Filippo Spiga, and Luca Brandt

    Marco Crialesi-Esposito, Nicol `o Scapin, Andreas D. Demou, Marco Edoardo Rosti, Pedro Costa, Filippo Spiga, and Luca Brandt. Flutas: A gpu-accelerated finite difference code for multiphase flows. Computer Physics Communications, 284:108602, 2023. ISSN 0010-4655

    work page 2023

Show all 43 references
  1. [9]

    Enhancing performance of high-speed engineering flow computations: The uranos case study

    Francesco De Vanna. Enhancing performance of high-speed engineering flow computations: The uranos case study. Procedia Computer Science, 255:23–32, 2025. ISSN 1877-0509. Proceedings of the Second EuroHPC user day

    2025

  2. [10]

    Improved locality of the phase-field lattice-boltzmann model for immiscible fluids at high density ratios

    Abbas Fakhari, Travis Mitchell, Christopher Leonardi, and Diogo Bolster. Improved locality of the phase-field lattice-boltzmann model for immiscible fluids at high density ratios. Physical Review E, 96(5):053301, 2017

    2017

  3. [11]

    Hybrid recursive regularized lattice boltzmann simulation of humid air with application to meteorological flows

    Yongliang Feng, Pierre Boivin, J ´erˆome Jacob, and Pierre Sagaut. Hybrid recursive regularized lattice boltzmann simulation of humid air with application to meteorological flows. Physical Review E, 100(2):023304, 2019

    2019

  4. [12]

    Hybrid recursive regularized thermal lattice boltzmann model for high subsonic compressible flows

    Yongliang Feng, Pierre Boivin, J ´erˆome Jacob, and Pierre Sagaut. Hybrid recursive regularized thermal lattice boltzmann model for high subsonic compressible flows. Journal of Computational Physics, 394:82–99, 2019

    2019

  5. [13]

    Turbulence: the legacy of AN Kolmogorov

    Uriel Frisch and Andrei Nikolaevich Kolmogorov. Turbulence: the legacy of AN Kolmogorov. Cambridge university press, 1995

    1995

  6. [14]

    Growth of cloud droplets in a turbulent environment

    Wojciech W Grabowski and Lian-Ping Wang. Growth of cloud droplets in a turbulent environment. Annual review of fluid mechanics, 45(1): 293–324, 2013

    2013

  7. [15]

    Parallel computing and domain decomposition

    William D Gropp. Parallel computing and domain decomposition. In Fifth International Symposium on Domain Decomposition Methods for Partial Differential Equations, pages 349–361. Publ by Soc for Industrial and Applied Mathematics Publ, 1992

    1992

  8. [16]

    Discrete lattice e ffects on the forcing term in the lattice boltzmann method

    Zhaoli Guo, Chuguang Zheng, and Baochang Shi. Discrete lattice e ffects on the forcing term in the lattice boltzmann method. Physical review E, 65(4):046308, 2002

    2002

  9. [17]

    Computer architecture: a quantitative approach

    John L Hennessy and David A Patterson. Computer architecture: a quantitative approach. Elsevier, 2011

    2011

  10. [18]

    Study of high–reynolds number isotropic turbulence by direct numerical simulation

    Takashi Ishihara, Toshiyuki Gotoh, and Yukio Kaneda. Study of high–reynolds number isotropic turbulence by direct numerical simulation. Annual Review of Fluid Mechanics, 41:165–180, 2009. doi: 10.1146 /annurev.fluid.010908.165203

    2009

  11. [19]

    E fficient implementation of weighted eno schemes

    Guang-Shan Jiang and Chi-Wang Shu. E fficient implementation of weighted eno schemes. Journal of Computational Physics , 126(1): 202–228, 1996. ISSN 0021-9991. doi: https: //doi.org/10.1006/jcph.1996.0130. URL https://www.sciencedirect.com/science/ article/pii/S0021999196901308

    1996

  12. [20]

    J ¨ahne and H

    B. J ¨ahne and H. Haußecker. Air-water gas exchange. Annual Review of Fluid Mechanics , 30(V olume 30, 1998):443–468, 1998. ISSN 1545-4479

    1998

  13. [21]

    A continuous surface tension force formulation for di ffuse-interface models

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

    2005

  14. [22]

    The lattice boltzmann method

    Timm Kr ¨uger, Halim Kusumaatmaja, Alexandr Kuzmin, Orest Shardt, Goncalo Silva, and Erlend Magnus Viggen. The lattice boltzmann method. Springer International Publishing, 10(978-3):4–15, 2017. 11

    2017

  15. [23]

    Thread-safe multiphase lattice boltzmann model for droplet and bubble dynamics at high density and viscosity contrasts.arXiv preprint arXiv:2501.00846, 2025

    Marco Lauricella, Adriano Tiribocchi, Sauro Succi, Luca Brandt, Aritra Mukherjee, Michele La Rocca, and Andrea Montessori. Thread-safe multiphase lattice boltzmann model for droplet and bubble dynamics at high density and viscosity contrasts.arXiv preprint arXiv:2501.00846, 2025

    2025 arXiv

  16. [24]

    Cans-fizzy: A gpu-accelerated finite di fference solver for turbulent two-phase flows

    Giandomenico Lupo, Peter Wellens, and Pedro Costa. Cans-fizzy: A gpu-accelerated finite di fference solver for turbulent two-phase flows. arXiv preprint arXiv:2502.04189, 2025

    2025 arXiv

  17. [25]

    Increasing stability and accuracy of the lattice boltzmann scheme: recursivity and regularization

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

    2015 arXiv

  18. [26]

    Theory of turbulence

    WD McComb. Theory of turbulence. Reports on Progress in Physics, 58(10):1117, 1995

    1995

  19. [27]

    On bubble clustering and energy spectra in pseudo-turbulence

    Julian Martinez Mercado, Daniel Chehata Gomez, Dennis Van Gils, Chao Sun, and Detlef Lohse. On bubble clustering and energy spectra in pseudo-turbulence. Journal of fluid mechanics, 650:287–306, 2010

    2010

  20. [28]

    Lattice Boltzmann modeling of complex flows for engineering applications

    Andrea Montessori and Giacomo Falcucci. Lattice Boltzmann modeling of complex flows for engineering applications. Morgan and Claypool Publishers, 2018

    2018

  21. [29]

    Thread-safe lattice boltzmann for high-performance computing on gpus

    Andrea Montessori, Marco Lauricella, Adriano Tiribocchi, Mihir Durve, Michele La Rocca, Giorgio Amati, Fabio Bonaccorso, and Sauro Succi. Thread-safe lattice boltzmann for high-performance computing on gpus. Journal of Computational Science, 74:102165, 2023

    2023

  22. [30]

    High-order thread-safe lattice boltzmann model for high performance computing turbulent flow simulations

    Andrea Montessori, Michele La Rocca, Giorgio Amati, Marco Lauricella, Adriano Tiribocchi, and Sauro Succi. High-order thread-safe lattice boltzmann model for high performance computing turbulent flow simulations. Physics of Fluids, 36(3), 2024

    2024

  23. [31]

    A thread-safe lattice boltzmann model for multicomponent turbulent jet simula- tions

    Andrea Montessori, Luiz A Hegele, and Marco Lauricella. A thread-safe lattice boltzmann model for multicomponent turbulent jet simula- tions. AIAA Journal, 63(3):1005–1012, 2025

    2025

  24. [32]

    Tlbfind: a thermal lattice boltzmann code for concentrated emulsions with finite-size droplets

    Francesca Pelusi, Matteo Lulli, Mauro Sbragaglia, and Massimo Bernaschi. Tlbfind: a thermal lattice boltzmann code for concentrated emulsions with finite-size droplets. Computer Physics Communications, 273:108259, 2022

    2022

  25. [33]

    Analysis of the heat transfer fluctuations in the rayleigh– b´enard convection of concentrated emulsions with finite-size droplets

    Francesca Pelusi, Stefano Ascione, Mauro Sbragaglia, and Massimo Bernaschi. Analysis of the heat transfer fluctuations in the rayleigh– b´enard convection of concentrated emulsions with finite-size droplets. Soft Matter, 19(37):7192–7201, 2023

    2023

  26. [34]

    Podvigina and A

    O. Podvigina and A. Pouquet. On the non-linear stability of the 1:1:1 abc flow. Physica D: Nonlinear Phenomena , 75(4):471–508, 1994. ISSN 0167-2789. doi: https: //doi.org/10.1016/0167-2789(94)00031-X. URL https://www.sciencedirect.com/science/article/ pii/016727899400031X

    1994

  27. [35]

    Stephen B. Pope. Turbulent Flows. Cambridge University Press, 2000

    2000

  28. [36]

    Theoretical model for k- 3 spectra in dispersed multiphase flows

    Fr ´ed´eric Risso. Theoretical model for k- 3 spectra in dispersed multiphase flows. Physics of fluids, 23(1), 2011

    2011

  29. [37]

    A gpu-ready pseudo-spectral method for direct numerical simulations of multiphase turbulence

    Alessio Roccon. A gpu-ready pseudo-spectral method for direct numerical simulations of multiphase turbulence. Procedia Computer Science, 240:17–30, 2024. ISSN 1877-0509. Proceedings of the First EuroHPC user day

    2024

  30. [38]

    Continuous growth of droplet size variance due to condensation in turbulent clouds

    Gaetano Sardina, Francesco Picano, Luca Brandt, and Rodrigo Caballero. Continuous growth of droplet size variance due to condensation in turbulent clouds. Phys. Rev. Lett., 115:184501, Oct 2015

    2015

  31. [39]

    Schramm, Elaine N

    Laurier L. Schramm, Elaine N. Stasiuk, and D. Gerrard Marangoni. Surfactants and their applications. Annu. Rep. Prog. Chem., Sect. C: Phys. Chem., 99:3–48, 2003

    2003

  32. [40]

    S. Succi. The Lattice Boltzmann equation: For complex states of flowing matter. Oxford University Press, 2018

    2018

  33. [41]

    Succi, G

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

    2019 doi

  34. [42]

    Pressure spectra in turbulent flows in the inertial and the dissipation ranges

    Sipei Zhao, Eva Cheng, Xiaojun Qiu, Ian Burnett, and Jacob Chia chun Liu. Pressure spectra in turbulent flows in the inertial and the dissipation ranges. The Journal of the Acoustical Society of America, 140:4178, 2016. doi: 10.1121 /1.4968881

    2016

  35. [43]

    Xiaojue Zhu, Everett Phillips, Vamsi Spandan, John Donners, Gregory Ruetsch, Joshua Romero, Rodolfo Ostilla-M ´onico, Yantao Yang, Detlef Lohse, Roberto Verzicco, Massimiliano Fatica, and Richard J.A.M. Stevens. Afid-gpu: A versatile navier–stokes solver for wall- bounded turb...

    2018

Pith tools

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