Buonasera,
Il codice che devo sviluppare tratta l'andamento del fouling factor in uno scambiatore di calore.
Il fouling factor viene considerato come una resistenza che cresce nel tempo: quindi ho scelto di valutare questa equazione differenziale con ode45. Il problema sorge quando devo valutare la decrescita dello ione Ca2+ dipendente dalla resistenza dovuta al fouling.
[Ca2+]f = k*Rf, k costante.
[Ca2+]i = [Ca2+]-[Ca2+]f
Con il pedice "i" si intende la concentrazione attuale, con "f" quella consumata mentre quella senza pedice è la concentrazione all'instante precedente.
Tutti i parametri che dipendono da [Ca2+]i si dovrebbero aggiornare per calcolare la nuovo resistenza.
L'unica idea che mi viene in mente è di integrare la ODE con un loop, ma non riesco a capire come far lavorare la ODE con i parametri appena aggiornati.
Molte grazie per l'aiuto.