Salve, sono nuovo e mi sono appena presentato. Chiedo cortesemente il vostro aiuto per la risoluzione della seguente equazione differenziale alle derivate parziali (corredata delle condizioni a contorno):
$(delc)/(delt)$=D/r^2*$(del)/(delr)$(r^2*$(delc)/(delr)$)
-Il primo layer va da 0 a un raggio Ri, in cui il coefficiente di diffusione è D=Df e le condizioni a contorno sono le seguenti:
$\{(t=0, c=cf),(r=0, (∂cf)/(∂r)=0),(r=Ri, (∂cf)/(∂r)=(∂cs)/(∂r))}$
L'ultima condizione è quella all'interfaccia che può anche essere scritta come: r=Ri, cf=10*cs
-Il secondo layer va dal raggio Ri al raggio Rs, in cui il coefficiente di diffusione è D=Ds e le condizioni a contorno sono le seguenti:
$\{(t=0, c=0),(r=Ri, (∂cf)/(∂r)=(∂cs)/(∂r)),(r=Rs, (∂cs)/(∂r)=0)}$
Dispongo di Matlab, e vorrei realizzare un qualche script con il metodo delle differenze finite ma non riesco proprio (altri metodi sono ben accetti ). Grazie anticipatamente.