Nella notazione indiciale/di Einstein, la divergenza di un tensore, sarà un vettore e viene scritta come:
[2u11) + (2u12) + (2u13)]
[2u21) + (2u22) + (2u23)]
[2u31) + (2u32) + (2u33)]
+ 2 [2u21) + 2u22) + 2 (2u23)]
Riprendendo la scrittura di energia cinetica fatta in precedenza, e notando le 4 divergenze, possiamo scrivere:
+ 2kt k uj = 1 / p - 2 [k| ρuj uxi+ Lj
+ 2k
2 k [2uj - 2uxiuj + fu.i
[2uj - 2u
Lj fu i con:
j, = -ku - pfu,
÷ rappresentano i fenomeni di trasporto dell'energia cinetica
ψ flusso di energia cinetica
linee di livello
livello alto di energia cinetica
spostamento di energia cinetica, da livelli alti a livelli bassi
Nel caso del 2o termine (equazione iniziale dell'energia cinetica), sarà trasportata dal gradiente della velocità, più è grande il gradiente più sarà grande il trasporto.
Nella notazione indiciale di Einstein, la divergenza di un tensore, sarà un vettore e viene scritta come:
∂/∂xj [
| (∂u1/∂x1), (∂u1/∂x2), (∂u1/∂x3) |
| (∂u2/∂x1), (∂u2/∂x2), (∂u2/∂x3) |
| (∂u3/∂x1), (∂u3/∂x2), (∂u3/∂x3) |
]
Riprendendo la scrittura di energia cinetica fatta in precedenza, e notando le 4 divergenze, possiamo scrivere:
∂k/∂t + ∂/∂xj kuj = - 1/⍴ ∂p/∂xi - ⍴ ∂k/∂xj - ∂ui/∂xj ∂ui/∂xj + fi ui
→ ∂k/∂t = ∂Ψj/∂xj - ⍴ ∂ui/∂xj ∂ui/∂xj + fi ui
con: Ψj = -kuj -⍴ pui + ⍴∂k/∂xj
Flusso di energia cinetica
Nel caso del 2° termine (equazione iniziale dell' energia cinetica), sarà trasportata dal gradiente della velocità, più è grande il gradiente più sarà grande il trasporto.
Nel caso del 3° termine sarà trasportata dal gradiente della velocità più la pressione.Per il 4° termine sarà trasportata dalla viscosità.
Andiamo ad integrare l'equazione iniziale, perché così andiamo a vedere il contenuto dell'insieme turbolento:
∫V ∂k/∂t dV = -∫V ∂k ui/∂xi dV - ∫V 1/ρ ∂p ui/∂xi dV + ∫V ν ∂2k/∂xj∂xj dV - ∫V ν ∂ui uj/∂xj dV + ∫V fi ∂ui dV
1° termine:
Posso portare fuori la variabile di integrazione parziale non essendo una variabile di integrazione:
∫V ∂k/∂t dV = d/dt ∫V k dV
2°/3°/4° termini:
Tutti e 3 sono divergenze di un flusso, utilizzando il Teorema di Gauss della divergenza:
∫V ∂/∂xj dV = ∫∂V ui n̂j ds
Nel nostro caso:
-∫V ∂k ui/∂xi dV - ∫V 1/ρ ∂p ui/∂xi dV + ∫V ν ∂2k/∂xj∂xj dV = ∫∂V -(k ui - 1/ρ p ui δij + ν ∂k/∂xi) n̂j ds
∂V → SUPERFICIE CHE CONTORNA IL VOLUME V
PER LA DIVERGENZA DI xi DOBBIAMO INSERIRE IL DELTA
DI KRONECHER PERCHE SAREBBE:
∫V 1/ρ ∂ρxi / ∂xi dV = -∫∂V 1/ρ ρxi ni dS
MA NEGLI ALTRI TERMINI HO L'INDICE j, QUINDI INSERISCO IL
DELTA DI KRONECHER: