Sostituisce il parametro lunghezza_conduzione_y_mm con un termine
q_y = -(h_esterno + h_interno)/spessore * (T - T_amb) ricavato dal
bilancio sull'intero volume cilindrico, eliminando un parametro libero.
Riduce anche l'inerzia del sensore IR a 0.01 s (pirometro fast-response).
Co-Authored-By: Claude Sonnet 5 <noreply@anthropic.com>
- Convezione anche sui bordi laterali x = 0 e x = lunghezza (nuovo
coefficiente h_bordi_W_m2K in ARIA)
- Termine di conduzione lungo y verso il resto della fascetta, assunto a
temperatura ambiente: q_y = -k*(T - T_amb)/L_y^2 con L_y configurabile
(lunghezza_conduzione_y_mm in FASCETTA)
- La temperatura iniziale del campo e' la temperatura ambiente del run
(rimosso temperatura_iniziale_C)
- Fisica estratta in prepara_stato_termico e passo_implicito, condivise
tra simulate.py e plot_animazione.py per evitare duplicazione
- Entrambi i termini sono lineari e restano nella matrice fattorizzata
una volta per run
Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
Nuovo script plot_animazione.py che riproduce la fisica di run_0001
(stesso seed di simulate.py) e anima tre pannelli con asse x allineato:
profilo di flusso q(x) con le sorgenti in transito e il verso di marcia,
campo di temperatura T(x, z) nella sezione con il sensore IR e la sua
linea di vista, temperatura nel punto osservato (vera e con inerzia del
sensore) con cursore temporale. Con backend non interattivo salva una GIF
in dataset/ invece di aprire la finestra.
Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
davide