Discussion :

[ginafort]

bonjour a tous,

j'aimerais écrire sous fortran la dérivée seconde de la vitesse avec un pas non régulier.

j'ai écrit cela mais j'aimerais savoir si c'est correct ou pas?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
au centre: d2udy2=2.*(u(i-1,j)-2*u(i,j)+u(i+1,j))/        
     /(((y(j+1)-y(j))**2 + (y(j)-x(j+1))**2)    
 
à droite: d2udy2=(u(i,j)-2*u(i+1,j)+u(i+2,j))
  /(y(j+1)-y(j))**2
 
à gauche: d2udy2=(u(i,j)-2*u(i-1,j)+u(i-2,j))
/(y(j)-y(j-1))**2

je vous remercie

[edfed]

par définition, la dérivée, c'est la difference entre deux échantillons consecutifs, ici, on, bosse en numérique, donc, il s'agit d'echantillons.
et donc, la derivée seconde, c'est la difference entre deux echantillons de la derivée.

donc, il y a besoin de 3 valeurs consecutives pour determiner la derivée seconde.

donc:
d2=(X(n+1)-X(n))-(X(n)-X(n-1))

en fortran, ça doit pas être bien eloigné

evidement, il est possible de simplifier grandement la formule.
d2=X(n+1)+X(n-1)-2X(n)
et si l'on utilise les vraies capacités du numérique, l'opération ne prend alors que quelques instructions.
le coefficient 2 etant remplacé par un simple décalage logique à gauche.

[Ehouarn]

Bonjour,

par définition, la dérivée, c'est la difference entre deux échantillons consecutifs, ici, on, bosse en numérique, donc, il s'agit d'echantillons.
et donc, la derivée seconde, c'est la difference entre deux echantillons de la derivée.

donc, il y a besoin de 3 valeurs consecutives pour determiner la derivée seconde.

donc:
d2=(X(n+1)-X(n))-(X(n)-X(n-1))

en fortran, ça doit pas être bien eloigné

evidement, il est possible de simplifier grandement la formule.
d2=X(n+1)+X(n-1)-2X(n)
et si l'on utilise les vraies capacités du numérique, l'opération ne prend alors que quelques instructions.
le coefficient 2 etant remplacé par un simple décalage logique à gauche.

Sur ce coup, je dirai que tu mérites un point de moins que ton pseudo...
(comment peux-tu même imaginer que la position des points n'intervienne pas dans l'évaluation de la dérivée?)

j'ai ecrie cela mais j'aimerais savoir si c'est correct ou pas?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
 
 
au centre: d2udy2=2.*(u(i-1,j)-2*u(i,j)+u(i+1,j))/        
     /(((y(j+1)-y(j))**2 + (y(j)-x(j+1))**2)    
 
à droite: d2udy2=(u(i,j)-2*u(i+1,j)+u(i+2,j))
  /(y(j+1)-y(j))**2
 
à gauche: d2udy2=(u(i,j)-2*u(i-1,j)+u(i-2,j))
/(y(j)-y(j-1))**2

Ca ne me semble pas correct; comment as-tu obtenu ces formules?
Pour déterminer les stencils de dérivation, on cherche le polynôme d'interpolation passant par les points, puis on évalue sa/ses dérivées; C'est simplement assez fastidieux si on considère le cas général de points non-équirépartis.
Mais ca se trouve dans les livres de méthodes numériques, et aussi sur le net, par exemple là: http://www.lavrovo.fr/Arithmurgistan.../lagrange.html

Ehouarn

[edfed]

Sur ce coup, je dirai que tu mérites un point de moins que ton pseudo...
(comment peux-tu même imaginer que la position des points n'intervienne pas dans l'évaluation de la dérivée?)

es tu bien sur de savoir ce qu'est le calcul differentiel?
la derivée ne peut fournir aucune information sur l'offset de depart d'une courbe.

il s'agit que d'une difference entre deux echantillons concecutifs.

peut être que ma notation t'as induit en erreur.
mais si tu veu, je peu te fournir du code assembleur pour demontrer le calcul differentiel (qui n'a rien de compliqué) tel que je l'ai expliqué.

peut etre que j'aurais du remplacer n par t :/

[FR119492]

Salut!
La meilleure solution me semble être de faire un spline cubique, ce qui revient à remplacer ta fonction par une suite de polynômes du deuxième degré, que tu peux ensuite facilement dériver deux fois.
Jean-Marc Blanc

[ginafort]

merci a tous
voila pour ecrire ma dérivée seconde je me suis référé a:

méthode de calcul numérique d
J.P.Noigier
3ème édition
pp 123-133

l'auteur écrit pour un pas constant et il fais le schéma général de son développement, donc j'ai suivi sa méthode pour l'ordre 1 sauf qu'au lieu de prendre h qui est la meme distance entre les 2pt i-1, i et les 2pt i, i+1 j'ai pris les coordonnées de ces points !

car j'ai un maillage tres fin au niveau d'une parois puis il commence a devenir grossier et j'aimerais calculer la dérivée seconde de la vitesse le long de ces points.

donc est ce que c'est faux ce que j'ai fait ?

je vous remercie