Discussion :

[azzeddin]

salut tous le monde
le probléme
dans le maillage d'une carré je veux donner les valeurs de la vitesse au bord
Vx et Vy à chaque coté Vx =Vy = 0 sauf Vx pour coté supérieur égal à 1
nl : nombre des noeuds de ligne
nc :nombre des noeuds de collone

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
ntotal =nl*nc
j'essaye avec ça mais je n'arrive pas :(
 do i = 1, nl
	  vx( i) = 0.
	  vy( i) = 0. 
 enddo
 
   do i=nl,nl,nc*nl-nc
	vx( i) = 0.
        vy( i) = 0.
    enddo
 
do i=nl,nl,nc*nl
	vx( i) = 0.
	vy( i) = 0.
enddo
 
 
do i=nc*nl-nl,nc*nl
	vx( i) = 1.
	vy( i) = 0.
endddo

salutations

[plegat]

Salut

j'essaye avec ça mais je n'arrive pas

Tu essayes de faire quoi exactement? Parce qu'à essayer de lire ton bout de code, on comprend que tu n'as pas compris ce que tu devais faire...


De plus, ça veut dire quoi:

pour coté supérieur égal à 1

?

Essaye de comprendre le problème avant de chercher sa solution...

[FR119492]

Salut!
Contrairement à notre ami Plegat, je fais semblant d'avoir compris: tu cherches à intégrer l'équation de Navier-Stokes en 2d sur un domaine carré. Ton fluide est-il un liquide ou un gaz?
Jean-Marc Blanc

[azzeddin]

Salut



Tu essayes de faire quoi exactement? Parce qu'à essayer de lire ton bout de code, on comprend que tu n'as pas compris ce que tu devais faire...


De plus, ça veut dire quoi:

?

Essaye de comprendre le problème avant de chercher sa solution...

oui monsieur j'ai bien compris le problème je veux donner les conditions de la vitesse au bord d'une carré ;la vitesse a 2 composantes Vx et Vy
Vx=Vy= 0 dans tous le bord du carré sauf le la coté supérieur la 1ère composante Vx=1 et merci

[azzeddin]

Salut!
Contrairement à notre ami Plegat, je fais semblant d'avoir compris: tu cherches à intégrer l'équation de Navier-Stokes en 2d sur un domaine carré. Ton fluide est-il un liquide ou un gaz?
Jean-Marc Blanc

salut oui justement c'est ça le problème
le fluide incompressible visqueux "liquide par exemple "
merci

[plegat]

oui monsieur j'ai bien compris le problème je veux donner les conditions de la vitesse au bord d'une carré ;la vitesse a 2 composantes Vx et Vy
Vx=Vy= 0 dans tous le bord du carré sauf le la coté supérieur la 1ère composante Vx=1 et merci

Je n'ai pas dit que tu n'avais pas compris le problème, mais ce que tu devais en faire...

Je présume, vu que les tableaux Vx et Vy n'ont qu'un indice, que tu utilises des tableaux 1D?

Que signifie la syntaxe:

Code :

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

do i=nl,nl,nc*nl-nc

?

Pour i de nl à nc*nl-nc par pas de nl? C'est ça?
(c'est quel langage au passage?)

Dans ce cas, les deux portions de code suivantes:

Code :

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

1
2
3
4
5
6
7
8
9
do i=nl,nl,nc*nl-nc
	vx( i) = 0.
        vy( i) = 0.
enddo
 
do i=nl,nl,nc*nl
	vx( i) = 0.
	vy( i) = 0.
enddo

travaillent partiellement sur les mêmes indices, non?

Dans tous les cas, ça pourrait nous aider que tu commentes un peu ton code... histoire que l'on sache ce que tu as voulu faire... en détails... pas en répétant le sujet...

salut oui justement c'est ça le problème
le fluide incompressible visqueux "liquide par exemple "
merci

Oserais-je dire que les implémentations de ce problème "pullulent" sur le net?

[FR119492]

Salut!

salut oui justement c'est ça le problème

Alors, il fallait le dire tout de suite!

Maintenant, avec l'équation de Navier-Stokes, il y a un problème: en réalité, tu étudies ce qui se passe dans un domaine 3d que tu réduis en 2d. Ton domaine réduit est limité par

Code :

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

1
2
0 <  x < Xmax
0 < y < Ymax

mais, pour le domaine réel, est-ce que tu as

Code :

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

-infini < z < +infini

ou

Code :

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

z = 0

Jean-Marc Blanc

[azzeddin]

Salut!

Alors, il fallait le dire tout de suite!

Maintenant, avec l'équation de Navier-Stokes, il y a un problème: en réalité, tu étudies ce qui se passe dans un domaine 3d que tu réduis en 2d. Ton domaine réduit est limité par

Code :

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

1
2
0 <  x < Xmax
0 < y < Ymax

mais, pour le domaine réel, est-ce que tu as

Code :

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

-infini < z < +infini

ou

Code :

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

z = 0

Jean-Marc Blanc

oui , je travaille avec Z=0

[azzeddin]

Je n'ai pas dit que tu n'avais pas compris le problème, mais ce que tu devais en faire...

Je présume, vu que les tableaux Vx et Vy n'ont qu'un indice, que tu utilises des tableaux 1D?

Que signifie la syntaxe:

Code :

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

do i=nl,nl,nc*nl-nc

?

Pour i de nl à nc*nl-nc par pas de nl? C'est ça?
(c'est quel langage au passage?)

Dans ce cas, les deux portions de code suivantes:

Code :

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

1
2
3
4
5
6
7
8
9
do i=nl,nl,nc*nl-nc
	vx( i) = 0.
        vy( i) = 0.
enddo
 
do i=nl,nl,nc*nl
	vx( i) = 0.
	vy( i) = 0.
enddo

travaillent partiellement sur les mêmes indices, non?

Dans tous les cas, ça pourrait nous aider que tu commentes un peu ton code... histoire que l'on sache ce que tu as voulu faire... en détails... pas en répétant le sujet...




Oserais-je dire que les implémentations de ce problème "pullulent" sur le net?

ah oui il faut bien comprendre le problème
pour "do i=nl,nl,nc*nl-nc" des boucles pour circuler les cotés du carré
je trouve que une seule indice suffit mais je ne suis pas sur ;
à chaque coté du carré je fais une boucle pour circuler tous les noeuds au bord
pour la coté inf du carré une boucle pour i = 1 à nl je donne à Vx et Vy 0
pour la coté gauche aussi une boucle pour i=nl à nc*nl-nc avec pas nl aussi Vx et Vy = 0
pour la coté droite pour i = nl à nl*nc avec pas nl Vx=Vy=0
et pour la coté supérieur du carré pour i=nl*nc-nl à nc*nl Vx=1 et Vy=0
et merci beaucoup

[FR119492]

Salut!

oui , je travaille avec Z=0

Donc, si je t'ai bien compris, tu étudies l'écoulement d'un liquide visqueux entre deux plaques infiniment rapprochées.
Jean-Marc Blanc

[plegat]

pour la coté inf du carré une boucle pour i = 1 à nl je donne à Vx et Vy 0

ok... donc nl c'est le nombre de colonnes, et pas de lignes (vu que les lignes verticales sont plutôt appelées "colonnes"...)

pour la coté gauche aussi une boucle pour i=nl à nc*nl-nc avec pas nl aussi Vx et Vy = 0
pour la coté droite pour i = nl à nl*nc avec pas nl Vx=Vy=0

Là... tu parcours deux fois le côté... droit! Il faudrait donc en corriger une pour faire le côté gauche.
De plus la borne supérieure serait plutôt nc*nl-nl si tu veux rester cohérent avec tes notations.

et pour la coté supérieur du carré pour i=nl*nc-nl à nc*nl Vx=1 et Vy=0
et merci beaucoup

nl*nc-nl est le point juste en-dessous de nc*nl. Donc à ne pas définir pour le côté supérieur.

Fais-toi un schéma de la position de tes indices, tu prend des petites valeurs pour nl et nc, ça sera très simple de voir si ce que tu fais est correct... ou pas!

[azzeddin]

Salut!

Donc, si je t'ai bien compris, tu étudies l'écoulement d'un liquide visqueux entre deux plaques infiniment rapprochées.
Jean-Marc Blanc

l'écoulement dans une cavité carré

[azzeddin]

ok... donc nl c'est le nombre de colonnes, et pas de lignes (vu que les lignes verticales sont plutôt appelées "colonnes"...)



Là... tu parcours deux fois le côté... droit! Il faudrait donc en corriger une pour faire le côté gauche.
De plus la borne supérieure serait plutôt nc*nl-nl si tu veux rester cohérent avec tes notations.




nl*nc-nl est le point juste en-dessous de nc*nl. Donc à ne pas définir pour le côté supérieur.

Fais-toi un schéma de la position de tes indices, tu prend des petites valeurs pour nl et nc, ça sera très simple de voir si ce que tu fais est correct... ou pas!

je fais un petit shéma et voila mon code mais je n'arrive pas
nl=np nombre de poit de ligne
nc=mp nombre de point de colonne

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Vx(1, )=Vx
Vx(2, )=Vy
coté inf 
do i=1,np
	vx(1, i) = 0.
	vx(2, i) = 0.
enddo
coté gauche
do i=np+1,np,mp*np-np+1             // avec pas np
	vx(1, i) = 0.
	vx(2, i) = 0.
enddo
coté droite 
 
do i=np,np,mp*np                       //avec pas np 
	vx(1, i) = 0.
	vx(2, i) = 0.
enddo
coté sup 
do i=mp*np-np+1,mp*np
	vx(1, i) = 1.
	vx(2, i) = 0.
enddo

et merci beaucoup

[plegat]

mais je n'arrive pas

ok, mais... tu n'arrives pas à quoi???

Ca serait quand même un grand pas pour l'humanité si tu arrivais à nous dire quel est ton soucis exactement...

"Je n'arrive pas", ça ne nous dit rien...

[FR119492]

Salut!

l'écoulement dans une cavité carré

Alors, c'est tout faux, parce que ton problème est tridimensionnel et non bidimensionnel: ta cavité a un fond. La réduction de 3d à 2d n'a de sens que dans deux cas très particuliers:

1. La profondeur est négligeable: Zmax-Zmin est très petit par rapport à Xmax-Xmin et Ymax-Ymin.
2. La profondeur est "infinie": Zmax-Zmin est très grand par rapport à Xmax-Xmin et Ymax-Ymin.

Jean-Marc Blanc

[azzeddin]

Salut!

Alors, c'est tout faux, parce que ton problème est tridimensionnel et non bidimensionnel: ta cavité a un fond. La réduction de 3d à 2d n'a de sens que dans deux cas très particuliers:

1. La profondeur est négligeable: Zmax-Zmin est très petit par rapport à Xmax-Xmin et Ymax-Ymin.
2. La profondeur est "infinie": Zmax-Zmin est très grand par rapport à Xmax-Xmin et Ymax-Ymin.

Jean-Marc Blanc

bonjour
oui dans notre problème on a :
La profondeur est négligeable: Zmax-Zmin est très petit par rapport à Xmax-Xmin et Ymax-Ymin.

[azzeddin]

ok, mais... tu n'arrives pas à quoi???

Ca serait quand même un grand pas pour l'humanité si tu arrivais à nous dire quel est ton soucis exactement...

"Je n'arrive pas", ça ne nous dit rien...

bonjour
je n'arrive pas à donner les valeurs de vitesses au bord avec l'algo en haut ;
je crois qu'il y une erreur au niveau des boucle pour numéroter les noeuds

[plegat]

bonjour
je n'arrive pas à donner les valeurs de vitesses au bord avec l'algo en haut ;
je crois qu'il y une erreur au niveau des boucle pour numéroter les noeuds

Il n'y a pas besoin de venir sur un forum pour cela, un simple débuggage sur un carré 3x3 permet de le vérifier et de corriger en conséquence...

Dans tous les cas, mélanger des notations np, nc, mp, nl, ça ne facilite pas la correction! Mais tes boucles ont l'air correctes, avec les remarques suivantes:

coté gauche
do i=np+1,np,mp*np-np+1 // avec pas np
vx(1, i) = 0.
vx(2, i) = 0.
enddo

le point mp*np-np+1 est le coin supérieur gauche. Voir si il appartient au côté gauche, ou au bord supérieur.

coté droite

do i=np,np,mp*np //avec pas np
vx(1, i) = 0.
vx(2, i) = 0.
enddo

idem pour mp*np qui est le coin supérieur droit

coté sup
do i=mp*np-np+1,mp*np
vx(1, i) = 1.
vx(2, i) = 0.
enddo

Là tu réaffectes les points mp*np-np+1 et mp*np que tu as définis auparavant. Quelles valeurs sont les bonnes???

Un conseil. Si tu as du mal avec les tableaux 1D, passe sur des tableaux 2D, ça sera plus simple à gérer.

Mais l'erreur n'est pas là. De toute façon on ne sait toujours pas quel est le problème en fait, sauf que ça ne marche pas!