Méthodes des Fonctions

de Base Radiale

Pour la Solution D`équations Aux Dérivées Partielles

Établissements

À propos de..

Le projet de recherche: Méthodes de fonctions de base radiale pour la solution de EDPs, a pour objectif premier d'étudier, rechercher et développer de nouveaux algorithmes numériques basés sur la théorie des fonctions de base radiales pour la solution de Èquations aux Dérivées Partielles (EDP). Ces techniques, développées au cours de cette dernière décennie, constituent un outil puissant pour résoudre des problèmes réels à grande échelle. De par sa nature, ce projet est interdisciplinaire puisque son objectif inclut également la modélisation mathématiques de problèmes abordant différents domaines de connaissances.

Un autre objectif important de ce projet est d'élaborer une bibliothèque de codes informatiques permettant d'implémenter ces méthodes. Il ne s’agit pas seulement d'utiliser des “boîtes noires” pour résoudre des problèmes d'EDP, mais de renforcer la recherche et le développement de systèmes numériques et informatiques pour ces nouvelles méthodes.

Ce portail contient du matériel qui couvre trois objectifs principaux divisés en : recherche, enseignement et développement de logiciels:

  • 01 Recherche:

    Dans ce domaine, nous avons l’intention de réaliser des recherches tant dans le domaine de méthodes numériques appliquées que dans des domaines appartenant à différentes disciplines.

  • Nous avons inclu des documents et des codes informatiques introductifs ainsi que des explications, qui visent à aider autant les professeurs que les étudiants qui étudient et/ou enseignent dans ce domaine.

  • Le troisième objectif est celui de construire une bibliothèque de programmes libres dont la flexibilité permet à l'utilisateur d'intégrer de nouvelles routines ou des algorithmes.

 

Numerical Methods for the Solution of Partial Differential Equations: A Review

(Méthodes numériques pour la solution d'équations aux dérivées partielles : une critique)

Recherche

Fonctions de bass radiale pour la solution d'EDP:

  • Commentaires Généraux

    Bien que les méthodes des fonctions radiales de base radiales servant à la résolution d'EDP aient été développées dans les années 1990, leur impact est resté dans les communautés des ingénieurs et des mathématiciens (voir [I, II, III, IV, V , VI] et les références qui y sont contenues). L’absence de mailles pour l’élaboration de ces schémas numériques explique en grande partie le succès de ces techniques. En simulations réelles, il est fréquent que 70% du temps de traitement soit dû à la génération de mailles. Dans ce contexte, nous remarquons les éléments suivants qui font que ces méthodes sont plus puissantes que les méthodes classiques :

    Aucune de ces méthodes ne requiert une construction de mailles, élément qui réduit considérablement la complexité algorithmique de ces méthodes. Il existe des stratégies simples pour affiner le nombre de nœuds dans les régions où la fonction présente des oscillations et des gradients élevés. En comparaison, les techniques de nœuds adaptés sont beaucoup moins coûteuses que les algorithmes de mailles adaptées.

  • Le traitement des problèmes de dimensions supérieures à deux et avec des frontières complexes a des résultats plus simples et efficaces qu'avec des méthodes dépendants de mailles. Il existe plusieurs façons de classer les méthodes numériques de fonctions radiales de base pour la solution d'équations différentielles partielles. En particulier, elles peuvent être regroupées en deux catégories :
    • Méthodes de base radiale:
    → Méthodes globales
    → Méthodes locales
    • Méthodes des moindres carrés mobiles

  • Ce qui distingue ces deux approches est que, dans le cas des méthodes de base radiale, l'ordre de convergence du schéma dépend du type de noyau de base radial, alors que pour les méthodes des moindres carrés mobiles, même quand le noyau est radial, l'ordre de convergence est déterminé par la technique des moindres carrés. En d'autres termes, les méthodes de base radiales héritent de l'ordre de convergence du noyau radial qui le détermine, qui peut être de type exponentiel, alors que les méthodes de moindres carrés mobiles ont un ordre de convergence algébrique indépendant du noyau radial utilisé.

    Une autre différence importante est que les méthodes de base radial sont invariantes par rapport à la dimension de l'espace, c'est-à-dire que la nature radiale de la méthode fait que le système algébrique de Gramm correspond à la méthode ayant la même complexité numérique en une, deux ou plus dimensions. En revanche, les méthodes de moindres carrés mobiles, dont la caractérisation dépend d'un réglage numérique, ont une complexité numérique dépendant de la dimension de l'espace euclidien.

    D'un autre côté, les méthodes globales de type radial dérivent d'un système algébrique unique d'équations incluant tous les nœuds du domaine, alors que les méthodes des moindes carrés mobiles nous conduisent à un ensemble de systèmes algébriques relativement petits par rapport au nombre total de nœuds. Simultanément, alors que les méthodes de moindres carrés mobiles sont consistants à l'EDP, les méthodes basés en noyaux radiaux en reproduisent pas nécessairement des polynômes et en sont donc pas consistants à l'EDP.

  • Dans le cas des méthodes globales de fonctions de base radiale, un des principaux problèmes est que le conditionnement de la matrice du système augmente lorsque le nombre de nœuds augmente. Cela en s eproduit pas dans le cas des méthodes de MLS, étant donné que les systèmes locaux en permettent pas de construire des matrices bandes et à diagonales dominantes. Les méthodes locales, cependant, ne peuvent pas gémérer des schémas convergeants exponentiellement, tandis que les globales le peuvent. De plus, il a été développé, pour les méthodes globales, des techniques de décomposition de domaine et de préconditionnement, qui permettent de résoudre des problèmes pour un grand nombre de nœuds.

  • Les deux approches ont eu un impact important dans différentes disciplines et sont efficacité pour résoudre différents problèmes est une question qui fait l'objet de la recherche actuelle. Parmi les méthodes de ces deux groupes, les plus importantes sont :
    Méthodes de base radiale:
    • Emplacement symétrique et asymétrique
    • Quadrature différentielle
    • Méthodes de fonctions radiale à support compact
    • Méthodes de solutions fondamentales

    Méthodes de moindres carrés mobiles:
    • Petro Galerkin sans maille (MLPG)
    • Elements finis par partition de l'unité (PUFEM)
    • Noyaux reproduits par partition de l'unité (RKPM)
    • Hydrodynamique de particules lisses (SPH)

  • Dans un travail récent, Fasshauer [VII] a introduit une formulation abstraite unifiant les deux approches – pour le cas des noyaux radiaux – dans un cadre commun. En particulier, l'auteur a démontré que les méthodes locales obtenues par les moindres carrés peuvent se transformer en méthodes globales lorsque le nombre de nœuds du support des techniques locales entourent tous les nœuds du domaine.


    [I] S. N. Atluri and S. Shen, The Meshless Local Petrov-Galerkin (MLPG) Method, Tech Science Press, Encino, CA, 2002.
    [II] M. D. Buhmann, Radial Basis Functions : Theory and Implementations, Cambridge University Press, Cambridge, 2003.
    [II] E. W. Cheney and W. A. Light, A Course in Approximation Theory, Brooks/Cole, Pacific Grove, CA, 1999.
    [IV] G. R. Liu, Mesh Free Methods: Moving beyond the Finite Element Method, CRC Press, Boca Raton, FL, 2002.
    [V] Wendland, H., Scattered Data Approximation, Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge, UK., 2005
    [VI] G. E. Fasshauer, Meshfree Approximation Methods with Matlab Interdisciplinary Mathematical Sciences - Vol. 6 World Scientific Publishers, Singapore, 2007
    [VII] G. E. Fasshauer, Dual Bases and Discrete Reproducing Kernels: A Unified Framework for RBF and MLS Approximation, Engineering Analysis with Boundary Elements, Volume 29, Issue 4, April 2005, pages 313-325

  • II Apliccation des méthodes de RBF dans différentes disciplines:
    1. Mathématiques financières:
    L'équation différentielle partielle de Black-Scholes a été fréquemment utilisée pour la modélisation des prix des contrats financiers. Les modèles d'option de tarification en une dimension pour les options américaines et européennes ont été explorées par des méthodes de RBF, par Hon et al. [1,2] et en deux dimensions par Fasshauer et al [3] et par Marcozzi et al. [4]. Ils ont formulé une méthode basée sur des fonctions quasi-radiales en une dimension [5]. Récemment, Petterson et al. [6] ont résolu le problème multi-dimensionnel pour l'option européenne basket call par une méthode d'emplacement asymétrique avec noyau multiquadratique et noeuds adaptés, reportant que ses résultats sont de 20 à 40 fois plus rapides que les résultats obtenus avec des différences finies adaptées.

    [1] Y.-C. Hon, X.-Z. Mao, A radial basis function method for solving options pricing models, J. Financial Engineering 8 (1999) 31–49.
    [2] Z.Wu, Y.-C. Hon, Convergence error estimate in solving free boundary diffusion problem by radial basis functions method, Engrg. Anal. Bound. Elem. 27 (2003) 73–79.
    [3] G. E. Fasshauer, A. Q. M. Khaliq, D. A. Voss, Using meshfree approximation for multi-asset American option problems, J. Chinese Institute Engineers 27 (2004) 563–571.
    [4] M. D. Marcozzi, S. Choi, C. S. Chen, On the use of boundary conditions for variational formulations arising in financial mathematics, Appl. Math. Comput. 124 (2001) 197–214.
    [5] Y. C. Hon, A quasi-radial basis functions method for American options pricing, Comput. Math. Appl. 43 (3) (2002) 513–524.
    [6] Ulrika Pettersson, Elisabeth Larsson, Gunnar Marcusson, and Jonas Persson, Improved Radial Basis Function Methods for Multi-Dimensional Option Pricing, Technical Report 2006-028, Uppsala University, 2006.

  • Dans [7], N. Flyer et al, le problème de la modélisation du mouvement des flux géophysiques se déplaçcant dans une sphère dominée par des limites connectées est posé. Ceci est le premier ouvrage dans lequel sont résolues des équations différentielles partielles hyperboliques grâce aux fonctions de base radiale. Les auteurs rapportent que le nombre de noeuds requis est inférieur que celui utilisé par des méthodes spectrales. De plus, la méthode RBF permet de passer plus de temps pour obtenir la même précision que celles utilisés par les méthodes spectrales classiques.
    Dans Jichun Li [8] , Hon [9] et Wong et al. [10], les équations pour les dépôts en eau peu profonde (shallow water equations) sont résolues par la méthode de placement asymétrique avec noyau multiquadratique, En particulier, dans [20], une étude est réalisée pour le port de Tolo et HongKong, qui est comparée avec des techniques d'éléments finis et des données de terrain, reportées par des supérieurs sur la technique des EF et une excellente coïncidence avec ces données.
    [7] N. Flyer and G.B. Wright, Transport schemes on a sphere using radial basis functions. J. Comp. Phys., to appear (2007).
    [8] Jichun Li, C.S. Chen, Darrel Pepper, Yitung Chen, "Mesh-free method for groundwater modeling", Boundary Elements XXIV, eds. C.A. Brebbia, A. Tadeu and V. Popov, WIT Press, Southampton, Boston, pp. 115-124 (2002).
    [9] Y. C. Hon, K. F. Cheung, X. Z. Mao, and E. J. Kansa, .A Multiquadric Solution for Shallow Water Equations", ASCE J. Hydraulic Engineering, 125, (5), 524-533 (1999).
    [10] Wong, S. M., Hon, Y. C., and Golberg, M. A. 2002. Compactly supported radial basis functions for shallow water equations. Appl. Math. Comput. 127, 1 (Mar. 2002), 79-101
    La modélisation des écoulements diphasiques et triphasiques est abordée par Hon dans [11] et [12] en utilisant une approche de placement asymétrique eulérien utilisant des multiquadriques. A. Iske et coll. [13] [14] abordent le même problème en développant une méthode semi-lagrangienne qui utilise des noyaux de plaques minces pour la solution des équations de Buckley-Leverett, afin de modéliser des problèmes en deux phases applicables à la récupération primaire des puits en l'industrie pétrolière.
    [11] Y. C. Hon, M. W. Lu, W. M. Xue, and X. Zhou, "Multiquadric Method for the Numerical Solution of a Biphasic Model", Internat. J. Appl. Sci. Comput., 88, 153-175 (1997).
    [12] Y. C. Hon , M. Lu , M. W. Xue and X. Zhou, "Numerical Algorithm for Triphasic Model of Charged and Hydrated Soft Tissues", Computacional Mechanics, 29(1), , pp 1-15 (2002)
    [13] A. Iske and M. Käser: Two-Phase Flow Simulation by AMMoC, an Adaptive Meshfree Method of Characteristics. Computer Modeling in Engineering & Sciences (CMES) 7(2), 2005, 133-148.
    [14] J. Behrens, A. Iske, and M. Käser: Adaptive Meshfree Method of Backward Characteristics for Nonlinear Transport Equations, in Meshfree Methods for Partial Differential Equations, M. Griebel and M. A. Schweitzer (eds.), Springer-Verlag, Heidelberg, 2002, 21-36.
    Chen et coll. [15] étudient la modélisation du transport des contaminants pour les réservoirs d'eau à l'aide de techniques de placement RBF. L'étude Le numérique comprend plusieurs cas : diffusion pure, advection et dispersion pour une source continue, advection et dispersion pour une source instantanée.
    [15] J. Li, Y. Chen and D. Pepper, Radial basis function method for 1-D and 2-D groundwater contaminant transport modeling, Computational Mechanics, (32) 10-15, 2003

  • La solution de l'équation de Navier Strokes, tant stationnaire que dépendante du temps, a été abordé par Shu, et al. [16] [17] [18], où se résoud le problème de la cavité carrée, Ghia, grâce à la méthode de quadrature différentielle pour les nombre de Reinolds de 10^5.
    [16] Shu, C., Khoo, B. C., and Yeo, K. S. 1994. Numerical solutions of incompressible Navier—Stokes equations by generalized differential quadrature. Finite Elem. Anal. Des. 18, 1-3 (Dec. 1994), 83-97.
    [17] C. Shu, H. Ding and K. S. Yeo (2005), 'Computation of incompressible Navier-Stokes equations by local RBF-based differential quadrature method', CMES-Computer Modeling in Engineering and Sciences, 7, 195-205.
    [18] H. Ding, C. Shu, K. S. Yeo and D. Xu (2006), 'Numerical Computation of Three-dimensional Incompressible Viscous Flows in the Primitive Variable Form by Local Multiquadric Differential Quadrature Method', Computer Methods in Applied Mechanics and Engineering, 195, 516-533.
    Chinchapatnam et al. [19] ont recherché l'équation non stationnaire de convection diffusion en une ou deux dimensions pour différents nombres de Peclet par un placement symétrique et asymétrique. Le schéma numérique utilise des fonctions radiales globales, telles que les multiquadratiques, multiquadratiques inverses, splines en plaque mince et d'ordre 5. Dans cette étude comparative, les auteurs ont conclu que la méthode symétrique est meilleure que l'asymétrique, et que pour des grands nombres de Peclet il recquièrent une forte densité de noeuds pour obtenir une bonne approximation.
    [19] Chinchapatnam, P.P., Djidjeli, K and Nair, P.B. (2006). Unsymmetric and symmetric meshless schemes for the unsteady convection–difusión equation. Computer Methods in Applied Mechanics and Engineering, 195, (19-22), 2432-2453.

  • Il existe peu de références, aussi loin que nous le savons, sur les algorithmes basés sur les fonctions radiales globales qui résolvent des problèmes d'analyse structurale. , le travail de Ferreira et al. [20] qui aborde le problème des efforts de poutres et plaques déformables, ainsi que le travail de Zhang et al. [21] sur l'élasticité plane. Récemment, Tiago et al. [22] utilisent la méthode de placement asymétrique avec des noyaux radiaux distincts pour résoudre et analyser les problèmes à une dimension comme des problèmes d'élasticité linéaire statique, de vibration libre et analyse linéaire de stabilité, pour des problèmes physiqes non-linéaires de modèles endommagés.
    [20] A. J. M. Ferreira, C. M. C. Roque, and P. A. L. S. Martins. Radial basis functions and higher-order shear deformation theories in the analysis of laminated composite beams and plates. Composite Structures, 66(1–4):287–293, 2004.
    [21] X. Zhang, K. Z. Song, M. W. Lu, and X. Liu. Meshless methods based on collocation with radial basis functions. Computational Mechanics, 26(4):333–343, 2000.
    [22] C. Tiago and V. Leitão, Application of radial basis functions to linear and non-linear structural analysis problems, Computers & Mathematics with Applications, 51(8), 1311-1334, 2006.

Un état de l'art dans les méthodes de fonctions radiales de base :

  • -Article de revue:

    Méthodes numériques pour résoudre les EDP par les RBF

  • • Placement symétrique et asymétrique
    • Quadrature différentielle
    • Galerkin.
    • Méthodes de préconditionnement
    • Décomposition des domaines
    → Algorithme avec traslepe
    → Algorithme sans traslepe
    • Algorithmes de nœuds adaptés
    → Algorithmes locaux
    → Algorithmes globaux
    • Visualisation des solutions

  • Flux de phase non miscible. Equation de Buckley Leverett
    Modélisation du transport des contaminants de l’eau d'un réservoir

 
Flux de phase non miscible. Equation de Buckley Leverett
Bibliothèque C++
Fonctions de base radiale - UNAM
Bibliothèque pour les problèmes de fonction de base radiale
RBF++
Documentation version 1.0 mai 2008
UNAM

Comme nous l'avons mentionné dans la page principale, un objectif important de ce projet est de construire une bibliothèque de logiciels permettant d'implémenter des algorithmes sur les fonctions de base radiales pour la solution d'EDP. Un objectif majeur de ce point-là est de permettre à l'utilisateur d'inclure d'une manière facile de nouveaux modules et procédures par l'intermédiaire d'une stratégie de programmation par objet en C++.
La bibliothèque est actuellement au stade initial, nous invitons les développeurs de logiciels intéressés par ce domaine à ajouter des modules orientés objet suivant l'approche du logiciel libre (licence GPL). Pour cela, contactez Daniel A, Cervantes pour toute question liée à cette bibliothèque,

Documentation
Classe
Dossiers
Un logiciel différent peut ètre trouvé dans:
Templates Units for Numerical Applications: Radial Basis Functions.
Reference Manual
TUNA::RBF
Class List
Radial Basis Function Library for PDEs Problems
RBF++
Documentation Ver. 1.0 May 2008
UNAM
Class List
LA::Matrix Base class for matrices
LA::Point Base class for Points
LA::Point2dT Class for two dimension points
LA::Point3dT Class for three dimension points
RBF::GramMatrix Gram matrix class
RBF::MQ Multicuadric kernel
RBF::TPS2 Second order thin plate spline kernel
RBF::TPS4 Fourth order thin plate spline kernel
RBF::TPS_2DX Second order partial derivative of thin plate spline (with respect to )
RBF::TPS_2DY Second order partial derivative of thin plate spline (with respect to )
RBF::TPSInterp Thin plate spline interpolator class
Lista de Archivo
Radial Basis Function Library for PDEs Problems
RBF++
Documentation Ver. 1.0 May 2008
UNAM
File List
    
RBF++-0.1/examples/linear_alg/linear_alg.cpp [code]	 
RBF++-0.1/examples/tps_interpolation/tps_interpolation.cpp [ code]	 
RBF++-0.1/examples/tps_poisson/tps_poisson.cpp [code]	 
RBF++-0.1/include/LA/Matrix.h [code]	 
RBF++-0.1/include/LA/Point.h [code]	 
RBF++-0.1/include/LA/Point2dT.h [code]	 
RBF++-0.1/include/LA/Point3dT.h [code]	 
RBF++-0.1/include/RBF/GrammMatrix.h [code]	 
RBF++-0.1/include/RBF/RBF.h [code]	 
RBF++-0.1/include/RBF/TPSInterp.h [code]	 
RBF++-0.1/src/LA/Matrix.cpp [code]	 
RBF++-0.1/src/LA/Point.cpp [code]	 
RBF++-0.1/src/RBF/GrammMatrix.cpp [code]	 
RBF++-0.1/src/RBF/TPSInterp.cpp [code] 


Meteorología y Geofísica:

Dans [7], N. Flyer et al. Le problème de la modélisation du mouvement des écoulements géophysiques dans une sphère dominée par des termes connectifs est posé. Il s'agit de Le premier travail dans lequel des équations aux dérivées partielles hyperboliques sont résolues à l'aide de fonctions de base radiales. Les auteurs rapportent que le nombre de nœuds requis est bien inférieur à celui utilisé par les méthodes spectrales. De plus, la méthode RBF permet des pas de temps beaucoup plus longs pour obtenir la même précision que ceux utilisés par les méthodes spectrales classiques.

Dans Jichun Li [8], Hon [9] et Wong et al. [10] les équations pour les réservoirs d'eau peu profonde (équations d'eau peu profonde) sont résolues à l'aide de la méthode placement asymétrique avec noyaux multiquadriques. En particulier, dans [20] une étude est réalisée pour les ports de Tolo et HonKong, qui est comparée à techniques d'éléments finis et de données de terrain, rapportant une supériorité sur la technique FE et une excellente coïncidence avec les données de terrain.
[7] N. Flyer and G.B. Wright, Transport schemes on a sphere using radial basis functions. J. Comp. Phys., to appear (2007).
[8] Jichun Li, C.S. Chen, Darrel Pepper, Yitung Chen, "Mesh-free method for groundwater modeling", Boundary Elements XXIV, eds. C.A. Brebbia, A. Tadeu and V. Popov, WIT Press, Southampton, Boston, pp. 115-124 (2002).
[9] Y. C. Hon, K. F. Cheung, X. Z. Mao, and E. J. Kansa, .A Multiquadric Solution for Shallow Water Equations", ASCE J. Hydraulic Engineering, 125, (5), 524-533 (1999).
[10] Wong, S. M., Hon, Y. C., and Golberg, M. A. 2002. Compactly supported radial basis functions for shallow water equations. Appl. Math. Comput. 127, 1 (Mar. 2002), 79-101
La modélisation des écoulements diphasiques et triphasiques est abordée par Hon dans [11] et [12] en utilisant une approche de placement asymétrique eulérien utilisant des multiquadriques. A. Iske et coll. [13] [14] abordent le même problème en développant une méthode semi-lagrangienne qui utilise des noyaux de plaques minces pour la solution des équations de Buckley-Leverett, afin de modéliser des problèmes en deux phases applicables à la récupération primaire des puits en l'industrie pétrolière.

[11] Y. C. Hon, M. W. Lu, W. M. Xue, and X. Zhou, "Multiquadric Method for the Numerical Solution of a Biphasic Model", Internat. J. Appl. Sci. Comput., 88, 153-175 (1997).
[12] Y. C. Hon , M. Lu , M. W. Xue and X. Zhou, "Numerical Algorithm for Triphasic Model of Charged and Hydrated Soft Tissues", Computacional Mechanics, 29(1), , pp 1-15 (2002)
[13] A. Iske and M. Käser: Two-Phase Flow Simulation by AMMoC, an Adaptive Meshfree Method of Characteristics. Computer Modeling in Engineering & Sciences (CMES) 7(2), 2005, 133-148.
[14] J. Behrens, A. Iske, and M. Käser: Adaptive Meshfree Method of Backward Characteristics for Nonlinear Transport Equations, in Meshfree Methods for Partial Differential Equations, M. Griebel and M. A. Schweitzer (eds.), Springer-Verlag, Heidelberg, 2002, 21-36.
Chen et coll. [15] étudient la modélisation du transport des contaminants pour les réservoirs d'eau à l'aide de techniques de placement RBF. L'étude numérique comprend plusieurs cas : diffusion pure, advection et dispersion pour une source continue, advection et dispersion pour une source instantanée.
[15] J. Li, Y. Chen and D. Pepper, Radial basis function method for 1-D and 2-D groundwater contaminant transport modeling, Computational Mechanics, (32) 10-15, 2003


Dynamique des fluides:

La solution de l'équation de Navier Stokes, à la fois stationnaire et dépendante du temps, a été abordée par Shu et al. [16] [17] [18], où résout le problème de la cavité carrée, Ghia, en utilisant la méthode de quadrature différentielle pour les nombres de Reinolds de 10^5.
[16] Shu, C., Khoo, B. C., and Yeo, K. S. 1994. Numerical solutions of incompressible Navier—Stokes equations by generalized differential quadrature. Finite Elem. Anal. Des. 18, 1-3 (Dec. 1994), 83-97.
[17] C. Shu, H. Ding and K. S. Yeo (2005), 'Computation of incompressible Navier-Stokes equations by local RBF-based differential quadrature method', CMES-Computer Modeling in Engineering and Sciences, 7, 195-205.
[18] H. Ding, C. Shu, K. S. Yeo and D. Xu (2006), 'Numerical Computation of Three-dimensional Incompressible Viscous Flows in the Primitive Variable Form by Local Multiquadric Differential Quadrature Method', Computer Methods in Applied Mechanics and Engineering, 195, 516-533.
Chinchapatnam et coll. [19] étudient l'équation non stationnaire de la diffusion par convection en une et deux dimensions pour différents nombres de Péclet par placement symétrique et asymétrique. Le schéma numérique utilise des fonctions radiales globales, telles que les multiquadratiques, les multiquadratiques inverses, les splines à plaques minces et les quintiques. Dans cette étude comparative, les auteurs concluent que la méthode symétrique n’est que légèrement meilleure que la méthode asymétrique. et que pour les grands nombres de Péclet, une forte densité de nœuds est nécessaire pour obtenir une bonne approximation.
[19] Chinchapatnam, P.P., Djidjeli, K and Nair, P.B. (2006). Unsymmetric and symmetric meshless schemes for the unsteady convection–difusión equation. Computer Methods in Applied Mechanics and Engineering, 195, (19-22), 2432-2453.


Élasticité et analyse structurelle

Il existe peu de références, à notre connaissance, sur des algorithmes basés sur des fonctions radiales globales qui résolvent des problèmes d'analyse structurelle. Notamment, les travaux de Ferreira et al. [20] qui traite du problème de contrainte des poutres et plaques composites déformables, ainsi que les travaux de Zhang et al. [21] en élasticité plane. Récemment, Tiago et al. [22] utilisent la méthode de placement asymétrique utilisant différents noyaux radiaux pour résoudre et analyser des problèmes dans une dimension allant de l'élasticité linéaire statique, des vibrations libres et de l'analyse de la stabilité linéaire, aux problèmes physiques non linéaires de modèles endommagés.
[20] A. J. M. Ferreira, C. M. C. Roque, and P. A. L. S. Martins. Radial basis functions and higher-order shear deformation theories in the analysis of laminated composite beams and plates. Composite Structures, 66(1–4):287–293, 2004.
[21] X. Zhang, K. Z. Song, M. W. Lu, and X. Liu. Meshless methods based on collocation with radial basis functions. Computational Mechanics, 26(4):333–343, 2000.
[22] C. Tiago and V. Leitão, Application of radial basis functions to linear and non-linear structural analysis problems, Computers & Mathematics with Applications, 51(8), 1311-1334, 2006.


ENSEIGNEMENT



Codes en matlab et C++
Programmes d'interpolation de surface.

Interpolation de surface à partir de données aléatoires

Algoritmo
Interpolation utilisant des fonctions radiales.
Interpolation à partir de plaques minces : codes dans Matlab

Programme principal
Programme principal et routines auxiliaires (compressées)

Interpolation à partir de plaques minces : codes en C++
Programme principal
Programme principal et routines auxiliaires (compressées)

Programmes pour résoudre l'équation de diffusion par convection

Diffusion par convection
Convection 2D pure, utilisant un placement non symétrique.
Matlab
Routine principale

Routine principale et auxiliaire (compressée)

C++
Ut + Ux + Uy =0

Principales routines
Méthode explicite, TPS+Pol grade 6

Méthode implicite TPS+Pol niveau 6

Méthode explicite Noyau multiquadratique

Routine principale et auxiliaire (compressée)

Convection - Diffusion 2D, utilisant un placement non symétrique.
Matlab
Routine principale

Routine principale et auxiliaire (compressée)

C++
Ut + mu*(Ux + Uy) = beta*grad(U)

Principales routines
Méthode implicite, TPS + Pol grade 6

Méthode de noyau multiquadrique implicite

Routine principale et auxiliaire (compressée)

Avis aux étudiants:

Annonce

Les étudiants de premier cycle et de troisième cycle sont invités à obtenir leur diplôme en résolvant des problèmes décrits en termes d'EDP en utilisant des techniques de calcul parallèle et des fonctions à base radiale, ainsi qu'en utilisant la bibliothèque développée à partir de ce projet.

Pour les sujets de thèse de premier cycle et de troisième cycle, veuillez contacter le Dr Pedro González Casanova (Instituto de Matemáticas, UNAM).
TEL: 56 22 47 47, OU ENVOYEZ-MOI UN EMAIL: casanova@matem.unam.mx

Lectures recommandées :
→ Wendland, H., Scattered Data Approximation, Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge, UK., 2005
→ G. E. Fasshauer, Meshfree Approximation Methods with Matlab Interdisciplinary Mathematical Sciences - Vol. 6 World Scientific Publishers, Singapore, 2007




Membres

Universitaires:

Dr. Pedro González Casanova Henríquez

Coordinador del proyecto

Instituto de Matemáticas, UNAM

Dr. Dominique Apprato

(project Ecos Nord)

Université de Pau et des Pays de l'Adour | UPPA · Mathématiques Appliquées

Dr. Mohammed-Najib Benbourhim,

(Research-gate)

(Paul Sabatier University - Toulouse III, Toulouse, France).

Dr. Daniel Cervantes Cabrera

INFOTEC CDMX

Dr. Luis Miguel de la Cruz Salas,

Instituto de Geofísica, UNAM

Dr. Christian Gout,

(Laboratoire de Mathématiques de l'INSA de Rouen, EA 3226, France)

Dr. Jose Antonio Muñoz Gómez,

CUCSUR, Universidad de Guadalajara

Dr. Gustavo Rodríguez Gómez,

INAOE

Étudiants:

Doctorat

  • Tesis doctoral de Daniel Cervantes Cabrera, (UNAM), marzo de 2017.
  • M. en C. Jorge Zavaleta Sánchez,(UNAM)

Publications

D. Cervantes, P. Gonzalez Casanova, C. Gout, L. H. Juarez, R. Resendiz, Vector field approximation using radial basis functions, Journal of Computational and Applied Mathematics, Ed. Elsevier, Volume 240, 1 March 2013, Pages 163–173. R. Reséndiz, L. H. Juárez, P. González-Casanova, D. A. Cervantes, C. Gout, Numerical Study of Wind Field Adjustment with Radial Basis Functions, Environmental Science and Engineering, Ed. Springer, 2013, Part 4, 371-378. P. González-Casanova, J. Antonio Muñoz-Gómez, G. Rodríguez-Gómez, Node Adaptive Domain Decomposition Method by Radial Basis Functions. Journal of Numerical Methods for Partial Differential Equation , Volume 25 Issue 6, Pages 1482 - 1501, 2009. J. Antonio Muñoz-Gómez, P. González-Casanova, G. Rodríguez-Gómez, "Exponential Convergence of Multiquadric Collocation Method: a Numerical Study". Revista de Ingeniería, Investigacion y Tecnologia, Issue 3, pages 197-206, 2009. V. Pereyra, G. Scherer, P. Gonzalez Casanova, Radial function collocation solution of partial differential equations in irregular domains , International Journal of Computing Sciences and Mathematics, Vol. 1, No. 1, 2007. D. Cervantes, P. Gonzalez-Casanova, C. Gout, Seafloor surfaces approximation from rapidly varying bathymetric data using pre-processing, Pages 1742 - 1744, IEEE Transaction on Geosciences and Remote sensing, 2007. J. Antonio Muñoz-Gómez, P. González-Casanova, G. Rodríguez-Gómez, Domain Decomposition by Radial Basis Functions for Time Dependent Partial Differential Equations, Proccedings of the IASTED: International Conference on Advances in Computer Science and Technology, pp. 105-109, (2006). J. Antonio Muñoz-Gómez, P. González Casanova, G. Rodríguez-Gómez, Adaptive Node Refinement Collocation Method for Partial Differential Equations, ENC-06 Sixth Mexican International Conference on Computer Sciences, published by IEEE, pp 70-77, (2006). P. González-Casanova, C. Gout, Approximation Near the Boundary of a Radial Quasi interpolant in Polygonal Domains , Int. J. of Comp. Math. 81 (12), pp. 1145-1455 (2004). González-Casanova, Pedro. A Hilbertian approach to quasi-interpolation methods. Neural Parallel Sci. Computations, 5 (1997), no. 1-2, 257--274. Benbourhim, M. N.; González-Casanova, P. Generalized variational quasi interpolants in $(H\sp m\sb O(\Omega))\sp n$. Bol. Soc. Mat. Mexicana, 3 (1997) , no. 2, 347--354. Benbourhim and P. González Casanova,. Local Interpolation of non smooth functions , Proceedings of the Thirth International Conference on Approximation and Optimization, Puebla, México 1995, European Matematical Society, Electronic Library of Mathematics, en Julio de (1997). González-Casanova, Pedro; Wolf, Kurt Bernardo. Interpolation for solutions of the Helmholtz equation. Numer. Methods Partial Differential Equations ,11 (1995), no. 1, 77--91. Atteia M., Benbourhim M. N., González Casanova P. Quasi-interpolant Elastic Manifolds. Reporte de Investigación IIMAS-UNAM, Vol. 3, No. 25, (1993) . González Casanova P. and Álvarez R. Splines in Geophysics. Golden Anniversary Issue, Geophysics. Vol. 50 No. 12, December 1985. pp. 2831-2848.


Second Latin American SCAT Workchop,

Island Cozumel, México, November 2007.

Seminario de Modelación Matemática y Computacional

Instituto de Geofísica de la UNAM

Seminarios internos:

TUNA::RBF. Dr. Luís Miguel de la Cruz, 28 de agosto 2008. Lugar: DGSCA.


El uso de patrones para el diseño de software científico
Dr. Gustavo Rodríguez Gómez, 4 de septiembre 2008. Lugar: DGSCA.