2025-11-16T02:07:11.622279

Evaluation of Statistical Consistency in Synthetic Turbulence under Wavenumber Bounds

Lin, Liu, Wang et al.
The random Fourier method (RFM) is widely employed for synthetic turbulence due to its mathematical clarity and simplicity. However, deviations remain between prescribed inputs and synthetic results, and the origin of these errors has not been fully evaluated. This study aims to systematically evaluate the effects of spectral coefficient calibration, grid constraints, and wavenumber bounds on the accuracy of RFM-generated turbulence. The results show that the recalibration of spectral coefficients is essential to ensure consistency in turbulent kinetic energy. The upper wavenumber bound, determined by grid resolution, controls the overall turbulent kinetic energy level, whereas the lower bound, set by computational domain size, governs the fidelity of the energy spectrum in the low-wavenumber range. Moreover, extending wavenumber bounds that exceed the grid-constrained bounds may improve turbulent kinetic energy accuracy but simultaneously amplifies spectral deviations. These findings clarify the distinct roles of coefficient calibration and wavenumber bounds in ensuring the statistical consistency of synthetic turbulence, providing practical guidance for selecting parameters in computational fluid dynamics and computational aeroacoustics applications.
academic

Évaluation de la Cohérence Statistique de la Turbulence Synthétique sous Contraintes de Nombre d'Onde

Informations Fondamentales

  • ID de l'article: 2510.13089
  • Titre: Evaluation of Statistical Consistency in Synthetic Turbulence under Wavenumber Bounds
  • Auteurs: Hongyuan Lin, Yi Liu, Shizhao Wang, Chun-Hian Lee
  • Classification: physics.flu-dyn
  • Date de publication: 15 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2510.13089

Résumé

La méthode de Fourier aléatoire (RFM) est largement utilisée pour synthétiser la turbulence en raison de sa clarté mathématique et de sa simplicité. Cependant, des écarts subsistent entre les entrées prédéfinies et les résultats synthétisés, et les sources de ces erreurs n'ont pas été suffisamment évaluées. Cette étude vise à évaluer systématiquement l'impact de l'étalonnage des coefficients spectraux, des contraintes de grille et des limites de nombre d'onde sur la précision de la turbulence générée par RFM. Les résultats montrent que le réétallonnage des coefficients spectraux est crucial pour assurer la cohérence de l'énergie cinétique turbulente. La limite supérieure du nombre d'onde, déterminée par la résolution de grille, contrôle le niveau global d'énergie cinétique turbulente, tandis que la limite inférieure, fixée par la taille du domaine de calcul, contrôle la fidélité du spectre énergétique dans la plage des bas nombres d'onde. De plus, l'extension de la plage de nombres d'onde au-delà des limites imposées par la grille peut améliorer la précision de l'énergie cinétique turbulente, mais amplifie simultanément les écarts spectraux.

Contexte et Motivation de la Recherche

Définition du Problème

La turbulence synthétique revêt une importance fondamentale en dynamique des fluides numérique (CFD) et en aéroacoustique numérique (CAA), avec des applications incluant l'étude du bruit à large bande, la génération d'écoulement turbulent d'entrée pour les simulations haute fidélité, et l'étude de la diffusion des particules. Une précision insuffisante de la turbulence synthétique peut introduire du bruit artificiel ou déformer les mécanismes physiques étudiés.

Limitations des Méthodes Existantes

Bien que la méthode de Fourier aléatoire (RFM) soit largement adoptée en raison de sa simplicité conceptuelle et de sa commodité, elle présente deux défauts majeurs en matière de cohérence statistique :

  1. Erreurs au niveau énergétique : manifestées par des écarts dans l'énergie cinétique turbulente totale (TKE)
  2. Erreurs au niveau spectral : manifestées par la distorsion du spectre énergétique turbulent E(κ) dans la plage de nombres d'onde discrétisés

Motivation de la Recherche

Bien que des recherches antérieures aient abordé la détermination des coefficients ou proposé des critères empiriques de nombres d'onde, peu se sont concentrées sur les sources d'erreurs, ce qui peut entraîner des écarts inattendus dans les applications pratiques. Cette étude évalue systématiquement l'impact de l'étalonnage des coefficients spectraux, des contraintes de grille et des limites de nombre d'onde sur la cohérence statistique de la turbulence générée par RFM.

Contributions Fondamentales

  1. Analyse systématique des erreurs : première évaluation systématique de l'impact de l'étalonnage des coefficients spectraux, des contraintes de grille et des limites de nombre d'onde sur la précision de la turbulence synthétique
  2. Identification des sources d'erreurs : clarification du fait que les erreurs TKE proviennent principalement de la troncature aux hauts nombres d'onde, tandis que les erreurs E(κ) proviennent d'une couverture insuffisante des bas nombres d'onde
  3. Perfectionnement du cadre théorique : établissement d'un cadre théorique selon lequel la sélection des limites de nombre d'onde doit tenir compte des contraintes de grille
  4. Principes directeurs pratiques : fourniture de conseils pratiques pour la sélection des paramètres dans les applications CFD et CAA

Détails de la Méthode

Définition de la Tâche

Construire un champ de vitesse fluctuante satisfaisant les contraintes statistiques prédéfinies, de manière à reproduire avec haute fidélité le spectre énergétique cible et l'énergie cinétique turbulente.

Architecture du Modèle

Fondamentaux de la Méthode de Fourier Aléatoire

La composante de vitesse fluctuante u' à une position spatiale arbitraire x peut être exprimée comme :

u'(x) = 2∑(m=1 to M) u_m cos(k_m·x + ψ_m)σ_m

Où :

  • M : nombre total de modes
  • u_m : amplitude du m-ième mode
  • k_m : vecteur d'onde du m-ième mode
  • ψ_m : angle de phase aléatoire
  • σ_m : vecteur unitaire de direction

Spectre de von Kármán-Pao

La forme de spectre de von Kármán-Pao largement utilisée est adoptée :

E(κ) = αu'²/κ_e · (κ/κ_e)⁴/[1+(κ/κ_e)²]^(17/6) · exp[-2(κ/κ_η)²]

Où α est le coefficient spectral, κ_e est le nombre d'onde contenant l'énergie, et κ_η est le nombre d'onde de Kolmogorov.

Étalonnage du Coefficient Spectral

Pour tenir compte des effets du nombre de Reynolds fini, un nombre de Reynolds local est introduit :

Re_L = TKE²/(εν)

Le coefficient spectral étalonné α et le nombre d'onde contenant l'énergie κ_e sont déterminés en résolvant une équation impliquant des fonctions hypergéométriques confluentes.

Points d'Innovation Technique

  1. Cadre d'analyse systématique : établissement d'un cadre d'analyse complet tenant compte de l'étalonnage des coefficients spectraux et des limites de nombre d'onde
  2. Théorie des contraintes de grille : clarification des limites de nombre d'onde imposées par la grille : κ_min = 2π/L, κ_max = π/Δx
  3. Mécanisme de séparation des erreurs : séparation des erreurs TKE et spectrales, analyse respective de leurs facteurs d'influence

Configuration Expérimentale

Paramètres de Référence

  • Fluctuation de vitesse : u' = 0,25 m/s
  • Taux de dissipation : ε = 0,5402 m²/s³
  • Énergie cinétique turbulente prédéfinie : TKE_pre = 0,09375 m²/s²

Configuration de Grille

  • Grille de référence : 128³ points
  • Longueur du domaine de calcul : 0,18π m
  • Nombre de modes : 10³ (distribution linéaire)

Indicateurs d'Évaluation

Erreur Relative TKE

Err_TKE = |TKE_syn - TKE_pre|/TKE_pre

Erreur Relative Spectrale

Err_E(κ_j) = |E_syn(κ_j) - E_pre(κ_j)|/E_pre(κ_j)

Résultats Expérimentaux

Efficacité de l'Étalonnage du Coefficient Spectral

Combinaison de CoefficientsTKE_synErr_TKEErr_E,meanErr_E,rms
Premier ensemble de coefficients0,0775517,28%6,143%6,717%
Deuxième ensemble de coefficients0,0842410,14%6,606%7,200%

L'étalonnage du coefficient spectral réduit l'erreur TKE de 17,28% à 10,14%, validant l'importance de l'étalonnage des coefficients.

Effet de l'Extension des Limites de Nombre d'Onde

Lorsque la limite supérieure est étendue à 10κ_max :

  • L'erreur TKE diminue significativement à 1,239%
  • Mais l'erreur spectrale augmente fortement à 28,50% (moyenne) et 41,15% (RMS)
  • Indiquant que l'extension au-delà des limites imposées par la grille provoque des effets de repliement

Impact de la Limite Supérieure du Nombre d'Onde

Avec une résolution de grille augmentant de 32³ à 512³ :

  • Err_TKE diminue de 32,8% à 0,64%
  • L'erreur spectrale reste essentiellement dans la plage 5-8%
  • Confirmant que la limite supérieure contrôle principalement la cohérence TKE

Impact de la Limite Inférieure du Nombre d'Onde

Effet de l'augmentation du domaine de calcul (réduction de κ_min) :

  • Err_TKE diminue de 36,12% à environ 8%
  • L'erreur spectrale diminue de 17% à moins de 5%
  • Confirmant que la limite inférieure contrôle principalement la fidélité spectrale aux bas nombres d'onde

Travaux Connexes

Développement Historique

  • Kraichnan (1970) : première proposition de la méthode de Fourier aléatoire
  • Karweit et al. (1991) : sélection de la plage de nombres d'onde basée sur des données expérimentales
  • Béchara et al. (1994) : introduction de l'échelle intégrale et du nombre d'onde de Kolmogorov
  • Bailly & Juvé (1999) : connexion entre le domaine de calcul et la résolution de grille
  • Saad et al. (2017) : association de la plage de nombres d'onde à la théorie de l'échantillonnage

Comparaison avec les Contributions de cet Article

Par rapport aux recherches existantes, cet article analyse systématiquement pour la première fois les sources d'erreurs plutôt que de fournir simplement des critères empiriques, offrant une base théorique pour la sélection des paramètres.

Conclusions et Discussion

Conclusions Principales

  1. L'étalonnage du coefficient spectral est une condition préalable à la cohérence TKE : dans une plage de nombres d'onde non restreinte, les coefficients réétallonnés peuvent éliminer les écarts entre les valeurs calculées et la TKE prédéfinie
  2. La limite supérieure du nombre d'onde contrôle principalement la cohérence TKE : l'amélioration de la résolution de grille et donc de la limite théorique supérieure peut réduire significativement l'erreur TKE
  3. La limite inférieure du nombre d'onde contrôle la fidélité spectrale aux bas nombres d'onde : l'expansion du domaine de calcul et donc la réduction de la limite inférieure peuvent améliorer la cohérence E(κ) dans la plage des bas nombres d'onde
  4. La sélection de la plage de nombres d'onde doit tenir compte des contraintes de grille : lorsque les limites dépassent les contraintes imposées par la grille, des effets de repliement se produisent, entraînant une augmentation significative de l'erreur E(κ)

Limitations

  1. Forme spectrale unique : l'étude est principalement basée sur le spectre de von Kármán-Pao ; l'applicabilité d'autres formes spectrales nécessite une vérification supplémentaire
  2. Hypothèse d'isotropie : l'étude suppose une turbulence isotrope ; les cas d'anisotropie dans les applications pratiques ne sont pas abordés
  3. Distribution de modes fixe : une distribution de modes linéaire est adoptée ; l'impact d'autres distributions de modes reste à étudier

Directions Futures

  1. Optimisation du nombre et de la distribution des modes : étude de l'impact du nombre et de la distribution des modes dans des limites de nombres d'onde fixes
  2. Extension à la turbulence anisotrope : extension de la méthode aux scénarios de turbulence anisotrope
  3. Validation dans les applications pratiques : vérification des principes directeurs de sélection des paramètres dans les applications CFD et CAA spécifiques

Évaluation Approfondie

Points Forts

  1. Force systématique : première séparation et analyse systématique de l'impact de différents facteurs sur la précision de la turbulence synthétique
  2. Fondations théoriques solides : établissement d'un cadre théorique complet reliant les sources d'erreurs aux mécanismes physiques
  3. Valeur pratique élevée : fourniture de principes directeurs clairs pour la sélection des paramètres, d'importance majeure pour les applications d'ingénierie
  4. Conception expérimentale rationnelle : analyse progressive de l'impact de chaque facteur par la méthode du contrôle des variables

Insuffisances

  1. Analyse des coûts de calcul manquante : absence de discussion sur l'impact de différentes sélections de paramètres sur l'efficacité computationnelle
  2. Cas de vérification limités : analyse basée sur un seul cas de référence, manque de vérification diversifiée
  3. Comparaison insuffisante avec la turbulence réelle : absence de comparaison directe avec les données DNS ou expérimentales

Impact

  1. Contribution théorique : fourniture d'aperçus théoriques importants pour le domaine de la turbulence synthétique
  2. Valeur d'application : fourniture d'outils pratiques de sélection de paramètres pour les ingénieurs CFD et CAA
  3. Signification méthodologique : le cadre d'analyse établi peut être généralisé à d'autres méthodes de synthèse de turbulence

Scénarios d'Application

  1. Aéroacoustique numérique : prédiction du bruit à large bande et conception acoustique
  2. Génération d'écoulement turbulent d'entrée : définition des conditions aux limites pour les simulations haute fidélité
  3. Étude de la diffusion des particules : modélisation de la diffusion dans les applications environnementales et industrielles

Références Bibliographiques

Cet article cite 21 références importantes couvrant le développement historique et les progrès récents des méthodes de synthèse de turbulence, fournissant une base théorique solide pour la recherche. Les références clés incluent les travaux fondateurs de Kraichnan, l'analyse de la théorie de l'échantillonnage par Saad et al., ainsi que les méthodes récentes d'étalonnage des coefficients par Guglielmi et al.