Researcher profile

Christophe Coreixas

Christophe Coreixas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2026arXiv

Multi-GPU Acceleration of PALABOS Fluid Solver using C++ Standard Parallelism

This article presents the principles, software architecture, and performance analysis of the GPU port of the lattice Boltzmann software library Palabos (J. Latt et al., "Palabos: Parallel lattice Boltzmann solver", Comput. Math. Appl. 81, 334-350, (2021)). A hybrid CPU-GPU execution model is adopted, in which numerical components are selectively assigned to either the CPU or the GPU, depending on considerations of performance or convenience. This design enables a progressive porting strategy, allowing most features of the original CPU-based codebase to be gradually and seamlessly adapted to GPU execution. The new architecture builds upon two complementary paradigms: a classical object-oriented structure for CPU execution, and a data-oriented counterpart for GPUs, which reproduces the modularity of the original code while eliminating object-oriented overhead detrimental to GPU performance. Central to this approach is the use of modern C++, including standard parallel algorithms and template metaprogramming techniques, which permit the generation of hardware-agnostic computational kernels. This facilitates the development of user-defined, GPU-accelerated components such as collision operators or boundary conditions, while preserving compatibility with the existing codebase and avoiding the need for external libraries or non-standard language extensions. The correctness and performance of the GPU-enabled Palabos are demonstrated through a series of three-dimensional multiphysics benchmarks, including the laminar-turbulent transition in a Taylor-Green vortex, lid-driven cavity flow, and pore-scale flow in Berea sandstone. Despite the high-level abstraction of the implementation, the single-GPU performance is similar to CUDA-native solvers, and multi-GPU tests exhibit good weak and strong scaling across all test cases.

preprint2025arXiv

GPU-based compressible lattice Boltzmann simulations on non-uniform grids using standard C++ parallelism: From best practices to aerodynamics, aeroacoustics and supersonic flow simulations

Despite decades of research, creating accurate, robust, and efficient lattice Boltzmann methods (LBM) on non-uniform grids with seamless GPU acceleration remains challenging. This work introduces a novel strategy to address this challenge by integrating simple yet effective components: (1) parallel algorithms in modern C++, (2) conservative cell-centered grid refinement, (3) local boundary conditions, and (4) robust collision models. Our framework supports multiple lattices (D2Q9, D2Q13, D2Q21, D2Q37, D3Q27, etc) tailored to various flow conditions. It includes collision models with polynomial and numerical equilibria, a second distribution for polyatomic behavior, a Jameson-like shock sensor, and generalizes Rohde's refinement strategy. The framework's accuracy and robustness is validated across diverse benchmarks, including lid-driven cavity flows, Aeolian noise, 30P30N airfoil aerodynamics, inviscid Riemann problems, and viscous flows past a NACA airfoil in transonic and supersonic regimes. Modern C++ further enables our framework to reach GPU-native performance, while ensuring high portability, modularity, and ease of implementation. Notably, weakly compressible LBMs achieve state-of-the-art GPU efficiency on non-uniform grids, while fully compressible LBMs benefit from acceleration equivalent to thousands of CPU cores in the most compute-intensive cases. Our advanced performance models incorporate neighbor-list and asynchronous time-stepping effects, providing new insights into the performance decomposition of LB simulations on non-uniform grids. Overall, this study sets a new standard for portable, tree-based LBMs, demonstrating that a combination of well-chosen components can achieve high performance, accuracy, and robustness across various flow conditions. As a final proof-of-concept, adaptive mesh refinement is proposed for subsonic and supersonic applications.

preprint2022arXiv

Local mesh refinement sensor for the lattice Boltzmann method

A novel mesh refinement sensor is proposed for lattice Boltzmann methods (LBMs) applicable to either static or dynamic mesh refinement algorithms. The sensor exploits the kinetic nature of LBMs by evaluating the departure of distribution functions from their local equilibrium state. This sensor is first compared, in a qualitative manner, to three state-of-the-art sensors: (1) the vorticity norm, (2) the Q-criterion, and (3) spatial derivatives of the vorticity. This comparison shows that our kinetic sensor is the most adequate candidate to propose tailored mesh structures across a wide range of physical phenomena: incompressible, compressible subsonic/supersonic single phase, and weakly compressible multiphase flows. As a more quantitative validation, the sensor is then used to produce the computational mesh for two existing open-source LB solvers based on inhomogeneous, block-structured meshes with static and dynamic refinement algorithms, implemented in the Palabos and AMROC-LBM software, respectively. The sensor is first used to generate a static mesh to simulate the turbulent 3D lid-driven cavity flow using Palabos. AMROC-LBM is then adopted to confirm the ability of our sensor to dynamically adapt the mesh to reach the steady state of the 2D lid-driven cavity flow. Both configurations show that our sensor successfully produces meshes of high quality and allows to save computational time.

preprint2021arXiv

PalaCell2D: A framework for detailed tissue morphogenesis

In silico, cell based approaches for modeling biological morphogenesis are used to test and validate our understanding of the biological and mechanical process that are at work during the growth and the organization of multi-cell tissues. As compared to in vivo experiments, computer based frameworks dedicated to tissue modeling allow us to easily test different hypotheses, and to quantify the impact of various biophysically relevant parameters. Here, we propose a formalism based on a detailed, yet simple, description of cells that accounts for intra-, inter- and extra-cellular mechanisms. More precisely, the cell growth and division is described through the space and time evolution of the membrane vertices. These vertices follow a Newtonian dynamics, meaning that their evolution is controlled by different types of forces: a membrane force (spring and bending), an adherence force (inter cellular spring), external and internal pressure forces. In addition to the cells dynamics, our formalism further relies on a lattice Boltzmann method, using the Palabos library, to simulate the diffusion of chemical signals. The latter aims at driving the growth and migration of a tissue by simply changing the state of the cells. All of this leads to an accurate description of the growth and division of cells, with realistic cell shapes and where membranes can have different properties. While this work is of methodological nature, we also propose to validate our framework through simple, yet biologically relevant benchmark tests at both single-cell and full tissue scales. This includes free and chemically controlled cell tissue growth in an unbounded domain. The ability of our framework to simulate cell migration, cell compression and morphogenesis under external constraints is also investigated in a qualitative manner.

preprint2020arXiv

Efficient supersonic flows through high-order guided equilibrium with lattice Boltzmann

A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on an extended equilibrium state that is constructed to reproduce the first 13 moments of the Maxwell-Boltzmann distribution exactly. This extends the validity of the standard 5-constraint (mass, momentum and energy) approach and leads to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high-Reynolds number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the flow past a NACA0012 airfoil (Mach number $\mathrm{Ma}=1.5$, Reynolds number $\mathrm{Re=10^4}$) confirms the excellent behavior of this model in a low-viscosity and supersonic regime. The proposed model is substantially more efficient than the previous 5-moment D3Q343 DDF-LBM and opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach.

preprint2020arXiv

Three-dimensional central-moments-based lattice Boltzmann method with external forcing: A consistent, concise and universal formulation

The cascaded or central-moments-based lattice Boltzmann method (CM-LBM) is a robust alternative to the more conventional BGK-LBM for the simulation of high-Reynolds number flows. Unfortunately, its original formulation makes its extension to a broader range of physics quite difficult. To tackle this issue, a recent work [A. De Rosis, Phys. Rev. E 95, 013310 (2017)] proposed a more generic way to derive concise and efficient three-dimensional CM-LBMs. Knowing the original model also relies on central moments that are derived in an adhoc manner, i.e., by mimicking those of the Maxwell-Boltzmann distribution to ensure their Galilean invariance a posteriori, a very recent effort [A. De Rosis and K. H. Luo, Phys. Rev. E 99, 013301 (2019)] was proposed to further generalize their derivation. The latter has shown that one could derive Galilean invariant CMs in a systematic and a priori manner by taking into account high-order Hermite polynomials in the derivation of the discrete equilibrium state. Combining these two approaches, a compact and mathematically sound formulation of the CM-LBM with external forcing is proposed. More specifically, the proposed formalism fully takes advantage of the D3Q27 discretization by relying on the corresponding set of 27 Hermite polynomials (up to the sixth order) for the derivation of both the discrete equilibrium state and the forcing term. The present methodology is more consistent than previous approaches, as it properly explains how to derive Galilean invariant CMs of the forcing term in an a priori manner. Furthermore, while keeping the numerical properties of the original CM-LBM, the present work leads to a compact and simple algorithm, representing a universal methodology based on CMs and external forcing within the lattice Boltzmann framework.