A x=0 la normale uscente è -x, quindi la condizione Robin corretta è
dT/dx - (h/k)(T-T_amb) = 0, speculare a x=L dove vale dT/dx + (h/k)(T-T_amb) = 0.
Il segno errato causava T(0,t) sotto T_amb (~12°C) con sorgente attiva.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
Rende la normalizzazione degli input simmetrica: x/L e t/T_END
entrambi in [0,1], indipendentemente dal valore di L in config.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
Il termine convettivo H_CONV/K*(T-T_AMB) domina il residuo Robin di un
fattore H*L/K=10 rispetto al gradiente, rendendo L_bc ~100x sovrastimata
rispetto a L_pde. Sostituisce grad_char=(Q/K)² con _bc_scale=max(Q/K,
H*T_char/K)² (~2.25e6) per bilanciare correttamente i termini della loss.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
Estrae la logica interattiva in cleanup_dir() e aggiunge un menu di primo
livello con le opzioni: FDM, PINN, Entrambi, Annulla.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
- forward(): divide t per T_END prima di passarlo alla rete, evita saturazione
di Tanh per t∈[0,10] e migliora la sensibilità temporale del modello
- _pde_scale: include il picco gaussiano della sorgente come denominatore;
con GAUSS_SIGMA=0.01 il picco (~60 °C/s) supera T_char/T_END (15), rendendo
la loss PDE non normalizzata senza questa correzione
- PATIENCE 100→500, SCHED_PATIENCE 30→150: il training ha ora spazio per
convergere prima che l'early stopping o lo scheduler blocchino l'ottimizzatore
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
W_BC=1 causava temperatura sub-ambiente (11.93°C < T_AMB=20°C) perché
i bordi erano poco vincolati; 5 è bilanciamento tra evitare trivial solution
(W_BC=10) e rispettare le BC di Robin (W_BC=1).
N_F aumentato per coprire meglio il dominio temporale tardo dove l'errore
cresceva (max 12.87°C a t=10s).
GAUSS_SIGMA ridotto per avvicinarsi alla sorgente puntuale del FDM e
ridurre il mismatch fisico che causava L2=11.6%.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
Il training collassava alla soluzione banale T=T_AMB perché W_BC=10 spingeva
Adam a soddisfare le Robin BC (trivialmente, con gradiente zero) sacrificando
la PDE. Fix: W_PDE=10, W_BC=1 così la PDE domina il gradiente fin dal primo
epoch. LBFGS_STEPS: 20→200 perché L-BFGS era l'unico ottimizzatore a fare
progressi reali. forward(): rimossa moltiplicazione t_norm che causava
vanishing gradient su dT/dx e d²T/dx².
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
La normalizzazione introdotta (x,t in [0,1]²) rendeva il minimo banale
net=0 (T=T_AMB ovunque) troppo accessibile, causando il collasso del training.
Soluzione: T = T_AMB + T_char * (t/T_END) * net(x_norm, t_norm).
Così T(x,0) = T_AMB per costruzione (vincolo hard) e la rete deve trovare
soluzioni non banali per t>0. La loss IC resta ma è sempre 0.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
- config.py: aggiunge GAUSS_SIGMA = 0.02 nella sezione parametri fisici
- model.py: T_char, grad_char, pde_scale diventano costanti di modulo (_T_char,
_grad_char, _pde_scale) calcolate una sola volta all'import
- engine.py: closure L-BFGS definita una volta sola fuori dal loop
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
- forward(): divide (x,t) per (L, T_END) prima di passare alla rete,
così le due dimensioni hanno la stessa scala indipendentemente da T_END
- heat_pinn_loss: calcola dT_dt e dT_dx in un singolo backward pass
usando autograd.grad con lista [t_f, x_f]
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
- visualizer.py: sostituisce animations/ con results/pinn/TIMESTAMP/,
nomi fissi (heatmap.html, animation.html, comparison.html) come FDM
- config.py: aggiunge sezioni architettura, sampling, Adam, L-BFGS, loss weights
- model.py: costruisce HeatPINN dinamicamente da HIDDEN_SIZE/N_HIDDEN_LAYERS;
heat_pinn_loss legge pesi W_PDE/W_IC/W_BC da config
- engine.py: tutti i parametri di training letti da config
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
- model.py: aggiunge termine sorgente Gaussiana (σ=0.02) nella PDE loss
per approssimare δ(x − X_SRC); sostituisce BC Neumann a x=0 con Robin
- engine.py: clustering collocation vicino X_SRC anziché x=0;
downsample FDM su entrambi gli assi spaziale e temporale in evaluate_model()
- visualizer.py: downsample FDM su entrambi gli assi prima del plot
- app.py: aggiorna header con fisica corrente
- CLAUDE.md: aggiorna PDE, BC e note architetturali
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
Sostituisce la BC Neumann ghost-cell a x=0 con BC Robin su entrambi
i bordi. Q(t) viene iniettato come termine sorgente puntuale al nodo
più vicino a X_SRC, dopo le BCs per non essere sovrascritto.
Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
davide