Brevet FR3021405A1 — publié le 20151127

Contexte de l'invention L'invention concerne un procédé de mesure d'une grandeur physique d'un écoulement de fluide. Le domaine de l'invention est la mécanique des fluides. Un domaine plus particulier de l'invention concerne la détermination d'une grandeur physique d'un écoulement de fluide autour d'un profil, tel que, par exemple, la pression du fluide et/ou la vitesse du fluide. Un domaine de l'invention concerne également l'optimisation des profils en vue de diminuer les traînées, ainsi que l'amélioration des algorithmes de fonctionnement des simulateurs de vol. Un domaine de l'invention est également la représentation des écoulements de fluides autour d'objets mobiles. On peut citer des profils d'ailes d'avion pour lesquels on souhaite évaluer la traînée et la portance. On connaît des méthodes de calcul en deux dimensions par points singuliers. On connaît également les simulations en trois dimensions par éléments finis. L'une des difficultés en mécanique des fluides est d'acquérir la valeur de la grandeur physique de l'écoulement dans un espace à trois dimensions. Une autre difficulté est de modéliser en mécanique des fluides des profils ou objets complexes en trois dimensions. 2 0 En effet, le fait de devoir mesurer une grandeur physique d'un écoulement de fluide en trois dimensions alourdit considérablement les temps de calcul. L'invention vise à obtenir un dispositif et un procédé, qui permettent de mesurer au moins une grandeur physique d'au moins un écoulement de fluide dans un espace tridimensionnel, qui pallient les inconvénients de l'état de la technique, en 2 5 étant fiables, rapides et directement applicables à la mécanique des fluides. Résumé de l'invention A cet effet, un premier objet de l'invention est un dispositif de mesure d'au moins une première grandeur physique d'au moins un écoulement de fluide dans un espace tridimensionnel ayant au moins une interface prédéterminée, située entre au moins deux milieux, caractérisé en ce que le dispositif comporte : au moins un calculateur ayant des moyens de prescription d'au moins une condition aux limites concernant la première grandeur physique de l'écoulement de 5 fluide prise au moins en un premier point test prédéterminé de l'interface, le calculateur ayant des moyens de calcul configurés pour calculer par méthode DPSM à partir de la condition aux limites, la première grandeur physique de l'écoulement de fluide en au moins un deuxième point de l'espace, différent du premier point test. 10 Un deuxième objet de l'invention est un procédé de mesure d'au moins une première grandeur physique d'au moins un écoulement de fluide dans un espace tridimensionnel ayant au moins une interface prédéterminée, située entre au moins deux milieux, caractérisé en ce que au cours d'une première étape, on prescrit par au moins un calculateur au 15 moins une condition aux limites concernant la première grandeur physique de l'écoulement de fluide prise en au moins un premier point test prédéterminé de l'interface, au cours d'une deuxième étape, on calcule par le calculateur par méthode DPSM à partir de la condition aux limites la première grandeur physique de 20 l'écoulement de fluide en au moins un deuxième point de l'espace, différent du premier point test. Un troisième objet de l'invention est un programme d'ordinateur, comportant des instructions pour la mise en oeuvre du procédé de mesure tel que décrit ci-dessus, lorsqu'il est mis en oeuvre sur un calculateur. 25 Liste des figures des dessins L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif en référence aux dessins annexés, sur lesquels : - les figures 1 A et 1B représentent une interface baignant dans un fluide, dont on cherche à déterminer au moins une grandeur physique par le procédé et le dispositif de mesure suivant un mode de réalisation de l'invention, - la figure 2 représente un synoptique modulaire d'un dispositif de mesure 5 permettant la mise en oeuvre du procédé de mesure suivant un mode de réalisation de l'invention, - la figure 3 représente un organigramme du procédé de mesure suivant un mode de réalisation de l'invention, - les figures 4A, 4B, 36B représentent une interface baignant dans un fluide, 10 dont on cherche à déterminer au moins une grandeur physique par le procédé de mesure et le dispositif suivant un mode de réalisation de l'invention, - les figures 5A, 5B, 6A, 6B, 7A, 7B, 8, 9, 10, 11A, 11B, 12A, 12B, 12C, 12D, 12E, 12F, 12G, 13A, 13B, 14A, 14B, 15A, 15B, 15C, 16A, 16B, 17A, 17B, 17C, 18A, 18B, 19N, 19P, 190, 19Q, 19R, 19S, 19T, 19U, 19V, 36C, 36D, 37A, 15 37B, 37C, 37D, 37E, 37F, 37G, 40A, 40B, 40C, 41A, 41B, 41C, 44A, 44b, 44C, 45A, 45B, 45C représentent des grandeurs physiques d'un écoulement de fluide, ayant été calculées par des modes de réalisation suivant l'invention, - les figures 15D, 15E, 15F, 16C, 16D, 16E, 16F, 16G, 16H, 18C, 18D, 18E, 18F, 18G, 18H, 19B, 19C, 19D, 19E, 19F, 19G, 19H, 191, 19J, 19K, 19L, 19M, 20 36A, 42A, 42B, 42C, 42D, 42E, 42F, 42G, 42H, 43A, 43B, 43C, 43D, 43E, 43F, 43G, 43H, 46A, 46B, 46C, 46D, 46E, 46F, 46G, 46H, 47A, 47B, 47C, 47D, 47E, 47F, 47G, 47H, 48A, 48B, 48C, 48D, 49A, 49B, 49C, 49D, représentent la valeur des sources selon différents exemples, - les figures 19A, 20, 21, 22, 23, 24, 25, 26, 27A, 27B, 27C, 28A, 28B, 28C, 25 29, 30, 31, 32, 33, 34, 38A, 38B, 38C représentent des modes de réalisation du procédé de mesure et du dispositif de mesure suivant l'invention, - la figure 35 représente un organigramme récapitulant 13 cas particuliers qui seront traités dans la suite de la description. 30 Description détaillée de l'invention Aux figures 1A, 1B, 4A et 4B est représentée une interface I située entre un premier milieu M1 et un deuxième milieu M2, dont au moins l'un M2 est soumis à l'écoulement de fluide F. Par exemple, le milieu M1 n'est pas soumis à l'écoulement de fluide F. Par exemple, l'interface I peut être ou comprendre une surface SUR imperméable au fluide F. L'interface peut comprendre au moins une surface d'au moins un solide ou plus généralement d'un objet, en étant appelée profil dans ces cas. Par exemple, le premier milieu M1 est constitué par l'intérieur du solide ou de l'objet, délimité par la surface SUR, tandis que le deuxième milieu M2 est constitué par l'extérieur du solide ou de l'objet. Suivant un mode de réalisation, l'interface I sépare le premier milieu Ml par rapport au deuxième milieu M2. Par exemple, la surface SUR est feuil& autour du premier milieu M1. Par ailleurs, l'objet ou solide ou interface I peut être mobile, comme par exemple une aile d'avion, représentée à titre d'exemple aux figures 1B, 4B et 36B. Le fluide F peut être par exemple l'air, l'eau ou autres. Il peut être prévu plusieurs interfaces I et I' distinctes. Bien entendu, dans ce qui suit, l'interface I et/ou la surface SUR peut être ou comporter une aile d'avion. On décrit ci-dessous un dispositif et un procédé de mesure d'au moins une grandeur physique d'au moins un écoulement de fluide F dans l'espace tridimensionnel, dans lequel se trouve l'interface I prédéterminée. 2 0 On définit par exemple un repère orthonormé x, y, z de l'espace tridimensionnel, ou x et y sont deux directions horizontales et z est une direction verticale ascendante. Le dispositif de mesure comporte au moins un calculateur ayant des moyens CAL1 de prescription d'au moins une condition L aux limites concernant la 25 première grandeur physique de l'écoulement de fluide F prise au moins en un premier point test prédéterminé P1 de l'interface I. Le calculateur CAL a des moyens CAL2 de calcul configurés pour calculer par méthode DPSM à partir de la condition L aux limites, la première grandeur physique de l'écoulement de fluide F en au moins un deuxième point P de l'espace, 3 0 différent du premier point test Pl . L'invention concerne également un dispositif de mesure d'au moins une première grandeur physique (V, P) d'au moins un écoulement de fluide (F) dans un espace tridimensionnel ayant au moins une interface (I) prédéterminée, située entre au moins deux milieux (M1, M2), caractérisé en ce que le dispositif comporte : au moins un calculateur ayant des moyens (E) de prescription d'au moins une condition (L) aux limites concernant au moins une grandeur physique de l'écoulement de fluide (F) prise au moins en un premier point test prédéterminé (P1) de l'interface (I), le calculateur (CAL) ayant des moyens de calcul configurés pour calculer 10 par méthode DPSM à partir de la condition (L) aux limites, au moins une grandeur physique de l'écoulement de fluide (F) en au moins un deuxième point (P) de l'espace, différent du premier point test (Pl). Etapes du procédé de mesure 15 Au cours d'une première étape El, on prescrit par au moins un calculateur au moins une condition L aux limites concernant la première grandeur physique de l'écoulement de fluide F prise en au moins un premier point test prédéterminé P1 de l'interface I. Au cours d'une deuxième étape E2, on calcule par le calculateur CAL par 2 0 méthode DPSM à partir de la condition L aux limites la première grandeur physique de l'écoulement de fluide F en au moins un deuxième point P de l'espace, différent du premier point test Pl. Suivant un mode de réalisation, au cours d'une étape préalable, on détermine les coordonnées de l'interface I dans l'espace. 25 Suivant un mode de réalisation, la première étape El peut comprendre par exemple l'étape préalable, dans laquelle on entre sur des moyens ME d'entrée de données du calculateur CAL, les coordonnées du premier point test prédéterminé P1 de l'interface I, en lequel sont situées les conditions L aux limites prédéterminées. Suivant un mode de réalisation, le calculateur CAL peut être un ordinateur 3 0 ayant comme moyens ME d'entrée de données, par exemple d'un clavier ou de tout autre accès d'entrée de données, et comme moyens SD de sortie de données un écran d'affichage ou autre. Le calculateur ou ordinateur CAL comporte un ou plusieurs programmes informatiques prévus pour mettre en oeuvre les étapes et moyens décrits. Suivant un mode de réalisation, on détell line les coordonnées de l'interface I, et/ou on entre les coordonnées de l'interface I par les moyens ME d'entrée de données du calculateur CAL, et/ou on entre sur les moyens ME d'entrée de données du calculateur CAL les conditions aux limites prédéterminées L en plusieurs points prédéterminés Pl de l'interface I ou en toute l'interface I. Ainsi, le procédé de mesure peut comporter une phase de prétraitement consistant à décrire la 10 topologie tridimensionnelle du problème et en particulier à en dessiner chacun des objets (interface I, surface SUR). Suivant un autre mode de réalisation, la ou les condition(s) L aux limites et/ou le ou les premier(s) point(s) test sont prescrits dans le calculateur CAL, par exemple par un programme informatique, ou en étant préenregistrés dans une 15 mémoire ou en provenant d'un autre calculateur. La méthode DPSM est la méthode des sources ponctuelles réparties, qui sera décrite ci-après. La première grandeur physique peut être par exemple une grandeur macroscopique du fluide. La première grandeur physique peut être par exemple l'une parmi la vitesse 2 0 V de l'écoulement de fluide F, la pression P du fluide F et/ou autre. La grandeur physique est par exemple calculée en un ou plusieurs points P du deuxième milieu M2 ou dans tout le milieu M2. La première grandeur physique peut notamment être le vecteur vitesse de l'écoulement de fluide F dans l'espace en un ou plusieurs points P du milieu M2 ou dans tout le milieu M2. Dans ce qui suit, on 2 5 considère la première grandeur V, P. Suivant un mode de réalisation, le calculateur CAL comporte des deuxièmes moyens CAL3 de calcul, à partir de la première grandeur physique ayant été calculée, d'une propriété de l'interface I dans l'écoulement de fluide F, différente de la première grandeur physique. 3 0 Suivant un mode de réalisation, la propriété est la traînée induite par l'interface I dans l'écoulement de fluide F et/ou la portance de l'interface I dans l'écoulement de fluide F et/ou une force de frottement de l'interface I dans l'écoulement de fluide F. Suivant un mode de réalisation, le calculateur CAL comporte des moyens SD de sortie pour fournir la première grandeur physique ayant été calculée et/ou la propriété ayant été calculée. Suivant un mode de réalisation, au cours d'une troisième étape E3 du procédé, on fournit, sur les moyens SD de sortie de données du calculateur CAL, la grandeur physique ayant été calculée par le calculateur CAL et/ou la propriété ayant été calculée. Suivant un mode de réalisation, le calculateur CAL comporte des troisièmes 10 moyens CAL4 de calcul d'une deuxième grandeur de commande d'un actionneur à partir au moins de la première grandeur physique ayant été calculée et/ou à partir au moins de la propriété ayant été calculée. Suivant un mode de réalisation, les moyens CALI et/ou CAL2 et/ou CAL3 et/ou CAL4 sont automatiques, en étant par exemple chacun mis en oeuvre par un 15 programme informatique. Suivant un autre mode de réalisation, le calculateur envoie sur ses moyens SD de sortie la première grandeur ayant été calculée à un autre calculateur chargé de calculer une grandeur de commande d'un actionneur en fonction au moins de cette première grandeur. 2 0 Le dispositif suivant l'invention comporte des moyens correspondants pour la mise en oeuvre des étapes du procédé. Les moyens ME, CALI, CAL2, CAL3, CAL4 peuvent être répartis sur des ordinateurs ou calculateurs ou processeurs différents, ou être réalisés par un même ordinateur ou calculateur ou processeur. 2 5 Suivant un mode de réalisation, lors de la deuxième étape E2, on prescrit la position d'au moins un troisième point P2 distinct de l'interface I, en lequel se trouve au moins une source S de fluide. On calcule par méthode DPSM au moins un paramètre de la source S de fluide à partir de la condition L aux limites concernant la première grandeur 3 0 physique de l'écoulement de fluide F prise au premier point test Pl de l'interface I. On calcule par méthode DPSM la première grandeur physique de l'écoulement de fluide F située au deuxième point P de l'espace, à partir du paramètre de la source S de fluide, ayant été calculé. Ainsi, l'on positionne au moins une source S de fluide en au moins un 5 deuxième point P2 distinct de l'interface I. Suivant un mode de réalisation, les moyens CAL2 de calcul sont configurés pour calculer par méthode DPSM : - au moins un paramètre d'au moins une source S de fluide modélisée, située en au moins un troisième point P2 de position prescrite, distinct de l'interface I, à 10 partir de la condition L aux limites, - la première grandeur physique de l'écoulement de fluide F située au deuxième point P de l'espace, à partir du paramètre de la source S de fluide, ayant été calculé. Chaque source est ponctuelle. 15 Le ou les paramètres peuvent correspondre par exemple à X ou à a dans ce qui suit. Suivant un mode de réalisation, la source S de fluide est séparée du premier point test prédéterminé Pl de l'interface I. Il peut être prévu une pluralité de sources S de fluide en respectivement une 2 0 pluralité de deuxièmes points P2 distincts de l'interface I. Dans ce cas, les sources ont des positions spatiales distinctes les unes des autres aux points P2. Dans le cas de plusieurs sources dans un groupe, les sources de ce groupe ont des positions spatiales distinctes les unes des autres aux points P2. Suivant un mode de réalisation, l'interface I comprend au moins une surface 25 SUR d'au moins un solide, imperméable au fluide F, au moins une source S de fluide associée au premier point test PI étant située à distance du premier point test P1 en dessous de la surface SUR du solide, de l'autre côté de l'écoulement de fluide F. Il peut être prévu plusieurs sources S de fluide associée au premier point test Pl et situées à distance du premier point test P1 en dessous de la surface SUR du 3 0 solide, de l'autre côté de l'écoulement de fluide F. Il peut être prévu, en association avec une pluralité de premiers points tests P1, une pluralité de groupes comprenant chacun une ou plusieurs sources S de fluide, situées à distance du premier point test Pl en dessous de la surface SUR du solide, de l'autre côté de l'écoulement de fluide F. Dans ces modes de réalisation, la ou les sources S de fluide associées au premier point test P1 sont situées dans le premier milieu MI non soumis à l'écoulement de fluide F. Dans ce cas, la source S est modélisée pour envoyer du fluide (ou une force) à travers l'interface I et la surface SUR dans le deuxième milieu M2 soumis à l'écoulement de fluide F. Par exemple, cette ou ces sources S peuvent être à une distance minimale en dessous de la surface SUR. La distance entre cette ou ces sources et la surface SUR peut être différente entre les sources. Ce mode de réalisation peut être pris en combinaison avec tous les autres modes de réalisations décrits. Ces modes de réalisation sont appelés ci-après premier mode de réalisation sous-jacent. Suivant un mode de réalisation, il est prévu, en plus du ou des premiers points tests de l'interface I, un ou plusieurs autres points tests P5 non situés à 15 l'interface I et situés dans l'écoulement de fluide F. Suivant un mode de réalisation, les moyens de prescription étant prévus pour prescrire au moins une condition aux limites concernant la première grandeur physique de l'écoulement de fluide F prise en ce ou ces autres points tests P5. Suivant un mode de réalisation, les moyens CAL1 de prescription sont 2 0 prévus pour prescrire une condition aux limites comportant au moins l'une parmi : - au premier point test PI à l'interface I une composante noiniale de vitesse du fluide F nulle, - au premier point test Pl à l'interface I une composante tangentielle de vitesse du fluide F nulle, 2 5 - au premier point test Pl à l'interface I une composante tangentielle de vitesse du fluide F prescrite, pouvant être non nulle, - à l'autre point test P5 à la surface délimitant la couche limite CL de l'écoulement de fluide F une composante normale de vitesse du fluide F nulle, - à l'autre point test P5 à la surface délimitant la couche limite CL de 30 l'écoulement de fluide F une composante normale de vitesse du fluide F conservée, - à l'autre point test P5 à la surface délimitant la couche limite CL de l'écoulement de fluide F une composante normale de vitesse du fluide F conservée modulo la masse volumique du fluide, - à l'autre point test P5 à la surface délimitant la couche limite CL de 5 l'écoulement de fluide F une continuité de la composante tangentielle de la vitesse du fluide F, par exemple égale à la vitesse initiale du fluide, - pression du fluide, - continuité de la pression. Par exemple, la vitesse est nulle à la surface des objets imperméables. 10 Suivant un mode de réalisation, il est prévu au moins une source S de fluide non associée au premier point test Pl et située dans l'écoulement de fluide F. Il peut être prévu au moins une source S de fluide en association avec au moins un autre point test P5 prescrit non situé à l'interface I et situé dans l'écoulement de fluide F. Il peut être prévu plusieurs sources S de fluide associées à un autre point test P5 non 15 situé à l'interface I et situé dans l'écoulement de fluide F. Il peut être prévu une pluralité de ces autres points tests P5 non situés à l'interface I et situés dans l'écoulement de fluide F. Il peut être prévu, en association avec une pluralité d'autres points tests P5 non situés à l'interface I et situés dans l'écoulement de fluide F une pluralité de groupes comprenant chacun une ou plusieurs sources S de 20 fluide. Ce ou ces autres points tests P5 peut être à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F. Par exemple, cette ou ces sources S fluide non associée au premier point test P1, et/ou située dans l'écoulement de fluide F et/ou associée avec au moins un autre point test P5, peuvent être à une 25 distance minimale de la surface SCL. La distance entre cette ou ces sources et la surface SCL peut être différente entre les sources. Des moyens sont prévus pour déterminer les coordonnées de la surface SCL de la couche limite. On définit une première zone Z1, ou zone intérieure Z1, située entre la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et 30 l'interface I, et une deuxième zone Z2, ou zone extérieure Z2, située au-dessus de la surface SCL et située de l'autre côté par rapport à l'interface I. La ou les sources S fluide non associée au premier point test Pl peuvent être soit dans la première zone Z1, soit dans la deuxième zone Z2, et à distance de la surface SCL délimitant la couche limite CL. Le fluide est visqueux dans la couche limite, c'est-à-dire dans la première zone intérieure Z1, tandis que le fluide est parfait au-dessus de la couche limite, c'est-à-dire dans la deuxième zone extérieure Z2. La couche limite CL est délimitéé par la surface SCL de glissement. La couche limite CL se prolonge par un sillage. Dans le cas d'une aile d'avion en conditions de vol, la surface SCL de la couche limite CL est à quelques millimètres de la surface SUR de l'aile d'avion, formée dans ce cas par l'interface I. Ainsi que cela est représenté aux figures 1B, 4B et 20, dans les modes de réalisation prévoyant un groupe BL de sources, pouvant être dans ce qui suit, un doublet de sources S1, S2, trois sources S1, S2, S3, appelées également triplet de sources S1, S2, S3 ou autres, ce groupe BL de sources peut être associé à un point test (premier point test Pl situé à l'interface I ou autre point test P5 non situé à l'interface I). Chaque source SI, S2, S3 du groupe BL se trouve à une position P2 différente, d'une source à l'autre du groupe B. Dans ce qui suit, chaque groupe BL ou BL' est également appelé bloc BL ou BL'. Suivant une possibilité de disposition représentée à la figure 20, les sources Sl, S2, S3 du groupe BL peuvent être situées autour d'un point central SB ou ST. Le point central SB ou ST peut se trouver sur la normale N à l'élément EM de maillage de la surface SUR de l'interface I ou de la surface SCL délimitant la couche limite CL. Le point SB est sur la normale N dans le sens rentrant dans l'élément EM de maillage, tandis que le point ST est sur la normale N dans le sens sortant de l'élément EM de maillage. Selon la modélisation, on peut utiliser soit une 2 5 source au point central SB et/ou ST, soit un doublet de sources autour du point central SB et/ou ST, soit un triplet de sources autour du point central SB et/ou ST, ou plusieurs de ces derniers. Bien entendu, cette possibilité de disposition n'est pas limitative et une autre disposition des sources peut être prévue. 30 Calcul des sources et de la grandeur physique de l'écoulement de fluide par méthode DPSM : On décrit ci-dessous plus en détails des modes de réalisation de ce calcul. Suivant un mode de réalisation, dans les moyens CAL2 de calcul est enregistrée une matrice globale de résolution M, comportant au moins un coefficient Mid dépendant à la fois d'une valeur prescrite caractérisant le fluide et de la distance 5 Rif entre une source S de fluide et un quatrième point P de l'espace, le produit de la matrice globale de résolution M, prise aux premiers points test Pl comme quatrième point P de l'espace, multipliée par un premier vecteur J des paramètres des sources S de fluide, étant égal à un deuxième vecteur C des conditions aux limites concernant la première grandeur physique V de l'écoulement 10 de fluide prise aux premiers points test Pl, selon l'équation C = M * J, où * désigne la multiplication. Suivant un mode de réalisation, le coefficient Mid représente la quantité de fluide envoyé de la j-ième source Sj sur le i-ième point test Pi. Mij peut être un réel ou lui-même une matrice (par exemple 3x3), selon le cas. 15 Suivant un mode de réalisation, pour chaque couple de surfaces, on a une matrice de couplage. Suivant un mode de réalisation, la matrice M globale de résolution est composée de matrices de couplage, elles-mêmes composées de coefficients Suivant un mode de réalisation, la valeur prescrite caractérisant le fluide est 2 0 la masse volumique p du fluide et/ou la viscosité cinématique p. du fluide. Suivant un mode de réalisation, il est prévu parmi les moyens CAL2 de calcul - un moyen d'inversion de la matrice globale de résolution M, prise aux premiers points test P1 comme quatrième point P de l'espace, pour calculer la 2 5 matrice inverse M'1, et - un moyen de calcul du premier vecteur J des paramètres des sources S de fluide par multiplication de la matrice inverse M-1 par le deuxième vecteur C des conditions aux limites de la première grandeur physique V de l'écoulement de fluide prise aux premiers points test Pl, selon l'équation 30 J = M-1 * C. Suivant un mode de réalisation, il est prévu parmi les moyens CAL2 de calcul : - un moyen de calcul de la première grandeur physique de l'écoulement de fluide F au deuxième point P de l'espace ayant une position déterminée, en multipliant la matrice globale de résolution M, calculée à ladite position déterminée du deuxième point P comme quatrième point P, par le premier vecteur J des paramètres des sources S de fluide. Il est par exemple prévu une pluralité de quatrièmes points P. La méthode DPSM permet de modéliser par une ou plusieurs sources S de 10 fluide ou de force, placées en des positions prédéterminées, la grandeur physique de l'écoulement de fluide devant être calculée. Le ou les paramètres définissant la ou les sources S de fluide ou de force sont obtenus tout d'abord par calcul à partir des conditions aux limites L de l'interface I, et à partir des positions prédéterminées des sources S de fluide ou de 15 force. Une fois les sources S de fluide ou de force calculées, on obtient à partir de celles-ci la valeur de la grandeur physique de l'écoulement de fluide. Suivant un mode de réalisation, la méthode DPSM comporte de mailler régulièrement les surfaces SUR des objets (ou interfaces I), et éventuellement de la 2 0 la surface SCL de la couche limite CL, afin d'y créér un ensemble de premiers points test P1 (ou P5 pour SCL). Le maillage comporte une pluralité d'éléments EM de maillage de la surface SUR ou SCL, chaque élément EM de maillage étant une portion de plan, par exemple triangulaire ou d'une autre forme, dans lequel se trouve par exemple le premier point test Pl ou l'autre point test P5. 25 Suivant un mode de réalisation, d'un seul côté ou de part et d'autre de chaque point test, on dispose un groupe BL contenant une ou plusieurs sources élémentaires S (appelées également « singularités »), l'ensemble de ces groupes étant destiné in fine à synthétiser les grandeurs physiques dans les milieux respectifs Ml et M2 adjacents à ces surfaces. 3 0 L'interaction entre les sources S et les points tests P 1 appartenant au même objet ou à la même interface I peut ensuite s'écrire sous forme d'une première matrice d'autocouplage. L'interaction entre les sources S et les points tests P1 appartenant à des objets différents va quant à elle s'écrire sous la forme d'une deuxième matrice d'intercouplage. Suivant un mode de réalisation, ces matrices élémentaires satisfaisant l'ensemble des conditions aux limites L entre les différents milieux produisent une troisième matrice globale M, dont l'inversion donne accès à la valeur des sources S. La matrice globale M de résolution et le vecteur C de conditions aux limites permettent, par inversion de cette matrice M, de trouver la valeur de chaque source élémentaire S. Suivant un mode de réalisation, la matrice globale M de résolution est carrée. Suivant un mode de réalisation, un nombre Ns de sources S de fluide et un nombre Np de points tests (premiers points test Pl et éventuellement autres points tests P5) distincts les uns des autres sont prévus, de telle sorte que la matrice globale 15 de résolution M soit carrée. Suivant un mode de réalisation, le nombre Np de points tests (premiers points test Pl et éventuellement autres points tests P5) des conditions aux limites L est supérieur ou égal au nombre Ns de sources S de fluide. Suivant un mode de réalisation, le nombre Np de points tests (premiers 2 0 points test P1 et éventuellement autres points tests P5) des conditions aux limites L est égal au nombre Ns de sources S de fluide. Une étape préalable de détermination de la matrice globale de résolution est prévue, qui est réalisés par les moyens CAL2 de calcul du calculateur ainsi que cela est décrit ci-dessus. 2 5 Suivant un mode de réalisation, le post-traitement consite ensuite à exploiter le fluide (ou force) envoyé par l'ensemble des sources S pour calculer les grandeurs physiques dans tout l'espace (à l'intérieur comme à l'extérieur des objets ou interface I). Par superposition des sources S ayant été calculées, on synthétise par ces 3 0 sources la première grandeur physique. Suivant un mode de réalisation, la méthode DPSM peut être appliquée à un écoulement de fluide F laminaire. La méthode DPSM permet de modéliser en mécanique des fluides des profils ou objets I complexes (par exemple avion ou aile d'avion) en trois 5 dimensions (au lieu de deux dimensions avec les méthodes par points singuliers exploitant uniquement le potentiel scalaire). La méthode DPSM permet de réduire notablement les temps de calculs par la réduction du nombre d'élements (maillage des surfaces et non pas de l'espace). De plus, la méthode DPSM permet de résoudre des problèmes dans lesquels plusieurs objets I sont en interaction (approche d'un 10 avion vers le sol impliquant un effet de sol par exemple). La méthode DPSM conduit à une équation semi-analytique du type C=M*J, où M est la matrice globale de résolution, qui contient l'ensemble des sous-matrices de couplage entre les différents objets (ou interfaces I) du problème, C le vecteur conditions aux limites, et J la valeur des sources, obtenues par inversion de la 15 relation précédente: J=inv(M)*C. Les notations suivantes sont adoptées ci-dessous en ce qui concerne la matrice M, qui peut comprendre d'autres matrices indiquées ci-dessous : - V : le vecteur vitesse de l'écoulement, V = (V', VY, V') - V sN = 2.1V : la composante normale de la vitesse de l'écoulement sur une 2 0 surface S (ou SUR) d'un objet ou d'une interface, - La surface S (ou SUR) a ns points tests P1 , - MBA : matrice de couplage entre les points tests appartenant à l'objet B (ou à l'interface B) et les sources associées à l'objet A (ou à une autre interface A), 25 Pour la surface de l'objet A ou de l'interface A: OEA ou : vecteur des paramètres des sources associées à l'objet A ou à l'interface A, de dimension nA, .)\,A1 : vecteur des paramètres des sources S1 associées à l'objet A ou à l'interface A dans le cas des sources S1 des triplets de sources S1, S2, S3, de 3 0 dimension nA 1 Al : paramètre du vecteur km (c'est-à-dire la composante de dimension lx1 de ce vecteur km ayant ces nm composantes), - Distance entre les points i et j : Ru= \i(x,- +(y,- y i)2 +(zi-z On décrit ci-dessous un mode de réalisation de cette méthode DPSM dans le cas d'une source ponctuelle scalaire au point P2 qui émet un flux ai homogène à un débit en m3/s est donné ci-dessous. Bien entendu, ce mode de réalisation n'est pas limitatif. L'équation de Navier-Stokes est : Di; p. dt = p.k> - P + 2 L'accélération d'une particule de fluide est donnée par : Di; ai; ± v dt at On a donc : av av Y+ av, av, 1 ap rila2V a2Vx 1, 02 r V '' +V - +Vz = + + + x - at ax Y ay az p ax p ax2 aY2 az2 ) av avy av av 1 ap ,7(a2v a2V a2V Y Y Y Y Y -1-1; +V +V = + + + at - 0x 3- ay z az p ay P , axe ay2 0z2 " av z av av av_ 1 ap n / a 2vz 02v 02v z +V +V z +V =-g + ' + ' + 2 0 at 1 z ax v ay az p az p ax ay ay2 az2 L'équation d'Euler, pour un fluide parfait et une viscosité négligeable ri = 0 , Di> p. = dt p.g -VP ce qui donne : 2 5 rav, av ' av, av' 1 ND +vt. +v +v = at ' ax y 0y z az p ax av avy av y av 1 ap Y +, +V +V at 1/ ax Y 0y z az p ay av, av av av z 1 ap at ±v' axz ±vY a; + Vz aZ = g Paz 30 L'équation d'Euler en régime stationnaire donne : av y t +v av av +v av 1 aP X = x ax Y aY z saz p ax av avav. 1 ND Y Y Y Vx + V + V ax Y ay z az p ay v avz +v avz +v avz - g 1 aP x ax Y aY z .. CZ p az est : Pour un écoulement incompressible de fluide : dive =. , ,- ..i5 =av x +av U avz =--- , Y + Pour un écoulement a ay az roi 17 => x -V> irrotationnel de fluide : avz avy ô ay az avx avz az ax avy ÔVx ax ay Cette dernière condition implique que la vitesse du fluide dérive d'un potentiel scalaire : Pour un écoulement incompressible et irrotationnel de fluide : V28 0 Pour un écoulement incompressible et irrotationnel de fluide, les équations d'Euler sont découplées en potentiel 8 et en vitesse v. Dans le mode de réalisation décrit ci-dessous, on suppose dans les équations 2 0 ci-dessous que la méthode DPSM est mise en oeuvre dans le cas d'un écoulement stationnaire, incompressible et irrotationnel de fluide. La solution de l'équation de Poisson, calculée en un point P ou P1=Pi de coordonnées x' yi, z, pour une source S=Sj ponctuelle isotrope (ou scalaire) au point P2 de coordonnées xj, yj, zj, qui émet le flux a ou 2^,, est le potentiel scalaire 8: 25 V28 = 0 = 47z-pR,, avec R,, = \/(x, - )2 + (yi yi )2 + zi )2 Pour la source scalaire de fluide, le calcul de la vitesse V de l'écoulement au point P ou P1=Pi de coordonnées xi, yi, z, dérive du potentiel selon l'équation 3 0 suivante : = (x. - x .) ( (x,- x1) v . .23 4.7z-pR 4;z-pRi3 v =(.Yi-1) .2 (Y 2 ;) Y 4.7z-pRi 4 7i-pR (z -z;) 47-role (z - z;) v = .23 z 47cpR r v r vY (xi - xi) r x .2 vz qi; qi); 47-rpRii3. (yi - 321) = -grad = (zi - z i) Dans ce mode de réalisation, les conditions aux limites L peuvent comprendre des conditions aux limites intrinsèques (CLI), qui prévoient qu'à l'interface entre deux milieux définis par leur densité respectives (pi, p2), on considère la continuité du potentiel et la vitesse normale (modulo rho) à l'interface en chaque point de l'interface. On écrit : I81 = 82 = /02-V2N Dans ce mode de réalisation, les conditions aux limites L peuvent comprendre des conditions aux limites utilisateurs (CLU) qui permettent de fixer 2 0 des conditions aux limites sur les objets : vitesse normale nulle et/ou vitesse nulle, et/ou autres. Pour une source S ponctuelle isotrope au point P2 de coordonnées xj, y j, zi' constituée d'un triplet de source (ou groupe de sources ou bloc de sources, le triplet de source se situant autour du point P2 aux coordonnées (xik,Yik,zik) k=1,2,3) qui 2 5 émet le flux a ou 23) (homogène à un débit en kg/s), la vitesse au point P ou P1=Pi de coordonnées xi, yi, zi est : r xi - Xjixi - X12 - X13 4;TR 112 U3 / \ ( ai al a2 = [Tij 1 a2 a3 / \ 3 / \ 3 r 3 0 ( 2 1/3 - Yi- .1212 Yi- y13 V = V Y 477Ê 47TRi3j24TCRü3 V / \ Z Zi Zji Zi -z Zi -2 10 15 20 25 30 Ru situées KJ.2. y, zi et émettant respectivement un P1=Pi (d'indices i, i K1/3 . 1=Pi de coordonnées xi, yi, zi est : (d'indices respectivement flux de fluide de coordonnées -x1 )2 (yi )2 +(z1 2.11 )2 = \i(xi - xj2 )2 + (yi .1,12 )2 + (zi zi2)2 Pour Ns (2j z13)2 = V jet-source )2 + (yi - y )2 ) (Zi xi, Qobjet-source x, sources S=Sj ponctuelles P ou ab' - points P2 de coordonnées 4 et pour Np points tests yi, zi, la vitesse au point P ou P Qobjet aux ai ou xi, 4ix')1 7):BSIA "S' i A + KSI B "51.8 ± WESis 'S1S - V ; 'SIE ± QTz ESIA 'S I A + QTz ESIR 'SIB 4- QTz nSis 'SIS =VBz Vs?' i sur l'entrée d'air donne 3*ne équations : V \ -SIE 4- QTx ES, A 'S1 A 4- QTx ESIR 'S1B + QTx ESis 'S1S = V Ex \\ S ) ) 'SIE + Q;ESI A -S1 A 4- Q-7çSIB 'SIB ± QiiESis 'SIS = V EY' S +Oz S =V; 'S I E ± QTz ESIA 'SI A 4- QTz ESIB ' 1B ,--TESis ' 1S sur l'entrée d'air donne 3*ns équations : 'SIE ± QTxsSi A 'SI A ± QTA s. SIR "S 18 + QTx sSis v 'S1S ' Sx SIE ± Q-;:SSi A 'SI A + KSI n .SI B ± Q-73"SSis 'SIS - VSv 'SIE + QT,z r.S1 A 'SI A + KSI ' 'SI B + Oss,s.S1s = VSz rM globale de résolution est alors : , Qx ( QTx ES,'. ) OES" ( TES, E ) ( T sSI ( QY ) ( QsYiJS2 ) T ES, E ) Q-73':,, Y.s,, ( TEs,' ( OFst, ) QTz EsIA ( QTESI' ) ( QI; IJSz ) ( QTAs' ) e,,s,' ( QTnAsi' ) ( QT"Asis ) VQx ( Qx T 'Sis ) TBSIE ) KSI il ( TA n ) ( QY ( QY ) ( TESis ) TA E ) Q7'13SI A (T'An ( 0-'s,,. ) Qrz 's,,, ( TBs11, ) ( QTBSIS ) QTx ( ex sS1E ) Q7:x A A ( Q sS ) -7.x AB ) ( TsSIE ) Q*7Y'sS, A ( 0';'Y QY ( QY ---, sSin ) ( TsSis) ( Qz z ) ( Qrz ssis ) \ TsSi E ) QTz A A TSS113 Les conditions aux limites utilisées pour la résolution sont : VB Vs -= Vo. uy VE = VO. Uy Dans cet exemple, les sources sont représentées aux figures 19D à 19M, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source. On a ensuite calculé et visualisé la première grandeur (vitesse V) dans les les 2 plans XOY et YOZ passant par le centre de la sphère, aux figures 19N à 19V. On décrit ci-dessous un mode de réalisation dans le cas d'une source rotationnelle de fluide. Soit une source S ponctuelle rotationnelle au point P2=Pj(xj,yj,zj) qui émet un flux k (X'.,4,,À.,) homogène à un débit massique (kg/s) dans un milieu défini par 15 sa densité p. On calcule au point P ou P1=Pi de coordonnées (xi,yi,z,) le potentiel vecteur A , sa dérivée normale à l'interface I et la vitesse de l'écoulement y. 3 0 respectivement aux points P2 de coordonnées xj, yj, zj et émettant respectivement un Ar A= A Az , , j, j={1,...,N,}) situées Pour Ns sources S=Sj ponctuelles (d'indices 47rpRii 2 - z,) (y, - Y,) (zi - z1) 47rpRii -(Yi- .Y1) 47rpRij. r X 2y -= tif. 11,), 2 2 z \ z 1 47zpRii 1 r 20 25 /1 y =wii. 21' z z = -rot A = 4-77pRii -(xi- x1) 47z-pRij flux de fluide ai ou 4 et pour Np points testsP ou P I =Pi (d'indices i, i ={1,..,Np}) de coordonnées xi, yi, zi, la matrice M, pour différencier les cas sera notée W, T ou Y W= T tu Li Et on considère un vecteur de source de flux = (2x 2Y 2z) P Dans ce mode de réalisation, les conditions aux limites L peuvent comprendre des conditions aux limites intrinsèques (CLI), qui prévoient qu'à l'interface entre deux milieux définis par leur densité respectives (pi, p2), on considère la continuité du potentiel vecteur et sa dérivée calculée par rapport à la normale à l'interface en chaque point de l'interface. On écrit : = A2 aA, a A2 Pl- P2 - - a n an Dans ce mode de réalisation, les conditions aux limites L peuvent comprendre des conditions aux limites utilisateurs (CLU) qui peimettent de fixer des conditions aux limites des objets : vitesse normale nulle et/ou vitesse nulle, et/ou autres. 20 On décrit ci-dessous un mode de réalisation dans le cas d'une source de force. Une source de force est également appelée Stokeslet. Soit une source S ponctuelle de force au point P2=Pj(xj,yi,z,) qui émet un flux I, (Xx,X),,Xz) homogène à un débit massique (kg/s) dans un milieu défini par sa 25 viscosité II et sa viscosité cinématique p. On calcule au point P ou Pl-Pi de coordonnées (x,,y'z,) la vitesse de l'écoulement v et la pression p. R7r,L 1 87rt 0 30 v-1/ 0 0 4 (Yi-yJ). -y1) e 1 ( (2 ), (2 x (xi -xj) (yi - y ) -z ) = Pif v z z P = RuR u 4.7c Pour Ns sources S=Sj ponctuelles (d'indices j, j={1,...,N,}) situées respectivement aux points P2 de coordonnées xi, yi, zi et émettant respectivement un flux de fluide ai ou a.; et pour Np points tests P ou Pl-Pi (d'indices i, i ={1,..,Np{) de coordonnées xi, yi, zi, la matrice M sera notée K pour les stokeslets ou P pour la pression Ku (kis (Pi.1)1 3 4iz F 3 Y - Y3' 3 z - z3' 15 Dans cette dernière équation, les ( a a2 47-ri ; -F3'3 3 composantes longitudinales et transverses 47r 3 Y Y2' 3 47r z - z2' i; sont possibles (mais l'écoulement reste irrotationnel). / est calculée en r pour un triplet de singularités situé en r' 1, r'3 (sources scalaires). Les 3 sources ne doivent pas être confondues. Selon un deuxième exemple, il est prévu un dipôle de sources scalaires 2 0 centré en rl ' et orienté selon l'axe X: x - (x, r+Ax) x - (x '-Ax) 47r i; 471- 4x)3 - Y Y Y,' 4,r Y -(Y'±A43 4,r -(ïr - A.413 Z--Z,' z - z,' .a 4 7 rli; -('rAx)1> 4 21. - 2 2 5 Dans cette dernière équation, les composantes longitudinales et transverses sont possibles (mais l'écoulement reste irrotationnel). / est calculée en r pour un doublet de singularités situé en fl+Ax, r' 1-Ax (sources scalaires). Selon un troisième exemple, il est prévu 3 dipôles centrés en rl r2', r3' et 3 0 orientés respectivement selon X, Y et Z : V(r)=- 1 x -(x,'+ 4x) x - x- x, X - X - X, rai a, 4771? '+4x)r 47r -(7--;'--43 Y-Yi 47r r -(F2'+Ay) 3 - (Y2 4 A))) 47r r -(r,'-4x' Y - (y2 '-AY) 47r -(F3'-Az)3 Y Y - Y3' Y - Y3' 47z1F '+4413 47r -(ï'; z - z,' 47r r -(F2'+Ay)13 z- z,' 4>TF '-Ay) 3 472-1-/). Az) 3 47-/- '-Az) 3 z z,' z - z,' z-(z3'+4z) z-(z3'-3.z) 47t -Ci', '+3.43 42z 4;11 ï; 47c1i2 '-.403 4,r r -(r3'+4x3 4711? -(P3'-Az)i3 Dans cette dernière équation, les composantes longitudinales et transverses 5 sont possibles (mais l'écoulement reste irrotationnel). V est calculée en r pour pour trois doublets de singularités situés respectivement en r" 1 +Ax, r" 1 -Ax, r'2+Ay, r'2-Dy , r'3+Az, r'3-Az (sources scalaires). Si chacun des trois dipôles a un même centre situé en r', on peut enlever les 10 indices 1, 2, et 3 dans la matrice: x- x x -- x' x- x' ra V (7) = x - (x' a +Ax) x-(x'-4x) x-- x' a, rai 47tilF -(T-'+4x13 Y - y' 4,TV -(F'-Ax 3 Y - Y' 47r F -(F'+Ay)3 Y - (Y'+3Y) 47z1F -(P-Ay)13 Y - (1»-AY) 47tili-- -(?'+Az)i3 Y - Y' 47r F -(F'-4z)3 Y - y' 47-c F -(F'+Ax)s z- z' 47r F -(F'-4x)1' z- z' 47-cli; --(F'+Ay)3 z - z' 47r F -(i---"-Ay3 z- z' 4,r F -(/''+4413 z -(z'+Az) 47L17- -(ii»-A4 3 z-(z'-Az.) 421-1i; -(P+Ax)13 4;z-F -(r'-Ax){3 4;zi-(F4-3,43 -(F-Ay)13 47z1T- 47z- r -(r'-Az~ 15 Lorsque la viscosité n'est plus négligeable, aux nombres de Reynolds faibles, la pression et la vitesse ont pour expression : P -7= et V = 1 ( 1 + an 3 .;1) OÙ I est l'identité et a 4.7r.p r - r' 3 r -r r -r 2 0 un point tenseur. d'indice j Doù l'expression point d'indice i, pour une singularité situé au 2 5 rvx` de p calculée au , - i\ x i ).(y, - y, ) (x, - xi, ) (x,- x i).(z,- z ( six\ À iy riz / Vy yz, : 2iy + r (x,-xi) (z,-z V: ) (x,- et le 4.7r.p tenseur 1 ( (x,- R,31 - d'ordre 2 s'exprime pour xj) (x,- xj).(x,- xj) (x, - xi, R + Ru ,, Ry + R; R,, R,' , (Y,- Y,) +(Y,-Y,).(z,- z J) (Y,-Y,) (Y,-Y,).(x,-x,) (Y,-Y,) (Y,-Y,).(Y,-Y,) ± + 8.7r.p Ry R, R, R;, R,, R,i (z,-zi) (z,-zi).(x,-x.,) (z,--zi) (z,-zi).(y,-y,) (z,-z,) (z,-zi).(z,-z,) + + Ru R;: Ry ± R 1 R,, R,' Les conditions aux limites appliquées aux surfaces et contenues dans le vecteur C, peuvent être vectorielles et différentes en chaque point test Pl de la 3 0 surface SUR. Ceci ouvre la possibilité d'immerger un objet I dans un champ de vecteur de topologie irrégulière, d'une part, mais également de conférer à chaque point de l'objet un vecteur vitesse particulier. Ceci a été testé, par exemple avec comme interface I un cylindre en rotation (chaque point test P3 de sa surface SUR est associé à un vecteur vitesse tangentiel) immergé dans un fluide F avec gradient de vent. On décrit ci-dessous d'autres modes de réalisation de la matrice M globale de résolution. Suivant un mode de réalisation, la matrice globale M de résolution comporte plusieurs coefficients, dont au moins un est égal à (x,-xj)/(4.7c. p .Ru3), ou à (y,-yi)/(4. TC . p.R,i3) ou à (zi-z1)/(4.n.p.R,j3), ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur J des paramètres des sources Sj de fluide, où Ru est la distance entre la source Sj de fluide située au troisième point P2 de position prescrite xi, y, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point P de l'espace ayant des coordonnées xl, yi, z, selon les trois 2 0 directions x, y et z. Ce mode de réalisation correspond par exemple à au moins une source scalaire, pouvant être l'un des cas 1, 2 et 3 décrits ci-dessous. Suivant un mode de réalisation, la matrice globale M de résolution comporte plusieurs coefficients, dont au moins un est égal à 1 /(4.n.p.Rij), 2 5 ou à celui-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur J des paramètres des sources S) de fluide, où Ru est la distance entre la source Sj de fluide située au troisième point P2 de position prescrite et le quatrième point P de l'espace. Ce mode de réalisation correspond par exemple au cas 4 décrit ci-dessous. 3 0 Suivant un mode de réalisation, la matrice globale M de résolution comporte plusieurs coefficients, dont au moins un est égal à ±(x,-xj)/(4.7c.p.Ru), ou à ±(yi-yi)/(4.7c.p.Ru) ou à ±(zi-zj)/(4.n. p ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur J des paramètres des sources Sj de fluide, où Ru est la distance entre la source Sj de fluide située au troisième point P2 de position prescrite xj, y, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point P de l'espace ayant des coordonnées x' yi, zi selon les trois directions x, y et z. Ce mode de réalisation correspond par exemple à au moins une source rotationnelle de fluide pouvant être le cas 8 décrit ci-dessous. Suivant un mode de réalisation, la matrice globale M de résolution comporte plusieurs coefficients, dont au moins un est égal à ±(Xi-Xj)/(4.7I.p.R,J3), ou à +(yi-yi)/(4.71.p.Rij3) ou à ±(zi-zi)/(4.7c.p.R1J3), 2 0 ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur J des paramètres des sources Sj de fluide, où Ru est la distance entre la source Sj de fluide située au troisième point P2 de position prescrite xj, yj, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point P de l'espace ayant des coordonnées xi, yi, z, selon les trois 2 5 directions x, y et z. Ce mode de réalisation correspond par exemple à au moins une source scalaire de fluide avec au moins une source rotationnelle de fluide, pouvant être le cas 9 ou 10 décrit ci-dessous. Suivant un mode de réalisation, la matrice globale M de résolution comporte plusieurs coefficients, dont au moins un est égal à l'un parmi 30 (xi-xj)/(8.7c.u.R,j), (yryi)/(8.n.j.i.Ru), (xi-xj)2/(8 u.Ru), (yi-yi)2/(8.7-c.i.t.Ru),(zi-zi)2/(8 (xi-xj).(y,-y,)/(8.n.i.t.R,J3),(x,-xj).(zrzy(8.7c.i.t.Rii3), (zrzi).(yryi)/(8.7c.i.t.Rii3), ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur J des paramètres des sources Sj de fluide, où Ru est la distance entre la source Sj de fluide située au troisième point P2 de position 5 prescrite xJ, y, zi selon trois directions non coplanaires x, y et z de l'espace et le quatrième point P de l'espace ayant des coordonnées x' yi, z, selon les trois directions x, y et z. Ce mode de réalisation correspond par exemple à au moins une source de force, pouvant être l'un des cas 11, 12, 13 décrits ci-dessous. 10 On décrit ci-dessous des sélections et des positionnements particuliers de la source ou des sources. Sources Suivant un mode de réalisation, la source S de fluide est choisie parmi une 15 source S scalaire de fluide, une source S rotationnelle de fluide et une source de force ou toutes combinaisons de ces sources. Suivant un mode de réalisation, la source S de fluide comprend au moins une source S scalaire de fluide et au moins une source S rotationnelle de fluide. Suivant un mode de réalisation, la source S de fluide comprend au moins au 2 0 moins une source S rotationnelle de fluide et/ou au moins une source de force. Suivant un mode de réalisation, la source S de fluide comprend au moins quatre sources choisies parmi au moins une source S scalaire de fluide, au moins une source S rotationnelle de fluide et au moins une source S de force. Suivant un mode de réalisation, la source S scalaire de fluide est une source 2 5 de débit massique isotrope de fluide. Suivant un mode de réalisation, la source S rotationnelle de fluide est une source provoquant une rotation du fluide autour d'une direction déterminée. On peut ainsi prévoir deux sources rotationnelles et une source scalaire, ou quatre sources, comme par exemple trois sources rotationnelles et une source 3 0 scalaire. Suivant un mode de réalisation, l'invention peut également être généralisée jusqu'à six sources (trois sources rotationnelles et trois sources scalaires) avec des conditions aux limites correspondantes. Différents exemples de cette caractéristique de sources mixtes sont décrits ci-dessous. Suivant un mode de réalisation, la présence d'au moins une source rotationnelle permet de calculer la portance et/ou la traînée et/ou de tenir compte d'un écoulement réversible et/ou de tenir compte de pertes d'énergie autour de l'interface I. Dans tous les cas qui suivent, le fluide F est considéré incompressible. Le nombre de Reynolds est le rapport entre les forces d'inerties et les forces visqueuses. On considère des cas de nombres de Reynolds élevés, c'est-à-dire supérieurs à 3000, pour lesquels le fluide est considéré comme non visqueux. Dans ces cas, on néglige les tenues de viscosité. Il y a découplage entre les équations de pression et celles de la vitesse. L'équation de Poisson qui résulte de ces approximations permet de calculer la vitesse et admet pour solutions les fonctions harmoniques (en 1/R et leurs dérivées). La pression est calculée par l'équation de Bemouilli. On considère également d'autres cas de nombres de Reynolds inférieurs à 2000 pour lesquels les forces de viscosité sont prépondérantes. Entre les 2 se situe 2 0 une zone intermédiaire où les 2 phénomènes sont mélangés. Suivant un mode de réalisation, les sources élémentaires sont des sources scalaires (flux de fluide radial) utilisées seules ou en triplets, et/ou des sources vectorielles (flux de fluide transverse, ou rotationnel), dont l'unité est homogène à un débit massique (Kg /s). 25 Cas 1 Suivant un mode de réalisation, représenté à la figure 21, il est prévu comme source de fluide une seule source S scalaire de débit massique de fluide associée à chaque premier point test Pl, pour la condition L aux limites ayant au premier point 3 0 test Pl une composante normale de vitesse du fluide au premier point test Pl , qui est nulle à l'interface I, B. Ce mode de réalisation prévoit donc une seule source par groupe. Chacun des modes de réalisation et des variantes du cas 1 est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S située en dessous de la surface imperméable SUR du solide de l'interface et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. Ce mode de réalisation est par exemple utilisé pour un fluide incompressible, non visqueux (nombre de Reynolds élevé) et non rotationnel, sans tenir compte d'une couche limite, l'objet étant plongé dans un champ Vo de fluide. Dans ce cas, la solution est l'équation de Poisson, la vitesse dérive d'un potentiel scalaire. Ce cas 1 ainsi que les cas comportant uniquement des sources scalaires 10 permet de connaître les valeurs de vitesse ou de pression loin de la surface SUR de l'interface I, pouvant être une aile d'avion, le fluide étant considéré comme parfait. Dans la matrice M globale de résolution, le coefficient est : 4.7c.p.R,J3 15 On a : Les conditions aux limites sont, en projetant sur les normales à la surface B (ou SUR): C Ai = Vi'N B----(QB1 ,A'NB)*/18j = QN1 BA*2.8 Vo.NB La matrice M est : 2 0N avec VBN = -r70ÀTB V ( 2 r " - - - 8,1 le13 VN (-)N,N3 ,13 ) B,n B ) \\ 8,11B Rn c'est-à-dire en écriture condensée : - "0" = La solution complète (V) est la somme de la solution initiale (Vo) et de la 25 solution calculée (Vc). Pour la normale N à la surface SUR, on a: V.NB = (Vo + Vc)..A7, = Vo- .N,3-VO.NB= 0 Des exemples non limitatifs sont : pour l'interface I solide, un cylindre dans une soufflerie selon les figures 19A, 19B et 19C (figure 5A illustrant la valeur 30 absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et figure 5B représentant des lignes de courant de cet écoulement de fluide, ayant été calculées, pour un cylindre CYL ayant un axe parallèle à la direction x), deux cylindres (figure 6A illustrant la valeur absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et figure 6B représentant des lignes de courant de cet écoulement de fluide, ayant été calculée, pour un cylindre CYL1 et un autre cylindre CYL2 distant du cylindre CYL1, qui ont chacun un axe parallèle à la direction x), une sphère ou un profil NACA dans un écoulement, avec prise en compte de l'effet de sol pour illustrer le fait que le procédé permet une analyse des interactions entre objets. Le ou les cylindres sont des obstacles dans l'écoulement de fluide d'une soufflerie. Un autre exemple non limitatif est une aile d'avion selon la figure 36A, pour laquelle la figure 36B illustre la surface imperméable SUR tridimensionnelle de l'interface I, la figure 36C illustre le module de la vitesse de l'écoulement de fluide ayant été calculée en des points P situés dans le plan des axes y et z et la figure 36D représente des lignes de courant de cet écoulement de fluide, ayant été calculées en des points P situés dans le plan des axes y et z. Il a été prévu dans cet exemple un triplet de sources scalaires S' par groupe BL' selon le cas 2 ci-dessous et une source scalaire par groupe BL. La condition aux limites sur la surface SUR a été vent = Vy = 100 m/s (y compris sur le sol) et vitesse normale nulle sur l'aile. La méthode DPSM est semi-analytique. 2 0 Les résultats dans l'espace sont plus ou moins précis en fonction - du nombre de blocs sources utilisés (inversion de la matrice pour obtenir le vecteur condition), - du nombre de points sur lesquels les grandeurs sont calculées (multiplication d'une matrice par le vecteur condition précédent), 2 5 Les figures 37A, 37B, 37C, 37D, 37E, 37F montrent les différences de précision sur le module de la vitesse ayant été calculée en des points P situés dans le plan des axes y et z (figures 37A, 37C, 37E) et sur les lignes de courant de l'écoulement ayant été calculée en des points P situés dans le plan des axes y et z (figures 37B, 37D, 37F) autour de l'aile lorsque l'on change seulement le nombre de points de 3 0 l'espace : - figures 37A, 37B : 10000 points de l'espace, - figures 37C, 37D : 86000 points de l'espace, - figures 37E, 37F : 240000 points de l'espace. Les calculs sont tous fait avec une aile ayant environ 2000 points tests Pl. Bien entendu, les calculs peuvent être effectués dans plusieurs plans, ainsi 5 que cela est représenté par exemple à la figure 37G montrant les lignes de courant de l'écoulement ayant été calculée en des points P situés dans plusieurs plans parallèles aux axes y et z, écartés suivant l'axe x. Le calcul est effectué en trois dimensions, les résultats sont présentés le plus souvent dans un plan pour plus de lisibilité. 10 Cas 2 et 3: Suivant un mode de réalisation, représenté à la figure 22, il est prévu comme source de fluide trois sources S scalaires de débit massique de fluide associées à chaque premier point test Pl, pour la condition L aux limites ayant une composante 15 normale de vitesse du fluide et une composante tangentielle de vitesse du fluide au premier point test Pl, qui sont nulles à l'interface I. Ce mode de réalisation prévoit donc un triplet de sources par groupe. Chacun des modes de réalisation et des variantes des cas 2 et 3 est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S située en dessous de la surface imperméable 2 0 SUR du solide de l'interface et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. Cas 2 Suivant un mode de réalisation, représenté à la figure 23, il est prévu comme 2 5 source de fluide : - une seule source S scalaire de débit massique de fluide associée à chaque premier point test Pl d'une première interface I, pour la condition L aux limites ayant une composante normale de vitesse du fluide au premier point test Pl, qui est nulle à l'interface I, 3 0 - trois sources S' scalaires de débit massique de fluide associées à chaque premier point test P1' d'une deuxième interface I' distincte de la première interface I, pour la condition L aux limites ayant une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide au premier point test P1', qui sont nulles à l'interface I'. Ce mode de réalisation prévoit donc un triplet de sources S'par groupe BL'. Ce mode de réalisation est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A. Par exemple, l'interface I' peut être l'entrée E d'une soufflerie, dans le cas 10 de l'interface I se trouvant dans un écoulement de fluide provoqué par une soufflerie. Bien entendu, cet exemple n'est pas limitatif. Par exemple : Ru - la matrice M comporte les coefficients 4.7r.p.R,J3 - le vecteur vitesse est : V = Q11./1, 15 - Les conditions aux limites sont, en projetant sur les normales à la surface SUR: C Al =Vol = ei 1,A./1 0 + -- A1 Ai 2, A'i 1 A2+ Q Ai 3,A-2 A3 B,A'''''B i - la matrice M est, en écriture condensée : (7 x x ry r') Vox n x AA1 QAA2 nAA3 !AB /''Al V - Oz Z vAA1 Z V AA 2 AA3 ilz QnAY! 2A2 VO) ()Y AA1 Q AY A2 QAYA3 - VN n N n N /IN AB 2 A3 ,--)N B ° ) \BA1 BA2 BA3 BB i\ÀB 2 0 c'est-à-dire Vo rn AA QBA /L'A \V: -° \QANB QBNB )\i1B ) - la dimension de la matrice est : 3nA x 3nA 3nA x nB, ) on + n ) dim = ( -(311A+n8 x A e nsx3nA nsxns) 2 5 Les interfaces I et I' peuvent être des objets imperméables au fluide, qui sont en interaction, comme par exemple un avion I à l'approche du sol I'. Il peut en outre être prévu des conditions aux limites vectorielles (indiquées ci-dessus pour les trois sources scalaires) sur les interfaces (par exemple pour un gradient de vent ou autre). Cas 3 Suivant un mode de réalisation, représenté à la figure 24, il est prévu comme source de fluide : - trois sources S scalaires de débit massique de fluide associées à chaque premier point test Pl d'une première interface I, pour la condition L aux limites ayant une composante normale de vitesse du fluide au premier point test P I , qui est 10 nulle à l'interface I, - trois sources S' scalaires de débit massique de fluide associées à chaque premier point test P1' d'une deuxième interface I' distincte de la première interface I, pour la condition L aux limites ayant une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide au premier point test Pl ', qui 15 sont nulles à l'interface I'. Ce mode de réalisation prévoit donc un triplet de sources S par groupe BL et un triplet de sources S'par groupe BL'. Ce mode de réalisation est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas 20 ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. Dans ce mode de réalisation, il pourrait être prévu en variantes une seule source S scalaire de débit massique de fluide associée à chaque premier point test P 1 d'une première interface I, ou deux scalaires de débit massique de fluide associées à chaque premier point test Pl d'une première interface I, avec 25 une seule source S' scalaire de débit massique de fluide associée à chaque premier point test P1' d'une deuxième interface I' ou deux sources S' scalaires de débit massique de fluide associées à chaque premier point test Pl d'une deuxième interface I'. Par exemple, l'interface I' peut être l'entrée E d'une soufflerie, dans le cas 30 de l'interface I se trouvant dans un écoulement de fluide provoqué par une soufflerie. Bien entendu, cet exemple n'est pas limitatif. Par exemple, dans le cas de 3 sources scalaires S et 3 sources scalaires S' Ry - la matrice M comporte les coefficients Q' = 4.7-c.p.Ry3 le vecteur vitesse est : V = Qu - Les conditions aux limites sont, en projetant sur les normales à la surface SUR: CA, =V. V, = Q +A1 ,,A Ai 2,A.2A2 A /1 A' A3 BI,A' B2.A- 82 B3 , A' B3l Q +Ql +Q CA =VO C = 0 - la matrice M est, en écriture condensée : vox VY voz V; = 0 y; = 0 T7,; = 0) 1n n x x ry nx AA1 AA 2 ecil 3 QA x 81 QAB2 A83 Al - - - (y) 2A, (z) - - - 2A3 y BAI r 8,42 Qj x Be A3 Q BB1 QBB2 883 281 - - - (y) 2B2 (z) - - - )\2,63 c'est-à-dire r ( n n VO vAA V BA l'A -.=. j/ - B = Q \Q AB QBB i 2B i - la dimension de la matrice est : r3nAx3nA 3nAx3nB = (3nA+3n8)x(3nA+3nB) \3nB x 3nA 3nB x 3nB Les modes de réalisation des cas 2 et 3 sont par exemple utilisés pour un fluide incompressible, non visqueux (nombre de Reynolds élevé), non rotationnel, avec restitution du phénomène de couche limite CL par annulation de la composante tangentielle de vitesse à l'interface I. Des exemples sont : cylindre dans un 2 0 écoulement avec effet de sol, et zoom au voisinage des parois des objets pour illustrer le fait que le procédé permet des conditions aux limites complexes en chacun des points de la surface des objets. Pour cet exemple dans le cas 3, la figure 7A représente la valeur absolue de la composante Vy de la vitesse de l'écoulement du fluide ayant été calculée selon l'axe y et la figure 7B représente des lignes de 2 5 courant de cet écoulement de fluide, ayant été calculées, et ce pour une interface I dim15 solide, formée par un cylindre CYL ayant un axe parallèle à la direction x. Les figures 8, 9 et 10 illustrent des courbes Cl, C2, C3, C4 de la vitesse Vy calculée de l'écoulement de fluide F selon l'axe y sur un axe vertical, en fonction de l'altitude z en abscisses, pour respectivement les conditions aux limites suivantes du cylindre CYL de la figure 7 et du sol : - Cl : Vsoin = O et Vcyin = 0, - C2 : Vsoi =0 et Veyin = 0, C3 : Vns01=0 et Veyi = 0, - C4 : Vsoi =0 et Vey1= 0, où l'exposant n désigne la composante normale de la vitesse V et l'absence de l'exposant n désigne le vecteur vitesse. Suivant un mode de réalisation, chaque groupe BL contient trois sources S scalaires de débit massique de fluide, pour les conditions L aux limites ayant une composante nonnale de vitesse du fluide, qui est nulle à l'interface I et une composante tangentielle de vitesse du fluide, qui est imposée. La figure 11A représente la valeur absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et la figure 11B représente des lignes de courant de cet écoulement de fluide, ayant été calculées, et ce pour une interface I 2 0 solide, formé par un cylindre CYL ayant un axe parallèle à la direction x, où les conditions aux limites prévoient que la composante tangentielle du vecteur vitesse de l'écoulement de fluide est égale à une valeur prescrite du vent (par exemple ici 100 m.s-1) sur le cylindre. 2 5 Cas 4 Suivant un mode de réalisation, représenté à la figure 25, il est prévu comme source S de fluide au moins une première source S scalaire de débit massique de fluide en association à chaque premier point test Pl de l'interface I, au moins une deuxième source S20 scalaire de débit massique de fluide située au-dessus d'une 30 surface SCL délimitant une couche limite CL de l'écoulement de fluide F par rapport à l'interface I et au moins une troisième source S30 scalaire de débit massique de fluide située en dessous de la surface SCL délimitant la couche limite CL, pour la condition L aux limites ayant au premier point test P 1 de l'interface I une composante normale de vitesse du fluide, nulle à l'interface I, et une composante nottnale de vitesse conservée en un autre point test prescrit P5 de la surface SCL délimitant la couche limite CL, modulo la masse volumique du fluide. La au moins une première source S scalaire de débit massique de fluide est donc située dans la première zone intérieure Z 1 . La au moins une deuxième source S20 scalaire de débit massique de fluide est donc située dans la deuxième zone extérieure Z2. Il peut être prévu une seule source scalaire S de débit massique de fluide en association à chaque premier point test P1 de l'interface I ou trois sources scalaires S de débit massique de fluide en association à chaque premier point test Pl de l'interface I. Ce mode de réalisation prévoit donc une source S ou un triplet de sources S par groupe BL et deux sources S20, S30 par groupe BL". L'interface I est par exemple un objet B. Ce mode de réalisation est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. C 2 0 désigne la surface SCL. Ce cas apparait par exemple pour un nombre de Reynolds élevé en fluide irrotationnel, des milieux M2 fluides stratifiés dans lesquels sont immergés des objets solides. Les conditions aux limites sont par exemple des conditions aux limites utilisateurs (CLU) aux surfaces imperméables et des conditions aux limites 2 5 intrinsèques (CLI) à la surface SCL fluide-fluide. Il peut être prévu un objet isolé dans un écoulement irrotationnel, entouré de couches de différentes propriétés physiques, avec deux conditions aux limites par point test de la surface SCL fluide/fluide (composante normale de la vitesse conservée modulo rho et potentiel scalaire) (CLI) et les mêmes conditions aux limites que précédemment à la surface 3 0 SUR imperméables. Le coefficient de la matrice M en potentiel est par exemple : M = 1 47zpRii avec 6= MY .a On impose, comme conditions limites intrinsèques - l'égalité du potentiel de part et d'autre de la surface SCL : C A = 84,Top - A,Bottonz = M1:14,A '2 A - MA'2 B = - La vitesse est tangente à la surface SCL, c'est-à-dire que la vitesse normale est nulle Dans une forme condensée de la matrice : vx ( QAY 42 gAA3 0 VY (Y) 0 vz (z) 0 Pi -P2.vcAl2 =0 P2e A P2ee P2 Rec, 6i -6'2 =0 MCI Mc4 Mc4 MCCjbp v: =0 r' 'A VACB', 0 ' \ r 2 - \ /Ll (Y) 0 2A2 (z) 0 2A3 - P2' .', Pl' /1cTop M CCB', MCB /140, 0 QIBVB \a13 i ) VO r Q AC 'o, 0 \(2A\ Pl*V N1 P2 J/CN QAA 1V A C 2 = 0 /1 N 81-82=0 nN P1-VCC 'o, P2-/ICB MCB VB = 0 P1*.CA P2-CCrop- M CA M CCrop CC Bo, ./1.'C Bot 0 BN C7' 0 QBB 2B ) La dimension de la matrice est : ( 3n Ax3n A 3n Ax ne 3n Ax ne 3nA X nB 2 0 dim = r i c x3n A ncx ne ne x ne nC x nB = (3n A + 2nc + n B)x (3n A+ 2n c + n B) r i c x3n A ncx ne ne xn c nC xnB \, nB x3nA n B X nC nB X nC nB X nB J Cas 5 Suivant un mode de réalisation, représenté à la figure 26, il est prévu comme 2 5 source S de fluide au moins deux sources S scalaires de débit massique de fluide en association à chaque premier point test P1 de l'interface I, pour la condition L aux limites ayant au premier point test Pl de l'interface I une composante normale de vitesse du fluide, nulle à l'interface I, et une composante normale de vitesse nulle en un autre point test P5 prescrit d'une surface SCL délimitant une couche limite CL 30 de l'écoulement de fluide par rapport à l'interface I. Ce mode de réalisation prévoit donc deux source S par groupe BL et deux sources S20, S30 par groupe BL". En variante, il pourrait être prévu une seule source scalaire S de débit massique de fluide en association à chaque premier point test Pl de l'interface I ou trois sources scalaires S de débit massique de fluide en association à chaque premier point test Pl de l'interface I. L'interface I est par exemple un objet B. Chacun des modes de 5 réalisation et des variantes du cas 5 est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. C désigne la 10 surface SCL. Ce cas apparait par exemple pour un nombre de Reynolds élevé, pour un objet isolé dans un écoulement irrotationnel, entouré de couches de même propriétés physiques. On fixe des conditions aux limites utilisateur (CLU) sur la surface SCL fluide/fluide (par exemple composante normale de la vitesse nulle). Ce cas fait 15 intervenir des Conditions aux Limites Stratifiées » (CLS) : il n'y a pas besoin de blocs sources de part et d'autres des points test de la surface SCL, les sources situées dans les objets environnants I sont calculées pour prendre ces CLU en compte. Il y a les mêmes conditions que précédemment sur les surfaces imperméables : les CLU s'appliquent. Il peut être des sources complémentaires pour 2 0 prendre en compte les CLS qu'on vient de rajouter. Le cas de deux sources scalaires S est décrit ci-dessous. La matrice M est : I-'\ V x = ( nx Q Ax A2 n x Q AxBI Q Ax B2 \ A V AAl CA 2 QBA2 QAA3 i.--) y QCN B2 ( 2 Al V Y QC AI CAI (y) CB1 QBN B2 2A2 2A3 2.81 \,2 B2 y A \.,nN (Z) QBN BI V.: BA1 (I N VcAr = 0 CA3 V = B nN QBA3 c'est-à-dire : Vo QN N N 1 Q n \ 1 \ Vc =0 ,n nN -,1/1 AB1 QAB2 A V B N = 0 ,A1 Q. QCB2 2BI BA1 .881 QBB2 y X282 La dimension de la matrice est : dim (3nAx3nA 3n xn 3n xn = (3n A +2n8)x (3n A +2118) n8x3nA A B A B n8 x 3nA n xn n xn B B B B n xn n xn B B B B La surface SCL n'est pas une surface de séparation, elle porte juste les points tests P5. Les sources S' de A envoient du fluide à travers B. Cas 6 Suivant un mode de réalisation, représenté aux figures 27A, 27B et 27C, au cours d'une première itération des première et deuxième étapes El, E2, on prescrit une première condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide F prise au premier point test prédéterminé Pl de l'interface I et un premier nombre déterminé de sources 5100 de fluide, supérieur ou égal à un, pour calculer par méthode DPSM une première itération de la première grandeur physique (V, P) de l'écoulement de fluide F située au deuxième point P de l'espace à partir du au moins un paramètre du premier nombre déterminé de sources S de fluide, ayant été calculé, puis au cours d'au moins une deuxième itération des première et deuxième étapes El, E2, postérieure à la première itération, on prescrit une deuxième condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au premier point test prédéterminé (P1) de l'interface (I), la deuxième condition aux limites étant différente de la première 2 5 condition aux limites, on ajoute au premier nombre déterminé de sources S100 de fluide un deuxième nombre déterminé de sources S200 de fluide, supérieur ou égal à un, pour calculer par méthode DPSM au moins une deuxième itération de la première grandeur physique (V, P) de l'écoulement de fluide F située au deuxième point P de l'espace à partir du au moins un paramètre des premier et deuxième nombres déterminés de sources S100, 5200 de fluide, ayant été calculé. Le calculateur comporte un cinquième moyen CALS de calcul configuré pour calculer la première itération, un sixième moyen CAL6 de calcul configuré 5 pour calculer la deuxième itération. Un exemple de première itération est représenté à la figure 27A (premier condition aux limites étant vitesse normale du fluide nulle, vitesse tangentielle du fluide non nulle). Un exemple de deuxième itération est représenté à la figure 27B (deuxième condition aux limites étant vitesse tangentielle du fluide nulle). Bien 10 entendu, une ou plusieurs autres itérations peuvent être effectuées, comme par exemple une troisième itération représentée à la figure 27C (par exemple circulation ou rotation de la surface SUR). Ce cas est appelé calcul séquentiel intégré. Dans ce calcul séquentiel intégré, les résultats obtenus dans l'étape n-1 sont intégralement intégrés dans l'étape n. 15 Ainsi, suivant un mode de réalisation, au cours d'une première itération de la deuxième étape E2 de calcul, chaque bloc BL contient une première source S scalaire de débit de fluide pour les conditions L aux limites ayant une composante normale de vitesse du fluide, qui est nulle à l'interface I, et on met en oeuvre la méthode DPSM pour calculer le débit de fluide de la première source S scalaire à 2 0 partir des conditions L aux limites. Puis, au cours d'une deuxième itération de la deuxième étape E2 de calcul, chaque bloc BL contient la première source S scalaire de débit massique de fluide, une deuxième source S de débit massique de fluide et deux sources S rotationnelles de débit massique de fluide pour les conditions L aux limites ayant une composante normale de vitesse du fluide, qui est nulle à l'interface 2 5 I, et une composante tangentielle de vitesse du fluide, qui est prescrite. Ce mode de réalisation est par exemple utilisé pour un fluide incompressible, non visqueux (nombre de Reynolds élevé) avec prise en compte de la vorticité par itérations successives. Cela permet par exemple de tenir compte d'un objet en rotation ou pas, pouvant être une sphère immobile ou en rotation dans un écoulement avec effet de 3 0 sol. Cette seconde méthode peut être couplée avec un calcul de la couche limite à l'issue de la première étape. Par exemple, une surface de glissement SG représentative de la couche limite entourant l'objet est définie géométriquement, et maillée de points tests. Cette surface peut être soit une surface de glissement issue du premier calcul, soit la surface d'une couche limite calculée avec un modèle adapté, par exemple un modèle de Prandtl, qui définit une épaisseur de couche limite égale à 8(l) = ,/ (où / est la distance parcourue selon une ligne de courant R, par rapport au point d'arrêt amont). Dans un mode de réalisation de l'invention, les blocs de sources BL disposés à l'intérieur du profil sont agencés pour envoyer du fluide simultanément dans la couche limite et dans le milieu extérieur à la couche limite selon des conditions aux limites fixées conjointement aux points test de l'interface I et de la surface de glissement SG. Les sources envoient du fluide dans ces deux milieux, supposés identiques de chaque côté de la couche limite. Ce cas est illustré dans le cas de 3 sources scalaires et 1 source rotationnelle. Dans une seconde réalisation de l'invention, on prend en compte la viscosité dans la couche limite, et la surface de glissement SG devient une interface séparant 2 milieux de propriétés différentes : proche de la surface avec nombre de Reynolds faible, au-delà de cette surface avec nombre de Reynolds élevé. Ce cas sera traité 2 0 plus loin. Ce mode de réalisation permet par exemple une superposition des solutions scalaires et vectorielles par itérations successives : intégration des résultats de l'itération n dans l'itération n+1. Dans ce cas, on calcule d'abord une solution scalaire à une source par bloc dans le profil, donc on obtient la forme générale des 2 5 lignes de l'écoulement. Puis on intègre ce résultat dans le calcul suivant qui, lui, introduit d'autres conditions sur l'objet (annulation de la composante tangentielle, objet en rotation, etc..). Par exemple, la première itération donne la première matrice M suivante : 30 Al\ BB Q BNA 2 Q BNA 3 BA 3 nn QAXB QAXAI QAXAI QAAI Al Q B Q A 2 Q Ay A 2 Q A- Y A 2 A2 Z Z Z QAB QAA3 QAA3 QZ A \, A3/\11A3 46 c'est-à-dire V N = 0v r,Ar nN 3 \ \ VO ) QBB BA A AB Q AA )\ dont la dimension est : dim r nB X nB nB x 3nA = \3nAX nB 3nA x 3nA 7 = (nB + 3n A)x (n B + 3nA) A la deuxième itération, on augmente par exemple la première matrice M, en ajoutant les 3 sources de B et en changeant les conditions aux limites : on annule la vitesse tangentielle calculée à la première itération, pour obtenir la deuxième 10 matrice M suivante : /1-1B ( Id 0 0 r/1.1 VT QBB Q88123 QBT A Â2B B - -VBretape' AB Q AB123 Q A71 ,7\,2 A ) VO ) dont la dimension est : 15 n B X n8 11B X3n B 11B X 3nA dim - 3n B x n B 3n B x 3n B 3n B x3n A = (4% + 3n A)x (4n B + 3nA) 371.AxnB 3nA x 3n B 3nA x 3nA 20 {Première itération : il, B Deuxième itération : 2,8 Troisième itération : 238 A la troisième itération, on obtient la troisième matrice M suivante : MO 0000 0000Vil1B 0 Id 0 0 0 0 0 0 0 Ci ç12/31 00/d00 0 0000/12B2 000/d0 0 0000 ,Â2B3 (B B BB (Id 082 083 «A Q712 QBAI (iI3B1 (Y) 113)32 (z) - - ° AB3 \ QAB «BI «Le Qiic BI «B2 aB3 0,4 QvA,2 el/13 (Y) - - - 24 (z) A3i) I 3 ( Id 0 0 0`0.18` '18 228 0 Id 0 0 228 VB Vrotation Q BBI Q882 Q883 Q BA 23B VO 7 Q ABI Q AB2 Q AB3 Q AA ) \il A 7 c'est-à-dire dont la dimension est : n x3n 8 3nBx3nB 3nBx3nB 3nAx3n B nB x3nA 3nBx3nA 3nBx3nA 3nAx3nA) =(7nB+3nA)x(7nB+3nA) dim = ( nBxnB nBx3n 8 3nBxnB 3nBx3nB 3nBxnB 3n8x3nB 3nAx118 3nAx3n B La vitesse en un point M est : aB V = M81 QMB2 QM83 Q MA , 238 ) Cas 7 Suivant un mode de réalisation, représenté aux figures 28A, 28B et 28C, au cours d'une première itération des première et deuxième étapes El, E2, on prescrit une première condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide F prise au premier point test prédéterminé Pl de l'interface I et un premier nombre déterminé de sources 5100 de fluide, supérieur ou 2 0 égal à un, pour calculer par méthode DPSM une première itération de la première grandeur physique (V, P) de l'écoulement de fluide F située au deuxième point P de l'espace à partir du au moins un paramètre du premier nombre déterminé de sources S100 de fluide, ayant été calculé, puis au cours d'au moins une deuxième itération des première et deuxième 25 étapes El, E2, postérieure à la première itération, on prescrit une deuxième condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide F prise au premier point test prédéterminé Pl de l'interface I et un deuxième nombre déterminé de sources 5200 de fluide, supérieur ou égal à un, pour calculer par méthode DPSM une deuxième itération de la première grandeur 30 physique (V, P) de l'écoulement de fluide F située au deuxième point P de l'espace à partir du au moins un paramètre du deuxième nombre déterminé de sources 5200 de fluide, ayant été calculé, la deuxième condition aux limites étant différente de la première condition aux limites, puis on additionne la première itération de la première grandeur physique (V, P) de l'écoulement de fluide F et la deuxième itération de la première grandeur 5 physique (V, P) de l'écoulement de fluide F. Le calculateur comporte un cinquième moyen CALS de calcul configuré pour calculer la première itération, un sixième moyen CAL6 de calcul configuré pour calculer la deuxième itération et un septième moyen CAL7 de calcul configuré pour additionner la première itération et la deuxième itération. 10 Un exemple de première itération est représenté à la figure 28A (première condition aux limites étant vitesse normale du fluide nulle, vitesse tangentielle du fluide calculée). Un exemple de deuxième itération et d'addition de la première itération et de la deuxième itération est représenté à la figure 28B (deuxième condition aux limites étant vitesse tangentielle imposée, par l'utilisateur ou 15 dépendant du premier calcul). Bien entendu, une ou plusieurs autres itérations peuvent être effectuées, comme par exemple une troisième itération représentée à la figure 28C (par exemple circulation ou rotation de la surface SUR). Ce cas est appelé calcul séquentiel sommé. Dans ce calcul séquentiel sommé, les résultats obtenus dans l'étape n-1 peuvent être partiellement injectés 2 0 dans l'étape n sous forme de conditions aux limites. Ce mode de réalisation permet par exemple une superposition des solutions scalaires et vectorielles par itérations successives : sommation des résultats des différentes itérations pour l'obtention du résultat final. Ce cas peut être identique au précédent, mais chaque solution est élaborée successivement. Ceci conduit à la 2 5 notion de sources additionnelles. Par exemple : on calcule d'abord une solution scalaire à une source par bloc dans le profil, puis on rajoute des sources hybrides pour annuler la composante tangentielle de la vitesse sur le profil : il en résulte des jeux de sources additionnelles pour tous les objets en interaction dans notre problème modélisé. Puis on impose une rotation à l'objet, etc. La solution sera la 3 0 résultante de chacune des solutions élémentaires obtenues à chaque étape. Ce cas permet de calculer la géométrie d'une couche limite CL et de la surface SCL la délimitant, à l'issue de la première itération et de la prendre en compte au cours de la deuxième itération. Ainsi, par exemple, on impose des conditions aux limites sur l'aile à la 5 première itération, puis on impose de nouvelles conditions aux limites à la deuxième itération. Les nouvelles conditions aux limites dépendent de la première itération ou sont intrinsèques, comme dans le cas de la recirculation constante. Par exemple, le résultat de la première itération pour une seule source scalaire S par bloc BL ou une seule source scalaire S associée à chaque premier 10 point test P 1 (cas 1 ci-dessous) a donné un vecteur vitesse à la surface SUR de l'interface I (en l'occurrence une aile I, SUR d'avion) selon la figure 38A (vecteur vitesse allant du bord avant BAV d'attaque de l'aile I, SUR au bord arrière de fuite BAR de l'aile I, SUR à la fois du côté de l'extrados EXTR de celle-ci et de l'intrados INTR de celle-ci). Lors de cette première itération, on impose comme 15 conditions aux limites la nullité de la composante normale de la vitesse de l'écoulement de fluide à la surface SUR de l'aile I, SUR. Lors de la deuxième itération, on peut imposer à la surface SUR de l'aile I, SUR des conditions aux limites étant une vitesse opposée à la vitesse calculée à la première itération selon l'exemple de la figure 38B (vecteur vitesse allant du bord 20 arrière de fuite BAR de l'aile I, SUR au bord avant BAV d'attaque de l'aile I, SUR à la fois du côté de l'extrados EXTR de celle-ci et de l'intrados INTR de celle-ci). Ainsi, en additionnant la première itération et la deuxième itération, on obtient une vitesse nulle de l'écoulement de fluide à la surface SUR de l'aile I, SUR. Ou, lors de la deuxième itération, on peut imposer à la surface SUR de l'aile 25 I, SUR des conditions aux limites étant une vitesse constante tangente selon l'exemple de la figure 38C (recirculation constante autour de l'aile I, SUR, à savoir vecteur vitesse allant du bord avant BAV d'attaque de l'aile I, SUR au bord arrière de fuite BAR de l'aile I, SUR du côté de l'extrados EXTR de celle-ci, et vecteur vitesse allant du bord arrière de fuite BAR de l'aile I, SUR au bord avant BAV 30 d'attaque de l'aile I, SUR du côté de l'intrados INTR de celle-ci). Un premier exemple des figures 38A et 38B est illustré ci-dessous. La figure 39 représente un exemple de vecteur conditions aux limites sur l'aile I, SUR. Les figures 40A et 41A représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à la première itération selon le cas 1 d'une source scalaire S par bloc BL ou une source scalaire S associée à chaque premier point test Pl. Ces figures correspondent à la figure 38A et ont été établies pour l'aile I, SUR de la figure 39. Les figures 40B et 41B représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à la deuxième itération selon le cas de trois sources scalaires S par bloc BL ou trois sources scalaires S associées à chaque premier point test P 1. Les conditions aux limites de cette deuxième itération comprennent une vitesse à la surface SUR, opposée à la vitesse calculée à la première itération. En outre, les conditions aux limites de cette deuxième itération comprennent une vitesse nulle à l'entrée et à la sortie de la soufflerie et au sol. Ces figures correspondent à la figure 38B. Les figures 40C et 41C représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération des figures 40B et 41B et de respectivement la première itération des figures 40A et 41A. On remarque que les vitesses à la deuxième itération sont beaucoup plus petites qu' à la première itération. Les figures 42A, 42B, 42C, 42D, 42E, 42F, 42G et 42H représentent les sources de la première itération des figures 40A et 41A, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de 25 la source. Les figures 43A, 43B, 43C, 43D, 43E, 43F, 43G et 43H représentent les sources de la deuxième itération des figures 40B et 41B, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source. On remarque que l'influence des sources de la soufflerie diminue, alors 30 que l'influence des sources de l'aile I, SUR croît. Un deuxième exemple des figures 38A et 38B est illustré ci-dessous. La figure 39 représente un exemple de vecteur conditions aux limites sur l'aile I, SUR. Les figures 44A et 45A représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à la première itération selon le cas 1 d'une source scalaire S par bloc BL ou une source scalaire S associée à chaque premier point test P1. Ces figures correspondent à la figure 38A et ont été établies pour l'aile I, SUR de la figure 39. Les figures 44B et 45B représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à la deuxième itération selon le cas d'une source scalaire S et de deux sources rotationnelles SR par bloc BL ou d'une source scalaire S et de deux sources rotationnelles SR associées à chaque premier point test Pl (les deux sources rotationnelles SR étant par exemple une source rotationnelle autour d'une direction parallèle à l'axe x et une source rotationnelle autour d'une autre direction parallèle à l'axe y). On appelle cela une source S, SR dite hybride. Les conditions aux limites de cette deuxième itération comprennent une recirculation constante autour de l'aile I, SUR selon la figure 38C, par exemple avec une vitesse tangente de 20 m/s. Les figures 44C et 45C représentent respectivement les lignes de champ de 20 l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération des figures 44B et 44B et de respectivement la première itération des figures 44A et 45A. Les figures 46A, 46B, 46C, 46D, 46E, 46F, 46G et 46H représentent les sources de la première itération des figures 44A et 45A, où l'abscisse est l'indice du 25 point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source. Les figures 47A, 47B, 47C, 47D, 47E, 47F, 47G et 47H représentent les sources de la deuxième itération des figures 44B et 45B, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux 30 de la source. Les figures 48A et 49A représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération et de la première itération dans le deuxième exemple précité des figures 38A et 38B pour un angle d'incidence de l'aile I, SUR de 0° et à une première distance déterminée, éloignée du sol. Les figures 48B et 49B représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération et de la première itération dans le deuxième exemple précité des figures 38A et 38B pour un angle d'incidence de l'aile I, SUR de 10° et à la première distance déterminée de l'aile I, SUR, éloignée du sol. Les figures 48C et 49C représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération et de la première itération dans le deuxième exemple précité des figures 38A et 38B pour un angle d'incidence de l'aile de 0° et à une deuxième distance déterminée de l'aile I, SUR, proche du sol. La deuxième distance déterminée est inférieure à la première distance déteiminée. Les figures 48C et 49C représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération et de la première itération dans le deuxième exemple précité des figures 38A et 38B pour un angle d'incidence de l'aile I, SUR de 0° et à une deuxième distance déterminée de l'aile I, SUR, proche du sol. Les figures 48D et 49D représentent respectivement les lignes de champ de l'écoulement et le module de la vitesse de l'écoulement et correspondent à l'addition de la deuxième itération et de la première itération dans le deuxième exemple précité des figures 38A et 38B pour un angle d'incidence de l'aile I, SUR de 10° et à la deuxième distance déteiniinée de l'aile I, SUR, proche du sol. Ci-dessus, les figures 36A à 49D correspondent à une aile de type NACA. Dans une variante, lors de la deuxième itération, la vitesse imposée sur la surface SUR de l'aile I, SUR est calculée à partir de la première itération. Bien entendu, ce qui est décrit pour l'aile I, SUR, pour laquelle le fluide est l'air, est valable pour toute surface SUR, toute interface I et tout fluide. Aux figures 41B, 45B, la couleur noire ou gris foncée correspond à une valeur située dans le bas de l'échelle de gris représentée sur la droite. Aux autres figures, Vyl, Vy2, Vy3, Vy4, Vy5 désignent des domaines de valeurs de la vitesse Vy, indiqués dans l'échelle de gris représentée sur la droite. Par exemple, le déroulement de ce mode de réalisation peut être le suivant : Etape 1 (première itération): on génère la matrice pour une source scalaire S100. (source : 1) Conditions aux limites CL 1: vitesse sur A = vitesse soufflerie, vitesse notinale sur B =0. Le calcul nous donne une vitesse tangentielle sur B Etape 2 (deuxième itération): On superpose des conditions aux limites (nécessite de nouvelles sources 5200 à la place des sources SI00), Conditions aux limites CL2: Vitesse sur A nulle Vitesse sur B opposée à la vitesse calculée à l'étapel précédemment. En sommant les 2 étapes, on obtient Vitesse = VitesseEtapel + Vitesse Etape2 VA = VA Etapel 2 A Etape2 - V soelerie et VB = VBEtapel V BEtape2 ° Etape 3: On superpose de nouvelles conditions aux limites sur l'objet (par exemple : circulation autour d'un profil ou objet en rotation donné ici pour 2 5 exemple), en prenant de nouvelles sources 5300. Par exemple, la matrice M à la première itération est : N (N /a B = 0 n BB ïBA '1B VO 2 &AB QAA )Y11 A ) de dimension : dim _ ( 118 x n B n x 3nA = (n B +3n A)x (n B +3n A) \AX nB 3nA X 3nA ) Par exemple, la matrice M à la deuxième itération est : V A= 0 (VT -VT ( T T \( 2 \QAB2 QAA iY12A B Betapel .VBB2 BA "2B de dimension : dim z3nBx3nB 3nBx3nA\ (3nB+3nA)x(3nB +3nA) 3nAx3nB 3nAx3nA) Par exemple, la matrice M à la troisième itération est (V -V ( v \ B - rotation BB3 QBA '3B V A O ) AB3 QAA A3 A y de dimension : dim = (3nB x 3nB 3nBx3nA = (3nB + 3nA) x (3n8 +3n A) 3nAx3nB 3nAx3nA En additionnant les itérations, la vitesse en un point M est : ( ( \ VA/I =(QMB1 QMA ( +(QMB2 QMA (QMB3 QMA B B 10 / \22A \3A Suivant un mode de réalisation dans les cas 6 et 7, la résolution est effectuée par au moins deux itérations successives ou plus. Par exemple, on résout le problème avec des sources scalaires dans la première itération, puis on ajoute des sources rotationnelles dans la deuxième itération pour satisfaire éventuellement le 15 cas d'un fluide rotationnel. Un exemple est celui d'une sphère dans un écoulement, qui peut être modélisé en DPSM par des sources scalaires, mais si cette sphère est en rotation (effet Magnus), il est nécessaire d'ajouter à la solution précédente une solution rotationnelle (vortex). Chacun des modes de réalisation et variante des cas 6 et 7 est par exemple 2 0 combiné avec le premier mode de réalisation sous-jacent précité (source S et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. C désigne la surface SCL. 25 Cas 8 Suivant un mode de réalisation, représenté à la figure 29, il est prévu comme source de fluide au moins une source SR rotationnelle de fluide en association à chaque premier point test Pl, pour les conditions L aux limites ayant une composante normale de vitesse du fluide au premier point test Pl, qui est nulle à l'interface I, et une composante tangentielle de vitesse du fluide au premier point test Pl, qui est prescrite. Dans ce mode de réalisation, il peut être prévu une seule source SR rotationnelle de fluide en association à chaque premier point test P1 de l'interface I ou deux sources SR rotationnelles de fluide en association à chaque premier point test P1 de l'interface I ou trois sources SR rotationnelles de fluide en association à chaque premier point test Pl de l'interface I. Ce mode de réalisation prévoit donc 1, 2 ou 3 sources rotationnelles SR par groupe BL. L'interface I est par exemple un objet B. Chacun des modes de réalisation et variante du cas 8 est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source SR et/ou S' située en dessous de la surface imperméable SUR du solide de l'interface I et/ou I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. C désigne la surface SCL. 2 0 Ce cas est par exemple utilisé pour un nombre de Reynolds élevé en fluide rotationnel. Ce cas s'applique par exemple à des conditions aux limites introduisant ou nécessitant de la vorticité (objet en mouvement, en rotation, ou dans un écoulement quelconque). Les géométries traitées sont identiques aux cas précédents. En variante à ce mode de réalisation, il est prévu comme source de fluide : 2 5 - au moins une sources SR rotationnelle de fluide associée à chaque premier point test P1 d'une première interface I, pour les conditions L aux limites ayant une composante normale de vitesse du fluide au premier point test P1, qui est nulle à l'interface I, et une composante tangentielle de vitesse du fluide au premier point test P 1, qui est prescrite, 3 0 - au moins une source S' scalaire de débit massique de fluide associée à chaque premier point test P1' d'une deuxième interface I' distincte de la première interface I. Il pourrait être prévu en variantes une seule source SR rotationnelle de débit massique de fluide associée à chaque premier point test Pl d'une première interface I, ou deux sources SR rotationnelles de fluide associées à chaque premier point test P1 d'une première interface I ou trois sources SR rotationnelles de fluide associées à chaque premier point test Pl d'une première interface I, avec une seule sources S' scalaire de débit massique de fluide associée à chaque premier point test Pi' d'une deuxième interface I' ou deux sources S' scalaires de débit massique de fluide associées à chaque premier point test P1' d'une deuxième interface ou trois sources S' scalaires de débit massique de fluide associées à chaque premier point test Pl". Ce cas 8 est décrit ci-dessous pour trois sources SR rotationnelles de fluide en association à chaque premier point test Pl de l'interface I avec trois sources S' scalaires de débit massique de fluide associées à chaque premier point test Pl '. La matrice de couplage entre les points tests P5---Pi de A et les sources S'=Sj de B s" écrit -(z1 -z1) 2, 2y >1,z ( z 0 (z - z j) 0 4;z- pRu (x, - xi) - (x - x .) 47-cpRii 41-cpRu - (Y » 0 42-cpRii 47z-pRij avec T BA p La matrice M globale de résolution comprenant cette matrice de couplage /. Vo AA T f' A'\ (-) VB - ) BA TBB Ce cas 8, ainsi que les cas 9 et 10 ci-dessous permet de modéliser un fluide parfait mais avec des conditions aux limites nulles à la surface SUR de l'interface I, pouvant être par exemple une aile d'avion. 3 0 Cas 9 et 10 est :25 Suivant un mode de réalisation, il est prévu comme source de fluide au moins une source S scalaire de débit massique de fluide et au moins une source SR rotationnelle de fluide en association à chaque premier point test P1, pour la condition L aux limites ayant au moins une composante normale de vitesse du fluide au premier point test Pl, qui est nulle à l'interface I, et une composante tangentielle de vitesse du fluide au premier point test Pl, qui est prescrite. La figure 30 est une illustration de ce mode de réalisation. Ce mode de réalisation est dit hybride ou à sources S, SR hybrides. L'interface I est par exemple un objet B. Chacun des modes de réalisation et des variantes des cas 9 et 10 est par 10 exemple combiné avec le premier mode de réalisation sous-jacent précité (source S, SR située en dessous de la surface imperméable SUR du solide de l'interface I et de l'autre côté de l'écoulement de fluide F et/ou source S' située en dessous de l'interface I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un 15 objet A ou E précité. Dans ce mode de réalisation, il peut être prévu une seule, deux ou trois sources S scalaires de débit massique de fluide et une seule, deux ou trois sources SR rotationnelles de fluide en association avec chaque premier point test P 1 . Il peut être prévu en plus dans ce mode de réalisation au moins une (par 2 0 exemple une seule ou deux ou trois) source S' scalaire de débit massique de fluide associée à chaque premier point test P1' d'une deuxième interface I' distincte de la première interface I. Cas 9 2 5 Suivant un mode de réalisation, il est prévu comme source de fluide une source S scalaire de débit massique de fluide et deux sources SR rotationnelles de fluide en association à chaque premier point test P1, pour les conditions L aux limites ayant une composante normale de vitesse du fluide au premier point test Pl, qui est nulle à l'interface I, et une composante tangentielle de vitesse du fluide au 3 0 premier point test P 1, qui est prescrite. Ce cas 9 est décrit ci-dessous pour trois sources S' scalaires de débit massique de fluide associées à chaque premier point test P1' d'une deuxième interface I' distincte de la première interface I. Par exemple, la matrice de couplage entre les points tests Pl-Pi de A et les Yi - Y 42-cpRu3 47cpRu3 427-pRii3 47rpRij.3 zi-z yi-yj xi-xj 0 0 0 1 0 4 n-pRu 3 47-cpRii3 4 irpRij3 Par exemple, on introduit la matrice de couplage hybride H3 suivante, qui est, avec 2 sources rotationnelles Jx et Jy et une source scalaire : r Zi - Zj x - Xj 3 -3 RU 0 Yi-Y, R3 u xi-xi zi-zj R..3 R..3 / 0 : 2 1 zi - z ..3 Yi-Y D 3 47rp La matrice M globale de résolution, comprenant la matrice de couplage, est : 15 ( Vo ( Q A A IPA; (A VN =0 HB3 B y1.8 B ) \QBA de dimension : dim = (3n A +3nB)x(3n A +3n8) Ce mode de réalisation est par exemple utilisé pour un fluide incompressible, non visqueux (nombre de Reynolds élevé) avec prise en compte de 2 0 la vorticité en une seule phase de résolution. Des exemples sont : objet ou cylindre en rotation dans un écoulement avec effet de sol, et éventuellement gradient de vent pour illustrer le fait que le procédé permet des conditions aux limites complexes en chacun des points de la surface des objets. Dans ce mode de réalisation, la figure 12A représente la valeur absolue de la 2 5 composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et la figure 12B représente des lignes de courant de cet écoulement de fluide, ayant été calculées, et ce pour une interface I solide, formée par un cylindre CYL en sources S'=Sj de B s'écrit : r xi - xï Yi yi 47-cpRij 3 o o (),"\ 2 0 zi z 47rpRij 3 0 zi-z xi- xi r 2x Y rotation dans le sens horaire autour d'un axe parallèle à la direction x, où les conditions aux limites prévoient que la composante tangentielle du vecteur vitesse de l'écoulement de fluide est égale à une valeur prescrite du vent (par exemple ici 100 m.s-1) sur le cylindre. Le ou les cylindres sont des obstacles dans l'écoulement de fluide d'une soufflerie. Cela permet de modéliser l'effet de sol et les effets couplés avec le cylindre. Les figures 12C à 12G représentent les sources (figures 12C et 12D : sources rotationnelles représentées en un vecteur ligne de 3x1600 environ selon X et Y ; figures 12E, 12F, 12G : sources scalaires selon 1 tranche du cylindre en commençant en haut, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source ; les nuages denses des figures 12C et 12D représentent les même courbes condensées des figures 12E, 12F, 12G). Dans ce mode de réalisation, la figure 13A représente la valeur absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et la figure 13B représente des lignes de courant de cet écoulement de fluide, ayant été calculées, et ce pour une interface I solide, formée par un cylindre CYL en rotation dans le sens horaire autour d'un axe parallèle à la direction x, où les conditions aux limites prévoient que la composante tangentielle du vecteur vitesse de l'écoulement de fluide est égale à cinq fois une valeur prescrite du vent (cette 2 0 valeur prescrite du vent étant par exemple ici 100 m.s-1) sur le cylindre. Dans ce mode de réalisation, la figure 14A représente la valeur absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et la figure 14B représente des lignes de courant de cet écoulement de fluide, ayant été calculées, et ce pour une interface I solide, formée par un cylindre CYL en 2 5 rotation dans le sens horaire autour d'un axe parallèle à la direction x avec un gradient de vent, où les conditions aux limites prévoient que la composante tangentielle du vecteur vitesse de l'écoulement de fluide est égale à un quart d'une valeur prescrite du vent (cette valeur prescrite du vent étant par exemple ici 100 m.s 1) sur le cylindre. 30 Cas 10 Suivant un mode de réalisation, il est prévu comme source de fluide une source S scalaire de débit massique de fluide et trois sources SR rotationnelles de fluide en association à chaque premier point test Pl, pour la condition L aux limites ayant au premier point test Pl de l'interface I, une composante normale de vitesse 5 du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface I, et une composante normale nulle de vitesse d'un autre point test P5 prescrit d'une surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I. Ce mode de réalisation est décrit ci-dessous pour trois sources S' scalaires de débit massique de fluide associées à chaque premier 10 point test Pl' d'une deuxième interface I' distincte de la première interface I. Ce mode de réalisation est illustré à la figure 31. Par exemple, ce mode de réalisation a des conditions aux limites stratifiées à la surface SCL. Par exemple, la matrice M globale de résolution comporte par exemple la 15 matrice H4 ci-dessous : ( 0 zi- z j yi-yj xi- xj Y Hii4 .2 - zi /1z 4 npRii3 471-pRii3 47rpRii3 0 X/ Yi Y/ Yi 4 npRii3 Yi - 4 IrpRii3 471-pRii3 0 Zi - Zi 4;z-pRi j3 4 zpRii3 471-pRu3 de dimension : dim (3nAx3nB 3nAxnB)=(3n,x4nB) La matrice M globale de résolution est par exemple : ( \ rn TI 4 VO AA " AB 4 Vci-v = = QCA H CB Bh 2 5 " 4 V = 0 \,Q BA HBB \ B de dimension : ( 3nA r 3n A x 3n A 3n A x 4nB ( 3n 4nB / nC nB nc X3n A ncx4n8 3nB I 3n B x3n A 3nB x 4nB c est-à-dire de dimension dim = (3n A + 4n B)x (3n A + 4n B) 30 20 L'une des variantes des cas 9 et 10 ci-dessus est illustrée ci-dessous. Cette variante prévoit trois sources S scalaires de débit massique de fluide et une source SR rotationnelle de débit massique de fluide en association avec chaque premier point test Pl, pour les conditions L aux limites ayant au premier point test Pl une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface I, et une composante normale nulle de vitesse d'une couche limite CL de l'écoulement de fluide F. Cette variante est illustrée à la figure 15F. Cette variante est par exemple utilisée pour un fluide incompressible, non visqueux (nombre de Reynolds élevé), rotationnel, avec des conditions L aux limites imposées en une frontière connue de la couche limite CL, dont la topologie aura été préalablement établie. Des exemples sont : cylindre dans un écoulement avec effet de sol, et zoom au voisinage des parois, pour illustrer le fait que le procédé permet de fixer des conditions aux limites complexes en chacun des points de la surface des objets et en la frontière de la couche limite. Les figures 15D et 15E représentent un exemple de sources SR rotationnelles associées aux points tests Pl de l'interface I, pour une interface I solide, formée par un cylindre CYL ayant un axe parallèle à la direction x, les sources scalaires étant présentes et n'étant pas représentées. Dans cet exemple, la figure 15A représente la valeur absolue de la 2 0 composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y, et les figures 15B et 15C représentent des lignes de courant de cet écoulement de fluide, ayant été calculées, et ce pour l'interface I solide, formée par le cylindre CYL ayant un axe parallèle à la direction x, pour les conditions L aux limites ayant une composante normale de vitesse du fluide et une composante tangentielle de vitesse 2 5 du fluide, qui sont nulles à l'interface I. Dans cette variante, la figure 16A illustre des courbes R1, R2, R3, R4, R5 de la vitesse normale Vri calculée de l'écoulement de fluide F des figures 15A, 15B et 15C, et la figure 16B illustre des courbes R6, R7, R8, R9, R10 de la vitesse tangentielle Vt calculée de l'écoulement de fluide F des figures 15A, 15B et 15C, en 30 fonction de l'angle 'y autour de l'axe central du cylindre (parallèle à x) en abscisses, sur un cercle centré sur cet axe ayant respectivement le rayon suivant, pour une vitesse normale nulle et une vitesse tangentielle nulle aux points tests du cylindre et une vitesse normale nulle aux points tests de la couche limite: - R1 et R6 : cercle de rayon 20 m, qui est le rayon du cylindre CYL, - R2 et R7: cercle de rayon 21.25 m, - R3 et R8 : cercle de rayon 22.5 m, - R4 et R9 : cercle de rayon 23.75 m, R5 et R10 : cercle de rayon 25 m. Les figures 16C à 16H représentent les sources (figures 16C et 16D : sources rotationnelles selon X et Y; figures 12E, 16F, 16G, 16H: sources scalaires selon 1 10 tranche du cylindre en commençant en haut, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source). Dans cette variante, la figure 17A représente la valeur absolue de la composante Vy de la vitesse de l'écoulement de fluide ayant été calculée selon l'axe y et les figures 17B et 17C représentent des lignes de courant de cet écoulement de 15 fluide, ayant été calculées, et ce pour une interface I solide, formée par un cylindre CYL en rotation dans le sens horaire autour d'un axe parallèle à la direction x, où les conditions aux limites prévoient que la composante tangentielle du vecteur vitesse de l'écoulement de fluide est égale à une valeur prescrite du vent (par exemple ici 100 m.s-I) sur le cylindre. 2 0 Dans cette variante, les figures 18A et 18B représentent des courbes R11, R12, R13, R14, R15, R16, R17, R18, R19, R20 analogues aux courbes R1, R2, R3, R4, R5, R6, R7, R8, R9, R10 des figures 16A et 16B, pour une vitesse normale nulle, une vitesse tangentielle égale à une valeur prescrite du vent (par exemple ici 100 m.s-1) aux points tests du cylindre et une vitesse normale nulle aux points tests 2 5 de la couche limite. Les figures 18C à 18H représentent les sources (figures 18C et 18D : sources rotationnelles selon X et Y; figures 18E, 18F, 18G, 18H: sources scalaires selon 1 tranche du cylindre en commençant en haut, où l'abscisse est l'indice du point dans le maillage de point sur l'objet correspondant et l'ordonnée est le flux de la source). 30 Dans une autre variante des cas 9 est 10, il est prévu comme source S de fluide, en association avec chaque premier point test Pl trois sources S scalaires de fluide et une source SR rotationnelle de de fluide, pour les conditions L aux limites ayant au premier point test Pl une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface I, et en au moins un autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I une composante normale nulle de vitesse à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F. Cas 11, 12, 13 Suivant un mode de réalisation, représenté aux figures 32, 33, 34, il est prévu comme source S de fluide au moins une septième source S7 de force en association à chaque premier point test Pl, pour la condition L aux limites ayant à chaque premier point test Pl une composante noiniale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface I. Chacun des modes de réalisation et des variantes des cas 11, 12 et 13 est par exemple combiné avec le premier mode de réalisation sous-jacent précité (source S7 et/ou S8 et/ou S9 située en dessous de la surface imperméable SUR du solide de 2 0 l'interface I et de l'autre côté de l'écoulement de fluide F et/ou source S' située en dessous de l'interface I' et de l'autre côté de l'écoulement de fluide F), ce qui est le cas ci-dessous. L'interface I est par exemple un objet B. L'interface I' est par exemple un objet A ou E précité. Dans chacun des modes de réalisation des cas 11, 12 et 13, il peut être prévu une seule, deux ou trois sources de force (S7 et/ou S8 2 5 et/ou S9) en association avec chaque premier point test Pl. Il peut être prévu en plus dans chacun ce mode de réalisation au moins une (par exemple une seule ou deux ou trois) source S' scalaire de débit massique de fluide associée à chaque premier point test Pl' d'une deuxième interface I' distincte de la première interface I. 30 Ces cas 11, 12 et 13 peuvent être utilisés pour un fluide incompressible pour un nombre de Reynolds faible, c'est-à-dire inférieur ou égal à 2000. On prend en compte les termes de viscosité, ce qui aboutit à l'équation de Stokes, dont les fonctions de Green sont les `Stokeslet". Ces fonctions donnent les équations de pression et celles de la vitesse, qui apparaissent couplées de façon linéaire dans l'équation de Stokes. Cas 11 Suivant un mode de réalisation, représenté à la figure 32, il est prévu comme source S de fluide: - en association avec chaque premier point test Pl de l'interface I au moins 10 trois septième, huitième et neuvième sources S7, S8, S9 de force, situées en des deuxième points P20, P21, P22, qui sont distincts de l'interface I, - en association avec au moins un autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, en dessous de la surface SCL délimitant la couche limite CL du côté de 15 l'interface I, trois première, deuxième et troisième sources Sl, S2, S3 de force, et - en association avec le au moins autre point test prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, au-dessus de la surface SCL délimitant la couche limite CL trois quatrième, cinquième et sixième sources S4, S5, S6 de force, 20 pour les conditions L aux limites ayant au premier point test Pl une composante normale de vitesse du fluide F et une composante tangentielle de vitesse du fluide F, qui sont nulles à l'interface I, et à l'autre point test P5 une composante normale de vitesse du fluide F nulle à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante 25 tangentielle de la vitesse du fluide F à la surface délimitant la couche limite CL. La première source Si de force et/ou la deuxième source S2 de force et/ou la troisième source S3 de force est donc située dans la première zone intérieure Z1. La quatrième source S4 de force et/ou la cinquième source S5 de force et/ou la sixième source S6 de force est donc située dans la deuxième zone extérieure Z2 . 30 Des exemples sont : cylindre en rotation ou fixe dans un écoulement, avec prise en compte de l'effet de sol pour illuster le fait que le procédé permet une analyse des interactions entre objets. Les sources élémentaires S7, S8, S9 sont des triplets vectoriels dont l'unité est homogène à une force (en N), ce qui est logique puisque apparaît la notion de pertes par frottement, donc de traînée. Par exemple, à nombre de Reynolds faibles, les équations de Stokes ont pour 5 solution les Stokeslets, à savoir : o o av =ki; V= Vv V 10 19 = 4.7t- ( (x, -xj) ° -Y,) o o (z, ) ij 1 '(x; -x;) (yi - y .) (xi -x]).(zi (y-yi).(zi-z» «.; Par exempel, sur chaque point test P5 de la surface SCL, on pose 2 conditions (Egalité de la pression et de la normale). 15 On a par exemple 3 points test P5 pour chaque bloc de 6 sources Si, S2, S3, S4, S5, S6 de la surface SCL de la couche limite. La matrice M globale de résolution exprimant la pression est P en stokeslets et M en scalaire est : Vo ( 0K 0 r A Csi N N nN P2VCext Pl VCi int P2CA P1KCCse P2KCCs1 P1KCB kse 2 0 P2- 191= 0 McA PCCse CCsi PCB 2csi Va = 0 0 KBCse 0 Kas ) Â.A / Cas 12 Suivant un mode de réalisation, représenté à la figure 33, il est prévu comme 2 5 source S de fluide: - en association avec chaque premier point test Pl de l'interface I au moins une septième source S7 de force, - en association avec au moins un autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à 30 l'interface I, en dessous de la surface SCL délimitant la couche limite CL du côté de l'interface I, une première source Si scalaire de fluide, et - en association avec le au moins autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, au-dessus de la surface SCL délimitant la couche limite CL au moins une quatrième source S4 de force, pour les conditions L aux limites ayant au premier point test Pl , une composante normale de vitesse du fluide F et une composante tangentielle de vitesse du fluide F, qui sont nulles à l'interface I, et à l'autre point test P5 une composante normale de vitesse du fluide F nulle à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante 10 tangentielle de la vitesse du fluide F à la surface délimitant la couche limite CL. Suivant un mode de réalisation, représenté à la figure 33, il est prévu comme source S de fluide: - en association avec chaque premier point test Pl de l'interface I au moins trois septième, huitième et neuvième sources S7, S8, S9 de force, situées en des 15 deuxième points P20, P21, P22 qui sont distincts de l'interface I, - en association avec au moins un autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, en dessous de la surface SCL délimitant la couche limite CL du côté de l'interface I, une première source S1 scalaire de fluide, et 2 0 - en association avec le au moins autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, au-dessus de la surface SCL délimitant la couche limite CL des quatrième, cinquième et sixième sources S4, S5, S6 de force, pour les conditions L aux limites ayant au premier point test Pl, une 2 5 composante normale de vitesse du fluide F et une composante tangentielle de vitesse du fluide F, qui sont nulles à l'interface I, et à l'autre point test P5 une composante normale de vitesse du fluide F nulle à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante tangentielle de la vitesse du fluide F à la surface délimitant la couche limite CL. 3 0 La première source Si scalaire de fluide est donc située dans la première zone intérieure Zl. La quatrième source S4 de force et/ou, la cinquième source S5 de force et/ou la sixième source S6 de force sont donc situées dans la deuxième zone extérieure Z2 . Les cas 12 et 13 peuvent être utilisés pour un nombre de Reynolds faible dans la couche limite, c'est-à-dire dans la première zone Z1 située entre la surface SCL et l'interface I, où le nombre de Reynolds est inférieur ou égal à 2000, et un nombre de Reynolds élevé au-delà, c'est-à-dire dans la deuxième zone Z2 située au-dessus de la surface SCL et située de l'autre côté par rapport à l'interface I, où le nombre de Reynolds est supérieur ou égal à 3000. Dans la zone Z1, on peut négliger les termes de viscosité et on prend les fonctions hamioniques en 1/R et leurs dérivées. Dans la zone Z2, on prend en compte les termes de viscosité et on utilise comme fonctions de Green les Stokeslet. On considère par exemple un fluide incompressible. Les conditions de raccord à la frontière SCL de la couche limite sont traitées comme des interfaces classiques en DPSM, avec des séries de sources ponctuelles de part et d'autre de cette surface SCL ou bien comme cela a été réalisé dans les cas précédents, toutes les sources sont à l'intérieur de l'interface I ou profil I. La matrice M globale de résolution est : 2 0 Vo = . KCC Q AC1 KCB K BB ). 3 \ V4i3 V4/3 =0 -CA KBC QCCI 2Cs Cext C int 0 0 2c1 VB = 0 de dimension : ( 3n Ax3n A 3n A x3nc 3n A x n, 3ne x3n, dim = 4nc x 3n A 4n, x3n, 4nc xnc 4n, x3n, = (3n A +4n, +3n8)x (3n A +3n, + ne +3%) 3ne 3n4 3n, x3nc 3n, x nc 3n,x3n8 2 5 Cas 13 Suivant un mode de réalisation, représenté à la figure 34, il est prévu comme source S de fluide: - en association avec chaque premier point test Pl de l'interface I au moins une septième source S7 de force, 3 0 - en association avec au moins un autre point test prescrit P5 de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, en dessous de la surface SCL délimitant la couche limite CL du côté de l'interface I, au moins une première source S1 scalaire de fluide et au moins une deuxième source rotationnelle S2 de fluide, - en association avec le au moins autre point test prescrit P5 de la surface 5 SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, au-dessus de la surface SCL délimitant la couche limite CL des quatrième, cinquième et sixième sources S4, S5, S6 de force, pour les conditions L aux limites ayant au premier point test P1 une composante normale de vitesse du fluide F et une composante tangentielle de 10 vitesse du fluide F, qui sont nulles à l'interface I, et à l'autre point test P5 une composante normale de vitesse du fluide F nulle à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante tangentielle de la vitesse du fluide F à la surface SCL délimitant la couche limite CL. 15 Suivant un mode de réalisation, il est prévu comme source S de fluide: - en association avec chaque premier point test P1 de l'interface I au moins trois septième, huitième et neuvième sources S7, S8, S9 de force, situées en des deuxième points P20, P21, P22 qui sont distincts de l'interface I, - en association avec au moins un autre point test P5 prescrit de la surface 20 SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, en dessous de la surface SCL délimitant la couche limite CL du côté de l'interface I, une première source S1 scalaire de fluide, une deuxième source rotationnelle S2 de fluide et une troisième source S3 de fluide, scalaire ou rotationnelle, 2 5 - en association avec le au moins autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I, au-dessus de la surface SCL délimitant la couche limite CL, trois quatrième, cinquième et sixième sources S4, S5, S6 de force, pour les conditions L aux limites ayant au premier point test P1 une 3 0 composante normale de vitesse du fluide F et une composante tangentielle de vitesse du fluide F, qui sont nulles à l'interface I, et à l'autre point test P5 une composante normale de vitesse du fluide F nulle à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante tangentielle de la vitesse du fluide F à la surface SCL délimitant la couche limite CL. Par exemple, il est représenté à la figure 34 la troisième source S3 de fluide, scalaire. La première source Si scalaire de fluide et/ou la deuxième source rotationnelle S2 de fluide et/ou la troisième source S3 de fluide, scalaire ou rotationnelle, est donc située dans la première zone intérieure Zl. La quatrième source S4 de force et/ou la cinquième source S5 de force et/ou la sixième source S6 de force sont donc situées dans la deuxième zone extérieure Z2 . Ainsi, on superpose à un ensemble de conditions limites L pris en un maillage de points test (premier point P1) répartis à la surface SUR de l'interface I, un autre ensemble de conditions aux limites L pris en un maillage de points tests P5 répartis à la surface SCL de la couche limite CL. On appelle cette superposition "conditions aux limites stratifiées". L'introduction de conditions aux limites stratifiées permet de tenir compte du problème de la couche limite CL. Par exemple, au contact d'un profil d'aile, la vitesse normale est nulle (l'aile est une surface SUR imperméable au fluide F), mais la vitesse tangentielle du fluide F s'annule 2 0 également (à cause de la viscosité). Au-delà de la couche limite CL dans laquelle ce phénomène de viscosité se fait sentir, le fluide F se comporte comme un fluide parfait, et à la surface SCL délimitant la couche limite CL, seule la composante normale de la vitesse est nulle. Ces conditions aux limites L permettent de prendre des conditions aux limites vectorielles à la surface SUR de l'interface I, auxquelles 2 5 on superpose les conditions aux limites scalaires à la surface SCL de la couche limite CL (exemple : composante normale de la vitesse nulle à la surface SCL de la couche limite CL). En variante aux cas 11, 12 et 13, il est prévu comme source S de fluide, en association avec chaque premier point test Pl, une source S scalaire de fluide et 3 0 trois sources S7, S8, S9 de force, pour les conditions L aux limites ayant au premier point test P1 une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface I, et en au moins un autre point test P5 prescrit de la surface SCL délimitant une couche limite CL de l'écoulement de fluide par rapport à l'interface I une composante nounale nulle de vitesse du fluide à la surface SCL délimitant la couche limite CL de l'écoulement de fluide F et une continuité de la composante tangentielle de la vitesse du fluide à la surface SCL délimitant la couche limite CL. La matrice M globale de résolution est : Vo u = AA K CC H3 K CB ( v6/3 _17613 QCA K BC ACh K BB ) 2Cs 2Ch rCext -Vent 3 VB = HC Ch 0 Dans les cas 11, 12 et 13, la source S5 et/ou S6 et/ou S7 située dans la zone 10 extérieure Z2 envoie du fluide dans la zone intérieure Zl au travers de la surface SCL, tandis que la source S1 et/ou S2 et/ou S3 située dans la zone intérieure Z1 envoie du fluide dans la zone extérieure Z2 au travers de la surface SCL. Les cas 11, 12 et 13 et leurs variantes peuvent s'appliquer à des secteurs micro-fluidiques et des écoulements visqueux. Dans les cas 11, 12 et 13, les 15 conditions aux limites peuvent comprendre la conservation de la pression de fluide à l'autre point test P5 et la conservation, à l'autre point test P5, de la composante normale de la vitesse modulo la masse volumique. Les 13 cas précités sont récapitulés sur l'organigramme représenté sur la figure 35. 20 D'une manière générale ci-dessus, le nombre de points tests P1 et/ou P5 et/ou Pl' et/ou le nombre de conditions aux limites L est adapté pour que la matrice globale M de résolution soit carrée et pour qu'ainsi le nombre d'équations reste égal au nombre d'inconnues. Par exemple, à la figure 32, on choisit trois points tests P5 pour chaque groupe de sources Sl, S2, S3, S4, S5, S6. Par exemple, à la figure 33, il 25 y a un rapport de 4/3 entre le nombre de points tests P5 et le nombre de blocs de sources S5, S6, S7, Sl. Par exemple, à la figure 32, on choisit trois points tests P5 pour chaque groupe de sources Sl, S2, S3, S4, S5, S6. Par exemple, à la figure 34, il y a un rapport de 6/3 entre le nombre de points tests P5 et le nombre de blocs de sources S5, S6, S7, S I, S2, S3. Le caractère de non réversibilité de l'écoulement peut être obtenu par un calcul des pertes dans chaque section de la couche limite le long de la progression 5 de l'écoulement, la géométrie des sections suivantes étant élaborée à partir des conditions obtenues dans chacune des sections précédentes. En aéronautique, un cas très intéressant est celui dans lequel l'environnement aérologique ou la position de l'avion change brutalement (variation d'angle d'incidence avion, rafale latérale ou verticale, gradient de vent, etc). Suivant un 10 mode de réalisation, ce cas est représenté par une simple modification des conditions aux limites L (vecteur C), la matrice inverse de M ne change pas, et on obtient une nouvelle valeur des sources S quasi instantanément. C'est l'un des avantages majeurs que de disposer d'un modèle semi-analytique, comparativement à une simulation 3D éléments finis par exemple. Dans l'exemple précité de l'avion, 15 lorsque ce dernier s'incline, chaque point test P3 placé sur le maillage de l'avion se déplace avec un vecteur différent: la vitesse horizontale de l'avion et la vitesse en chaque point liée au mouvement de roulis deviennnent les nouvelles conditions aux limites de l'avion. On peut donc évaluer le comportement d'un objet en mouvement quelconque 2 0 dont la trajectoire se situe dans un fluide dont la topologie varie, instant après instant (ce qu'on peut appeler du 4D pour quatre dimensions, puisqu'on rajoute le temps). D'autre part, la méthode DPSM en mécanique des fluides permet également la prise en compte d'un environnement de propriétés inhomogènes par zones 2 5 (exemple: couches stratifiées de différentes densités). Chaque limite entre domaines devient une interface I qui est maillée comme il est décrit précédemment. Les idées directrices propres à cette méthodologie consistent donc à accéder beaucoup plus facilement qu'avec les méthodes actuelles (3D éléments finis notamment) à une représentation réaliste des écoulements de fluide autour de profils 3 0 I. Les applications sont nombreuses, même si on peut rester dans le cadre des écoulements non turbulents: - Modélisation de tout ou partie d'un avion afin de mieux restituer son comportement en simulateur de vol ou optimiser le fonctionnement des pilotes automatiques ; optimisation des profils aéronautiques en vue de diminuer les trâinées et l'amélioration des algorithmes de fonctionnement des simulateurs de vol et l'entraînement des pilotes ; - Dimension réduite du problème par un maillage des surfaces SUR au lieu d'un maillage de tout l'espace avec donc des temps de calculs nettement réduits ; - Prise en compte de l'influence mutuelle de plusieurs objets I comme par exemple la modélisation de l'effet de sol éventuellement en présence d'un gradient 10 de vent (en vitesse et en direction), du fuselage ou du couple aile empennage ; - Optimisation de profils par injection d'écoulements d'air à certains endroits précis, ce qui nécessite une modélisation complexe ; - Analyse fine des profils d'ailes pour diminuer la tramée et ainsi réduire la consommation de carburant des avions, cette diminution ayant un impact sur le coût 15 d'exploitation et sur le bilan carbone des vols, - hydrodynamique, domaine maritime. Ainsi, l'application de la méthode DPSM aux équations de Navier-Stokes, qui régissent le comportement de l'écoulement de fluide F, permet d'obtenir une représentation en trois dimensions la plus fidèle possible de l'écoulement autour de 20 l'interface I, pouvant être un objet de forme complexe et/ou mobile. Suivant un mode de réalisation, l'invention permet notamment de juxtaposer la résolution du problème en trois dimensions à l'extérieur de la couche limite, où les équations se simplifient alors un peu en fluide incompressible, avec la résolution du problème dans la couche limite. Pour cela, suivant ce mode de réalisation, des 25 sources rotationnelles (donc vectorielles ou transversales) et des sources scalaires (donc radiales) sont disposées en dessous de l'interface I, c'est-à-dire dans le premier milieu Ml, et éventuellement de part et d'autre du sillage au-delà de l'interface I (surface de glissement). Ainsi, l'invention, suivant un mode de réalisation, utilise des sources hybrides ou au moins une source hybride de fluide (laquelle comporte 30 au moins une source scalaire et au moins une source rotationnelle). Cela permet, suivant ce mode de réalisation, de calculer les grandeurs dans l'espace environnant l'interface I (deuxième milieu M2) dans la couche limite et hors de la couche limite. L'introduction de blocs de sources hybrides, contenant des sources scalaires et vectorielles permet de satisfaire les conditions imposées par la vorticité. L'invention permet dans certains modes de réalisation une résolution 5 itérative, en 2 ou 3 étapes, permettant entre chaque étape de déterminer des grandeurs qui sont réinjectées dans les étapes suivantes, comme les vitesses tangentielles ou la définition de la couche limite. Les grandeurs peuvent ensuite être calculées facilement dans tout l'espace de travail et les différents objets peuvent être placés en interaction. Par ailleurs, le fait 10 de faire coexister des sources de natures différentes (rotationnelles et scalaires) qui interagissent mutuellement permet de fournir une solution précise aux problèmes traités. Un autre avantage est d'associer des conditions limites supplémentaires, éloignées des points tests 131 positionnés à la surface SUR de l'interface I, ce qui permet de donner une bonne représentation des phénomènes également dans la zone 15 proche de l'interface I (induisant une possible restitution des phénomènes de frottement dans la couche limite CL). L'invention permet ainsi d'accéder beaucoup plus facilement qu'avec les méthodes actuelles (éléments finis 3D) à une représentation réaliste des écoulements de fluides autour de profils. 2 0 L'invention permet notamment de modéliser tout ou partie d'un avion, afin de mieux restituer son comportement en simulateur de vol ou optimiser le fonctionnement des pilotes automatiques. Le procédé de mesure et/ou dispositif de mesure peut ainsi être embarqué sur un aéronef, tel que par exemple un avion ou un hélicoptère, ou sur un simulateur 2 5 de vol, ou sur un navire, ou plus généralement un véhicule de transport automoteur, pour y être mis en oeuvre. Ainsi, un autre objet de l'invention peut être un procédé de commande et ou un dispositif de commande pour commander un aéronef, tel que par exemple un avion ou un hélicoptère, ou un simulateur de vol, ou un navire, ou plus 3 0 généralement un véhicule de transport automoteur, dans lequel il est prévu une étape de calcul, dans laquelle on calcule à partir de la grandeur physique ayant été calculée par le procédé de mesure et/ou dispositif de mesure, au moins une grandeur de commande par un calculateur embarqué sur l'aéronef, simulateur de vol, navire ou véhicule de transport automoteur. Dans ce qui précède l'actionneur peut être un de ceux de cet aéronef, avion, hélicoptère, simulateur de vol, navire, ou plus généralement d'un véhicule de transport automoteur. L'invention pennet d'avoir une dimension réduite du problème par un maillage des surfaces au lieu d'un maillage de tout l'espace et donc d'avoir des temps de calcul nettement réduits. L'invention permet de prendre en compte l'influence mutuelle de plusieurs 10 objets comme, par exemple, la modélisation de l'effet de sol éventuellement en présence d'un gradient de vent (en vitesse et en direction), du fuselage ou du couple aile/empennage. L'invention permet également d'optimiser des profils (surfaces SUR de l'interface I) par injection d'écoulement d'air à certains endroits précis, ce qui 15 nécessite une modélisation complexe. L'invention permet également une analyse fine des profils d'ailes pour diminuer la traînée et ainsi réduire la consommation de carburant des avions. Cette diminution a un impact sur le coût d'exploitation et sur le bilan carbone des vols. Par conséquent, l'invention apporte une nette diminution du temps de calcul 20 par rapport aux techniques connues de résolution en deux dimensions par des méthodes de points singuliers (résolution en potentiel complexe et lignes portantes de Prandt, par exemple), qui avaient un temps de calcul trop important, notamment pour modéliser un avion complet en mouvement dans un environnement variable (gradient de vent, effets de sol). 25 Bien entendu, les modes de réalisation décrits ci-dessus peuvent être combinés entre eux. REVENDICATIONS1. Dispositif de mesure d'au moins une première grandeur physique (V, P) d'au moins un écoulement de fluide (F) dans un espace tridimensionnel ayant au 5 moins une interface (I) prédéterminée, située entre au moins deux milieux (Ml, M2), caractérisé en ce que le dispositif comporte : au moins un calculateur ayant des moyens (E) de prescription d'au moins une condition (L) aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au moins en un premier point test prédéterminé (P1) 10 de l'interface (I), le calculateur (CAL) ayant des moyens de calcul configurés pour calculer par méthode DPSM à partir de la condition (L) aux limites, la première grandeur physique (V, P) de l'écoulement de fluide (F) en au moins un deuxième point (P) de l'espace, différent du premier point test (P1). 15 2. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que l'interface (I) comprend au moins une surface (SUR) imperméable au fluide (F). 3. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que l'interface (I) comprend au moins une surface 2 0 (SUR) d'au moins un solide, imperméable au fluide (F). 4. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que la première grandeur physique (V, P) est la vitesse et/ou la pression du fluide. 5. Dispositif de mesure suivant l'une quelconque des revendications 2 5 précédentes, caractérisé en ce que le calculateur (CAL) comporte des troisièmes moyens de calcul d'une deuxième grandeur de commande d'un actionneur à partir au moins de la première grandeur physique (V, P) ayant été calculée. 6.. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que les moyens de calcul (E2) sont configurés pour 3 0 calculer par méthode DPSM :- au moins un paramètre d'au moins une source (S) de fluide modélisée, située en au moins un troisième point (P2) de position prescrite, distinct de l'interface (I), à partir de la condition (L) aux limites, - la première grandeur physique (V, P) de l'écoulement de fluide (F) située 5 au deuxième point (P) de l'espace, à partir du paramètre de la source (S) de fluide, ayant été calculé. 7. Dispositif de mesure suivant la revendication 6, caractérisé en ce que la source (S) de fluide est séparée du premier point test prédéterminé (P1) de l'interface (I). 8. Dispositif de mesure suivant l'une quelconque des revendications 6 10 et 7, caractérisé en ce que l'interface (I) comprend au moins une surface (SUR) d'au moins un solide, imperméable au fluide (F), au moins une source (S) de fluide associée au premier point test (PI) étant située à distance du premier point test (Pl) en dessous de la surface (SUR) du solide, de l'autre côté de l'écoulement de fluide (F). 15 9. Dispositif de mesure suivant l'une quelconque des revendications 6 à 8, caractérisé en ce qu'il est prévu au moins une source (S) de fluide non associée au premier point test (P1) et située dans l'écoulement de fluide (F). 10. Dispositif de mesure suivant l'une quelconque des revendications 6 à 9, caractérisé en ce que dans les moyens de calcul est enregistrée une matrice globale 2 0 de résolution M, comportant au moins un coefficient (Mu) dépendant à la fois d'une valeur prescrite caractérisant le fluide (p, p.) et de la distance (Ru) entre une source (S) de fluide et un quatrième point (P) de l'espace, le produit de la matrice globale de résolution M, prise aux premiers points test (Pl) comme quatrième point (P) de l'espace, multipliée par un premier vecteur J 2 5 des paramètres des sources (S) de fluide, étant égal à un deuxième vecteur C des conditions aux limites concernant la première grandeur physique (V) de l'écoulement de fluide prise aux premiers points test (P1), selon l'équation C = M * J, où * désigne la multiplication. 11. Dispositif de mesure suivant la revendication 10, caractérisé en ce qu'il 3 0 est prévu parmi les moyens de calcul- un moyen d'inversion de la matrice globale de résolution M, prise aux premiers points test (P1) comme quatrième point (P) de l'espace, pour calculer la matrice inverse M-1, et - un moyen de calcul du premier vecteur J des paramètres des sources (S) de 5 fluide par multiplication de la matrice inverse M'1 par le deuxième vecteur C des conditions aux limites de la première grandeur physique (V) de l'écoulement de fluide prise aux premiers points test (P1), selon l'équation M-I * C. 12. Dispositif de mesure suivant la revendication 10 ou 11, caractérisé en ce 10 qu'il est prévu parmi les moyens de calcul : - un moyen de calcul de la première grandeur physique (V, P) de l'écoulement de fluide (F) au deuxième point (P) de l'espace ayant une position déterminée, en multipliant la matrice globale de résolution M, calculée à ladite position déterminée du deuxième point (P) comme quatrième point (P), par le 15 premier vecteur J des paramètres des sources (S) de fluide. 13. Dispositif de mesure suivant l'une quelconque des revendications 6 à 12, caractérisé en ce que la source (S) de fluide est choisie parmi une source (S) scalaire de fluide, une source (S) rotationnelle de fluide et une source de force. 14. Dispositif de mesure suivant l'une quelconque des revendications 6 à 13, 2 0 caractérisé en ce que la source (S) de fluide comprend au moins une source (S) scalaire de fluide et au moins une source (S) rotationnelle de fluide. 15. Dispositif de mesure suivant l'une quelconque des revendications 6 à 13, caractérisé en ce que la source (S) de fluide comprend au moins une source (S) rotationnelle de fluide et/ou au moins une source de force. 2 5 16. Dispositif de mesure suivant l'une quelconque des revendications 6 à 15, caractérisé en ce que la source (S) de fluide comprend au moins quatre sources choisies parmi au moins une source (S) scalaire de fluide, au moins une source (S) rotationnelle de fluide et au moins une source (S) de force. 17. Dispositif de mesure suivant l'une quelconque des revendications 6 à 14 3 0 et 16, caractérisé en ce que la source (S) scalaire de fluide est une source de débit massique isotrope de fluide. 18. Dispositif de mesure suivant l'une quelconque des revendications 6 à 17, caractérisé en ce que la source (S) rotationnelle de fluide est une source provoquant une rotation du fluide autour d'une direction déterminée. 19. Dispositif de mesure suivant l'une quelconque des revendications 10 à 12, caractérisé en ce que la valeur prescrite caractérisant le fluide est la masse volumique p du fluide et/ou la viscosité cinématique u du fluide. 20.. Dispositif de mesure suivant la revendication 19, caractérisé en ce que (CAS 1, 2, 3) la matrice globale (M) de résolution comporte plusieurs coefficients, dont au moins un est égal à (x,-xj)/(4.n.p.R,j3), ou à (yi-yi)/(4. TC p.R,j3) ou à (zrzi)/(4.n. p.R,j3), ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur (J) des paramètres des sources (Sj) de fluide, où R1 est la distance entre la source (Sj) de fluide située au troisième point (P2) de position prescrite xj, y, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point (P) de l'espace ayant des coordonnées x' y,, z, selon les trois directions x, y et z. 21. Dispositif de mesure suivant la revendication 19, caractérisé en ce que (CAS 4) la matrice globale (M) de résolution comporte plusieurs coefficients, dont au 2 5 moins un est égal à 1/(4.7r. p .R,j), ou à celui-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur (J) des paramètres des sources (Sj) de fluide, où Ru est la distance entre la source (Sj) de fluide située au troisième point (P2) de position prescrite et le 3 0 quatrième point (P) de l'espace. 22. Dispositif de mesure suivant la revendication 19, caractérisé en ce que (CAS 8) la matrice globale (M) de résolution comporte plusieurs coefficients, dont au moins un est égal à ±(xi-xj)/(4.7c.p.R,j), ou à ±(yryj)/(4.7c.p.R,i) ou à ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur (J) des paramètres des sources (Sj) de fluide, où R- est la distance entre la source (Sj) de fluide située au troisième point (P2) de position prescrite xj, yj, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point (P) de l'espace ayant des coordonnées xi, yi, z, selon les trois directions x, y et z. 23. Dispositif de mesure suivant la revendication 19, caractérisé en ce que (CAS 9, 10) la matrice globale (M) de résolution comporte plusieurs coefficients, dont au moins un est égal à 2 0 ±(x;-xj)/(4.n.p.R,j3), ou à ±(yryj)/(4.7r.p.Rii3) ou à 2 5 ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur (J) des paramètres des sources (Sj) de fluide, où est la distance entre la source (Sj) de fluide située au troisième point (P2) de position prescrite xi, yj, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point (P) de l'espace ayant des coordonnées xi, yi, zi selon les trois 3 0 directions x, y et z. 24. Dispositif de mesure suivant la revendication 19, caractérisé en ce que (CAS 11) la matrice globale (M) de résolution comporte plusieurs coefficients, dont au moins un est égal à l'un parmi (yj-yi)/(8.n.t.t.R,j), (zi-zi)/(8 . TC. p-Rii), (y,-yi)2/(8 .7C. p-Rii),(Zi-Zj)2/(8 .7C. 1-1. Rij) (Xi-Xj).(yryi)/(8.n.i.i.R,j3), ou à l'un de ceux-ci, multiplié par une constante prescrite, par laquelle est divisée le premier vecteur (J) des paramètres des sources (Sj) de fluide, où est la distance entre la source (Sj) de fluide située au troisième point (P2) de position prescrite xi, yj, zj selon trois directions non coplanaires x, y et z de l'espace et le quatrième point (P) de l'espace ayant des coordonnées x' y,, z, selon les trois directions x, y et z. 25. Dispositif de mesure suivant l'une quelconque des revendications 6 à 24, 15 caractérisé en ce que le nombre (Np) de points tests (P1) des conditions aux limites (L) est supérieur ou égal au nombre (Ns) de sources (S) de fluide. 26. Dispositif de mesure suivant l'une quelconque des revendications 6 à 25, caractérisé en ce qu'un nombre Ns de sources (S) de fluide et un nombre Np de points tests (Pl) distincts les uns des autres sont prévus, de telle sorte que la matrice 20 globale de résolution M soit carrée. 27. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 1) il est prévu comme source de fluide une seule source (S) scalaire de débit massique de fluide associée à chaque premier point test (P 1), pour la condition (L) aux limites ayant au premier point test (P1) une composante 25 normale de vitesse du fluide au premier point test (P1), qui est nulle à l'interface (I). 28. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 2 et 3) il est prévu comme source de fluide trois sources (S) scalaires de débit massique de fluide associées à chaque premier point test (P1), 30 pour la condition (L) aux limites ayant une composante normale de vitesse du fluideet une composante tangentielle de vitesse du fluide au premier point test (Pl), qui sont nulles à l'interface (I). 29. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 8) il est prévu comme source de fluide au moins une source (S) rotationnelle de fluide en association à chaque premier point test (P1), pour les conditions (L) aux limites ayant une composante normale de vitesse du fluide au premier point test (P1), qui est nulle à l'interface (I), et une composante tangentielle de vitesse du fluide au premier point test (P1), qui est prescrite. 30. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 9 et 10) il est prévu comme source (S) de fluide au moins une source (S) scalaire de débit massique de fluide et au moins une source (S) rotationnelle de fluide en association à chaque premier point test (Pl), pour la condition (L) aux limites ayant au moins une composante noiniale de vitesse du fluide au premier point test (Pl), qui est nulle à l'interface (I), et une composante tangentielle de vitesse du fluide au premier point test (Pl), qui est prescrite. 31. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 4) il est prévu comme source (S) de fluide au moins une première source (S) scalaire de débit massique de fluide en association à chaque premier point test (Pl) de l'interface (I), au moins une deuxième source (S20) 2 0 scalaire de débit massique de fluide située au-dessus d'une surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide (F) surface (SCL) par rapport à l'interface (I) et au moins une troisième source (S30) scalaire de débit massique de fluide située en dessous de la surface (SCL) délimitant la couche limite (CL), pour la condition (L) aux limites ayant au premier point test (Pl) de l'interface (I) une 2 5 composante normale de vitesse du fluide, nulle à l'interface (I), et une composante normale de vitesse conservée en un autre point test (P5) prescrit de la surface (SCL) délimitant la couche limite (CL), modulo la masse volumique du fluide. 32. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 5) il est prévu comme source (S) de fluide au moins 3 0 deux sources (S) scalaires de débit massique de fluide en association à chaque premier point test (Pl) de l'interface (I), pour la condition (L) aux limites ayant aupremier point test (P1) de l'interface (I) une composante normale de vitesse du fluide, nulle à l'interface (I), et une composante normale de vitesse nulle en un autre point test prescrit (P5) d'une surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I). 33. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 11, 12, 13) il est prévu comme source (S) de fluide au moins une source (S7) de force en association à chaque premier point test (P1), pour la condition (L) aux limites ayant à chaque premier point test (P1) une composante normale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface (I). 34. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 11) il est prévu comme source (S) de fluide: - en association avec chaque premier point test (P1) de l'interface (I) au moins une septième source (S7) de force, - en association avec au moins un autre point test prescrit de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I), en dessous de la surface (SCL) délimitant la couche limite (CL) du côté de l'interface (I), au moins une première source (S1) de force, et - en association avec le au moins autre point test prescrit (P5) de la surface 20 (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I), au-dessus de la surface (SCL) délimitant la couche limite (CL) au moins une quatrième source (S4) de force, pour les conditions (L) aux limites ayant au premier point test (P1) une composante normale de vitesse du fluide (F) et une composante tangentielle de 25 vitesse du fluide (F), qui sont nulles à l'interface (I), et à l'autre point test (P5) une composante normale de vitesse du fluide (F) nulle à la surface (SCL) délimitant la couche limite (CL) de l'écoulement de fluide (F) et une continuité de la composante tangentielle de la vitesse du fluide (F) à la surface délimitant la couche limite (CL). 35. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, 3 0 caractérisé en ce que (CAS 12) il est prévu comme source (S) de fluide:- en association avec chaque premier point test (Pl) de l'interface (I) au moins une septième source (S7) de force, - en association avec au moins un autre point test prescrit (P5) de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à 5 l'interface (I), en dessous de la surface (SCL) délimitant la couche limite (CL) du côté de l'interface (I), une première source (Si) scalaire de fluide, et - en association avec le au moins autre point test prescrit (P5) de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I), au-dessus de la surface (SCL) délimitant la couche limite (CL) au 10 moins une quatrième source (S4) de force, pour les conditions (L) aux limites ayant au premier point test (Pl), une composante normale de vitesse du fluide (F) et une composante tangentielle de vitesse du fluide (F), qui sont nulles à l'interface (I), et à l'autre point test (P5) une composante normale de vitesse du fluide (F) nulle à la surface (SCL) délimitant la 15 couche limite (CL) de l'écoulement de fluide (F) et une continuité de la composante tangentielle de la vitesse du fluide (F) à la surface délimitant la couche limite (CL). 36. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce que (CAS 13) il est prévu comme source (S) de fluide: - en association avec chaque premier point test (PI) de l'interface (I) au 2 0 moins une septième source (S7) de force, - en association avec au moins un autre point test prescrit (P5) de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I), en dessous de la surface (SCL) délimitant la couche limite (CL) du côté de l'interface (I), au moins une première source (Si) scalaire de fluide et au 2 5 moins une deuxième source rotationnelle de fluide, - en association avec le au moins autre point test prescrit (P5) de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I), au-dessus de la surface (SCL) délimitant la couche limite (CL) des quatrième, cinquième et sixième sources (S4, S5, S6) de force, 3 0 pour les conditions (L) aux limites ayant au premier point test (Pl) une composante normale de vitesse du fluide (F) et une composante tangentielle devitesse du fluide (F), qui sont nulles à l'interface (I), et à l'autre point test (P5) une composante normale de vitesse du fluide (F) nulle à la surface (SCL) délimitant la couche limite (CL) de l'écoulement de fluide (F) et une continuité de la composante tangentielle de la vitesse du fluide (F) à la surface délimitant la couche limite (CL). 37. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce qu'il est prévu (cas 10) comme source (S) de fluide, en association avec chaque premier point test (P1) trois sources (S) scalaires de fluide et une source (S) rotationnelle de de fluide, pour les conditions (L) aux limites ayant au premier point test (P1) une composante nolluale de vitesse du fluide et une composante tangentielle de vitesse du fluide, qui sont nulles à l'interface (I), et en au moins un autre point test (P5) prescrit de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I) une composante normale nulle de vitesse à la surface (SCL) délimitant la couche limite (CL) de l'écoulement de fluide (F). 38. Dispositif de mesure suivant l'une quelconque des revendications 6 à 26, caractérisé en ce qu'il est prévu (CAS 13) comme source (S) de fluide, en association avec chaque premier point test (P1) une source (S) scalaire de fluide et trois sources (S) de force, pour les conditions (L) aux limites ayant au premier point test (P1) une composante normale de vitesse du fluide et une composante 2 0 tangentielle de vitesse du fluide, qui sont nulles à l'interface (I), et en au moins un autre point test (P5) prescrit de la surface (SCL) délimitant une couche limite (CL) de l'écoulement de fluide par rapport à l'interface (I) une composante normale nulle de vitesse du fluide à la surface (SCL) délimitant la couche limite (CL) de l'écoulement de fluide (F) et une continuité de la composante tangentielle de la 2 5 vitesse du fluide à la surface (SCL) délimitant la couche limite (CL). 39. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que le calculateur (CAL) comporte des deuxièmes moyens de calcul, à partir de la première grandeur physique (V, P) ayant été calculée, d'une propriété de l'interface (I) dans l'écoulement de fluide (F), différente 3 0 de la première grandeur physique (V, P). 40. Dispositif de mesure suivant la revendication précédente, caractérisé en ce que la propriété est la traînée induite par l'interface (I) dans l'écoulement de fluide (F) et/ou la portance de l'interface (I) dans l'écoulement de fluide (F) et/ou une force de frottement de l'interface (I) dans l'écoulement de fluide (F). 41. Dispositif de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que le calculateur (CAL) comporte des moyens (SD) de sortie pour fournir la première grandeur physique (V, P) ayant été calculée et/ou la propriété ayant été calculée. 42. Procédé de mesure d'au moins une première grandeur physique (V, P) 10 d'au moins un écoulement de fluide (F) dans un espace tridimensionnel ayant au moins une interface (I) prédéterminée, située entre au moins deux milieux (M1, M2), caractérisé en ce que au cours d'une première étape (E 1), on prescrit par au moins un calculateur au moins une condition (L) aux limites concernant la première grandeur physique 15 (V, P) de l'écoulement de fluide (F) prise en au moins un premier point test prédéterminé (Pl) de l'interface (I), au cours d'une deuxième étape (E2), on calcule par le calculateur (CAL) par méthode DPSM à partir de la condition (L) aux limites la première grandeur physique (V, P) de l'écoulement de fluide (F) en au moins un deuxième point (P) de 20 l'espace, différent du premier point test (P1). 43. Procédé de mesure selon la revendication 42, caractérisé en ce que lors de la deuxième étape (E2), on prescrit la position d'au moins un troisième point (P2) distinct de l'interface (I), en lequel se trouve au moins une source (S) de fluide, 25 on calcule par méthode DPSM au moins un paramètre de la source (S) de fluide à partir de la condition (L) aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au premier point test (Pl) de l'interface (I), on calcule par méthode DPSM la première grandeur physique (V, P) de 3 0 l'écoulement de fluide (F) située au deuxième point (P) de l'espace, à partir du paramètre de la source (S) de fluide, ayant été calculé. 44. Procédé de mesure selon la revendication 43, caractérisé en ce que (CAS 6) au cours d'une première itération des première et deuxième étapes (El, E2), on prescrit une première condition aux limites concernant la première grandeur 5 physique (V, P) de l'écoulement de fluide (F) prise au premier point test prédéterminé (Pl) de l'interface (I) et un premier nombre détenniné de sources (S) de fluide, supérieur ou égal à un, pour calculer par méthode DPSM une première itération de la première grandeur physique (V, P) de l'écoulement de fluide (F) située au deuxième point (P) de l'espace à partir du au moins un paramètre du 10 premier nombre déterminé de sources (S) de fluide, ayant été calculé, puis au cours d'au moins une deuxième itération des première et deuxième étapes (El, E2), postérieure à la première itération, on prescrit une deuxième condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au premier point test prédéteuniné (Pl) de 15 l'interface (I), la deuxième condition aux limites étant différente de la première condition aux limites, on ajoute au premier nombre déterminé de sources (S) de fluide un deuxième nombre déterminé de sources (S) de fluide, supérieur ou égal à un, pour calculer par méthode DPSM au moins une deuxième itération de la première grandeur physique (V, P) de l'écoulement de fluide (F) située au deuxième 20 point (P) de l'espace à partir du au moins un paramètre des premier et deuxième nombres déterniinés de sources (S) de fluide, ayant été calculé. 45. Procédé de mesure selon la revendication 43, caractérisé en ce que (CAS 7) au cours d'une première itération des première et deuxième étapes (El, E2), 25 on prescrit une première condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au premier point test prédéterminé (Pl) de l'interface (I) et un premier nombre déterminé de sources (S) de fluide, supérieur ou égal à un, pour calculer par méthode DPSM une première itération de la première grandeur physique (V, P) de l'écoulement de fluide (F) 30 située au deuxième point (P) de l'espace à partir du au moins un paramètre du premier nombre déterminé de sources (S) de fluide, ayant été calculé,puis au cours d'au moins une deuxième itération des première et deuxième étapes (El, E2), postérieure à la première itération, on prescrit une deuxième condition aux limites concernant la première grandeur physique (V, P) de l'écoulement de fluide (F) prise au premier point test prédéterminé (P 1) de 5 l'interface (I) et un deuxième nombre déterminé de sources (S) de fluide, supérieur ou égal à un, pour calculer par méthode DPSM une deuxième itération de la première grandeur physique (V, P) de l'écoulement de fluide (F) située au deuxième point (P) de l'espace à partir du au moins un paramètre du deuxième nombre déterminé de sources (S) de fluide, ayant été calculé, la deuxième condition aux 10 limites étant différente de la première condition aux limites, puis on additionne la première itération de la première grandeur physique (V, P) de l'écoulement de fluide (F) et la deuxième itération de la première grandeur physique (V, P) de l'écoulement de fluide (F). 46. Programme d'ordinateur, comportant des instructions pour la mise en 15 oeuvre du procédé de mesure selon l'une quelconque des revendications 42 à 45, lorsqu'il est mis en oeuvre sur un calculateur. 20