<p align="justify">Ce travail consiste à étudier des systèmes constitués par un grand nombre de grains, auxquels de l’énergie cinétique est fournie, et à étudier leurs similarités et leurs différences avec des fluides traditionnels. Je me concentre principalement sur la nature de non-équilibre de ces fluides granulaires, en montrant que, même si les méthodes de méchanique statistique y sont applicables, leurs propriétés sont très différentes de celles de systèmes à l’équilibre ou proches de l’équilibre :</p> <ul><li>Les fluides granulaires présentent des phénomènes de transport qui n’ont pas d’équivalent dans des fluides moléculaires, tels qu’un couplage spécifique entre flux de chaleur et gradient de densité. <li>Leur distribution de vitesse est en général différente de la distribution de Maxwell-Boltzmann, et présente une surpopulation pour les grandes vitesses. <li>Dans le cas de mélanges, différentes espèces de grains sont en général caractérisées par des énergies cinétiques différentes, i.e. ces systèmes sont sujet à une non-equipartition de leur énergie. <li>Ces fluides ont tendance à former des inhomogénéités spatiales spontanément. Cette propriété est illustrée en étudiant l’expérience du Demon de Maxwell appliquée aux systèmes granulaires.</ul> <p align="justify">Chacune de ces particularités est discutée en détail dans des chapitres distincts, où l’on applique différentes méthodes de méchanique statistique (équation de Boltzmann, transition de phase, mean field models…) et où l’on vérifie les prédictions théoriques par simulations numériques (MD, Monte Carlo…).</p>

Some protein gels are products of the dairy industry and some are used as pH-sensitive gels for the controlled delivery of biologically active substances. To understand the dynamics of drug delivery it is very important to establish a mathematical model of protein gel swelling. This required the identification and integration of theory and equations from a wide range of topics. The aim of this research was to develop a mathematical model of transport in polyelectrolyte gels (using the example of β-lactoglobulin protein gels). A complete mathematical model of protein gel swelling was established. The swelling process of protein gels in this thesis involved multicomponent diffusion, chemical ionisation and mechanical deformation. Diffusion of electrolyte solutions through protein gels was modelled using the generalised Maxwell-Stefan (GMS) equation. The swelling pressure as a driving force in the GMS equation was described by rubber elasticity theory. Thermodynamic factors including the charged protein effect were considered in the GMS equation. The model included pH as a variable so it could be applied to both acidic and alkaline cases. The model yielded a set of partial differential equations with algebraic equations for which COMSOL was selected as the simulation software. Although it was found that COMSOL could not always solve the model equations, numerical solutions were obtained for several cases. The model predicted that the equilibrium swelling degree of the gel decreased with high concentration of salts in the bulk solution. The model also predicted that the non-ideal effects were not always small and they depended on the activity coefficients of the species. Satisfactory solutions could not be obtained for all cases using commercial software such as COMSOL Multiphysics. It was shown that COMSOL did not conserve mass but conservativeness was critical in this application because pH and hence the net protein charge is very sensitive to the mass of hydrogen present. In the future, research should be carried out to improve the pressure model in the GMS equation. Theoretical research on Manning condensation theory should be done to modify Manning’s model for more robust prediction of activities of water and ions with protein, and experiments should be done to validate the performance of the modified Manning model. Efforts should be made to write the programming code for a finite volume method to solve the system in three dimensions.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

Student Number :0008015M - MSc Dissertation - School of Electrical and Information Engineering - Faculty of Engineering and the Built Environment / Shaft voltages and currents are an unavoidable characteristic of rotating machines, though they are typically a nuisance this work shows that shaft signals can be used for fault diagnosis. This work focussed on shaft voltages present on synchronous generators. Measurements on a 4-pole generator found that the angular position of the shaft could be determined from the shaft voltage. An experimental 20 kVA 2-pole synchronous generator was designed and built which resembled a full-size 600 MVA turbo-generator. The effects of a static eccentricity on the shaft voltage were successfully determined firstly through FEM simulation and then verified with physical measurements. Shaft voltages can be used to diagnose static eccentricity, future work should investigate other faults. In addition, computer simulation was found to be effective and simulation and measurements of operating machines (such as turbo-generators) should be considered.

"This paper gives a brief account of the life of William Maxwell under such aspects as childhood, education, travels, writings and awards. A discussion of his novels presented in the order of publication will include critical reaction to Mr. Maxwell's skill as a writer and his themes. Following a summary and conclusions, an appendix will give a chronological list of Mr. Maxwell's books with sources of reviews, a list of his short stories and articles with annotations for those materials available, and a general bibliography of the major sources used"--Introduction. / Typescript. / "August, 1957." / "Submitted to the Graduate Council of Florida State University in partial fulfillment of the requirements for the degree of Master of Arts." / Advisor: Ruth H. Rockwood, Professor Directing Paper. / Includes bibliographical references (leaves 72-73).

Mak, Tsz Fan. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2013. / Includes bibliographical references (leaves 33-36). / Electronic reproduction. Hong Kong : Chinese University of Hong Kong, [2012] System requirements: Adobe Acrobat Reader. Available via World Wide Web. / Abstracts also in Chinese.

Zhang Kai. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2005. / Includes bibliographical references (leaves 99-102). / Abstracts in English and Chinese. / Chapter 1 --- Introduction --- p.6 / Chapter 1.1 --- The Generalized PML Theory --- p.6 / Chapter 1.1.1 --- Background --- p.6 / Chapter 1.1.2 --- Derivation --- p.8 / Chapter 1.1.3 --- Reflection Properties --- p.11 / Chapter 1.2 --- Unified Formulation --- p.12 / Chapter 1.2.1 --- "Face-, Edge- and Corner-PMLs" --- p.12 / Chapter 1.2.2 --- Unified PML Equations in 3D --- p.15 / Chapter 1.2.3 --- Unified PML Equations in 2D --- p.16 / Chapter 1.2.4 --- Examples of PML Formulations --- p.16 / Chapter 1.3 --- Inhomogeneous Initial Conditions --- p.23 / Chapter 2 --- Numerical Analysis of PMLs --- p.25 / Chapter 2.1 --- Continuous PMLs --- p.26 / Chapter 2.1.1 --- PMLs for Wave Equations --- p.27 / Chapter 2.1.2 --- Finite PMLs for Wave Equations --- p.31 / Chapter 2.1.3 --- Berenger's PMLs for Maxwell Equations --- p.33 / Chapter 2.1.4 --- Finite Berenger's PMLs for Maxwell Equations --- p.35 / Chapter 2.1.5 --- PMLs for Acoustic Equations --- p.38 / Chapter 2.1.6 --- Berenger's PMLs for Acoustic Equations --- p.39 / Chapter 2.1.7 --- PMLs for 1-D Hyperbolic Systems --- p.42 / Chapter 2.2 --- Discrete PMLs --- p.44 / Chapter 2.2.1 --- Discrete PMLs for Wave Equations --- p.44 / Chapter 2.2.2 --- Finite Discrete PMLs for Wave Equations --- p.51 / Chapter 2.2.3 --- Discrete Berenger's PMLs for Wave Equations --- p.53 / Chapter 2.2.4 --- Finite Discrete Berenger's PMLs for Wave Equations --- p.56 / Chapter 2.2.5 --- Discrete PMLs for 1-D Hyperbolic Systems --- p.58 / Chapter 2.3 --- Modified Yee schemes for PMLs --- p.59 / Chapter 2.3.1 --- Stability of the Yee Scheme for Wave Equation --- p.61 / Chapter 2.3.2 --- Decay of the Yee Scheme Solution to the Berenger's PMLs --- p.62 / Chapter 2.3.3 --- Stability and Convergence of the Yee Scheme for the Berenger's PMLs --- p.67 / Chapter 2.3.4 --- Decay of the Yee Scheme Solution to the Hagstrom's PMLs --- p.70 / Chapter 2.3.5 --- Stability and Convergence of the Yee Scheme for the Hagstrom's PMLs --- p.75 / Chapter 2.4 --- Modified Lax-Wendroff Scheme for PMLs --- p.80 / Chapter 2.4.1 --- Exponential Decays in Parabolic Equations --- p.80 / Chapter 2.4.2 --- Exponential Decays in Hyperbolic Equations --- p.82 / Chapter 2.4.3 --- Exponential Decays of Modified Lax-Wendroff Solutions --- p.86 / Chapter 3 --- Numerical Simulation --- p.93 / Bibliography --- p.99

在本文中，我們為了三維空間的馬克士威方程組(Maxwell’s equation)制定和分析了一套新種類的交錯間斷伽遼金(discontinuous Galerkin)方法，同時考慮了時間依賴性和時間諧波的馬克士威方程組。我們用了空間離散上交錯笛卡兒網格，這種方法具有許多良好的性質。首先，我們的方法所得出的數值解保留了電磁能量，並自動符合了高斯定律的離散版本。第二，質量矩陣是對角矩陣，從而時間推進是顯式和非常有效的。第三，我們的方法是高階準確，最佳收斂性在這裏會被嚴格地證明。第四，基於笛卡兒網格，它也很容易被執行，並可視為是典型的Yee’s Scheme的以及四邊形的邊有限元的推廣。最後，超收斂結果也會在這裏被證明。 / 在本文中，我們還提供了幾個數值結果驗證了理論的陳述。我們計算了時間依賴性和時間諧波的馬克士威方程組數值收斂結果。此外，我們計算時間諧波馬克士威方程組特徵值問題的數值特徵值，並與理論特徵值比較結果。最後，完美匹配層(Perfect Matching Layer)吸收邊界的問題也有實行其數值結果。 / We develop and analyze a new type of staggered discontinuous Galerkin methods for the three dimensional Maxwell’s equations in this paper. Both time-dependent and time-harmonic Maxwell’s equations are considered. The spatial discretization is based on staggered Cartesian grids which possess many good properties. First of all, our method has the advantages that the numerical solution preserves the electromagnetic energy and automatically fulfills a discrete version of the Gauss law. Second, the mass matrices are diagonal, thus time marching is explicit and is very efficient. Third, our method is high order accurate and the optimal order of convergence is rigorously proved. Fourth, it is also very easy to implement due to its Cartesian structure and can be regarded as a generalization of the classical Yee’s scheme as well as the quadrilateral edge finite elements. Lastly, a superconvergence result, that is the convergence rate is one order higher at interpolation nodes, is proved. / In this paper, we also provide several numerical results to verify the theoretical statements. We compute the numerical convergence order using L2-norm and discrete-norm respectively for both the time-dependent and time-harmonic Maxwell’s equations. Also, we compute the numerical eigenvalues for the time-harmonic eigenvalue problem and compare the result with the theoretical eigenvalues. Lastly, applications to problems in unbounded domains with the use of PML are also presented. / Detailed summary in vernacular field only. / Detailed summary in vernacular field only. / Yu, Tang Fei. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2012. / Includes bibliographical references (leaves 46-49). / Abstracts also in Chinese. / Chapter 1 --- Introduction and Model Problems --- p.1 / Chapter 2 --- Staggered DG Spaces --- p.4 / Chapter 2.1 --- Review on Gauss-Radau and Gaussisan points --- p.5 / Chapter 2.2 --- Basis functions --- p.6 / Chapter 2.3 --- Finite Elements space --- p.7 / Chapter 3 --- Method derivation --- p.14 / Chapter 3.1 --- Method --- p.14 / Chapter 3.2 --- Time discretization --- p.17 / Chapter 4 --- Energy conservation and Discrete Gauss law --- p.19 / Chapter 4.1 --- Energy conservation --- p.19 / Chapter 4.2 --- Discrete Gauss law --- p.22 / Chapter 5 --- Error analysis --- p.24 / Chapter 6 --- Numerical examples --- p.29 / Chapter 6.1 --- Convergence tests --- p.30 / Chapter 6.2 --- Diffraction by a perfectly conducting object --- p.30 / Chapter 6.3 --- Perfectly matched layers --- p.37 / Chapter 7 --- Time Harmonic Maxwell’s equations --- p.40 / Chapter 7.1 --- Model Problems --- p.40 / Chapter 7.2 --- Numerical examples --- p.40 / Chapter 7.2.1 --- Convergence tests --- p.41 / Chapter 7.2.2 --- Eigenvalues tests --- p.41 / Chapter 8 --- Conclusion --- p.45 / Bibliography --- p.46

In the first part, an efficient and reliable a posteriori error estimate is derived for solving three-dimensional static Maxwell's equations based on the lowest order edge elements of the first family. We propose an adaptive finite element method and establish convergence of the adaptive scheme in energy norm under a restriction on the initial mesh size. Any prescribed error tolerance is thus achieved in a finite number of steps. For discretization based on the lowest order edge elements of the second family, a similar adaptive method is designed which guarantees convergence without any initial mesh size restriction. The proofs rely mainly on error and oscillation reduction estimates as well as the Galerkin orthogonality of the edge element approximation. For time-dependent Maxwell's equations, we deduce an efficient and reliable a posteriori error estimate, upon which an adaptive finite element method is built. / In this thesis, we will address three typical problems with discontinuous coefficients in a general Lipschitz polyhedral domain, which are often encountered in numerical simulation of electromagnetism. / The second part deals with a saddle point problem arising from Maxwell's equations. We present an adaptive finite element method on the basis of the lowest order edge elements of the first family and prove its convergence. The main ingredients of the proof are a novel quasi-orthogonality, which replaces the usual Pythagoras relation, which fails in this case, all error reduction depending on an efficient and reliable a posteriori error estimate and an oscillation reduction. We show that this adaptive scheme is a contraction for the sum of some energy error plus the oscillation. Likewise, the above result is generalized to the discretization by the lowest order edge elements of the second family. / We introduce in the third part an adaptive finite element method for solving the eigenvalue problem of the Maxwell system based on an inverse iterative method. By modifying the exact inverse iteration algorithm involving an inner saddle point solver, we construct an adaptive inverse iteration finite element algorithm, which consists of an inexact inner adaptive procedure for a discrete mixed formulation in place of the original saddle point problem. An efficient and reliable a posteriori error estimate is obtained and the convergence of the inner adaptive method is proved. In addition, the important convergence property of the algorithm is studied, which ensures the errors between true solutions (eigenfunction and eigenvalue) and iterative ones to fall below any given tolerance within a finite number of iterations. / Xu, Yifeng. / "June 2007." / Adviser: Jun Zou. / Source: Dissertation Abstracts International, Volume: 69-01, Section: B, page: 0357. / Thesis (Ph.D.)--Chinese University of Hong Kong, 2007. / Includes bibliographical references (p. 166-175). / Electronic reproduction. Hong Kong : Chinese University of Hong Kong, [2012] System requirements: Adobe Acrobat Reader. Available via World Wide Web. / Electronic reproduction. [Ann Arbor, MI] : ProQuest Information and Learning, [200-] System requirements: Adobe Acrobat Reader. Available via World Wide Web. / Abstracts in English and Chinese. / School code: 1307.

Les méthodes intégrales pour résoudre des EDP, et en particulier le système de Maxwell, sont bien connues depuis environ vingt ans. Après discrétisation par éléments finis, un système linéaire plein apparaît, ce qui rend toute implémentation numérique difficile voire impossible. Pour les opérateurs d'ordre positif, quelques travaux ont été menés avec succès pour rendre creuse la matrice du système discret. Quelques difficultés restaient pour le problème de Maxwell: espace(s) d'énergie, présence d'un opérateur d'ordre négatif, et donc choix des ondelettes pour la résolution. Dans cette thèse, je donne une méthode pour ramener le système de Maxwell, issu d'un problème de diffraction en régime harmonique, à une étude sur des espaces de Sobolev classiques définis sur une surface, en utilisant des décompositions de Hodge. Je donne aussi une méthode de compression pourvu que les ondelettes vérifient certaines conditions (moments nuls, stabilité). La méthode de compression donnée fonctionne même avec des ondelettes formées à partir de polynômes de degré un, malgré la présence d'un opérateur d'ordre négatif, sans perturber des taux de convergence optimaux. L'analyse a été faite sur une surface fermée (sans bord) régulière simplement connexe, puis sur une partie à bord polygonal d'une telle surface (plaque ouverte). Les espaces fonctionnels et la compression de matrice, bien plus compliqués dans ce dernier cas, ont été étudiés en détail.