Questo lavoro esamina una discretizzazione nel tempo di tipo estrapolazione lineare di un modello di rilassamento di fase \( 2D \) con convezione e reazione dipendenti dalla temperatura. Il modello consiste in una equazione a derivate parziali di tipo avvezione-diffusione nell'incognita temperatura e in una disequazione variazionale ordinaria con doppio ostacolo \( \pm 1 \) per fase variabile. Sotto condizioni di stabilità si dimostra che questo schema converge con ordini ottimali \( \mathcal{O} (\tau | \log \tau |^{1/2}) \) per temperatura ed entalpia, ed ordini \( \mathcal{O} (\tau^{1/2} | \log \tau |^{1/2}) \) per flusso di calore quando il passo \( \tau \) della discretizzazione temporale tende a \( 0 \).