Cette thèse est consacrée à l’étude d’expansions de certaines structures algébriques et leur place dans la classification modèle-théorique des structures, initiée par Shelah. La première partie aborde de manière abstraite l’expansion d’une théorie par un prédicat aléatoire –ou générique– pour une sous-structure modèle d’un réduit de la théorie. Nous éla- borons un critère pour l’existence d’une telle expansion, qui est vérifié pour certaines théories de structures algébriques. En particulier, nous montrons l’existence de sous-groupes additifs génériques pour certaines théories de corps, ainsi que de sous-groupes multiplicatifs génériques pour les corps algébriquement clos en toute caractéristique. Nous étudions aussi la conservation de diverses notions de néostabilité, en particulier nous montrons que cette expansion préserve la propriété NSOP 1 , mais en général ne préserve pas la simplicité. Nous produisons par cette construction de nouveaux exemples de structures NSOP 1 non simples, et faisons une étude toute particulière de l’une d’entre elles : l’expansion d’un corps algébriquement clos de caractéristique positive par un sous-groupe additif générique. La deuxième partie étudie les expansions du groupe des entiers par des valuations p-adiques. Nous montrons l’élimination des quantificateurs dans un langage naturel et calculons le dp-rang d’une telle expansion : il est égal au nombre de valuations considérées. L’expansion du groupe des entiers par une seule valuation p-adique est donc une nouvelle expansion dp-minimale du groupe des entiers. Enfin, nous montrons que cette dernière n’admet pas de structures intermédiaires : tout ensemble définissable dans l’expansion est soit définissable dans le groupe des entiers, soit capable de “reconstruire” la valuation en utilisant seulement la structure additive / This thesis is concerned with the expansions of some algebraic structures and their fit in Shelah’s classification landscape. The first part deals with the expansion of a theory by a random –or generic– predicate for a substructure model of a reduct of the theory. We describe a setup allowing such an expansion to exist, which is suitable for several algebraic structures. In particular, we obtain the existence of additive generic subgroups of some theories of fields and multiplicative generic subgroups of algebraically closed fields in all characteristic. We also study the preservation of certain neostability notions, for instance, the NSOP 1 property is preserved but the simplicity is not in general. Thus, this construction produces new examples of NSOP 1 not simple theories, and we study in depth a particular example: the expansion of an algebraically closed field of positive characteristic by a generic additive subgroup. The second part studies expansions of the groups of integers by p-adic valuations. We prove quantifier elimination in a natural language and compute the dp-rank of these expansions: it equals the number of distinct p-adic valuations considered. Thus, the expansion of the integers by one p-adic valuation is a new dp-minimal expansion of the group of integers. Finally, we prove that the latter expansion does not admit intermediate structures: any definable set in the expansion is either definable in the group structure or is able to "reconstruct" the valuation using only the group operation

De nombreux problèmes de calcul numérique demandent parfois à effectuer des calculs très précis. L'étude desystèmes dynamiques chaotiques fournit des exemples très connus: la stabilité du système solaire ou l’itération à longterme de l'attracteur de Lorenz qui constitue un des premiers modèles de prédiction de l'évolution météorologique. Ons'intéresse aussi aux problèmes d'optimisation semi-définie positive mal-posés qui apparaissent dans la chimie oul'informatique quantique.Pour tenter de résoudre ces problèmes avec des ordinateurs, chaque opération arithmétique de base (addition,multiplication, division, racine carrée) demande une plus grande précision que celle offerte par les systèmes usuels(binary32 and binary64). Il existe des logiciels «multi-précision» qui permettent de manipuler des nombres avec unetrès grande précision, mais leur généralité (ils sont capables de manipuler des nombres de millions de chiffres) empêched’atteindre de hautes performances. L’objectif majeur de cette thèse a été de développer un nouveau logiciel à la foissuffisamment précis, rapide et sûr : on calcule avec quelques dizaines de chiffres (quelques centaines de bits) deprécision, sur des architectures hautement parallèles comme les processeurs graphiques et on démontre des bornesd'erreur afin d'être capables d’obtenir des résultats certains. / Many numerical problems require some very accurate computations. Examples can be found in the field ofdynamical systems, like the long-term stability of the solar system or the long-term iteration of the Lorenz attractor thatis one of the first models used for meteorological predictions. We are also interested in ill-posed semi-definite positiveoptimization problems that appear in quantum chemistry or quantum information.In order to tackle these problems using computers, every basic arithmetic operation (addition, multiplication,division, square root) requires more precision than the ones offered by common processors (binary32 and binary64).There exist multiple-precision libraries that allow the manipulation of very high precision numbers, but their generality(they are able to handle numbers with millions of digits) is quite a heavy alternative when high performance is needed.The major objective of this thesis was to design and develop a new arithmetic library that offers sufficient precision, isfast and also certified. We offer accuracy up to a few tens of digits (a few hundred bits) on both common CPU processorsand on highly parallel architectures, such as graphical cards (GPUs). We ensure the results obtained by providing thealgorithms with correctness and error bound proofs.

Human cancers exhibit the common phenotype of elevated glycolytic metabolism, which causes acidification of the tissue microenvironment and may facilitate tumour invasion. In this thesis, we use mathematical models to address a series of open problems underlying the glycolytic tumour phenotype and its attendant acidity. We first explore tissue-scale consequences of metabolically-derived acid. Incorporating more biological detail into a canonical model of acidity at the tumour-host interface, we extend the range of tumour behaviours captured by the modelling framework. We then carry out an asymptotic travelling wave analysis to express invasive tumour properties in terms of fundamental parameters, and find that interstitial gaps between an advancing tumour and retreating healthy tissue, characteristic of aggressive invasion and comprising a controversial feature of the original model, are less significant under our generalised formulation. Subsequently, we evaluate a potential role of lactate---historically assumed to be a passive byproduct of glycolytic metabolism---in a perfusion-dependent metabolic symbiosis that was recently proposed as a beneficial tumour behaviour. Upon developing a minimal model of dual glucose-lactate consumption in vivo and employing a multidimensional sensitivity analysis, we find that symbiosis may not be straightforwardly beneficial for our model tumour. Moreover, new in vitro experiments, carried out by an experimental collaborator, place U87 glioblastoma tumours in a weakly symbiotic parameter regime despite their clinical malignancy. These results suggest that intratumoural metabolic cooperation is unlikely to be an important role for lactate. Finally, we examine the complex pH regulation system that governs expulsion of metabolically derived acid loads across tumour cell membranes. This system differs from the healthy system by expression of only a few key proteins, yet its dynamics are non-intuitive in the crowded and poorly perfused in vivo environment. We systematically develop a model of tumour pH regulation, beginning with a single-cell scenario and progressing to a spheroid, within a Bayesian framework that incorporates information from in vitro data contributed by a second experimental collaborator. We predict that a net effect of pH regulation is a straightforward transmembrane pH gradient, but also that existing treatments are unable to disrupt the system strongly enough to cause tumour cell death. Taken together, our models help to elucidate previously unresolved features of glycolytic tumour metabolism, and illustrate the utility of a combined mathematical, statistical, and experimental approach for testing biological hypotheses. Opportunities for further investigation are discussed.

Large amplitude internal solitary waves in the coastal ocean are commonly modelled with the Korteweg-de Vries (KdV) equation or a closely related evolution equation. The characteristic feature of these models is the solitary wave solution, and it is well documented that these provide the basic paradigm for the interpretation of oceanic observations. However, often internal waves in the ocean survive for several inertial periods, and in that case, the KdV equation is supplemented with a linear non-local term representing the effects of background rotation, commonly called the Ostrovsky equation. This equation does not support solitary wave solutions, and instead a solitary-like initial condition collapses due to radiation of inertia-gravity waves, with instead the long-time outcome typically being an unsteady nonlinear wave packet. The KdV equation and the Ostrovsky equation are formulated on the assumption that only a single vertical mode is used. In this thesis we consider the situation when two vertical modes are used, due to a near-resonance between their respective linear long wave phase speeds. This phenomenon can be described by a pair of coupled Ostrovsky equations, which is derived asymptotically from the full set of Euler equations and solved numerically using a pseudo-spectral method. The derivation of a system of coupled Ostrovsky equations is an important extension of coupled KdV equations on the one hand, and a single Ostrovsky equation on the other hand. The analytic structure and dynamical behaviour of the system have been elucidated in two main cases. The first case is when there is no background shear flow, while the second case is when the background state contains current shear, and both cases lead to new solution types with rich dynamical behaviour. We demonstrate that solitary-like initial conditions typically collapse into two unsteady nonlinear wave packets, propagating with distinct speeds corresponding to the extremum value in the group velocities. However, a background shear flow allows for several types of dynamical behaviour, supporting both unsteady and steady nonlinear wave packets, propagating with the speeds which can be predicted from the linear dispersion relation. In addition, in some cases secondary wave packets are formed associated with certain resonances which also can be identified from the linear dispersion relation. Finally, as a by-product of this study it was shown that a background shear flow can lead to the anomalous version of the single Ostrovsky equation, which supports a steady wave packet.

There exists substantial applications motivating the study of nonlinear longitudinal wave propagation in layered (or laminated) elastic waveguides, in particular within areas related to non-destructive testing, where there is a demand to understand, reinforce, and improve deformation properties of such structures. It has been shown [76] that long longitudinal waves in such structures can be accurately modelled by coupled regularised Boussinesq (cRB) equations, provided the bonding between layers is sufficiently soft. The work in this thesis firstly examines the initial-value problem (IVP) for the system of cRB equations in [76] on the infinite line, for localised or sufficiently rapidly decaying initial conditions. Using asymptotic multiple-scales expansions, a nonsecular weakly nonlinear solution of the IVP is constructed, up to the accuracy of the problem formulation. The asymptotic theory is supported with numerical simulations of the cRB equations. The weakly nonlinear solution for the equivalent IVP for a single regularised Boussinesq equation is then constructed; constituting an extension of the classical d'Alembert's formula for the leading order wave equation. The initial conditions are also extended to allow one to separately specify an O(1) and O(ε) part. Large classes of solutions are derived and several particular examples are explicitly analysed with numerical simulations. The weakly nonlinear solution is then improved by considering the IVP for a single regularised Boussinesq-type equation, in order to further develop the higher order terms in the solution. More specifically, it enables one to now correctly specify the higher order term's time dependence. Numerical simulations of the IVP are compared with several examples to justify the improvement of the solution. Finally an asymptotic procedure is developed to describe the class of radiating solitary wave solutions which exist as solutions to cRB equations under particular regimes of the parameters. The validity of the analytical solution is examined with numerical simulations of the cRB equations. Numerical simulations throughout this work are derived and implemented via developments of several finite difference schemes and pseudo-spectral methods, explained in detail in the appendices.

The main motivation behind writing this thesis was to construct numerical methods to approximate solutions to differential equations driven by rough paths, where the solution is considered in the rough path-sense. Rough paths of inhomogeneous degree of smoothness as driving noise are considered. We also aimed to find applications of these numerical methods to stochastic differential equations. After sketching the core ideas of the Rough Paths Theory in Chapter 1, the versions of the core theorems corresponding to the inhomogeneous degree of smoothness case are stated and proved in Chapter 2 along with some auxiliary claims on the continuity of the solution in a certain sense, including an RDE-version of Gronwall's lemma. In Chapter 3, numerical schemes for approximating solutions to differential equations driven by rough paths of inhomogeneous degree of smoothness are constructed. We start with setting up some principles of approximations. Then a general class of local approximations is introduced. This class is used to construct global approximations by pasting together the local ones. A general sufficient condition on the local approximations implying global convergence is given and proved. The next step is to construct particular local approximations in finite dimensions based on solutions to ordinary differential equations derived locally and satisfying the sufficient condition for global convergence. These local approximations require strong conditions on the one-form defining the rough differential equation. Finally, we show that when the local ODE-based schemes are applied in combination with rough polynomial approximations, the conditions on the one-form can be weakened. In Chapter 4, the results of Gyurko & Lyons (2010) on path-wise approximation of solutions to stochastic differential equations are recalled and extended to the truncated signature level of the solution. Furthermore, some practical considerations related to the implementation of high order schemes are described. The effectiveness of the derived schemes is demonstrated on numerical examples. In Chapter 5, the background theory of the Kusuoka-Lyons-Victoir (KLV) family of weak approximations is recalled and linked to the results of Chapter 4. We highlight how the different versions of the KLV family are related. Finally, a numerical evaluation of the autonomous ODE-based versions of the family is carried out, focusing on SDEs in dimensions up to 4, using cubature formulas of different degrees and several high order numerical ODE solvers. We demonstrate the effectiveness and the occasional non-effectiveness of the numerical approximations in cases when the KLV family is used in its original version and also when used in combination with partial sampling methods (Monte-Carlo, TBBA) and Romberg extrapolation.

We study the convergence of a nonlinear approximation method introduced in the engineering literature for the numerical solution of a high-dimensional Fokker--Planck equation featuring in Navier--Stokes--Fokker--Planck systems that arise in kinetic models of dilute polymers. To do so, we build on the analysis carried out recently by Le~Bris, Leli\`evre and Maday (Const. Approx. 30: 621--651, 2009) in the case of Poisson's equation on a rectangular domain in $\mathbb{R}^2$, subject to a homogeneous Dirichlet boundary condition, where they exploited the connection of the approximation method with the greedy algorithms from nonlinear approximation theory explored, for example, by DeVore and Temlyakov (Adv. Comput. Math. 5:173--187, 1996). We extend the convergence analysis of the pure greedy and orthogonal greedy algorithms considered by Le~Bris, Leli\`evre and Maday to the technically more complicated situation of the elliptic Fokker--Planck equation, where the role of the Laplace operator is played out by a high-dimensional Ornstein--Uhlenbeck operator with unbounded drift, of the kind that appears in Fokker--Planck equations that arise in bead-spring chain type kinetic polymer models with finitely extensible nonlinear elastic potentials, posed on a high-dimensional Cartesian product configuration space $\mathsf{D} = D_1 \times \dotsm \times D_N$ contained in $\mathbb{R}^{N d}$, where each set $D_i$, $i=1, \dotsc, N$, is a bounded open ball in $\mathbb{R}^d$, $d = 2, 3$. We exploit detailed information on the spectral properties and elliptic regularity of the Ornstein--Uhlenbeck operator to give conditions on the true solution of the Fokker--Planck equation which guarantee certain rates of convergence of the greedy algorithms. We extend the analysis to discretized versions of the greedy algorithms.

In studies of the ocean it has become conventional to retain only the component of the Coriolis force associated with the radial component of the Earth’s rotation vector, the so-called “traditional approximation”. We investigate the role of the “non-traditional” component of the Coriolis force, corresponding to the non-radial component of the rotation vector, in transporting abyssal waters across the equator. We first derive a non-traditional generalisation of the multi-layer shallow water equations, which describe the flow of multiple superposed layers of inviscid, incompressible fluid with constant densities over prescribed topography in a rotating frame. We derive these equations both by averaging the three-dimensional governing equations over each layer, and via Hamilton’s principle. The latter derivation guarantees that conservation laws for mass, momentum, energy and potential vorticity are preserved. Within geophysically realistic parameters, including the complete Coriolis force modifies the domain of hyperbolicity of the multi-layer equations by no more than 5%. By contrast, long linear plane waves exhibit dramatic structural changes due to reconnection of the surface and internal wave modes in the long-wave limit. We use our non-traditional shallow water equations as an idealised model of an abyssal current flowing beneath a less dense upper ocean. We focus on the Antarctic Bottom Water, which crosses the equator in the western Atlantic ocean, where the bathymetry forms an almost-westward channel. Cross-equatorial flow is strongly constrained by potential vorticity conservation, which requires fluid to acquire a large relative vorticity in order to move between hemispheres. Including the complete Coriolis force accounts for the fact that fluid crossing the equator in an eastward/westward channel experiences a smaller change in angular momentum, and therefore acquires less relative vorticity. Our analytical and numerical solutions for shallow water flow over idealised channel topography show that the non-traditional component of the Coriolis force facilitates cross-equatorial flow through an almost-westward channel.

La première réalisation de ce travail est la construction unifiée et automatique d’un milieu continu équivalent à un treillis de poutres, dans le domaine élastique, en adoptant un modèle de poutres de Bernoulli. Une extension a été réalisée au domaine plastique, selon un algorithme de suivi de la loi de comportement après écrouissage. Suivant l’ordre des développements asymptotiques choisi, on obtient pour le comportement élastique un milieu continu classique ou micropolaire. On se restreint dans ce dernier cas aux treillis à cellules élémentaires centro-symétriques. Les codes de calculs obtenus fournissent de façon automatique les lois de comportement effectives et les modules mécaniques homogénéisés. Une grande variété de treillis, existants ou originaux, a été étudiée. Les résultats ont été systématiquement comparés aux données de la littérature et vérifiés par des simulations éléments finis avec une bonne concordance. La méthode utilisée montre également une capacité à prédire et comprendre le comportement atypique de certains treillis dits auxétiques présentant des coefficients de contraction négatifs. L’homogénéisation dans le domaine plastique a été limitée aux treillis à dominante extensionnelle. Le domaine de résistance élastique a été construit pour différents treillis, et un algorithme d’évolution du comportement avec écrouissage, de type retour-radial a été conçu et implémenté dans un code dédié. Un modèle de poutre élastoplastique à écrouissage isotrope est utilisé. L’application de l’algorithme à une simulation de charge-décharge montre une bonne concordance entre le treillis homogénéisé et les simulations éléments finis / The first achievement of this work is to construct a unified and effective continuum equivalent to a lattice of beams, in the elastic domain, using a Bernoulli beam model. An extension has been done to calculate the elastic domain resistance of such lattices and to build an algorithm for monitoring the constitutive law taking into account work hardening. The choice of the asymptotic expansions leads to a classical continuous or to a micropolar elastic continuum. We restrict in this last case our study to lattices with centro-symmetric unit cells. The numerical codes developed provide the stress-strain relationship and the effective mechanical moduli. A wide variety of trusses has been studied, either existing or original, including typical geometries of foams and various auxetic lattices, exhibiting negative contraction coefficients. The results were systematically compared with data from literature and verified by finite element simulations with a good agreement. The homogenization in the plastic range has been limited to stretching dominated lattices. The equilibrium equations of the discrete asymptotic homogenization have been used to automatically obtain the elastic resistance domain for several trusses, and a return-mapping algorithm for the follow up of the stress-strain relationship including hardening has been conceived and implemented in a dedicated code. An isotropic hardening elastoplastic model of the beam has been used. The application of the algorithm to the simulation of a loading-unloading cycle shows a good agreement between the homogenized lattice and finite element simulations

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