Bonjour,
Mon problème d'une façon générale est de résoudre l'équation de la chaleur avec fortran.
J'ai obtenu un fichier de résultats pour les positions x à chaque itération comme suit:
Ce qui signifie x(1)=0.000;x(2)=9.4868...e-3;...P.S:pour toutes les itérations on a les mêmes valeurs de x(i)Xstep adimensionné= 0.100000000000 Realposition selon xr, le xtep réel n'est pas constant
0.00000000000 9.486832980505E-03 1.341640786500E-02 1.643167672515E-02
1.897366596101E-02 2.121320343560E-02 2.323790007724E-02 2.509980079602E-02
2.683281573000E-02 2.846049894152E-02 3.000000000000E-02
0.00000000000 9.486832980505E-03 1.341640786500E-02 1.643167672515E-02
1.897366596101E-02 2.121320343560E-02 2.323790007724E-02 2.509980079602E-02
2.683281573000E-02 2.846049894152E-02 3.000000000000E-02
0.00000000000 9.486832980505E-03 1.341640786500E-02 1.643167672515E-02
1.897366596101E-02 2.121320343560E-02 2.323790007724E-02 2.509980079602E-02
2.683281573000E-02 2.846049894152E-02 3.000000000000E-02....
et un fichier de températures à chaque itération de temps comme suit:
ce qui signifie qu'à chaque itération, on affiche les différentes valeurs de températures à différentes positions. exemple: pour l'itéartion N°1; temps=4.5e-5s;T(x=0)=300; T(x=1)=293.000420000; T(x=2)=293.000000000...Timestep= 4.500000E-05température au cours du temps
iter= 1temps= 4.500000000000E-05 300.000000000 293.000420000
293.000000000 293.000000000 293.000000000 293.000000000
293.000000000 293.000000000 293.000000000 293.000000000
293.000000000
iter= 2temps= 9.000000000000E-05 300.000000000 293.000839933
293.000000042 293.000000000 293.000000000 293.000000000
293.000000000 293.000000000 293.000000000 293.000000000
293.000000000
iter= 3temps= 1.350000000000E-04 300.000000000 293.001259798
293.000000126 293.000000000 293.000000000 293.000000000
293.000000000 293.000000000 293.000000000 293.000000000
293.000000000
iter= 4temps= 1.800000000000E-04 300.000000000 293.001679597
293.000000252 293.000000000 293.000000000 293.000000000
293.000000000 293.000000000 293.000000000 293.000000000
293.000000000...
J'essaye de tracer ces résultats (c-à-d T(x,t) avec gnuplot. Pour ce faire, il faut que je présente un fichier de données sous la forme suivante:
comment procéder pour créer un fichier de ce type (avec les valeurs de x et T en colonnes pour chaque itération) ainsi de séparer ces blocs par de lignes vides.
Merci pour votre coopération,
Partager