Au cœur des moyens de transport, de transformation d'énergie, et de biens d'équipements, les machines tournantes peuvent avoir des comportements dynamiques complexes dus à de multiples sources de non linéarités liées aux paliers hydrodynamiques, à la présence de fissures, aux touches rotor-stator, ... Des phénomènes comme les décalages fréquentiels et donc de vitesses critiques, les cycles d'hystérésis avec sauts d'amplitudes, le changement brutal du contenu fréquentiel des réponses, sont des expressions de ces comportements. Résoudre les équations du mouvement induites par des modélisations avec des éléments finis de type poutre ou volumique, pour calculer les réponses à des sollicitations diverses (comme le balourd ou le poids propre), est réalisable avec des méthodes d'intégration pas à pas dans le temps mais au prix de temps de calcul prohibitifs. Cela devient particulièrement préjudiciable au stade du pré-dimensionnement où il est nécessaire de réaliser rapidement des études paramétriques. Aussi une alternative intéressante est de mettre en {\oe}uvre une méthode numérique, à la fois générale et efficace pour analyser la réponse non linéaire des rotors en régime stationnaire. La démarche proposée combine, dans un premier temps, la méthode de la balance harmonique (HBM) et la technique de bascule Temps-Fréquence (AFT) afin d'obtenir rapidement dans le domaine fréquentiel les réponses périodiques des rotors à grand nombre de degrés de liberté apportés par les éléments finis volumiques. Puis, l'association à la méthode de continuation par pseudo-longueur d'arc aboutit à établir continûment l'ensemble des solutions d'équilibre dynamique sur la plage de vitesse de rotation. Enfin la stabilité dynamique locale de la solution périodique est analysée grâce à des indicateurs de bifurcation basés sur l'évolution des exposants de Floquet. Ainsi sont détectées les bifurcations de branches de solutions périodiques de type point limite, point de branchement et notamment Neimark-Sacker. Leur localisation est déterminée précisément en résolvant un système augmenté constitué de l'équation du mouvement et d'une équation supplémentaire caractérisant le type de bifurcation considéré. En déclarant un paramètre du système (coefficient de frottement, jeu rotor/stator, amplitude de l'excitation,...) comme nouvelle variable, l'utilisation de la technique de continuation conjointement avec le système augmenté détermine directement le cheminement des bifurcations en fonction de ce paramètre sur la nappe des réponses non linéaires. Les suivis de bifurcations délimitent les zones de fonctionnement spécifiques, extraient efficacement l'essentiel du comportement dynamique et offrent ainsi une nouvelle approche pour dimensionner de façon efficace les systèmes notamment en rotation. Nombre des développements réalisés sont implantés dans le code de calcul Cast3M. / Generally speaking, the rotating systems utilized in the energy production have a small rotor-stator gap, are able to run during long periods, and are mounted on hydrodynamic bearings. Rotor-stator interactions in case of blade loss, crack propagation due to fatigue, and a variable stiffness due to the nonlinear restoring forces of the bearings can make the rotordynamics nonlinear and the responses complicated: significant amplitude and frequency shifts are introduced, sub- and super-harmonics appear, and hysteresis occurs. It is of great importance to understand, predict and control this complicated dynamics. Due to the large number of DOFs and the broad range of study frequency, the computation time for solving the equations of motion by a temporal integration method can be quite prohibitive. It becomes particularly disadvantageous at the design stage where a parametrical study need to be quickly performed. An alternative numerical method, which is general and effective at the same time, is proposed in order to analyse the nonlinear response of the rotors at steady state. Firstly, the periodic responses of nonlinear rotors are calculated in the frequency domain by combining harmonic balance method (HBM) and alternating frequency-time (AFT). With the help of continuation method, all dynamic equilibrium solutions of nonlinear systems are determined for the range of study frequency. Then, Floquet exponents which are the eigenvalues of Jacobian are sought for stability analysis of periodic solutions. Then the local stability of the periodic solution is analysed through the bifurcation indicators which are based on the evolution of Floquet exponents. The bifurcations of periodic solution branch, such as limit point, branch point, and Neimark-Sacker bifurcation, are thus detected. By declaring a system parameter (friction coefficient, rotor / stator gap, excitation amplitude, ...) as a new variable, applying once again the continuation method to the augmented system determines directly the bifurcation's evolution as a function of this parameter. Thus, parametric analysis of the nonlinear dynamic behaviour is achieved, the stability boundary or the regime change boundary is directly determined. Numerous developments are implemented in the calculation code Cast3M.

L’objectif de cette thèse est d’étudier, à l’aide de modèles mathématiques, le mécanisme de régulation qui permet au corps de maintenir une quantité optimale de plaquettes sanguines. Le premier chapitre présente le contexte biologique et mathématique. Dans un second chapitre, un modèle pour la mégacaryopoïèse est introduit qui suppose une régulation ponctuelle par le nombre de plaquettes du taux de différentiation des cellules souches vers la lignée mégacaryocytaire et du nombre de plaquettes produites par mégacaryocyte. Nous montrons que la dynamique de ce modèle est régie par une équation différentielle à retard x'(t) = -?x(t)+f(x(t))g(x(t-t)), et nous obtenons ensuite de nouvelles conditions suffisantes pour la stabilité et l’oscillation des solutions de cette équation. Dans le troisième chapitre, nous analysons un second modèle pour la mégacaryopoïèse qui considère cette fois-ci une régulation opérée en continu uniquement via la vitesse de maturation des mégacaryoblastes. L’analyse de stabilité nécessite d’adapter un cadre pré-existant aux cas où le paramètre de bifurcation n’est pas le retard, et permet de montrer que l’augmentation du taux de mort des mégacaryoblastes conduit à l’apparition de solutions périodiques, en accord avec les observations cliniques de la thrombopénie cyclique amégacaryocytaire. Le dernier chapitre est consacré l’analyse de stabilité d’une équation différentielle à deux retards qui apparait notamment dans le cadre de la mégacaryopoïèse lorsque l’on considère que les plaquettes ont une durée de vie limitée / The object of this thesis is the study, using mathematical models, of the regulation mechanism maintaining an optimal quantity of blood platelets. The first chapter presents the biological and mathematical context of the thesis. In a second chapter, we introduce a model for megakaryopoiesis assuming a regulation by the platelet quantity of both the differentiation rate of stem cells to the platelet cell line and the amount of platelets produced by each megakaryocyte. We show that the dynamic of this model corresponds to a delay differential equation x'(t) = -?x(t) + f(x(t))g(x(t - t)), and we obtain for this equation new sufficient conditions for stability and for the oscillation of solutions. In a third chapter, we analyze a second model for megakaryopoiesis in which the regulation is continuous through the maturation speed of megakaryocyte progenitors. The stability analysis requires to adapt a pre-existing framework to problems where the bifurcation parameter is not the delay, and allows to show that increasing the death rate of megakaryocyte progenitors leads to the onset of periodic solutions, in agreement with clinical observation of amegakaryocytic cyclical thrombocytopenia. The last chapter covers a differential equation with two delays that appears among others in a model of platelet production which considers that platelet death can both age-independent and age-dependent

Cette étude traite de la stabilité et la post-bifurcation des élastomères magnétorhéologiques isotropes (MRE). Les MRE sont des élastomères comprenant une fraction en volume fini de particules de fer magnétisables, réparties de façon aléatoire dans le volume. Plus précisément, un système de film/substrat magnéto-élastique non linéaire est exploité expérimentalement, numériquement et théoriquement pour obtenir un contrôle actif de la rugosité de la surface du film. L'interaction non-intuitive entre le champ magnétique et la déformation élastique est due au choix des matériaux et de la géométrie du système, à savoir un film composite de particules ferromagnétiques collé sur une fondation passive et compliante. La coopération de deux mécanismes qui sont par ailleurs indépendants, la pré-compression mécanique et le champ magnétique, permet de rapprocher la structure d'un état faiblement stable et puis de la rendre instable par des champs magnétiques ou mécaniques. Nous démontrons pour la première fois que le champ magnétique critique est une fonction décroissante de la pré-compression et vice versa. Les résultats expérimentaux sont ensuite sondés avec succès par des simulations à champs complets par éléments finis en grandes déformations et champs magnétiques. Une analyse théorique de bifurcation magnéto-mécanique sur un système magnéto-élastique infini est également utilisée pour explorer l'effet des propriétés combinées sur la réponse critique.Dans la perspective d'élargir l'activation de surface à de nouveaux motifs magnéto-mécaniques, nous étudions plus en détail la post-stabilité d'un bloc bi-couche entièrement magnétorhéologique. L'idée sous-jacente est de créer différents contrastes entre les couches de propriétés magnétiques/mécaniques et de déclencher une gamme de motifs de surface plus riche que celle déjà obtenue en utilisant un film MRE sur une fondation passive. Les calculs post-bifurcation des films MRE collés sur des substrats MRE permettent de mettre en évidence les modes morphologiques résultant de la (in)compatibilité des modes de champs indépendants. Le couplage magnéto-élastique permet le contrôle réversible marche/arrêt de la configuration de surface sous des champs magnétiques et mécaniques critiques ajustables et donc, cette étude constitue un premier pas vers des dispositifs haptiques et morphiques actifs. / The present work deals with the stability and post-bifurcation response of isotropic magnetorheological elastomers (MREs). MREs are elastomers comprising a finite volume fraction of magnetizable iron particles, distributed randomly in the volume. Specifically, a nonlinear magnetoelastic film/substrate system is experimentally, numerically and theoretically exploited to obtain active control of surface roughness. The non-intuitive interplay between magnetic field and elastic deformation owes to material and geometry selection, namely, a ferromagnetic particle composite film bonded on a compliant passive foundation. Cooperation of two otherwise independent loading mechanisms--mechanical pre-compression and magnetic field--allows to bring the structure near a marginally stable state and then destabilize it with either magnetic or mechanical fields. We demonstrate for the first time that the critical magnetic field is a decreasing function of pre-compression and vice versa. The experimental results are probed successfully with full-field finite element simulations at large strains and magnetic fields. A theoretical magnetomechanical bifurcation analysis on an infinite magnetoelastic system is further employed to explore the effect of the interlayer combined properties on the critical response and is compared with the available numerical results.  With the perspective of applying the principle of surface actuation to new magnetomechanically triggered patterns, we further investigate the post-bifurcation of an entirely magnetorheological bilayer block. The underlying idea is to create different interlayer contrasts of magnetic and mechanical properties allowing us to trigger a larger range of surface patterns than that already obtained when using a MRE film on a passive (magnetically insensitive) foundation. Post-bifurcation calculations of MRE films bonded on MRE substrates allow to reveal novel patterns that lead to significant curvature localisation and crinkling.  In all cases studied, the magnetoelastic coupling allows for the reversible on/off control of surface patterning under adjustable critical magnetic and mechanical fields for a single specimen and thus, this study constitutes a first step towards realistic active haptic and morphing devices.

Cette thèse s’inscrit dans le domaine mathématique de l’analyse des équations aux dérivées partielles non linéaires. Précisément, nous nous sommes intéressés à une classe de problèmes elliptiques et paraboliques avec coefficients singuliers. Ce manque de régularité pose un certain nombre de difficultés qui ne permettent pas d’utiliser directement les méthodes classiques de l’analyse non-linéaire fondées entre autres sur des résultats de compacité. Dans les démonstrations des principaux résultats, nous montrons comment pallier ces difficultés. Ceci suppose d’adapter certaines techniques bien connues mais aussi d’introduire de nouvelles méthodes. Dans ce contexte, une étape importante est l’estimation fine du comportement des solutions qui permet d’adapter le principe de comparaison faible, d’utiliser la régularité elliptique et parabolique et d’appliquer dans un nouveau contexte la théorie globale de la bifurcation analytique. La thèse se présente sous forme de deux parties indépendantes. 1- Dans la première partie (chapitre I de la thèse), nous avons étudié un problème quasi-linéaire parabolique fortement singulier faisant intervenir l’opérateur p-Laplacien. On a démontré l’existence locale et la régularité de solutions faibles. Ce résultat repose sur des estimations a priori obtenues via l’utilisation d’inégalités de type log-Sobolev combinées à des inégalités de Gagliardo-Nirenberg. On démontre l’unicité de la solution pour un intervalle de valeurs du paramètre de la singularité en utilisant un principe de comparaison faible fondé sur la monotonie d’un opérateur non linéaire adéquat. 2- Dans la deuxième partie (correspondant aux Chapitres II, III et IV de la thèse), nous sommes intéressés à l’étude de problèmes de bifurcation globale. On a établi pour ces problèmes l’existence de continuas non bornés de solutions qui admettent localement une paramétrisation analytique. Pour établir ces résultats, nous faisons appel à différents outils d’analyse non linéaire. Un outil important est la théorie analytique de la bifurcation globale qui a été introduite par Dancer (voir Chapitre II de la thèse). Pour un problème semi linéaire elliptique avec croissance critique en dimension 2, on montre que les solutions le long de la branche convergent vers une solution singulière (solution non bornée) lorsque la norme des solutions converge vers l’infini. Par ailleurs nous montrons que la branche admet une infinité dénombrable de "points de retournement" correspondant à un changement de l’indice de Morse des solutions qui tend vers l’infini le long de la branche. / This thesis is concerned with the mathematical study of nonlinear partial differential equations. Precisely, we have investigated a class of nonlinear elliptic and parabolic problems with singular coefficients. This lack of regularity involves some difficulties which prevent the straight-orward application of classical methods of nonlinear analysis based on compactness results. In the proofs of the main results, we show how to overcome these difficulties. Precisely we adapt some well-known techniques together with the use of new methods. In this framework, an important step is to estimate accurately the solutions in order to apply the weak comparison principle, to use the regularity theory of parabolic and elliptic equations and to develop in a new context the analytic theory of global bifurcation. The thesis presents two independent parts. 1- In the first part (corresponding to Chapter I), we are interested by a nonlinear and singular parabolic equation involving the p-Laplacian operator. We established for this problem that for any non-negative initial datum chosen in a certain Lebeque space, there exists a local positive weak solution. For that we use some a priori bounds based on logarithmic Sobolev inequalities to get ultracontractivity of the associated semi-group. Additionaly, for a range of values of the singular coefficient, we prove the uniqueness of the solution and further regularity results. 2- In the second part (corresponding to Chapters II, III and IV of the thesis), we are concerned with the study of global bifurcation problems involving singular nonlinearities. We establish the existence of a piecewise analytic global path of solutions to these problems. For that we use crucially the analytic bifurcation theory introduced by Dancer (described in Chapter II of the thesis). In the frame of a class of semilinear elliptic problems involving a critical nonlinearity in two dimensions, we further prove that the piecewise analytic path of solutions admits asymptotically a singular solution (i.e. an unbounded solution), whose Morse index is infinite. As a consequence, this path admits a countable infinitely many “turning points” where the Morse index is increasing.

Dans cette thèse, on s'est intéressé à la dynamique des poches de tourbillon pour des équations issues de la mécanique des fluides posées dans le plan. La thèse est composée de trois partie indépendantes. Un des objectifs est d'établir l'existence des tourbillons uniformément concentrés et rigides, c’est-à-dire, qui ne se déforment pas lors de l'évolution. Nous analysons deux configurations liées à la nature topologique du support: poches simplement et doublement connexes. Nos solutions sont obtenues via des techniques de bifurcations et d'analyse complexe. Le deuxième objectif est d'obtenir des précisions sur la structure globale du diagramme de bifurcation et sa réponse vis-à-vis des petites perturbations dans le modèle. Plus précisément, dans le deuxième chapitre on prouve l'existence de V-states doublement connexes dans un voisinage de l'anneau pour le modèle des surfaces quasi-géostrophique. On montre que l'on peut construire des branches de solutions qui sont des anneaux perturbés pour certaines valeurs explicites de vitesses angulaires qui sont liées aux fonctions hypergéométriques de Gauss et aux fonctions de Bessel. Le troisième chapitre porte sur l'étude de la structure du diagramme de bifurcation dans le cas doublement connexes pour l'équation d'Euler. Numériquement, près d'un cas dégénéré, les deux branches issues des deux vitesses angulaires possibles semblaient se rejoindre pour former un lacet. Nous avons prouvé analytiquement ce résultat. Le quatrième chapitre porte sur le modèle shallow water quasi-géostrophique. Dans une première partie, on prouve l'existence de V-states simplement connexes dans un voisinage du tourbillon de Rankine pour un nombre dénombrable de vitesses angulaires liées aux fonctions de Bessel modifiées. La deuxième partie porte sur la réponse du diagramme de bifurcation lorsque l'on fait varier un paramètre du modèle. On montre en particulier qu'une singularité présente lors d'un cas limite est éclatée. Notre étude analytique a été complétée par des simulations numériques portant sur les V-states limites pour les symétries deux et trois. / In this dissertation, we are concerned with the vortex dynamics for some equations arising in fluid mechanics. We distinguish three independent parts. One of the objectives is to prove the existence of uniformly concentrated rigid vortices, they do not change their shapes during the motion. We examine two configurations related to the topological nature of the support: simply and doubly connected vortex patches. Our solutions are obtained using bifurcation arguments and complex analysis tools. The second objective is to obtain some precisions on the global structure of the bifurcation diagram and its response to small perturbations. More precisely, in the second chapter we prove the existence of doubly connected V-states in a neighborhood of the annulus for the surface quasi-geostrophic model. We check that we can construct some branches of solutions which are perturbated annulus at some angular velocities related to hypergeometric Gauss functions and Bessel functions. The goal of the third chapter is to study the structure of the bifurcation diagram in the doubly connected case for Euler equations. Numerically, close to a degenerate case, the two branches of solutions come from the two angular velocities seems to merge to form a loop. We prove analytically this result. In the last chapter, we focus on the shallow quasi-geostrophic model. In the first part, we prove the existence of the simply V-states in a neighborhood of the Rankine Vortices for a countable number of angular velocities related to modified Bessel functions. In the second part, we study the reaction of the diagram bifurcation for small perturbations of the parameter. In particular, we prove that some singularities are broken due to a resonance phenomenon. Our analytical study is completed by numerical simulations on the limiting V-states for the two and three fold symetries.

Cette thèse s'appuie sur deux branches des systèmes dynamiques : la théorie du nombre de rotation des endomorphismes du cercle de degré un et des applications de l'anneau déviant la verticale, ainsi que la théorie des systèmes non-uniformément hyperboliques. Nous nous intéressons tout d'abord à une classe d'applications bimodales du cercle, dilatantes et affines par morceaux. Chaque application de cette famille possède un nombre de rotation presque sûr : c'est le nombre de rotation de presque tout point du cercle. Nous étudions sa régularité et montrons que le nombre de rotation presque sûr est irrationnel pour un ensemble de paramètres de mesure totale. Nous considérons ensuite les applications de l'anneau qui dévient la verticale et plus particulièrement les applications bimodales de la famille d'Arnol'd épaissie. Un rôle essentiel est joué par les orbites de torsion nulle. Elles permettent de montrer que l'ensemble des applications qui possèdent un nombre de rotation fixé, forme dans l'espace des paramètres une langue d'Arnol'd bordée par deux surfaces. La frontière des langues rationnelles est associée à des bifurcations selle-noeud et homoclines. Nous obtenons enfin des estimations sur la taille de l'ensemble de rotation et de l'attracteur de Birkhoff. L'appendice est consacré aux bifurcations selles-noeud d'ensembles hyperboliques localement maximaux dont la direction instable est de dimension un. Cette bifurcation préserve la décomposition géométrique de l'espace tangent en espaces stables et instables. En revanche, l'expansion dans la direction instable dégénère près d'une orbite périodique. Nous obtenons alors une bifurcation de codimension un.

L’ischémie aigüe (restriction de la perfusion suite à l’infarctus du myocarde) induit des changements majeurs des propriétés électrophysiologique du tissu ventriculaire. Dans la zone ischémique, on observe une augmentation du potassium extracellulaire qui provoque l’élévation du potentiel membranaire et induit un "courant de lésion" circulant entre la zone affectée et saine. Le manque d’oxygène modifie le métabolisme des cellules et diminue la production d’ATP, ce qui entraîne l’ouverture de canaux potassique ATP-dépendant. La tachycardie, la fibrillation ventriculaire et la mort subite sont des conséquences possibles de l’ischémie. Cependant les mécanismes responsables de ces complications ne sont pas clairement établis. La création de foyer ectopique (automaticité), constitue une hypothèse intéressante expliquant la création de ses arythmies. Nous étudions l’effet de l’ischémie sur l’automaticité à l’aide d’un modèle mathématique de la cellule ventriculaire humaine (Ten Tusscher, 2006) et d’une analyse exhaustive des bifurcations en fonction de trois paramètres : la concentration de potassium extracellulaire, le "courant de lésion" et l’ouverture de canaux potassiques ATP-dépendant. Dans ce modèle, nous trouvons que seule la présence du courant de lésion peut entrainer une activité automatique. Les changements de potassium extracellulaire et du courant potassique ATP-dépendant altèrent toutefois la structure de bifurcation. / Acute ischemia (restriction in blood supply to part of the heart which can result in myocardial infarction) induces major changes in the electrophysiological properties of the ventricular tissue. Extracellular potassium concentration increases in the ischemic zone, leading to an elevation of the resting membrane potential that creates an “injury current” between the infarcted and the healthy zone. In addition, the lack of oxygen impairs the metabolic activity of the myocytes and decreases ATP production, thereby affecting ATP-sensitive potassium channels. A frequent complication of myocardial infarction is tachycardia, fibrillation and sudden cardiac death, but the mechanisms underlying their initiation are still debated. One hypothesis is that these arrhythmias may be triggered by abnormal automaticity. We have investigated the effect of ischemia on myocyte automaticity by performing a comprehensive bifurcation analysi(fixed points, cycles and their stability) of a human ventricular myocyte model (Ten Tusscher, 2006) as a function of three ischemia-relevant parameters: extracellular potassium, injury current, and ATP-sensitive potassium current. In this single-cell model, we found that automatic activity was possible only in the presence of an injury current. Changes in extracellular potassium and ATP-sensitive potassium current significantly altered the bifurcation structure as a function of IS, including the occurrence of early-after-depolarization. The results provide a sound basis for studying higher-dimensional tissue structures representing an ischemic heart.

Les régions de montagne sont parfois le siège de phénomènes dévastateurs brutaux, tels que les éboulements, les glissements de terrains, ou encore les avalanches. Fondamentalement, ces phénomènes sont issus de la rupture d’un volume de matériau qui va alors s’écouler sous l’effet de son poids propre. Prévoir ces phénomènes impose donc de mieux décrire la phase initiale de rupture. Ce travail de thèse s’inscrit dans cette volonté, en considérant le cas très général de la rupture dans un milieu granulaire. Pour cela, des simulations numériques utilisant une méthode aux éléments discrets ont été mises en œuvre, afin de vérifier la pertinence des critères usuels de rupture et d’instabilité. Dans un second temps, des analyses plus fines basées sur le travail du second ordre calculé à partir des variables microscopiques ont été menées afin de cerner comment évolue la microstructure de l’assemblage granulaire juste avant et pendant la rupture. Finalement, l’influence du chemin de chargement et des paramètres de contrôle sur la nature (effective ou non effective) et le mode (diffus ou localisé) de rupture a été examinée à travers une série d’essais biaxiaux sur deux échantillons numériques bidimensionnels (dense et moyennement dense). / Landslides and avalanches are the most common brutal and devastating natural hazards in mountain regions. Basically, these phenomena are derived from the failure of a material volume which will then flow under the effect of its own weight. Therefore, predicting these phenomena requires a better description of the initial phase of failure. This work comes within this willingness, considering the very general case of failure in granular media. For this purpose, numerical simulations using a discrete element method were carried out in order to investigate the relevance of the common failure and instability criteria. In a second step, more detailed analyses based on the second order work computed from the microscopic variables were conducted in order to describe the evolution of the microstructure of a granular assembly just before and during the failure. Finally, the influence of the loading path and the control parameters on the nature (effective or non-effective) and the mode (diffuse or localized) of failure was discussed through a series of biaxial tests performed on two-dimensional numerical samples depicting respectively a dense and a medium dense dry sand.

La mémoire présente une étude numérique des trajectoires non-verticales d’objets en chute ou en ascension libre dans un fluide newtonien initialement au repos. Une méthode numérique originale combinant une discrétisation spatiale spectrale et la décomposition du domaine a été implémentée à cet effet. Le code obtenu a été exploité pour apporter de nouvelles connaissances sur des objets fixes et objets libres. Pour les objet fixes, ellipsoïdes et cylindres, l’étude très complète des divers états de la transition a permis d’établir un lien entre le scénario de transition de la sphère et du disque infiniment mince. La simulation numérique d’objets libres a apporté des résultats très complets sur la chute de disques minces et de cylindres de faible épaisseur. Plusieurs questions soulevées dans les travaux précédents ont trouvés des réponses claires. Une étude paramétrique exhaustive, jamais abordée précédemment, portant sur des ellipsoïdes est décrite dans le dernier chapitre du mémoire. / The manuscript presents a numerical study of non-vertical trajectories of bodies fallingor rising freely in a Newtonian fluid initially at rest. The original numerical methodcombining a spectral spatial discretization and the decomposition of the domain wasimplemented. The developed code was used to study fixed and freely moving objects. Forthe fixed bodies, spheroids and cylinders, a very complete study of various transitionalregimes was carried out and established a link between the transition scenario of a sphereand of a disc. The numerical simulation of free bodies brought about the results onpath instabilities of discs and cylinders of small thickness. Several questions raised inprevious studies were clearly answered. An exhaustive parametric study, that has never beentaken up before, of freely moving spheroids is described in the last chapter of the manuscript.

Ces travaux ont été menés dans le cadre du projet européen FACETS-ITN. Nous avons contribué à la simulation de cellules corticales grâce à des données expérimentales d'électrophysiologie comme référence et d'un circuit intégré neuromorphique comme simulateur. Les propriétés intrinsèques temps réel de nos circuits neuromorphiques à base de modèles à conductance, autorisent une exploration détaillée des différents types de neurones. L'aspect analogique des circuits intégrés permet le développement d'un simulateur matériel temps réel à l'échelle du réseau. Le deuxième objectif de cette thèse est donc de contribuer au développement d'une plate-forme mixte - matérielle et logicielle - dédiée à la simulation de réseaux de neurones impulsionnels. / This work has been supported by the European FACETS-ITN project. Within the frameworkof this project, we contribute to the simulation of cortical cell types (employingexperimental electrophysiological data of these cells as references), using a specific VLSIneural circuit to simulate, at the single cell level, the models studied as references in theFACETS project. The real-time intrinsic properties of the neuromorphic circuits, whichprecisely compute neuron conductance-based models, will allow a systematic and detailedexploration of the models, while the physical and analog aspect of the simulations, as opposedthe software simulation aspect, will provide inputs for the development of the neuralhardware at the network level. The second goal of this thesis is to contribute to the designof a mixed hardware-software platform (PAX), specifically designed to simulate spikingneural networks. The tasks performed during this thesis project included: 1) the methodsused to obtain the appropriate parameter sets of the cortical neuron models that can beimplemented in our analog neuromimetic chip (the parameter extraction steps was validatedusing a bifurcation analysis that shows that the simplified HH model implementedin our silicon neuron shares the dynamics of the HH model); 2) the fully customizablefitting method, in voltage-clamp mode, to tune our neuromimetic integrated circuits usinga metaheuristic algorithm; 3) the contribution to the development of the PAX systemin terms of software tools and a VHDL driver interface for neuron configuration in theplatform. Finally, it also addresses the issue of synaptic tuning for future SNN simulation.