Les vésicules sont des gouttes immergées dans un fluide externe visqueux, dont le rayon vautquelques dizaines de micromètres et entourées par une membrane imperméable constituée de lipides, dont l’épaisseur est approximativement 4 nm. La membrane d’une vésicule est un systèmeoriginal du point de vue mécanique : celle-ci présente à la fois des propriétés fluides (les lipidespeuvent s’écouler librement le long de la membrane, mais la surface est incompressible locale-ment) et des propriétés solides (la membrane résiste à la flexion). Les propriétés spécifiques de lamembrane rendent ce système à la fois très déformable et très contraint.Ce manuscrit s’intéresse à la modélisation d’une vésicule soumise à des efforts extérieurs d’o-rigine hydrodynamique, dans le régime de Stokes. Une attention particulière est consacrée à lasituation d’une vésicule qui sédimente. Cette situation est étudiée analytiquement dans le régimedes faibles déformations. Il est montré que plusieurs familles de solutions stationnaires non triv-iales existent, grâce aux propriétés spécifiques de la membrane. L’étude de la sédimentation d’unevésicule est poursuivie par le développement d’un code numérique capable de simuler de grandesdéformations. Pour cela, des méthodes numériques originales de calcul de prise en compte de laflexion et de l’incompressibilité surfacique sont développées. Ce code permet d’étudier la forma-tion d’un tube à l’arrière d’une vésicule en sédimentation. Ces tubes sont de fins (rapport d’aspecttypique longueur/rayon ∼ 100) cylindres connectés à la vésicule d’origine. Il est montré que cesformes tubes sont des formes stationnaires. Un modèle théorique est proposé et comparé auxsimulations numériques. Ce modèle met en lumière l’importance particulière de la tension dansces formes. Une modélisation mécanique basée sur un milieu de Cosserat surfacique courbé estégalement présentée, et permet d’identifier la contribution de la flexion au tenseur des contraintes.Cette contribution est un ingrédient indispensable pour comprendre les formes tubes. / Vesicles are drops of radius of a few tens micrometers, bounded by an impermeable lipidmembrane of approximately 4 nm thickness, and embedded in an external viscous fluid. Thevesicle membrane is an original system from the mechanical point of view : it presents bothincompressible fluid properties (the lipids can flow freely along the membrane, but membraneis incompressible locally) and solid properties (the membrane resists to bending). The specificproperties of the membrane make the system very deformable and very constrained at the sametime.This manuscript deals with the modelisation of a vesicle subjected to external stresses of hydrodynamical origin, in the Stokes regime. A particular attention is paid to the situation of asettling vesicle. This situation is studied analytically in the small deformation regime. It is foundthat several families of non-trivial stationnary shapes exist, owing to the specific properties ofthe membrane. The study of a settling vesicle is pursued by the development of a numerical codeable to deal with large deformations. Original numerical methods are developped to deal with thecomputation of the bending and with the surface incompressibility constraint. This code permitsto study the formation of tether at the rear of a settling vesicle. These tethers are thin (typicalaspect ratio : length/radius ∼ 100) cylinders of membrane connected to the original vesicle. Itis shown that these tethered shapes are stationary shapes. A theoretical model is proposed andcompared to numerical simulations. This model shows the particular importance of tension inthese shapes. A mechanical modelling based on a curved Cosserat surface is also presented, andpermits to identify the bending contribution to the stress tensor. This contribution is a salientingredient to understand tethered shapes.

A l’heure actuelle, la microfluidique est une science en plein développement ayant un besoin croissant de dispositifs permettant de générer des écoulements aux échelles micrométriques. Les phénomènes physiques mis en jeu lors du mouvement d’un fluide sont en effet majoritairement gouvernés par la viscosité (bas nombre de Reynolds) contrairement aux écoulements macroscopiques dominés par les effets inertiels.Dans cette thèse, les écoulements engendrés par le mouvement de parois mobiles ont été étudiés en vue d’une application aux micropompes, dispositifs essentiels en microfluidique.Dans une première partie, une étude analytique et numérique évalue la possibilité de générer un écoulement par un cylindre en rotation à proximité de parois mobiles.Les résultats obtenus du régime de Stokes (Re=0) jusqu’à un nombre de Reynolds Re=60 en régime stationnaire témoignent du potentiel notable d’intégration de cette géométrie dans les microsystèmes en tant que micropompes.Dans une seconde partie, une micropompe, basée sur un principe de fonctionnement novateur, est conçue par l’intermédiaire des techniques de microfabrication. Dans cette optique, le procédé de gravure RIE d’un élastomère est entièrement développé. Les performances de la micropompe en terme de pression et débit générés dépassent l’état de l’art des microsystèmes similaires et ceci en utilisant une technologie simple et bas-coût / Currently, microfluidic is a science field in constant development with an increasing need of devices able to generate flows at the micrometer order. At these length scales, physical phenomenons occurring in a moving fluid are mainly governed by its viscosity (low Reynolds number) contrary to macroscale flows dominated by inertial effects.In this thesis, a study on flows engendered by moving walls has been carried to fulfill to micropumps devices.In a first part, an analytical and a numerical study evaluates the possibility to generate a flow for a rotating cylinder close to moving boundaries.The results ranging from Stokes flows (Re=0) up to the low Reynolds number Re=60 in the stationary regime reveals the noticeable potential of integrating this device in microsystems as a micropump. In a second part, a new micropump, based on an innovative principle, is designed thanks to microfabrication technologies. In this perspective, the etching process of an elastomer called Silastic S is developed. Micropump performances in terms of pressure and flow rate are beyond the state of the art for similar microsystems and are achieved by using a simple and low-cost technology

Fluid mechanics of blood can be well approximated by a mixture model of a Newtonian fluid and deformable particles representing the red blood cells. Experimental and theoretical evidence suggests that the deformation and rheology of red blood cells is similar to that of phospholipid vesicles. Vesicles and red blood cells are both area preserving closed membranes that resist bending. Beyond red blood cells, vesicles can be used to investigate the behavior of cell membranes, intracellular organelles, and viral particles. Given the importance of vesicle flows, in this thesis we focus in efficient numerical methods for such problems: we present computationally scalable algorithms for the simulation of dilute suspension of deformable vesicles in two and three dimensions. Our method is based on the boundary integral formulation of Stokes flow. We present new schemes for simulating the three-dimensional hydrodynamic interactions of large number of vesicles with viscosity contrast. The algorithms incorporate a stable time-stepping scheme, high-order spatiotemporal discretizations, spectral preconditioners, and a reparametrization scheme capable of resolving extreme mesh distortions in dynamic simulations. The associated linear systems are solved in optimal time using spectral preconditioners. The highlights of our numerical scheme are that (i) the physics of vesicles is faithfully represented by using nonlinear solid mechanics to capture the deformations of each cell, (ii) the long-range, N-body, hydrodynamic interactions between vesicles are accurately resolved using the fast multipole method (FMM), and (iii) our time stepping scheme is unconditionally stable for the flow of single and multiple vesicles with viscosity contrast and its computational cost-per-simulation-unit-time is comparable to or less than that of an explicit scheme. We report scaling of our algorithms to simulations with millions of vesicles on thousands of computational cores.

Cette thèse propose d'étudier en détail une méthode générale de construction d'un réseaude pores interconnectés à partir d'une image 3D d'un matériau poreux afin de calculer laperméabilité absolue. Un squelette de l'espace poral est tout d'abord utilisé pour définir laposition des pores. Ce squelette est transformé en graphe, puis plusieurs étapes de fusiondes nœuds sont réalisées pour obtenir un réseau exploitable. Une ligne de partage des eauxprenant comme marqueurs les nœuds du graphe permet d'attribuer une géométrie auxpores. Le système d'équation linéaire à résoudre est construit en intégrant les équations deStokes sur des éléments de l'espace poral. Une méthode originale de calcul du coefficient deproportionnalité existant entre les pressions au centre de deux pores et le débit traversantla surface les séparant est proposée. Une application de l'intégralité de l'approche estégalement présentée sur un matériau réel. / A general method to build a pore network model from a 3D image of a porous material ispresented in this work in order to compute its permeability. A skeleton of the void spaceis first used to define pores position. This skeleton is converted to a first graph. Mergingsteps are necessary to obtain a relevant network. A volume is allocated to the pores using awatershed algorithm, starting with markers defined from the nodes of the graph. The linearsystem of equations to solve is deduced by integrating the Stokes equations on pore spacepartitioning elements. The proportionality coefficient linking the pressure drop betweentwo pore centers and the flow going through the surface separating them is calculated withan original method. Application of the complete process to a real material is presented.

Le mélange de solutés par les écoulements en milieux poreux contrôle les réactions chimiques dans un grand nombre d'applications souterraines, dont le transport et la remédiation des contaminants, le stockage et l'extraction souterrains d'énergie, et la séquestration du CO2. Nous étudions les mécanismes du mélange à l'échelle du pore et plus précisément comment la topologie de l'écoulement est reliée à la dynamique du mélange d'espèces conservatives; en particulier, l'émergence d'un mélange chaotique est-elle possible dans un milieu poreux tridimensionnel (3D) ? Nous calculons donc numériquement ou mesurons expérimentalement les vitesses d'écoulement et l'évolution temporelle des champs de concentration afin de caractériser la déformation et le mélange à l'échelle du pore. Une première étude, expérimentale, permet de caractériser le mélange dans un fluide s’écoulant à travers un milieu poreux bidimensionnel (2D). Nous mesurons les vitesses par suivi de microparticules solides (''PTV''). L’évolution temporelle de la distance séparant deux particules permet de caractériser la dynamique de la déformation lagrangienne. Des mesures de transport conservatif dans le même milieu fournissent l'évolution temporelle du gradient de concentration moyen (une mesure du mélange). À partir de ces résultats expérimentaux nous proposons la première validation expérimentale à l'échelle du pore de la théorie lamellaire du mélange, reliant les propriétés de la déformation du fluide à la dynamique du mélange. Dans une deuxième étude nous examinons les conditions d'apparition du mélange chaotique dans l’écoulement dans des milieux poreux 3D granulaires ordonnés. Nous effectuons des calculs numériques hautement résolus de d'écoulement de Stokes entre des sphères empilées selon une structure cristalline (cubique simple ou cubique centrée), périodique. La déformation lagrangienne, obtenue à partir des champs de vitesse à l'aide d'outils numériques développés spécifiquement, met en lumière une large variété de dynamiques de la déformation dans ces milieux 3D, selon l'orientation de l'écoulement. Quand la direction de l'écoulement n'est pas normale à l'un des plans de symétrie de réflection du cristal, l'évolution temporelle de la déformation est exponentielle, traduisant une advection chaotique. L’émergence (ou non) du chaos est contrôlée par un mécanisme similaire à la ''transformation du boulanger'': les particules fluides se déplaçant autour d'un grain solide se retrouvent séparées par une surface virtuelle (appelée “variété”) qui émerge de la surface du grain. De multiples variétés existent dans l’écoulement, et la façon dont elles s'intersectent contrôle la nature - chaotique ou non - du mélange, et l'intensité du chaos. En particulier, l'exposant de Lyapunov (une mesure du chaos), est contrôlé par la fréquence spatiale des intersections appropriées à la génération du chaos, nommées ''connections hétéroclinines'' entre variétés. L'image conventionnelle, 2D, des mécanismes du mélange, impose des contraintes topologiques qui ne permettent pas le développement de ces mécanismes 3D. Elle pourrait donc être inadaptée aux milieux poreux naturels. La troisième étude a deux objectifs: (i) fournir une preuve expérimentale de la nature chaotique de l'advection, par la visualisation des variétés et par l'obtention d'une mesure de l'exposant de Lyapunov; et (ii), évaluer si nos résultats numériques obtenus pour des milieux granulaires ordonnés peuvent être généralisés à des milieux désordonnés, plus proches des milieux naturels. L’expérience est fondée sur un empilement désordonné de sphères rendu transparent par l'ajustement optique du liquide avec les sphères. La fluorescence induite par laser (''LIF'') permet de détecter les variétés au sein de l'écoulement, et des techniques PTV de mesurer les vitesses d'écoulement et quantifier l'exposant de Lyapunov. Les premiers résultats expérimentaux sont prometteurs. / Solute mixing in porous media flows plays a central role in driving chemical reactions in a number of subsurface applications, including contaminant transport and remediation, subsurface energy storage and extraction, and CO2 sequestration. We study the mechanisms of solute mixing, in particular how the pore scale flow topology is related to the mixing dynamics of conservative solutes, with a particular emphasis on the possible emergence of chaotic mixing processes in three-dimensional (3D) porous media. To do so, we perform numerical computations or experimental measurements of the flow velocities and temporal evolution of the concentration fields, and characterize fluid deformation and mixing at the pore scale. This PhD work consists of three main studies. In the first study, we experimentally characterize mixing in a fluid flowing through a two-dimensional (2D) porous medium built by lithography. We measure the velocity distributions from Particle Tracking Velocimetry (PTV). The time evolution of the separation distance between two particles is analyzed to characterize the Lagrangian deformation dynamics. In parallel we perform conservative transport experiments with the same porous media, and quantify the temporal evolution of the mean concentration gradient, which is a measure of the mixing rate. From these experimental results we obtain the first experimental pore scale validation of the lamella mixing theory, which relates the fluid deformation properties to the mixing dynamics. In the second study, we investigate the conditions of emergence of chaotic mixing in the flow through 3D ordered granular porous media. In these periodic cubic crystalline packings (Simple Cubic - SC - and Body-Centered Cubic - BCC) of spheres, we are able to perform highly resolved computations of the 3D Stokes flow. Using custom-developed numerical tools to measure the Lagrangian deformation from the computed velocity fields, we uncover the existence of a rich array of Lagrangian deformation dynamics in these 3D media, depending on the flow orientation. When the flow direction is not normal to one of the reflection symmetry planes of the crystalline lattice, we find that the Lagrangian deformation dynamics follow an exponential law, which indicates chaotic advection. This chaotic behavior is controlled by a mechanism akin to the baker's transformation: fluid particles traveling around a solid grain along different paths end up either separated by, or on the same side of, a virtual surface projecting from the grain surface and called a manifold. Multiple such manifolds exist within the flow, and the way they intersect controls the nature of mixing (that is, either non-chaotic or chaotic), and the strength of chaos. We show in particular that the magnitude of the Lyapunov exponent (a measure of the vigor of chaos) is controlled by the spatial frequency of transverse connections between the manifolds (called heteroclinic intersections). We thus demonstrate that the conventional 2D picture of the mechanisms of mixing may not be adapted for natural porous media because that picture imposes topological constraints which cannot account for these important 3D mechanisms. The third study has two objectives: (i) provide experimental evidence of the chaotic nature of pore scale advection/mixing, both by visualizing the manifolds and by obtaining a quantitative estimate of the Lyapunov exponent; and (ii) assess if the results obtained numerically in ordered packings of spheres extend to random packings, which are closer to natural porous media. The experiment features a random packing of glass beads rendered transparent by optical index-matching between the fluid and solid grains. We use Laser Induced Fluorescence (LIF) to detect the manifolds, and PTV techniques to measure flow velocities and subsequently quantify Lyapunov exponent. The first experimental results are promising.

We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid with inertial forces. A new feature of this work is the implementation of the local inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. We introduce a new equation for the local Lagrange multiplier whose solution essentially provides a harmonic extension of the multiplier off the interface while maintaining the local inextensibility constraint near the interface. We also develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. To solve the equations, we use an adaptive finite element method with implicit coupling between the Navier-Stokes and the diffuse interface inextensibility equations. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and smaller differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to prevent the accumulation of stretching and compression errors very effectively. Simulations of two vesicles in an extensional flow show that local inextensibility plays an important role when vesicles are in close proximity by inhibiting fluid drainage in the near contact region.

De nombreux systèmes fluides dans les domaines de la biologie ou encore de la cosmétique sont limités par une interface dont les propriétés mécaniques régissent la stabilité. En particulier, les objets tels que des gouttes, vésicules ou polymersomes se déforment dans un écoulement simple et mènent à une grande richesse de dynamiques spatio-temporelles contrôlées par la nature des matériaux qui composent l'interface. Les travaux présentés concernent l'étude numérique de la déformation de ces objets dans un écoulement de Stokes, en particulier dans des situations où les viscosités de l'interface jouent un rôle important. Un code de calcul couplant intégrales de frontières et éléments finis a été utilisé afin de décrire la physique interfaciale et étudier leur comportement une fois plongés dans un écoulement. Ces travaux ont permis d'étudier l'influence des viscosités interfaciales sur la dynamique d'une goutte dans un écoulement extensionnel plan, leur influence sur sa dynamique de déformation et sur les conditions de rupture de celle-ci. Les études réalisées sur une vésicule fortement dégonflée et plongée dans un écoulement cisaillé ont caractérisé la bifurcation entre les deux familles de forme existantes dans ces conditions. Ces formes ayant une influence sur la dynamique de la vésicule dans l'écoulement, celle-ci a été étudiée dans le cadre d'un écoulement infini puis proche d'une paroi parallèle à l'écoulement. Enfin, de premiers résultats sur la dynamique d'un polymersome dans un écoulement cisaillé permettent de construire un diagramme de phase illustrant les différents comportement de cet objet en fonction de la viscosité de la membrane et du taux de cisaillement / There are many fluid systems in the biology, food industry, pharmacology or cosmestics fields that are bound by an interface which mechanical properties rule the system stability. Objects like droplets, vesicles or polymersomes change their shape in a simple flow which lead to a wealth of space and time dynamics. These properties are controlled by the nature of the interface material. The aim of this work is the numerical study of the deformation of droplets, vesicles and polymersomes in a Stokes flow, especially when the interfacial viscosities play an important role. A numerical computation code coupling boundary integrals and finite elements was used to describe the interfacial physics of these objects and study their behaviour when immerged in a flow. Multiple resolution strategies where developped to this end in order to optimize the numerical computation in the cas of an interface with viscosities.Using this work, the influence of interfacial viscosities on the dynamics of a droplet in an extensional flow is studied : in particular, their influence on the stretching dynamics of a droplet and its break up conditions was characterized. The study of a vesicle, droplet bounded by a lipid bilayer, strongly deflated and immerged in a shear flow detailed the bifurcation between two shape types existing for this system. These shapes have an influence on the vesicle dynamics under flow, which is studied for an unbounded flow and a near-wall flow. Finally, we show first results about the dynamics of a polymersome in a shear flow. We used them to build a phase diagram for the behaviour of this object depending on the membrane viscosity and the shear rate

Este trabalho propõe o estudo do projeto de dispositivos baseados em microcanais de fluido, tais como difusores, misturadores, válvulas, e trocadores de calor, através da aplicação do Método de Otimização Topológica (MOT). O MOT é um método computacional que permite obter um projeto otimizado de um sistema, através da distribuição de uma quantidade limitada de material num dado domínio de projeto. Neste caso, o MOT é aplicado a um domínio fluido, e permite obter a topologia otimizada (formato ótimo) dos microcanais, segundo uma determinada característica, seja esta, a minimização da perda de carga, ou a maximização da velocidade num dado ponto, ou ainda a maximização da troca de calor, no caso de trocadores de calor. Os canais utilizados nestas aplicações operam com baixo número de Reynolds, sendo um caso típico da aplicação das equações de escoamento de Stokes. A implementação do MOT é realizada sob a forma de rotinas computacionais, permitindo um projeto sistematizado dos canais. No processo de otimização, utiliza-se o Método dos Elementos Finitos (MEF) como método de análise dos fenômenos físicos envolvidos, e a Programação Linear Seqüencial (PLS) como algoritmo de otimização. Ao final, propõe-se um estudo multi-físico, aliando-se características otimizadas tanto do ponto de vista da eficiência do escoamento, quanto do ponto de vista da dissipação térmica no canal, combinando-os através de uma função multi-objetivo. Exemplos de projeto bidimensionais de dispositivos de fluido são apresentados para ilustrar o método. / This work proposes studying the design of micro channel devices, such as fluid diffusers, mixers, valves, and heat exchangers, through the application of the Topology Optimization Method (TOM). The TOM is a computational method that allows the distribution of a limited amount of material, inside a given design domain, in order to obtain an optimized system design. Herein, the TOM is applied to a fluidic domain, allowing the design of an optimized microchannel topology (optimal configuration), according to a given objective function, such as head loss minimization, maximum velocity in a given direction, or the heat transfer maximization, in a heat exchanger example. Especially this kind of channel devices, operates at low Reynolds number, thus, it can be modeled through Stokes flow equations. The optimization procedure applies the Finite Element Method (FEM) to perform the physical analysis, and Sequential Linear Programming (SLP) as the optimization algorithm. At the end, a multi-physics analysis is proposed, through a multi-objective cost function, that combines both flow and heat dissipation efficiency optimization. Two-dimensional designs of fluidic devices are presented as examples to illustrate the method.

In this thesis, simulations are performed to study the motion ofindividual cells in flow, focusing on the hydrodynamics of actively swimming cells likethe self-propelling microorganisms, and of passively advected objects like the red bloodcells. In particular, we develop numerical tools to address the locomotion ofmicroswimmers in viscoelastic fluids and complex geometries, as well as the motion ofdeformable capsules in micro-fluidic flows. For the active movement, the squirmer is used as our model microswimmer. The finiteelement method is employed to study the influence of the viscoelasticity of fluid on theperformance of locomotion. A boundary element method is implemented to study swimmingcells inside a tube. For the passive counterpart, the deformable capsule is chosen as the modelcell. An accelerated boundary integral method code is developed to solve thefluid-structure interaction, and a global spectral method is incorporated to handle theevolving cell surface and its corresponding membrane dynamics. We study the locomotion of a neutral squirmer with anemphasis on the change of swimming kinematics, energetics, and flowdisturbance from Newtonian to viscoelastic fluid. We also examine the dynamics of differentswimming gaits resulting in different patterns of polymer deformation, as well as theirinfluence on the swimming performance. We correlate the change of swimming speed withthe extensional viscosity and that of power consumption with the phase delay of viscoelasticfluids. Moreover, we utilise the boundary element method to simulate the swimming cells in astraight and torus-like bent tube, where the tube radius is a few times the cell radius. Weinvestigate the effect of tube confinement to the swimming speed and power consumption. Weanalyse the motions of squirmers with different gaits, which significantly affect thestability of the motion. Helical trajectories are produced for a neutralsquirmer swimming, in qualitative agreement with experimental observations, which can beexplained by hydrodynamic interactions alone. We perform simulations of a deformable capsule in micro-fluidic flows. We look atthe trajectory and deformation of a capsule through a channel/duct with a corner. Thevelocity of capsule displays an overshoot as passing around the corner, indicating apparentviscoelasticity induced by the interaction between the deformable membrane and viscousflow. A curved corner is found to deform the capsule less than the straight one. In addition, we propose a new cell sorting device based on the deformability of cells. Weintroduce carefully-designed geometric features into the flow to excite thehydrodynamic interactions between the cell and device. This interaction varies andclosely depends on the cell deformability, the resultant difference scatters the cellsonto different trajectories. Our high-fidelity computations show that the new strategy achievesa clear and robust separation of cells. We finally investigate the motion of capsule in awall-bounded oscillating shear flow, to understand the effect of physiological pulsation to thedeformation and lateral migration of cells. We observe the lateral migration velocity of a cellvaries non-monotonically with its deformability. / <p>QC 20140313</p>

The main concern of this Dissertation is focused on the derivation of novel integral formulation for simple problems. This alternative integral representations display a rapid decay as the complex parameter involved tends to infinity and are therefore suitable for numerical computations and for the study of the asymptotic properties of those solutions. There is also another important advantage attached to the novel formulae presented. These integral representations are useful for solving changing-type boundary value problems (such as Dirichlet data on part of the boundary and Neumann data on the complementary of the boundary). The following problems are analyzed: (a) The Laplacian operator in the interior of a Square, (b) the Laplacian operator in the interior and exterior of a Sphere and, (c) the Stokes' operator concerning the irrotational flow of an incompressible, viscous fluid. Moreover, the behaviour of the Gegenbauer functions of the first and second kind of general complex degree and order on the cut (-1, +1) are examined. / Με οδηγό μια νέα μεθοδολογία επίλυσης Μερικών Διαφορικών Εξισώσεων (ΜΔΕ), προβλήματα που σχετίζονται με την θεωρία Δυναμικού όπως επίσης και με την ροή Stokes, θα αναλυθούν. Απώτερος σκοπός αποτελεί η ανάπτυξη ολοκληρωτικών αναπαραστάσεων, η οποίες χαρακτηρίζονται από ταχεία σύγκλιση, με σκοπό να χρησιμοποιηθούν στην ασυμπτωτική μελέτη, στην αριθμητική ανάλυση όπως επίσης και στην επίλυση προβλημάτων μεικτών συνοριακών συνθηκών (π.χ. δεδομένα Dirichlet στο ένα κομμάτι του συνόρου και δεδομένα Neumann στο υπόλοιπο). Συγκεκριμένα, τα ακόλουθα προβλήματα αναλύονται: (α) Εξίσωση Laplace στο εσωτερικό ενός τετραγώνου, (β) εξίσωση Laplace στο εσωτερικό και εξωτερικό μιας σφαίρας και, (γ) εξίσωση αστρόβιλης ροής Stokes στο εσωτερικό ενός σφαιρικού κελύφους το οποίο στην συνέχεια καταλήγει, με οριακές διαδικασίες, στο εσωτερικό και εξωτερικό μιας σφαίρας. Τέλος, παρουσιάζονται αναπτύγματα και ασυμπτωτικές εκφράσεις των συναρτήσεων Gegenbauer.