Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Trias, F. Xavier, Álvarez-Farré, Xavier, Alsalti-Baldellou, Àdel, Gorobets, Andrey, Oliva, Assensi
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:https://arxiv.org/abs/2512.22994
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917173383921664
author Trias, F. Xavier
Álvarez-Farré, Xavier
Alsalti-Baldellou, Àdel
Gorobets, Andrey
Oliva, Assensi
author_facet Trias, F. Xavier
Álvarez-Farré, Xavier
Alsalti-Baldellou, Àdel
Gorobets, Andrey
Oliva, Assensi
contents Direct and large-eddy simulations of turbulence are often solved using explicit temporal schemes. However, this imposes very small time-steps because the eigenvalues of the (linearized) dynamical system, re-scaled by the time-step, must lie inside the stability region. In practice, fast and accurate estimations of the spectral radii of both the discrete convective and diffusive terms are therefore needed. This is virtually always done using the so-called CFL condition. On the other hand, the large heterogeneity and complexity of modern supercomputing systems are nowadays hindering the efficient cross-platform portability of CFD codes. In this regard, our leitmotiv reads: relying on a minimal set of (algebraic) kernels is crucial for code portability and maintenance! In this context, this work focuses on the computation of eigenbounds for the above-mentioned convective and diffusive matrices which are needed to determine the time-step à la CFL. To do so, a new inexpensive method, that does not require to re-construct these time-dependent matrices, is proposed and tested. It just relies on a sparse-matrix vector product where only vectors change on time. Hence, both implementation in existing codes and cross-platform portability are straightforward. The effectiveness and robustness of the method are demonstrated for different test cases on both structured Cartesian and unstructured meshes. Finally, the method is combined with a self-adaptive temporal scheme, leading to significantly larger time-steps compared with other more conventional CFL-based approaches.
format Preprint
id arxiv_https___arxiv_org_abs_2512_22994
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle An efficient eigenvalue bounding method: CFL condition revisited
Trias, F. Xavier
Álvarez-Farré, Xavier
Alsalti-Baldellou, Àdel
Gorobets, Andrey
Oliva, Assensi
Fluid Dynamics
Computational Physics
Direct and large-eddy simulations of turbulence are often solved using explicit temporal schemes. However, this imposes very small time-steps because the eigenvalues of the (linearized) dynamical system, re-scaled by the time-step, must lie inside the stability region. In practice, fast and accurate estimations of the spectral radii of both the discrete convective and diffusive terms are therefore needed. This is virtually always done using the so-called CFL condition. On the other hand, the large heterogeneity and complexity of modern supercomputing systems are nowadays hindering the efficient cross-platform portability of CFD codes. In this regard, our leitmotiv reads: relying on a minimal set of (algebraic) kernels is crucial for code portability and maintenance! In this context, this work focuses on the computation of eigenbounds for the above-mentioned convective and diffusive matrices which are needed to determine the time-step à la CFL. To do so, a new inexpensive method, that does not require to re-construct these time-dependent matrices, is proposed and tested. It just relies on a sparse-matrix vector product where only vectors change on time. Hence, both implementation in existing codes and cross-platform portability are straightforward. The effectiveness and robustness of the method are demonstrated for different test cases on both structured Cartesian and unstructured meshes. Finally, the method is combined with a self-adaptive temporal scheme, leading to significantly larger time-steps compared with other more conventional CFL-based approaches.
title An efficient eigenvalue bounding method: CFL condition revisited
topic Fluid Dynamics
Computational Physics
url https://arxiv.org/abs/2512.22994