Ocean submesoscales are the subject of increasing focus in the oceanographic literature; with instrumentation now more capable of observing them in situ and numerical models now able to reach the resolution required to more fully capture them. Submesoscales are typified by horizontal spatial scales of O(1 − 10) km, vertical scales O(100) m and time-scales of O(1) day and are known to be associated with regions of high vertical velocity and vorticity. Occurring most commonly at density fronts at the ocean surface they can control mixed layer restratification and provide an important control on fluxes between the atmosphere and the deep ocean. This thesis sets out to better understand the fundamental physical processes underpinning submesoscale instabilities using a number of idealised process models. Linear stability analysis complemented by non-linear, high-resolution simulations will be used initially to explore the ways in which submesoscale instabilities in the mixed layer may compete and interact with one another. In particular, we will investigate the way in which symmetric and ageostrophic baroclinic instabilities interact when simultaneously present in a flow, with focus on the growth rates and energetic pathways of previously unexplored dynamic instabilities that arise in this paradigm; three-dimensional, mixed symmetric-baroclinic instabilities. Further, these non-linear simulations will allow us to investigate the transition to dissipative scales that can occur in the classical Eady model via a multitude of small-scale secondary instabilities that result from primary submesoscale instabilities. Finally, observational data, taken aboard the SMILES project cruise to the Southern Ocean, helps to motivate the consideration of a new dynamical paradigm; the Eady model with superimposed high amplitude barotropic jet. Non-linear simulations investigate the extent to which the addition of such a jet is capable of damping submesoscale growth. The causes of this damping are then investigated using linear analysis. With this approach eventually demonstrated as being unable to fully explain growth rate reductions, we introduce a new framework combining potential vorticity mixing by submesoscale instabilities with geostrophic adjustment, which relaxes the flow back to a geostrophic balanced state. This framework will help to explain, conceptually, how non-linear eddies control the linear stability of the flow.

Cette thèse a pour objectif d'améliorer la compréhension de la dynamique d'une interaction entre une onde de choc droite ou oblique et une couche limite laminaire ou turbulente. En particulier, nous nous sommes intéressés aux mécanismes responsables de l'apparition d'oscillations auto-entretenues basses fréquences. Ce phénomène survient dans de nombreux cas applicatifs comme dans des entrées d'air d'avions supersoniques, autour d'un profil d'aile en régime transsonique et au sein de tuyère en régime de sur-détente. La première partie de ce mémoire traite des différentes études réalisées pour déterminer la phénoménologie de ce type de dynamique. Dans un deuxième temps, nous expliquons la stratégie retenue pour effectuer notre étude qui consiste à développer un outil d'étude des instabilités, adapté à des écoulements turbulents présentant une interaction entre une onde de choc et une couche limite. Le développement d'un outil CFD linéarisé couplé à une méthode de résolution d'un problème aux valeurs propres par une approche dite sans matrice ou de « time-stepping », a permis la réalisation d'une telle étude. Après une étape de validation de notre outil, nous avons étudié des cas d'écoulements présentant une interaction entre une onde de choc et une couche limite. Trois cas en particulier ont été traités. Le premier cas correspond à une interaction entre une onde de choc oblique impactant une couche limite laminaire se développant sur une plaque plane. Ce cas est généralement qualifié dans la littérature de cas de « réflexion de choc ». Nous montrons qu'un tel écoulement est globalement stable et que sa dynamique peut être caractérisée par des mécanismes de réceptivité et par la réponse de l'écoulement vis-à-vis de perturbations extérieures. Les deux autres cas abordés dans ce travail ont été le cas d'un écoulement transsonique autour d'un profil d'aile de type NACA0012 en régime d'entrée en tremblement aérodynamique et un cas de tuyère transsonique plane de type Sajben en régime de sur-détente. Dans le premier cas, l'analyse de stabilité nous permet de mettre en évidence le phénomène de « buffet » sur le profil NACA0012, ce qui montre que le phénomène est lié à une instabilité globale linéaire . Dans le second cas, l'analyse de stabilité ne permet pas d'expliquer le phénomène auto-entretenues basses fréquences, et montre que l'écoulement est linéairement globalement stable. Dans ce cas, la dynamique est convective, transitoire et pilotée par des mécanismes de réceptivité. / The general purpose of this study is to provide a better understanding of the dynamics of an interaction between a shock wave and a laminar or turbulent boundary layer. In particular, we were interested in mechanisms responsible for the emergence of low-frequency self-sustained oscillations. This phenomenon arises in numerous industrial cases as in air inlets of supersonic aircrafts, around a profile of wing in transonic regime and within over-extended nozzle. The first part of this report handles various studies carried out to determine the phenomenology of this kind of dynamics. Secondly, we explain the strategy adopted to make our study which consists in developing a tool of study of the instabilities adapted to turbulent flows including an interaction between a shock wave and a boundary layer. The development of a linearized CFD tool coupled with a method of resolution of a eigenvalue problem by a free-matrix approach ( " time-stepping " approach), allowed the realization of such a study. After a stage of validation of our tool, we studied cases of flows including an interaction between a shock wave and a boundary layer. Three cases in particular were handled. The first case corresponds to an interaction enter an oblique shock wave impacting on a laminar boundary layer developing on a flat plate. This case is generally qualified in the literature of case as " reflected shock wave". We show that such a flow is globally stable and that the dynamics of such a flow behaves as a selective noise amplifier, the dynamic is mainly driven by receptivity mechanisms and by the response of upstream white nose disturbance. Two other cases have been studied on this work, the case of a transonic flow around a profile wing of NACA0012 type around the onset of buffet phenomenon and the case of transonic nozzle of Sajben type on over-extended regime. In the first case, the global stability analysis allows us to highlight the buffet phenomenon of on the profile NACA0012, what shows that the phenomenon is linked to a linear global instability. In the second case, the analysis of stability does not allow to explain the self-sustained low frequencies phenomenon, and shows that the flow is linearly globally stable. In this case, the dynamics is convective, passing and piloted by receptivity mechanisms.

Dans cette thèse, on a étudié l’évolution de l’interface par électro-osmose entre deux couches de fluides dans un canal microfluidique. Les applications de ce problème concernent le mélange et le transport, sans contact avec des actionneurs, de fluides en micro-canal. De nombreuses questions restent toutefois posées lorsque le champ est oscillant en temps, notamment vis à vis de la stabilité de l'interface entre les deux fluides. Une analyse de stabilité linéaire basée sur une perturbation à l’interface a été réalisée pour un film mince d'électrolyte sous des champs électriques continus (constants) et alternatifs (dépendant du temps). Une analyse asymptotique avec une hypothèse de grande longueur d’onde des équations d'Orr-Sommerfeld a été appliquée afin de déterminer les seuils de stabilité paramétriques d'un film mince aqueux. L’accent a été mis sur les effets de la tension de surface, de la pression de disjonction pour l'interaction gaz-liquide-substrat, de l'amplitude et de la fréquence du champ électrique appliqué, ainsi que du potentiel zêta du substrat et de la surface libre. Une analyse comparative des profils de vitesse de l’état de base avec et sans contraintes de Maxwell à l’interface, a montré que les gradients de vitesse étaient importants à l'interface liquide-liquide avec les contraintes de Maxwell. De tels gradients sont essentiels à l'instabilité interfaciale sous l’action d’un champ électrique périodique car ils peuvent atténuer ou amplifier les ondes à l’interface. Parallèlement, un dispositif expérimental a été conçu et monté afin de caractériser l’écoulement électroosmotique dans un micro-canal rectangulaire. Avec l'aide d'une analyse PTV (« Particle Tracking Velocimetry »), les distributions de vitesse ont été obtenues et comparées aux prédictions théoriques. Cette comparaison a permis d’estimer le potentiel zêta du PDMS utilisé, valeur conforme à la valeur indiquée dans la littérature. / Since the past decade, use of electro-osmotic flow (EOF) as an alternative flow mechanism in microdevices is becoming more popular due to its less bulky and low maintenance system design. However, one of the biggest shortcomings for its usage in mainstream applications is that it requires the concerned liquid to be electrically conductive. One idea can be to use the flow of conductive fluids to transport non-conductive liquids passively via interfacial shear transfer. Such an idea can has numerous applications in a wide range of fields like bio-chemical processing (e.g. lab-on-a-chip reactors, mixers, etc.), to oil extraction from porous rock formations. One of the significant characteristics of micro-scale flows is high surface to volume ratio, which significantly highlights the role of multi-phase interfaces in such dynamics. The presence of a fluid-fluid interface in an EOF necessitates the characterization of the parameters responsible for hydrodynamic instability of such systems. The present work focuses on the role of steady and time-dependent electric stress (Maxwell stress), capillary force and disjoining pressure on fluid-fluid interfacial instability. A linear stability analysis of interfacial perturbation was performed for a thin film of electrolyte under DC and AC electric fields. Through long wave asymptotic analysis of the Orr-Sommerfeld equations, parametric stability thresholds of a thin aqueous film explored. Further, a set of experiments were performed in order to characterize the EOF in a rectangular microchannel. With the help of a Particle Tracking Velocimetry analysis, velocity distributions were obtained which agreed well to the theoretical values. This was further used to estimate PDMS zeta potential, which was found to be within the reported values in the existing literature. Liquid-liquid interfacial deformation was also explored under a time-periodic EOF and a wide range of the magnitudes of capillary force, and diffusive and convective transport.

The thesis is focused on the investigation of Rayleigh-Bénard problem in an extended setting approximating the conditions in the Earth's mantle. The aim is to evaluate the influence of depth- and temperature- dependent material parameters, dissipation, adiabatic heating/cooling and heat sources on the qualitative characteristics of thermal convection. We identify the critical values of dimensionless parameters that determine the onset of convection and characterize the dominating convection patterns in marginally supercritical states. These issues are addressed by the application of linear stability analysis and weakly non-linear analysis. It has been found that the character of convection differ substantially from the standard case of Rayleigh-Bénard convection. Powered by TCPDF (www.tcpdf.org)

<div>The growing need to improve the power density of compact thermal systems necessitates developing new techniques to modulate the convective heat transfer efficiently. In the present research, acoustic streaming is evaluated as a potential technology to achieve this objective. Numerical simulations using the linearized and fully non-linear Navier-Stokes equations are employed to characterize the physics underlying this process. The linearized Navier-Stokes equations accurately replicate the low-frequency flow unsteadiness, which is used to find the optimal control parameters. Local and global stability analysis tools were developed to identify the modes with a global and positive heat transfer effect.</div><div><br></div><div>High-fidelity numerical simulations are performed to evaluate the effect of the excitation at selected frequencies, directed by the linear stability analysis, on the heat and momentum transport in the flow. Results indicate that, under favorable conditions, superimposing an acoustic wave, traveling along with the flow, can <i>resonate</i> within the domain and lead to a significant heat transfer enhancement with minimal skin friction losses. Two main flow configurations are considered; at the fixed Reynolds number Re_b=3000, in the supersonic case, 10.1% heat transfer enhancement is achieved by an 8.4% skin friction increase; however, in the subsonic case, 10% enhancement in heat transfer only caused a 5.3% increase to the skin friction. The deviation between these two quantities suggests a violation of the Reynolds analogy. This study is extended to include a larger Reynolds number, namely Re_b=6000 at M_b=0.75 and a similar response is observed. The effect of excitation amplitude and frequency on the resonance, limit-cycle oscillations, heat transfer, and skin friction are also investigated here.</div><div><br></div><div>Applying acoustic waves normal to the flow in the spanwise direction disrupts the near-wall turbulent structures that are primarily responsible for heat and momentum transport near the solid boundary. Direct numerical simulations were employed to investigate this technique in a supersonic channel flow at M_b=1.5 and Re_b=3000. The external excitation is applied through a periodic body force in the spanwise direction, mimicking loudspeakers placed on both walls that are operating with a 180^o phase shift. By keeping the product of forcing amplitude A_f and pulsation period (<i>T</i>) constant, spanwise velocity perturbations are generated with a similar amplitude at different frequencies. Under this condition, spanwise pulsations at <i>T</i>=20 and <i>T</i>=10 show up to 8% reduction in Nusselt number as well as the skin friction coefficient. Excitation at higher or lower frequencies fails to achieve such high level of modulations in heat and momentum transport processes near the walls.<br> <br>In configurations involving a spatially-developing boundary layer, a computational setup that includes laminar, transitional, and turbulent regions inside the domain is considered and the impact of acoustic excitation on this flow configuration has been characterized. Large-eddy simulations with dynamic Smagorinsky sub-grid scale modeling has been implemented, due to the excessive computational cost of DNS calculations at high-Reynolds numbers. The optimal excitation frequency that resembles the mode chosen for the fully-developed case has been identified via global stability analysis. Fully non-linear simulations of the spatially-developing boundary layer subjected to the excitation at this frequency reveal an interaction between the <i>pulsations</i> and the perturbations originated from the tripping which creates a re-laminarization zone traveling downstream. Such technique can locally enhance or reduce the heat transfer along the walls.<br></div>

Reaction-diffusion models have been widely successful in explaining a large variety of patterning phenomena in biology ranging from embryonic development to cancer growth and angiogenesis. Firstly proposed by Alan Turing in 1952 and applied to a simple two-component system, reaction-diffusion models describe spontaneous spatial pattern formation, driven purely by interactions of the system components and their diffusion in space. Today, access to unprecedented amounts of quantitative biological data allows us to build and test biochemically accurate reaction-diffusion models of intracellular processes. However, any increase in model complexity increases the number of unknown parameters and thus the computational cost of model analysis. To efficiently characterize the behavior and robustness of models with many unknown parameters is, therefore, a key challenge in systems biology. Here, we propose a novel computational framework for efficient high-dimensional parameter space characterization of reaction-diffusion models. The method leverages the $L_p$-Adaptation algorithm, an adaptive-proposal statistical method for approximate high-dimensional design centering and robustness estimation. Our approach is based on an oracle function, which describes for each point in parameter space whether the corresponding model fulfills given specifications. We propose specific oracles to estimate four parameter-space characteristics: bistability, instability, capability of spontaneous pattern formation, and capability of pattern maintenance. We benchmark the method and demonstrate that it allows exploring the ability of a model to undergo pattern-forming instabilities and to quantify model robustness for model selection in polynomial time with dimensionality. We present an application of the framework to reconstituted membrane domains bearing the small GTPase Rab5 and propose molecular mechanisms that potentially drive pattern formation.

No description available.

Plusieurs études tant numériques qu’expérimentales font état de la présence d’instabilités thermiques dans des films liquides chauffés uniformément par le bas pour des conditions aux limites et d’écoulements particuliers. La présence de ces instabilités modifiera les transferts thermiques associés. Le sujet de ce travail de thèse consiste à étudier numériquement les instabilités thermoconvectives d’un écoulement laminaire tridimensionnel de convection mixte dans un canal horizontal à surface libre. Les variations de la tension de surface avec la température (effet Marangoni ou effet thermocapillaire) sont prises en compte. Bien que d’un intérêt certain pour de nombreuses applications industrielles, cette situation a été très peu étudiée d’un point de vue académique dans la configuration considérée ici. Dans cette de configuration plusieurs types de structures thermoconvectives sont susceptibles d’apparaître. Lorsque les forces induites par les courants de convection naturelle, forcée et thermocapillaire sont du même ordre de grandeur, les premiers résultats montrent un développement des instabilités sous forme de rouleaux convectifs longitudinaux stationnaires semblables à ceux rencontrés pour des écoulements de type Poiseuille-Rayleigh-Bénard. A notre connaissance, c’est la première fois que l’écoulement de convection de type Poiseuille-Rayleigh-Bénard associé aux effets Marangoni est étudié. Le nombre et la distribution spatiale des rouleaux convectifs le long du canal dépendent des conditions de l’écoulement. Nous proposons une étude numérique pour ces conditions particulières d’écoulement pouvant conduire à des instabilités avec une évaluation de leur effet sur les transferts de chaleur. Les équations de Navier-Stokes et de l’énergie sont résolues numériquement par la méthode de volumes finis en prenant en compte les effets thermocapillaires. Les résultats présentés concernent l’influence des paramètres contrôlant l’écoulement (nombres de Reynolds, de Rayleigh, de Biot, de Marangoni et le rapport de forme) sur les motifs de l’écoulement et les échanges thermiques. Dans une seconde partie du travail, complémentaire à la première, une analyse de stabilité linéaire de l’écoulement dans un canal ouvert à surface libre d’extension latérale infinie est réalisée en utilisant la méthode spectrale de type collocation Chebyshev pour résoudre un système aux valeurs propres. Les diagrammes de stabilité déterminant les seuils des paramètres conduisant à l’instabilité thermoconvective ont été obtenus et analysés, ainsi que les structures spatiales associées. / Several studies both numerical and experimental have reported the presence of thermal instabilities in liquid films uniformly heated from below for specific boundary conditions and flows. The presence of these instabilities modifies the associated heat transfer. The subject of this PhD thesis is to study numerically the instability of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The variations of the surface tension with the temperature (Marangoni effect or thermocapillary effect) are taken into account. Although of great interest for many industrial applications, this problem has received little attention from an academic point of view. In this configuration, several types of thermoconvective structures may appear. When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable, the results show the development of instabilities in the form of steady longitudinal convective rolls similar to those encountered in the Poiseuille-Rayleigh-Bénard flow. To our knowledge, this is the first time that the Poiseuille-Rayleigh-Bénard flow associated to the Marangoni effects has been investigated. The number and spatial distribution of the convective rolls along the channel depend on the flow conditions. We propose a numerical study on the flow conditions that could lead to thermal instabilities with an evaluation of their effect on the heat transfer. The coupled Navier-Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects. The results presented concern the influence of several control parameters (the Reynolds, Rayleigh, Biot and Marangoni numbers and the aspect ratio of the channel) on the flow patterns and heat transfer characteristics. In the second part of this work, complimentary to the first, a linear stability analysis of a horizontal liquid film flowing in an open channel, with infinite lateral extension and uniform heating from below, is carried out. An eigenvalue problem is obtained in the course of this analysis which is solved numerically using the Chebyshev collocation spectral method. The stability diagrams determining the threshold parameters leading to thermoconvective instabilities were obtained and analyzed as well as the associated spatial patterns.

Para o projeto de estruturas com perfis de aço formados a frio, é fundamental a compreensão dos fenômenos da instabilidade local e global, uma vez que estes apresentam alta esbeltez e baixa rigidez à torção. A determinação do carregamento crítico e a identificação do modo de instabilidade contribuem para o entendimento do comportamento dessas estruturas. Este trabalho avalia três metodologias para a análise linear de estabilidade de perfis de aço formados a frio isolados, com o objetivo de determinar os carregamentos críticos elásticos de bifurcação e os modos de instabilidade associados. Estritamente, analisa-se perfis de seção U enrijecido e Z enrijecido isolados, de diversos comprimentos e diferentes condições de vinculação e carregamento. Determinam-se os carregamentos críticos elásticos de bifurcação e os modos de instabilidade globais e locais por meio de: (i) análise com o Método das Faixas Finitas (MFF), através do uso do programa computacional CUFSM; (ii) análise com elementos finitos de barra baseados na Teoria Generalizada de Vigas (MEF-GBT), via uso do programa GBTUL; e (iii) análise com elementos finitos de casca (MEF-cascas) por meio do uso do programa ABAQUS. Algumas restrições e ressalvas com relação ao uso do MFF são apresentadas, assim como limitações da Teoria Generalizada de Viga e precauções a serem tomadas nos modelos de cascas. Analisa-se também a influência do grau de discretização da seção transversal. No entanto, não é feita avaliação em relação aos procedimentos normativos e tampouco análises não lineares, considerando as imperfeições geométricas iniciais, tensões residuais e o comportamento elastoplástico do material. / For the design of cold formed steel members, it is essential to understand the effects of local and global instability, since these members typically have a high slenderness and low torsion stiffness. The determination of critical loads and the associated buckling modes contribute to understand the behavior of these members. This work performs a evaluation of three methods for linear stability analysis of isolated cold-formed steel members in order to determine the elastic critical loads and the corresponding buckling modes. Specifically, Ue and Ze shape members were studied with various length, different boundary conditions and loads. The elastic critical loads and buckling modes are determined by means of: (i) analysis with the Finite Strip Method (FSM), by the computer program CUFSM, (ii) beam finite element analysis based on the Generalized Beam Theory (FEM-GBT), by GBTUL program, and (iii) Finite Element Method with shell analysis using ABAQUS program. Some restrictions and warnings regarding the use of the FSM are presented, as well as limitations of the Generalized Beam Theory and precautions to be taken in the shell models. It is also analyzed the influence of the degree of discretization of the cross section. In the present study, no evaluation was made with respect to normative procedures neither nonlinear analyses considering the initial geometric imperfections, residual stresses and elastoplastic behavior of the material.

Αντικείμενο της παρούσας εργασίας είναι η πρόβλεψη τόσο της κίνησης, αλληλεπίδρασης και παραμόρφωσης δύο φυσαλίδων λόγω μεταβολής της πίεσης στο περιβάλλον ιξώδες υγρό, όσο και της ανοδικής κίνησης μιας φυσαλίδας λόγω άνωσης σε ένα Νευτωνικό ή ιξωδοπλαστικό ρευστό. Για τη μοντελοποίηση των αλληλεπιδρώντων φυσαλίδων, αναπτύχθηκε μιας νέα ελλειπτική μεθόδος κατασκευής του υπολογιστικού πλέγματος προκειμένου να αντιμετωπιστούν επιτυχώς τα ιδιάζοντα σημεία (πόλοι) των φυσαλίδων και οι μεγάλες παραμορφώσεις των διεπιφανειών τους. Με τη μέθοδο αυτή η πύκνωση του πλέγματος περιορίζεται μόνο στις περιοχές που είναι αναγκαίο, μειώνοντας έτσι το υπολογιστικό κόστος και αυξάνοντας την ακρίβεια των υπολογισμών. Για την επίλυση των παρακάτω προβλημάτων χρησιμοποιήθηκε η μέθοδος των μικτών πεπερασμένων στοιχείων κατά Galerkin. Στην περίπτωση των αλληλεπιδρώντων φυσαλίδων έχει εξετασθεί η επίδραση του σχετικού μεγέθους τους, της συχνότητας και του εύρους μεταβολής της επιβαλλόμενης πίεσης και πότε οδηγούν σε έλξη ή άπωση των φυσαλίδων. Στην περίπτωση ελκτικής δύναμης, ακολουθείται η κίνηση και η παραμόρφωσή τους μέχρι του σημείου που έρχονται σε επαφή, όπου αυτό είναι εφικτό. Για τη μελέτη του προβλήματος της φυσαλίδας που ανέρχεται λόγω άνωσης, υποθέτουμε αξονική συμμετρία και μόνιμη κατάσταση. Σύγκριση των προβλέψεών μας για το σχήμα των φυσαλίδων και το πεδίο ροής γύρω τους με προηγούμενα θεωρητικά και πειραματικά αποτελέσματα για Νευτωνικά ρευστά έδειξε άριστη συμφωνία. Στην περίπτωση του ιξωδοπλαστικού ρευστού εξετάστηκαν λεπτομερώς οι παραμορφώσεις των φυσαλίδων σαν συνάρτηση των αριθμών Bingham, Bond και Αρχιμήδη και υπολογίσθηκαν οι συνθήκες υπό τις οποίες είναι δυνατή η παγίδευση της φυσαλίδας μέσα σε αυτό. / The present study deals with the numerical simulation of the motion, interaction and deformation of two bubbles due to variation of the pressure of the ambient Newtonian fluid, and the buoyancy-driven rise of a bubble in a Newtonian or a viscoplastic fluid. A new elliptic mesh generation method is developed in order to deal with the singular points (poles) of the bubbles and the large deformations of their surface. This method permits us to increase the mesh resolution only in the regions that is necessary, decreasing thus the computational cost and increasing the precision of our calculations. The following problems are solved using the mixed finite element/Galerkin method. In the case of the interacting bubbles the effect of their relative size, the frequency and the width of the imposed pressure is examined as well as the conditions that lead in attraction or repulsion of the bubbles. In the case that attractive forces exist, the motion and the deformation of the bubbles followed up to the point that they come in contact, whenever this is possible. In order to study the problem of the bubble that rises due to buoyancy, axial symmetry and steady flow is assumed. Our results for the shape of the bubbles and the flow around them are in very good agreement with previous theoretical and experimental results for Newtonian fluids. The deformations of the bubbles rising in a viscoplastic material are also examined for various values of the Bingham, Bond and Archimedes numbers and the conditions under which entrapment of a bubble is possible are determined.