Spelling suggestions: "subject:"stica ODEs"" "subject:"stick ODEs""
1 |
Exponential integrators: tensor structured problems and applicationsCassini, Fabio 21 April 2023 (has links)
The solution of stiff systems of Ordinary Differential Equations (ODEs), that typically arise after spatial discretization of many important evolutionary Partial Differential Equations (PDEs), constitutes a topic of wide interest in numerical analysis. A prominent way to numerically integrate such systems involves using exponential integrators. In general, these kinds of schemes do not require the solution of (non)linear systems but rather the action of the matrix exponential and of some specific exponential-like functions (known in the literature as φ-functions). In this PhD thesis we aim at presenting efficient tensor-based tools to approximate such actions, both from a theoretical and from a practical point of view, when the problem has an underlying Kronecker sum structure. Moreover, we investigate the application of exponential integrators to compute numerical solutions of important equations in various fields, such as plasma physics, mean-field optimal control and computational chemistry. In any case, we provide several numerical examples and we perform extensive simulations, eventually exploiting modern hardware architectures such as multi-core Central Processing Units (CPUs) and Graphic Processing Units (GPUs). The results globally show the effectiveness and the superiority of the different approaches proposed.
|
2 |
Numerical algorithms for differential equations with periodicityMontanelli, Hadrien January 2017 (has links)
This thesis presents new numerical methods for solving differential equations with periodicity. Spectral methods for solving linear and nonlinear ODEs, linear ODE eigenvalue problems and linear time-dependent PDEs on a periodic interval are reviewed, and a novel approach for computing multiplication matrices is presented. Choreographies, periodic solutions of the n-body problem that share a common orbit, are computed for the first time to high accuracy using an algorithm based on approximation by trigonometric polynomials and optimization techniques with exact gradient and exact Hessian matrix. New choreographies in spaces of constant curvature are found. Exponential integrators for solving periodic semilinear stiff PDEs in 1D, 2D and 3D periodic domains are reviewed, and 30 exponential integrators are compared on 11 PDEs. It is shown that the complicated fifth-, sixth- and seventh-order methods do not really outperform one of the simplest exponential integrators, the fourth-order ETDRK4 scheme of Cox and Matthews. Finally, algorithms for solving semilinear stiff PDEs on the sphere with spectral accuracy in space and fourth-order accuracy in time are proposed. These are based on a new variant of the double Fourier sphere method in coefficient space and standard implicit-explicit time-stepping schemes. A comparison is made against exponential integrators and it is found that implicit-explicit schemes perform better. The algorithms described in each chapter of this thesis have been implemented in MATLAB and made available as part of Chebfun.
|
3 |
Schémas d'ordre élevé pour des simulations réalistes en électrophysiologie cardiaque / High order schemes for realistic simulations in cardiac electrophysiologyDouanla Lontsi, Charlie 15 November 2017 (has links)
Les simulations numériques réalistes en électrophysiologie cardiaque ont un coût de calcul extrêmement élevé. Ce coût s’explique en grande partie par la raideur, à la fois en temps et en espace, d’une onde de « potentiel d’action » (PA). Par ailleurs, les phénomènes observés sont très instationnaires et s’étudient en temps long. Une description précise de la dynamique des PA est cruciale pour construire des modèles numériques pertinents d’un point de vue médical ou clinique. Cet aspect fondamental ne peut être contourné dans les études numériques réalistes.La raideur de l’onde de PA ne peut être captée numériquement qu’en ayant recours à des maillages très fins. Ces maillages très fins induisent un coût de calcul très important, et introduisent aussi des erreurs supplémentaires : les systèmes linéaires à résoudre deviennent très mal conditionnés. Au final, les erreurs numériques peuvent être particulièrement grandes dans les simulations alors que leur contrôle est évidemment essentiel pour assurer la fiabilité des résultats. Jusqu’à présent, très peu de résultats sont disponibles pour assurer cette fiabilité. Dans les faits, les erreurs sont la plupart du temps contrôlées par des procédés empiriques. Il existe quelques résultats théoriques étudiant la convergence et la stabilité des schémas numériques associés. En pratique, en plus d'avoir un contrôle de l'erreur sur le potentiel, il est aussi nécessaire d'avoir un contrôle de l’erreur sur des quantités macroscopiques décrivant la dynamique de l’onde de PA : temps d’activation, durée du PA, propriétés de restitution... Ces quantités ont en effet une interprétation physiologique qui permet de caractériser le caractère arythmogène des tissus.Les modèles sont des systèmes d’EDP de réaction-diffusion couplés avec des systèmes d’équations différentielles pouvant être très raides, les modèles ioniques. Ils sont actuellement discrétisés par éléments finis conforme (Lagrange) et par des schémas en temps d’ordre un ou deux. Dans ce travail, nous concevons et évaluons l’intérêt d'utiliser des méthodes d’ordre supérieure pour ces systèmes. Parallèlement nous introduisons d'une part une nouvelle classe de schémas appelé schémas exponentiel Adams Bashforth intégral (IEAB), et d'autre part des schémas Rush Larsen (RL) d'ordre élevé. Ces nouveaux schémas sont des schémas multipas de type exponentiels. Nous montrons qu'ils possèdent des bonnes propriétés de stabilité et permettent de faire face efficacement à la raideur des modèles ioniques. Les schémas que nous proposons sont comparés numériquement (en terme de précision, coût en temps de calcul et stabilité) à plusieurs schémas classiques, ainsi qu'aux schémas exponentiels (RL1, RL2) communément utilisés pour des simulations en électrophysiologie cardiaque. Nous proposons des techniques permettant de calculer avec précision les quantités d’intérêts cliniques (temps d’activation, de récupération, durée du potentiel d’action). Des résultats théoriques de convergence en temps et de convergence globale (espace et temps) sont énoncés et prouvés. Ces résultats sont ensuite illustrés numériquement à travers le modèle monodomaine et les modèles ioniques de Beeler Reuter, de Ten Tusscher et al. L’intérêt d'utiliser des schémas d'ordre élevés est aussi évalué sur des ondes spirales en 2D et 3D. / Realistic numerical simulations in cardiac electrophysiology have a computational cost of extremely high. This cost is largely explained by the stiffness both in time and space, of the action potential (AP) wave. Moreover, the observed phenomena are very unsteady and are studied in long time. A precise description of the dynamic of AP is crucial for constructing relevant numerical models, from a medical or clinical perspective. This fundamental aspect can not be circumvented in realistic numerical studies.The stiffness of AP wave can only be captured numerically, by using very fine meshes. In addition to the high computational cost, these very fine meshes also introduce additional errors : the linear systems to solve become very badly conditioned. In the end, the numerical errors can be particularly large whereas their control is obviously essential to ensure the reliability of the results. So far very few results are available to ensure this reliability. In practice, the errors are mostly controlled by empirical processes. In practice, in addition of having a control of the error on the potential, it is also necessary to have an error control on macroscopic quantities describing the dynamics of the AP wave : activation time, AP duration, properties of restitution ... These quantities have indeed a physiological interpretation which allows to characterize the arrhythmogenic character of the tissues.The models are systems of reaction diffusion PDE coupled with systems of differential equations that can be very stiffs (ionic models). They are currently discretized by conforming finite elements (Lagrange finite elements methods) and by schemes in time of order one or two. In this work, we design and evaluate the interest of using higher order methods for these systems. At the same time, we introduce on the one hand, a new class of schemes called Integral Exponential Adams Bashforth (IEAB) schemes and, on the other hand, high order Rush Larsen (RL) schemes. These new schemes are exponential time-stepping schemes. We show that they have good stability properties and can efficiently cope with the stiffness of ionic models. The schemes we propose are numerically compared (in terms of accuracy, CPU time and stability) with several classical schemes, as well as with the exponential schemes (RL1, RL2), commonly used for cardiac electrophysiology simulations. We propose good techniques for accurately calculating quantities of clinical interest (activation time, recovery time, duration of action potential). Theoretical results of convergence in time and global convergence (in space and time) are stated and proved. These results are then illustrated numerically through the monodomain model and the ionic models of Beeler Reuter, Ten Tusscher et al. The advantage of using high order schemes is also evaluated on spiral waves in 2D and 3D.
|
Page generated in 0.0269 seconds