Física Acoplada de la Capa Límite

Esta página detalla las ecuaciones físicas, las formulaciones de escala y los estudios de validación científica del esquema de acoplamiento de capa límite de Shardian Atmos (EML-SR).

1. Teoría de la Semejanza de Monin-Obukhov (MOST)

En la capa superficial (el 10% inferior de la capa límite planetaria), los flujos turbulentos de calor y momento se parametrizan utilizando las relaciones de semejanza de Monin-Obukhov:

\[u_* = \frac{\kappa u}{\ln\left(\frac{z}{z_{0m}}\right) - \psi_m\left(\frac{z}{L}\right)}\]

\[\theta_* = \frac{\kappa \Delta\theta}{\ln\left(\frac{z}{z_{0h}}\right) - \psi_h\left(\frac{z}{L}\right)}\]

Donde: * \(u_*\) es la velocidad de fricción (escala de momento). * \(\theta_*\) es la escala de temperatura (escala de flujo térmico). * \(\kappa \approx 0.4\) es la constante de von Kármán. * \(L\) es la longitud de Obukhov, que indica la altura en la que la producción de cizalladura es igual a la producción por flotabilidad de la energía cinética turbulenta. * \(\psi_m\) y \(\psi_h\) son las funciones de corrección de estabilidad para momento y calor.

2. El Modelo de Acoplamiento Shardian EML-SR

Los modelos tradicionales utilizan relaciones empíricas (como Dyer-Businger o Beljaars) para definir \(\psi_m\) y \(\psi_h\). Estas formulaciones suelen fallar en regímenes estables nocturnos (\(z/L > 0\)), provocando pérdidas excesivas de calor o un enfriamiento descontrolado.

Shardian Atmos reemplaza estos términos empíricos por un modelo simbólico no lineal descubierto por regresión simbólica (EML-SR). El resolvedor calcula de forma adaptativa los coeficientes de transferencia superficial (\(C_m\) y \(C_h\)):

\[C_m = \left(\frac{u_*}{u}\right)^2, \quad C_h = \frac{u_* \theta_*}{u \Delta\theta}\]

Clamp Adaptativo de Vegetación (Canopy Clamping)

En celdas de rejilla con una alta cobertura vegetal (\(\text{VEGFRA} > 0.7\)), la extracción de momento se regula mediante un factor de clamp de vegetación para evitar la desaceleración artificial del flujo cerca de la superficie:

\[z_{0m,\text{clamped}} = \max\left(z_{0m,\text{min}}, \; z_{0m} \cdot \exp\left(-\beta_{\text{canopy}} \cdot \text{VEGFRA}\right)\right)\]

Esto limita el esfuerzo cortante, permitiendo predecir velocidades de viento realistas a la altura de buje de las turbinas eólicas (80m a 120m).

3. Validación en Benchmarks Globales

Shardian Atmos ha sido validado utilizando cinco conjuntos de datos meteorológicos globales que representan diversos regímenes de capa límite:

graph TD
    BOMEX[BOMEX: Capa Límite Marina] --> A[Validación EML-SR]
    BUBBLE[BUBBLE: Capa Urbana de Alta Rugosidad] --> A
    COASTAL[COASTAL: Transición Tierra-Agua] --> A
    GABLS1[GABLS1: Capa Polar Altamente Estable] --> A
    WANGARA[WANGARA: Ciclo Diurno Semiárido] --> A

Principales Hallazgos Científicos:

GABLS1 (Capa Límite Estable): Los esquemas tradicionales desacoplan la capa límite, provocando un perfil de temperatura superficial plano. EML-SR mantiene una mezcla turbulenta débil pero continua, reduciendo el sesgo frío nocturno en un 51.3%.
BUBBLE (Canopy Urbano): Al ajustar las escalas de rugosidad térmica sobre topografía muy irregular, Shardian Atmos resuelve la isla de calor urbana con un RMSE un 36.5% menor que MOST.
WANGARA (Ciclo Diurno): Resuelve de manera robusta la transición de la capa convectiva diurna inestable a la estable nocturna sin generar oscilaciones numéricas.