Evaluación de la Consistencia Estadística en Turbulencia Sintética bajo Límites de Número de Onda

Información Básica

• ID del Artículo: 2510.13089
• Título: Evaluation of Statistical Consistency in Synthetic Turbulence under Wavenumber Bounds
• Autores: Hongyuan Lin, Yi Liu, Shizhao Wang, Chun-Hian Lee
• Clasificación: physics.flu-dyn
• Fecha de Publicación: 15 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.13089

Resumen

El método de Fourier aleatorio (RFM) se utiliza ampliamente para sintetizar turbulencia debido a su claridad matemática y simplicidad. Sin embargo, persisten desviaciones entre los parámetros de entrada preestablecidos y los resultados sintetizados, cuyas fuentes de error aún no han sido evaluadas adecuadamente. Este estudio tiene como objetivo evaluar sistemáticamente el impacto de la calibración de coeficientes espectrales, las restricciones de malla y los límites de número de onda en la precisión de la turbulencia generada por RFM. Los resultados demuestran que la recalibración de coeficientes espectrales es crucial para garantizar la consistencia de la energía cinética turbulenta. El límite de número de onda superior, determinado por la resolución de malla, controla el nivel general de energía cinética turbulenta, mientras que el límite inferior, establecido por el tamaño del dominio computacional, controla la fidelidad del espectro de energía en el rango de números de onda bajos. Además, extender el rango de números de onda más allá de los límites de restricción de malla puede mejorar la precisión de la energía cinética turbulenta, pero simultáneamente amplifica la desviación espectral.

Antecedentes de Investigación y Motivación

Definición del Problema

La turbulencia sintética es de importancia central en dinámica de fluidos computacional (CFD) y aeroacústica computacional (CAA), con aplicaciones que incluyen investigación de ruido de banda ancha, generación de entrada turbulenta para simulaciones de alta fidelidad e investigación de difusión de partículas. Si la precisión de la turbulencia sintética es insuficiente, puede introducir ruido artificial o distorsionar los mecanismos físicos bajo estudio.

Limitaciones de Métodos Existentes

Aunque el método de Fourier aleatorio (RFM) es ampliamente adoptado por su simplicidad conceptual y conveniencia, presenta dos deficiencias principales en consistencia estadística:

1. Errores a nivel energético: manifestados como desviaciones en la energía cinética turbulenta total (TKE)
2. Errores a nivel espectral: manifestados como distorsión del espectro de energía turbulenta E(κ) en el rango de números de onda discretizados

Motivación de la Investigación

Aunque investigaciones previas han abordado la determinación de coeficientes o propuesto criterios empíricos de números de onda, pocas se han enfocado en las fuentes de error, lo que puede resultar en desviaciones inesperadas en aplicaciones prácticas. Este estudio evalúa sistemáticamente el impacto de la calibración de coeficientes espectrales, restricciones de malla y límites de número de onda en la consistencia estadística de la turbulencia generada por RFM.

Contribuciones Principales

1. Análisis sistemático de errores: Primera evaluación sistemática del impacto de la calibración de coeficientes espectrales, restricciones de malla y límites de número de onda en la precisión de la turbulencia sintética
2. Identificación de fuentes de error: Clarificación de que los errores de TKE provienen principalmente del truncamiento de números de onda altos, mientras que los errores de E(κ) provienen de cobertura insuficiente de números de onda bajos
3. Perfeccionamiento del marco teórico: Establecimiento de un marco teórico donde la selección de límites de números de onda debe considerar restricciones de malla
4. Principios de orientación práctica: Provisión de orientación práctica para la selección de parámetros en aplicaciones de CFD y CAA

Explicación Detallada de Métodos

Definición de la Tarea

Construcción de campos de velocidad fluctuante que satisfagan restricciones estadísticas preestablecidas, logrando reproducción de alta fidelidad del espectro de energía objetivo y la energía cinética turbulenta.

Arquitectura del Modelo

Fundamentos del Método de Fourier Aleatorio

El componente de velocidad fluctuante u' en una posición espacial arbitraria x puede expresarse como:

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

Donde:

• M: número total de modos
• u_m: amplitud del m-ésimo modo
• k_m: vector de número de onda del m-ésimo modo
• ψ_m: ángulo de fase aleatorio
• σ_m: vector unitario de dirección

Espectro von Kármán-Pao

Se adopta la forma de espectro von Kármán-Pao ampliamente utilizada:

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

Donde α es el coeficiente espectral, κ_e es el número de onda que contiene energía, y κ_η es el número de onda de Kolmogorov.

Calibración de Coeficientes Espectrales

Para considerar efectos de número de Reynolds finito, se introduce el número de Reynolds local:

Re_L = TKE²/(εν)

El coeficiente espectral calibrado α y el número de onda de contención de energía κ_e se determinan resolviendo ecuaciones que contienen funciones hipergeométricas confluentes.

Puntos de Innovación Técnica

1. Marco de análisis sistemático: Establecimiento de un marco analítico completo que considera simultáneamente calibración de coeficientes espectrales y límites de números de onda
2. Teoría de restricciones de malla: Clarificación de que los límites de números de onda restringidos por malla son κ_min = 2π/L, κ_max = π/Δx
3. Mecanismo de separación de errores: Separación de errores de TKE y errores espectrales, analizando sus factores de influencia respectivamente

Configuración Experimental

Parámetros de Referencia

• Fluctuación de velocidad: u' = 0.25 m/s
• Tasa de disipación: ε = 0.5402 m²/s³
• Energía cinética turbulenta preestablecida: TKE_pre = 0.09375 m²/s²

Configuración de Malla

• Malla de referencia: 128³ puntos
• Longitud del dominio computacional: 0.18π m
• Número de modos: 10³ (distribución lineal)

Métricas de Evaluación

Error Relativo de TKE

Err_TKE = |TKE_syn - TKE_pre|/TKE_pre

Error Relativo Espectral

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

Resultados Experimentales

Efecto de Calibración de Coeficientes Espectrales

Los coeficientes espectrales calibrados reducen el error de TKE de 17.28% a 10.14%, validando la importancia de la calibración de coeficientes.

Efecto de Extensión de Límites de Números de Onda

Cuando el límite superior se extiende a 10κ_max:

• El error de TKE se reduce significativamente a 1.239%
• Pero el error espectral aumenta drásticamente a 28.50% (media) y 41.15% (RMS)
• Indicando que extender más allá de los límites de restricción de malla causa efectos de aliasing

Impacto del Límite de Número de Onda Superior

A medida que la resolución de malla aumenta de 32³ a 512³:

• Err_TKE disminuye de 32.8% a 0.64%
• El error espectral se mantiene esencialmente en el rango de 5-8%
• Confirmando que el límite superior controla principalmente la consistencia de TKE

Impacto del Límite de Número de Onda Inferior

Efecto de aumentar el dominio computacional (disminuyendo κ_min):

• Err_TKE disminuye de 36.12% a aproximadamente 8%
• El error espectral disminuye de 17% a menos del 5%
• Confirmando que el límite inferior controla principalmente la fidelidad espectral de números de onda bajos

Trabajo Relacionado

Desarrollo Histórico

• Kraichnan (1970): Primera propuesta del método de Fourier aleatorio
• Karweit et al. (1991): Selección de rango de números de onda basada en datos experimentales
• Béchara et al. (1994): Introducción de escala integral y número de onda de Kolmogorov
• Bailly & Juvé (1999): Conexión entre dominio computacional y resolución de malla
• Saad et al. (2017): Asociación de rango de números de onda con teoría de muestreo

Comparación de Contribuciones del Presente Trabajo

En comparación con investigaciones existentes, este trabajo analiza sistemáticamente las fuentes de error por primera vez, en lugar de solo proporcionar criterios empíricos, proporcionando una base teórica para la selección de parámetros.

Conclusiones y Discusión

Conclusiones Principales

1. La calibración de coeficientes espectrales es requisito previo para consistencia de TKE: En un rango de números de onda sin restricciones, los coeficientes recalibrados pueden eliminar desviaciones entre valores calculados y TKE preestablecido
2. El límite de número de onda superior controla principalmente la consistencia de TKE: Mejorar la resolución de malla y así elevar el límite superior teórico puede reducir significativamente el error de TKE
3. El límite de número de onda inferior controla la fidelidad espectral de números de onda bajos: Expandir el dominio computacional y así disminuir el límite inferior puede mejorar la consistencia de E(κ) en el rango de números de onda bajos
4. La selección de rango de números de onda debe considerar restricciones de malla: Cuando los límites exceden las restricciones de malla, ocurren efectos de aliasing, causando aumento drástico en errores de E(κ)

Limitaciones

1. Forma espectral única: La investigación se basa principalmente en el espectro von Kármán-Pao; la aplicabilidad de otras formas espectrales requiere verificación adicional
2. Suposición de isotropía: La investigación asume turbulencia isótropa; casos de anisotropía en aplicaciones prácticas no se abordan
3. Distribución de modos fija: Se adopta distribución de modos lineal; los efectos de otras distribuciones requieren investigación adicional

Direcciones Futuras

1. Optimización de cantidad y distribución de modos: Investigación del impacto de cantidad y distribución de modos dentro de límites de números de onda fijos
2. Extensión a turbulencia anisotrópica: Extensión del método a escenarios de campos de turbulencia anisotrópica
3. Validación en aplicaciones prácticas: Verificación de principios de orientación de selección de parámetros en aplicaciones específicas de CFD y CAA

Evaluación Profunda

Fortalezas

1. Sistematicidad fuerte: Primera separación y análisis sistemático del impacto de diferentes factores en la precisión de turbulencia sintética
2. Base teórica sólida: Establecimiento de marco teórico completo, vinculando fuentes de error con mecanismos físicos
3. Alto valor práctico: Provisión de principios de orientación claros para selección de parámetros, con importancia significativa para aplicaciones de ingeniería
4. Diseño experimental razonable: Análisis individual de factores mediante método de control de variables

Insuficiencias

1. Análisis de costo computacional ausente: Falta de discusión sobre impacto de diferentes selecciones de parámetros en eficiencia computacional
2. Casos de verificación limitados: Análisis basado solo en caso de referencia único, careciendo de verificación diversificada
3. Comparación insuficiente con turbulencia real: Falta de comparación directa con datos de DNS o experimentales

Impacto

1. Contribución teórica: Provisión de perspectivas teóricas importantes para el campo de turbulencia sintética
2. Valor de aplicación: Provisión de herramientas prácticas de selección de parámetros para ingenieros de CFD y CAA
3. Significado metodológico: El marco analítico establecido puede generalizarse a otros métodos de turbulencia sintética

Escenarios Aplicables

1. Aeroacústica computacional: Predicción de ruido de banda ancha y diseño acústico
2. Generación de entrada turbulenta: Establecimiento de condiciones de frontera para simulaciones de alta fidelidad
3. Investigación de difusión de partículas: Modelado de difusión en aplicaciones ambientales e industriales

Referencias Bibliográficas

Este artículo cita 21 referencias importantes que abarcan el desarrollo histórico y avances recientes en métodos de turbulencia sintética, proporcionando una base teórica sólida para la investigación. Las referencias clave incluyen el trabajo pionero de Kraichnan, análisis de teoría de muestreo de Saad et al., y métodos recientes de calibración de coeficientes de Guglielmi et al.