Ce travail de recherche est une contribution au développement de la méthode Décomposition Propre Généralisée (PGD) à la résolution de problèmes de diffusion-réaction raides dédiés à la mécanique des matériaux. Ce type d’équations est notamment rencontré lors de l’oxydation des matériaux polymères et il est donc nécessaire de mettre en place un outil pour simuler ce phénomène afin de prédire numériquement le vieillissement de certains matériaux composites à matrice organique utilisés dans l’aéronautique. La méthode PGD a été choisie dans cette thèse car elle permet un gain en temps de calcul notable par rapport à la méthode des éléments finis. Néanmoins cette famille d’équations n’a jamais été traitée avec cette méthode. Cette dernière se résume à la recherche de solutions d’Équations aux Dérivées Partielles sous forme séparée. Dans le cas d’un problème 1D transitoire, cela revient à chercher la solution sous la forme d’une représentation séparée espace-temps. Dans le cadre de cette thèse, un outil numérique a été mis en place permettant une flexibilité telle que différents algorithmes peuvent être testés. La diffusion Fickienne 1D est tout d’abord évaluée avec en particulier une discussion sur l’utilisation d’un schéma de type Euler ou Runge-Kutta à pas adaptatif pour la détermination des fonctions temporelles. Le schéma de Runge-Kutta permet de réduire notablement le temps de calcul des simulations.Ensuite, la mise en place de l’outil pour les systèmes d’équation de type diffusion-réaction nécessite des algorithmes de résolution de systèmes non linéaires, couplés et raides. Pour cela, différents algorithmes ont été implémentés et discutés.Dans le cas d’un système non linéaire, l’utilisation de la méthode de Newton-Raphson dans les itérations pour la recherche du nouveau mode permet de réduire le temps de calcul en limitant le nombre de modes à considérer pour une erreur donnée. En ce qui concerne les couplages, deux stratégies de résolution ont été évaluées. Le couplage fort mène aux mêmes conclusions que dans le cas non linéaire. Les systèmes raides mais linéaires ont ensuite été traités en implémentant l’algorithme de Rosenbrock pour la détermination des fonctions temporelles. Cet algorithme permet contrairement à Euler et à Runge-Kutta de construire une solution avec un temps de calcul raisonnable liée à l’adaptation du maillage temporel sous-jacent à l’utilisation de cette méthode. La résolution d’un système d’équations de diffusion-réaction raides non linéaires utilisée pour la prédiction de l’oxydation d’un composite issu de la littérature a été testée en utilisant les différents algorithmes mis en place. Néanmoins, les non linéarités et la raideur du système génèrent des équations différentielles intermédiaires à coefficients variables pour lesquelles la méthode de Rosenbrock montre ses limites. Il sera donc nécessaire de tester ou développer d’autres algorithmes pour lever ce verrou.Mots / This work presents the development of the Proper Generalized Decomposition (PGD) method for solving stiff reaction-diffusion equations in the framework of mechanics of materials. These equations are particularly encountered in the oxidation of polymers and it is therefore necessary to develop a tool to simulate this phenomenon for example for the ageing of organic matrix composites in aircraft application. The PGD method has been chosen in this work since it allows a large time saving compared to the finite element method. However this family of equations has never been dealt with this method. The PGD method consists in approximating a solution of a Partial Differential Equation with a separated representation. The solution is sought under a space-time separated representation for a 1D transient equation.In this work, a numerical tool has been developed allowing a flexibility to test different algorithms. The 1D Fickian diffusion is first evaluated and two numerical schemes, Euler and Runge-Kutta adaptive methods, are discussed for the determination of the time modes. The Runge-Kutta method allows a large time saving. The implementation of the numerical tool for reaction-diffusion equations requires the use of specific algorithms dedicated to nonlinearity, couplingand stiffness. For this reason, different algorithms have been implemented and discussed. For nonlinear systems, the use of the Newton-Raphson algorithm at the level of the iterations to compute the new mode allows time saving by decreasing the number of modes required for a given precision. Concerning the couplings, two strategies have been evaluated. The strong coupling leads to the same conclusions as the nonlinear case. The linear stiff systems are then studied by considering a dedicated method, the Rosenbrock method, for the determination of the time modes. This algorithm allows time saving compared to the Runge-Kutta method. The solution of a realistic nonlinear stiff reaction-diffusionsystem used for the prediction of the oxidation of a composite obtained from the literature has been tested by using the various implemented algorithms. However, the nonlinearities and the stiffness of the system generate differential equations with variable coefficients for which the Rosenbrock method is limited. It will be necessary to test or develop other algorithms to overcome this barrier.

The solving of stiff systems is still a contemporary sophisticated problem. The basic problem is the absence of precise definition of stiff systems. A question is also how to detect the stiffness in a given system of differential equations. Implicit numerical methods are commonly used for solving stiff systems. The stability domains of these methods are relatively large but the order of them is low.   The thesis deals with numerical solution of ordinary differential equations, especially numerical calculations using Taylor series methods. The source of stiffness is analyzed and the possibility how to reduce stiffness             in systems of ordinary differential equations (ODEs) is introduced. The possibility of detection stiff systems using explicit Taylor series terms is analyzed. The stability domains of explicit and implicit Taylor series are presented. The solutions of stiff systems using implicit Taylor series method are presented in many examples. The multiple arithmetic must be used in many cases. The new suitable parallel algorithm based on implicit Taylor series method with recurrent calculation of Taylor series terms and Newton iteration method (ITMRN) is proposed.