Simulating two-Dimensional decaying vortex structure with the Carleman-Lattice-Boltzmann method
Abstract
We perform a numerical study of the free decay of two-dimensional nonlinear vortex structures using the second-order-truncated Carleman-Lattice Boltzmann method (CLB2). To this purpose, we simulate the nonlinear vortex system on a 128×128 grid at two initial Reynolds numbers, Re≈15 and Re≈46, by resorting to a highly-tuned multi-GPU implementation of the CLB2 code. At Re ≈ 15, CLB2 shows significant accuracy all along the decay process, up to six diffusive times (corresponding to $950$ time steps), with structure-preserving errors below 10\% in amplitude, where by structure-preserving we imply that such errors can be almost entirely reabsorbed by a global time-dependent rescaling factor. At Re ≈ 46, the evolution is carried out over $2850$ time steps, also corresponding to six diffusive times. In this scenario, however, the errors grow in time well above 10\% in amplitude with a still relatively modest but appreciable effect on the spatial structure of the flow. These results indicate that CLB2 is capable of preserving the nonlinear evolution of coherent structures over sufficiently long time spans to reach an effective zero steady-state. However CLB2 shows a constant sub-dissipative trend (more energy than the corresponding Lattice Boltzmann configurations) with a deteriorating accuracy at increasing time and Reynolds number. The present results shed hope that further developments, such as suitable closures instead of truncation and/or various forms of regularization, may lead to a significant improvement of the CLB2 accuracy, prior or even without escalating to expensive higher-order truncations.
Related
- Residual Carleman lift for nonlinear lattice Boltzmann simulations
- Conserved and relaxing modes enable polynomial-success coherent Carleman lattice Boltzmann evolution
- Representation-split Carleman lattice Boltzmann with near-unit step success probability on quantum computers
- Quantum algorithm for high-dimensional nonlinear stochastic differential equations based on Carleman linearization
- Coupled dynamical Boltzmann transport equations with long-range electron-phonon and electron-electron interactions in two-dimensional materials
BibTeX
@misc{cappelli2026simulating,
title = {{Simulating two-Dimensional decaying vortex structure with the Carleman-Lattice-Boltzmann method}},
author = {Luca Cappelli and Massimo Berneschi and Sauro Succi},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.41653},
eprint = {2610.41653},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.41653},
}
RIS
TY - UNPB
TI - Simulating two-Dimensional decaying vortex structure with the Carleman-Lattice-Boltzmann method
AU - Luca Cappelli
AU - Massimo Berneschi
AU - Sauro Succi
PY - 2026
DA - 2026/10/30
PB - arXiv
DO - 10.48550/arxiv.2610.41653
UR - https://arxiv.org/abs/2610.41653
AB - We perform a numerical study of the free decay of two-dimensional nonlinear vortex structures using the second-order-truncated Carleman-Lattice Boltzmann method (CLB2). To this purpose, we simulate the nonlinear vortex system on a 128×128 grid at two initial Reynolds numbers, Re≈15 and Re≈46, by resorting to a highly-tuned multi-GPU implementation of the CLB2 code. At Re ≈ 15, CLB2 shows significant accuracy all along the decay process, up to six diffusive times (corresponding to $950$ time steps), with structure-preserving errors below 10\% in amplitude, where by structure-preserving we imply that such errors can be almost entirely reabsorbed by a global time-dependent rescaling factor. At Re ≈ 46, the evolution is carried out over $2850$ time steps, also corresponding to six diffusive times. In this scenario, however, the errors grow in time well above 10\% in amplitude with a still relatively modest but appreciable effect on the spatial structure of the flow. These results indicate that CLB2 is capable of preserving the nonlinear evolution of coherent structures over sufficiently long time spans to reach an effective zero steady-state. However CLB2 shows a constant sub-dissipative trend (more energy than the corresponding Lattice Boltzmann configurations) with a deteriorating accuracy at increasing time and Reynolds number. The present results shed hope that further developments, such as suitable closures instead of truncation and/or various forms of regularization, may lead to a significant improvement of the CLB2 accuracy, prior or even without escalating to expensive higher-order truncations.
ER -