no save
Assistance
Achat
News

Forum | programmation
résolution de l'equation de chaleur
mirinda , le ven. 02 mai 2008 à 15:06:20
salut ,les liens ce sont utile mais malheureusement on n'a pas encore fait des cours concernant ceci et moi je cherche des trucs trés détaillé.voilà le programme en scilab que j'ai fait mais je n'ai pas pu le terminer :

lg = 10. ; // lg = demie longueur du domaine
dx = 0.05 ; // dx = pas d'espace
nx = lg/dx ; // nx = demi nombre de mailles
cfl = 0.4 ; // cfl
dt = dx*dx*cfl ; // dt = pas de temps
nt = 500 ; // nt = nombre de pas de temps effectues
//
// initialisation
//
x=zeros(1,2*nx+1) ;
u0=zeros(1,2*nx+1) ;
for i=1:2*nx+1
x(i) = (i-nx-1)*dx ;
u0(i) = max(0.,1.-x(i)**2) ;
end
u=u0 ;
up=u0 ;
um=u0 ;
uexacte=u0 ;

tics=[4,10,4,10];
plotframe([-lg,-0.1,lg,1.1],tics);
plot2d(x,u0,1,"000")
xtitle ('donnee initiale' ,' ',' ') ;


halt() ;
////////////////////////////////////////////////////////////////
// schema explicite: cfl=0.4
////////////////////////////////////////////////////////////////
//
for n=1:nt
//
up = sh('+1',u) ;
um = sh('-1',u) ;

u=u + dt/(dx*dx)*(up+um-2*u) ;

plotframe([-lg,-0.1,lg,1.1],tics);
plot2d(x,u,[1,1],"100","schema explicite")
plot2d(x,u0,[2,2],"100","donnee initiale")
xtitle('schema explicite, cfl=0.4',' ',' ');
xpause(1000) ;
end

avec sh une fct qui décale un tableau d'element de +1 ou -1 .
ici j'ai traité le cas de shéma explicite , est ce que vous pouvez m'aider à le terminer ?merci et by.
PrécédentSacabouffe
mai 08
Sacabouffe
mai 08
Suivant
REPONSES
Sacabouffe
avr. 08
mirinda
avr. 08
Sacabouffe
avr. 08
mirinda
avr. 08
Sacabouffe
mai 08
mirinda
mai 08
Sacabouffe
mai 08
mirinda
mai 08
Sacabouffe
mai 08
Sacabouffe
mai 08
Version Web
Réalisé par RedShift
no save