PBL Turbolenza & Dispersione

Contenuto app

Statistiche campionarie: media, varianza, deviazione standard, skewness, kurtosis, covarianza, PDF/CDF.
Eddy-Covariance: covarianze, friction velocity, flussi turbolenti, z/L, L, T* e flusso di uno scalare.
MOST: profili medi U(z), Θ(z), Q(z) con funzioni Ψ convettive e stabili.
Metodo dei profili modificato: stima iterativa di u*, T*, q*, L, H0 e relative incertezze da profili multilivello.
Metodi indiretti: u* a un livello, Bowen ratio, Surface Renewal, SODAR/Weill, ceilometer, radiosondaggio.

Note operative

App pensata per hosting statico su Aruba: basta caricare i file così come sono.
Tutte le elaborazioni avvengono nel browser, senza backend.
I grafici si aggiornano in tempo reale dopo il calcolo o il caricamento di esempi/dati.
I moduli indicano sempre le equazioni PDF di riferimento per trasparenza.

Record attivonessuno

Campioni0

u*-

H0-

z/L-

h mixing-

Riepilogo formule implementate

Statistiche: [12.1], [12.3a-b], [12.4b-c], [12.5a], [12.6b-d]

EC: [13.14d-e], [13.16c], [13.17a-b], [13.18d], [13.19b]

MOST: [15.2a-c], [15.3a-b], [15.4a-c], [15.5], [15.6a-c]

Profili modificati: [15.38]-[15.40], [15.35a-b], [15.36b-c], [15.41]-[15.43]

Indirect / dispersione: [15.47a-h], [15.84a-c], [15.91a-c], [15.98b-e], [15.99c-f], [15.136], [15.135]

Grafico di controllo dataset

Input serie temporali

CSV time-series

Formato consigliato: intestazione con colonne tipo time,u,v,w,theta,T,q,chi_v,c.

Variabile statistica Coppia covarianza

Numero bin PDF Settori rosa/barre vento

Equazioni: media [12.1], varianza [12.3a-b], skewness [12.4b], kurtosis [12.4c], covarianza [12.5a], PDF/CDF [12.6b-d], conversione vento [12.13], velocità scalare media [12.14a].

KPI statistici

Media-

Varianza-

Std-

Skewness-

Kurtosis-

Cov(X,Y)-

V scal. media-

Serie temporale

PDF & CDF

Settori del vento

Anteprima dati

Eddy-Covariance

Colonna u Colonna v Colonna w Colonna θ sonica Colonna T reale Colonna q Colonna χv Colonna scalare c Quota z [m] ρd [kg/m³] ρa [kg/m³] Cp [J kg⁻¹ K⁻¹] k g [m s⁻²] Rotazione assi

Equazioni: u* [13.14d-e], H0 [13.16c], z/L [13.17a], T* [13.17b], correzione sonica [13.18d], flusso scalare [13.19b].

KPI EC

u*-

√(-u'w')-

H0v-

H0-

T*-

z/L-

L-

F scalare-

angolo α-

w'θ'-

w'T'-

w'c'-

Componenti vento

Scatter w vs θ

Varianze e covarianze

Log risultati

Profili MOST - forward

u* [m/s] T* [K] q* [kg/kg] L [m] z0m [m] z0h [m] z0q [m] d [m] Θ(z0h) [K] Q(z0q) [kg/kg] Heights z [m] ΨM stabile ΨH stabile Forma profilo

Equazioni: [15.2a-c], [15.3a-b], [15.4a-c], [15.5], [15.6a-c], [15.100a-c].

Tabella profili

U(z)

Θ(z)

Q(z)

Metodo dei profili modificato

Incolla colonne: z,U,Theta,Q

z0m [m] z0h [m] z0q [m] ρ [kg/m³] Cp [J kg⁻¹ K⁻¹] λ [J kg⁻¹] k g Iter max Tolleranza relativa L ΨM stabile ΨH stabile

Equazioni: [15.38]-[15.40], [15.35a-b], [15.36b-c], incertezze [15.41]-[15.43].

KPI stimati

u*-

δu*-

T*-

δT*-

q*-

δq*-

L-

H0-

δH0-

HE-

Iterazioni-

Stato-

Convergenza L

Fit velocità

Fit temperatura

Fit umidità

u* da un livello - metodo iterativo

U(z) [m/s] z [m] z0 [m] H0 [W/m²] T [K] ρ [kg/m³] Cp [J kg⁻¹ K⁻¹] k g ΨM stabile

Equazioni: iterativo [15.47a-c], diretto Hanna-Chang [15.47d-h].

KPI u* a un livello

u* neutro-

u* iterativo-

u* Hanna-Chang-

z/L-

L-

iterazioni-

Bowen ratio e bilancio energetico

RN [W/m²] G0 [W/m²] Δθ [K] Δq [kg/kg] Cp [J kg⁻¹ K⁻¹] λ [J kg⁻¹]

Equazioni: RN - G0 = H0 + HE [15.84a], B = H0/HE [15.84b], B = Cp/λ · Δθ/Δq [15.84c].

B-

H0-

HE-

Chiusura-

Surface Renewal

z [m] ρ [kg/m³] Cp [J kg⁻¹ K⁻¹] dT/dt [K/s]

Equazioni: ρCpV dT/dt = H0 A [15.91a], H0 = z ρCp dT/dt [15.91b].

H0-

zρCp-

SODAR / Metodo di Weill - convettivo

Incolla colonne: z,sigma_w

T [K] ρ [kg/m³] Cp [J kg⁻¹ K⁻¹] g [m s⁻²]

Equazioni: [15.98b-c], profilo lineare del flusso [15.98e].

w'θ'0-

H0-

h'-

h2-

SODAR - stabile

σw: colonne z,sigma_w. Profilo vento: colonne z,u.

T [K] k g ρ [kg/m³] Cp [J kg⁻¹ K⁻¹]

Equazioni: [15.99c-d], profilo vento [15.99f] per L e z0.

u*-

zi-

L-

z0-

H0-

Ceilometer / backscattering

Incolla colonne: z,b

Bm (opzionale) Bu (opzionale) s min [m] s max [m] Passo h [m] Passo s [m]

Equazione teorica del backscattering: B(z) = (Bm+Bu)/2 - (Bm-Bu)/2 erf((z-h)/s) [15.136]; EZT = 2.77 s.

EZT-

Bm-

Bu-

RMSE-

Radiosondaggio / altezza di rimescolamento

Incolla colonne: z,theta_v,U,V

Ri critico Soglia grad θv [K/m] g

Bulk Richardson: Ri(z₂) = g(z₂-z₁)/θv · [θv(z₂)-θv(z₁)] / ([U(z₂)-U(z₁)]²+[V(z₂)-V(z₁)]²) [15.135].

h grad(θv) max-

h Ri crit-

h stabile soglia-

Ri max-