Pattern generation by dissipative parametric instability

Nonlinear instabilities are responsible for spontaneous pattern formation in a vast number of natural and engineered systems, ranging from biology to galaxy buildup. We propose a new instability mechanism leading to pattern formation in spatially extended nonlinear systems, which is based on a periodic antiphase modulation of spectrally dependent losses arranged in a zigzag way: an effective filtering is imposed at symmetrically located wave numbers k and -k in alternating order. The properties of the dissipative parametric instability differ from the features of both key classical concepts of modulation instabilities, i.e., the Benjamin-Feir instability and the Faraday instabiltyity. We demonstrate how the dissipative parametric instability can lead to the formation of stable patterns in one- and two-dimensional systems. The proposed instability mechanism is generic and can naturally occur or can be implemented in various physical systems.

In-cavity pulse shaping and dissipative parametric instability mode-locking in fibre lasers

At the level of fundamental research, fibre lasers provide convenient and reproducible experimental settings for the study of a variety of nonlinear dynamical processes, while at the applied research level, pulses with different and optimised features – e.g., in terms of pulse duration, temporal and/or spectral shape, energy, repetition rate and emission bandwidth – are sought with the general constraint of developing efficient cavity architectures. In this work, we review our recent progress on the realisation of pulse shaping in passively- mode-locked fibre lasers by inclusion of an amplitude and phase spectral filter into the laser cavity. We present a fibre laser design in which pulse shaping occurs through filtering of a spectrally nonlinearly broadened pulse in the cavity. This strategy of pulse shaping is illustrated through the numerical demonstration of the laser operation in different pulse-generation regimes, including parabolic, flattop and triangular waveform generations, depending on the amplitude profile of the in-cavity spectral filter [1]. As an application of this general approach, we show that the use of an in-cavity flat-top spectral filter makes it possible to directly generate sinc-shaped Nyquist pulses of high quality and of a widely tunable bandwidth from the laser [2]. We also report on a recently-developed versatile erbium-doped fibre laser, in which conventional soliton, dispersion-managed soliton (stretched-pulse) and dissipative soliton mode-locking regimes can be selectively and reliably targeted by programming different group-velocity dispersion profiles and bandwidths on an in-cavity programmable filter [3]. Further, we report on our recent results on the passive mode locking of a Raman fibre laser by a recently predicted new type of parametric instability – the dissipative Faraday instability [4], where spatially periodic zig-zag modulation of spectrally dependent losses can lead to pattern formation in the temporal domain. High-order harmonic mode locking is achieved in a very simple experimental configuration, with the laser cavity including an optical fibre and two chirped fibre Bragg gratings, and no additional mode-locking elements. The results not only open up new possibilities for the design of mode-locked lasers, but extend beyond fibre optics to other fields of physics and engineering. References [1] S. Boscolo, C. Finot, H. Karakuzu, P. Petropoulos, “Pulse shaping in mode-locked fiber laser by in-cavity spectral filter,” Opt. Lett., vol. 39, pp. 438–441, 2014. [2] S. Boscolo, C. Finot, S. K. Turitsyn, “Bandwidth programmable optical Nyquist pulse generation in passively mode-locked fiber laser,” IEEE Photon. J., vol. 7, 7802008(8), 2015. [3] J. Peng, S. Boscolo, “Filter-based dispersion-managed versatile ultrafast fibre laser,” Sci. Rep., 2016, In press. [4] A. M. Perego, N. Tarasov, D. V. Churkin, S. K. Turitsyn, K. Staliunas, “Pattern generation by dissipative parametric instability,” Phys. Rev. Lett., vol. 116, 028701, 2016.

Self-pulsating nonlinear systems via dissipative parametric instability

The recently discovered dissipative parametric instability is presented in the framework of the universal complex Ginzburg-Landau equation. The pattern formation associated with the instability is discussed in connection to the relevant applications in nonlinear photonics especially as a new tool for pulsed lasers design.

Dissipative parametric modulation instability and pattern formation in nonlinear optical systems

We present the essential features of the dissipative parametric instability, in the universal complex Ginzburg- Landau equation. Dissipative parametric instability is excited through a parametric modulation of frequency dependent losses in a zig-zag fashion in the spectral domain. Such damping is introduced respectively for spectral components in the +ΔF and in the -ΔF region in alternating fashion, where F can represent wavenumber or temporal frequency depending on the applications. Such a spectral modulation can destabilize the homogeneous stationary solution of the system leading to growth of spectral sidebands and to the consequent pattern formation: both stable and unstable patterns in one- and in two-dimensional systems can be excited. The dissipative parametric instability provides an useful and interesting tool for the control of pattern formation in nonlinear optical systems with potentially interesting applications in technological applications, like the design of mode- locked lasers emitting pulse trains with tunable repetition rate; but it could also find realizations in nanophotonics circuits or in dissipative polaritonic Bose-Einstein condensates.

Dynamic stability of laminated FGM plates based on higher-order shear deformation theory

This paper conducts a dynamic stability analysis of symmetrically laminated FGM rectangular plates with general out-of-plane supporting conditions, subjected to a uniaxial periodic in-plane load and undergoing uniform temperature change. Theoretical formulations are based on Reddy's third-order shear deformation plate theory, and account for the temperature dependence of material properties. A semi-analytical Galerkin-differential quadrature approach is employed to convert the governing equations into a linear system of Mathieu-Hill equations from which the boundary points on the unstable regions are determined by Bolotin's method. Free vibration and bifurcation buckling are also discussed as subset problems. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with FGM layers made of silicon nitride and stainless steel. The influences of various parameters such as material composition, layer thickness ratio, temperature change, static load level, boundary constraints on the dynamic stability, buckling and vibration frequencies are examined in detail through parametric studies.

Fast magnetosonic wave excitation by an array of wires with time-modulated currents

The excitation of Fast Magnetosonic (FMS)waves by a cylindrical array of parallel tethers carrying timemodulated current is discussed. The tethers would fly vertical in the equatorial plane, which is perpendicular to the geomagnetic field when its tilt is ignored, and would be stabilized by the gravity gradient. The tether array would radiate a single FMS wave. In the time-dependent background made of geomagnetic field plus radiated wave, plasma FMS perturbations are excited in the array vicinity through a parametric instability. The growth rate is estimated by truncating the evolution equation for FMS perturbations to the two azimuthal modes of lowest order. Design parameters such as tether length and number, required power and mass are discussed for Low Earth Orbit conditions. The array-attached wave structure would have the radiated wave controlled by the intensity and modulation frequency of the currents, making an active experiment on non-linear low frequency waves possible in real space plasma conditions.

Magnetic pumping of whistler waves by tether current modulation

Magnetic excitation of whistlers by a square array of electrodynamic tethers is discussed. The array is made of perpendicular rows of tethers that carry equal, uniform, and time-modulated currents at equal frequency with a 90° phase shift. The array would fly vertical in the orbital equatorial plane, which is perpendicular to the geomagnetic field B0 when its tilt is ignored. The array radiates a whistler wave along B0. A parametric instability due to pumping by the background magnetic field through the radiated wave gives rise to two unstable coupled whistler perturbations. The growth rate is maximum for perturbations with wave vector at angles 38.36° and 75.93° from B0. For an experiment involving a wavefront that moves with the orbiting array, which might serve to study nonlinear wave interactions and turbulence in space plasmas, characteristic values of growth rate and parameters, such as the number of tethers and their dimensions and distances in the array, are discussed for low Earth orbit ambient conditions.

A model for turbulent combustion simulation of large scale hydrogen explosions

El agotamiento, la ausencia o, simplemente, la incertidumbre sobre la cantidad de las reservas de combustibles fósiles se añaden a la variabilidad de los precios y a la creciente inestabilidad en la cadena de aprovisionamiento para crear fuertes incentivos para el desarrollo de fuentes y vectores energéticos alternativos. El atractivo de hidrógeno como vector energético es muy alto en un contexto que abarca, además, fuertes inquietudes por parte de la población sobre la contaminación y las emisiones de gases de efecto invernadero. Debido a su excelente impacto ambiental, la aceptación pública del nuevo vector energético dependería, a priori, del control de los riesgos asociados su manipulación y almacenamiento. Entre estos, la existencia de un innegable riesgo de explosión aparece como el principal inconveniente de este combustible alternativo. Esta tesis investiga la modelización numérica de explosiones en grandes volúmenes, centrándose en la simulación de la combustión turbulenta en grandes dominios de cálculo en los que la resolución que es alcanzable está fuertemente limitada. En la introducción, se aborda una descripción general de los procesos de explosión. Se concluye que las restricciones en la resolución de los cálculos hacen necesario el modelado de los procesos de turbulencia y de combustión. Posteriormente, se realiza una revisión crítica de las metodologías disponibles tanto para turbulencia como para combustión, que se lleva a cabo señalando las fortalezas, deficiencias e idoneidad de cada una de las metodologías. Como conclusión de esta investigación, se obtiene que la única estrategia viable para el modelado de la combustión, teniendo en cuenta las limitaciones existentes, es la utilización de una expresión que describa la velocidad de combustión turbulenta en función de distintos parámetros. Este tipo de modelos se denominan Modelos de velocidad de llama turbulenta y permiten cerrar una ecuación de balance para la variable de progreso de combustión. Como conclusión también se ha obtenido, que la solución más adecuada para la simulación de la turbulencia es la utilización de diferentes metodologías para la simulación de la turbulencia, LES o RANS, en función de la geometría y de las restricciones en la resolución de cada problema particular. Sobre la base de estos hallazgos, el crea de un modelo de combustión en el marco de los modelos de velocidad de la llama turbulenta. La metodología propuesta es capaz de superar las deficiencias existentes en los modelos disponibles para aquellos problemas en los que se precisa realizar cálculos con una resolución moderada o baja. Particularmente, el modelo utiliza un algoritmo heurístico para impedir el crecimiento del espesor de la llama, una deficiencia que lastraba el célebre modelo de Zimont. Bajo este enfoque, el énfasis del análisis se centra en la determinación de la velocidad de combustión, tanto laminar como turbulenta. La velocidad de combustión laminar se determina a través de una nueva formulación capaz de tener en cuenta la influencia simultánea en la velocidad de combustión laminar de la relación de equivalencia, la temperatura, la presión y la dilución con vapor de agua. La formulación obtenida es válida para un dominio de temperaturas, presiones y dilución con vapor de agua más extenso de cualquiera de las formulaciones previamente disponibles. Por otra parte, el cálculo de la velocidad de combustión turbulenta puede ser abordado mediante el uso de correlaciones que permiten el la determinación de esta magnitud en función de distintos parámetros. Con el objetivo de seleccionar la formulación más adecuada, se ha realizado una comparación entre los resultados obtenidos con diversas expresiones y los resultados obtenidos en los experimentos. Se concluye que la ecuación debida a Schmidt es la más adecuada teniendo en cuenta las condiciones del estudio. A continuación, se analiza la importancia de las inestabilidades de la llama en la propagación de los frentes de combustión. Su relevancia resulta significativa para mezclas pobres en combustible en las que la intensidad de la turbulencia permanece moderada. Estas condiciones son importantes dado que son habituales en los accidentes que ocurren en las centrales nucleares. Por ello, se lleva a cabo la creación de un modelo que permita estimar el efecto de las inestabilidades, y en concreto de la inestabilidad acústica-paramétrica, en la velocidad de propagación de llama. El modelado incluye la derivación matemática de la formulación heurística de Bauwebs et al. para el cálculo de la incremento de la velocidad de combustión debido a las inestabilidades de la llama, así como el análisis de la estabilidad de las llamas con respecto a una perturbación cíclica. Por último, los resultados se combinan para concluir el modelado de la inestabilidad acústica-paramétrica. Tras finalizar esta fase, la investigación se centro en la aplicación del modelo desarrollado en varios problemas de importancia para la seguridad industrial y el posterior análisis de los resultados y la comparación de los mismos con los datos experimentales correspondientes. Concretamente, se abordo la simulación de explosiones en túneles y en contenedores, con y sin gradiente de concentración y ventilación. Como resultados generales, se logra validar el modelo confirmando su idoneidad para estos problemas. Como última tarea, se ha realizado un analisis en profundidad de la catástrofe de Fukushima-Daiichi. El objetivo del análisis es determinar la cantidad de hidrógeno que explotó en el reactor número uno, en contraste con los otros estudios sobre el tema que se han centrado en la determinación de la cantidad de hidrógeno generado durante el accidente. Como resultado de la investigación, se determinó que la cantidad más probable de hidrogeno que fue consumida durante la explosión fue de 130 kg. Es un hecho notable el que la combustión de una relativamente pequeña cantidad de hidrogeno pueda causar un daño tan significativo. Esta es una muestra de la importancia de este tipo de investigaciones. Las ramas de la industria para las que el modelo desarrollado será de interés abarca la totalidad de la futura economía de hidrógeno (pilas de combustible, vehículos, almacenamiento energético, etc) con un impacto especial en los sectores del transporte y la energía nuclear, tanto para las tecnologías de fisión y fusión. ABSTRACT The exhaustion, absolute absence or simply the uncertainty on the amount of the reserves of fossil fuels sources added to the variability of their prices and the increasing instability and difficulties on the supply chain are strong incentives for the development of alternative energy sources and carriers. The attractiveness of hydrogen in a context that additionally comprehends concerns on pollution and emissions is very high. Due to its excellent environmental impact, the public acceptance of the new energetic vector will depend on the risk associated to its handling and storage. Fromthese, the danger of a severe explosion appears as the major drawback of this alternative fuel. This thesis investigates the numerical modeling of large scale explosions, focusing on the simulation of turbulent combustion in large domains where the resolution achievable is forcefully limited. In the introduction, a general description of explosion process is undertaken. It is concluded that the restrictions of resolution makes necessary the modeling of the turbulence and combustion processes. Subsequently, a critical review of the available methodologies for both turbulence and combustion is carried out pointing out their strengths and deficiencies. As a conclusion of this investigation, it appears clear that the only viable methodology for combustion modeling is the utilization of an expression for the turbulent burning velocity to close a balance equation for the combustion progress variable, a model of the Turbulent flame velocity kind. Also, that depending on the particular resolution restriction of each problem and on its geometry the utilization of different simulation methodologies, LES or RANS, is the most adequate solution for modeling the turbulence. Based on these findings, the candidate undertakes the creation of a combustion model in the framework of turbulent flame speed methodology which is able to overcome the deficiencies of the available ones for low resolution problems. Particularly, the model utilizes a heuristic algorithm to maintain the thickness of the flame brush under control, a serious deficiency of the Zimont model. Under the approach utilized by the candidate, the emphasis of the analysis lays on the accurate determination of the burning velocity, both laminar and turbulent. On one side, the laminar burning velocity is determined through a newly developed correlation which is able to describe the simultaneous influence of the equivalence ratio, temperature, steam dilution and pressure on the laminar burning velocity. The formulation obtained is valid for a larger domain of temperature, steam dilution and pressure than any of the previously available formulations. On the other side, a certain number of turbulent burning velocity correlations are available in the literature. For the selection of the most suitable, they have been compared with experiments and ranked, with the outcome that the formulation due to Schmidt was the most adequate for the conditions studied. Subsequently, the role of the flame instabilities on the development of explosions is assessed. Their significance appears to be of importance for lean mixtures in which the turbulence intensity remains moderate. These are important conditions which are typical for accidents on Nuclear Power Plants. Therefore, the creation of a model to account for the instabilities, and concretely, the acoustic parametric instability is undertaken. This encloses the mathematical derivation of the heuristic formulation of Bauwebs et al. for the calculation of the burning velocity enhancement due to flame instabilities as well as the analysis of the stability of flames with respect to a cyclic velocity perturbation. The results are combined to build a model of the acoustic-parametric instability. The following task in this research has been to apply the model developed to several problems significant for the industrial safety and the subsequent analysis of the results and comparison with the corresponding experimental data was performed. As a part of such task simulations of explosions in a tunnel and explosions in large containers, with and without gradient of concentration and venting have been carried out. As a general outcome, the validation of the model is achieved, confirming its suitability for the problems addressed. As a last and final undertaking, a thorough study of the Fukushima-Daiichi catastrophe has been carried out. The analysis performed aims at the determination of the amount of hydrogen participating on the explosion that happened in the reactor one, in contrast with other analysis centered on the amount of hydrogen generated during the accident. As an outcome of the research, it was determined that the most probable amount of hydrogen exploding during the catastrophe was 130 kg. It is remarkable that the combustion of such a small quantity of material can cause tremendous damage. This is an indication of the importance of these types of investigations. The industrial branches that can benefit from the applications of the model developed in this thesis include the whole future hydrogen economy, as well as nuclear safety both in fusion and fission technology.

Temporal and spatio-temporal regimes of generation of Raman fibre lasers

Temporal dynamics of Raman fibre lasers tend to have very complex nature, owing to great cavity lengths and high nonlinearity, being stochastic on short time scales and quasi-continuous on longer time scales. Generally fibre laser intensity dynamics is represented by one-dimensional time-series, which in case of quasi-continuous wave generation in Raman fibre lasers gives little insight into the processes underlying the operation of a laser. New methods of analysis and data representation could help to uncover the underlying physical processes, understand the dynamics or improve the performance of the system. Using intrinsic periodicity of laser radiation, one dimensional intensity time series of a Raman fibre laser was analysed over fast and slow variation time. This allowed to experimentally observe various spatio-temporal regimes of generation, such as laminar, turbulent, partial mode-lock, as well as transitions between them and identify the mechanisms responsible for the transitions. Great cavity length and high nonlinearity also make it difficult to achieve stable high repetition rate mode-locking in Raman fibre lasers. Using Faraday parametric instability in extremely simple linear cavity experimental configuration, a very high order harmonic mode-locking was achieved in ò.ò kmlong Raman fibre laser. The maximum achieved pulse repetition rate was 12 GHz, with 7.3 ps long Gaussian shaped pulses. There is a new type of random lasers – random distributed feedback Raman fibre laser, which temporal properties cannot be controlled by conventionalmode-locking or Q-switch techniques and mechanisms. By adjusting the pump configuration, a very stable pulsed operation of random distributed feedback Raman fibre laser was achieved. Pulse duration varied in the range from 50 to 200 μs depending on the pump power and the cavity length. Pulse repetition rate scaling on the parameters of the system was experimentally identified.

Parametric resonances in a base-excited double pendulum

Two parametrically-induced phenomena are addressed in the context of a double pendulum subject to a vertical base excitation. First, the parametric resonances that cause the stable downward vertical equilibrium to bifurcate into large-amplitude periodic solutions are investigated extensively. Then the stabilization of the unstable upward equilibrium states through the parametric action of the high-frequency base motion is documented in the experiments and in the simulations. It is shown that there is a region in the plane of the excitation frequency and amplitude where all four unstable equilibrium states can be stabilized simultaneously in the double pendulum. The parametric resonances of the two modes of the base-excited double pendulum are studied both theoretically and experimentally. The transition curves (i.e., boundaries of the dynamic instability regions) are constructed asymptotically via the method of multiple scales including higher-order effects. The bifurcations characterizing the transitions from the trivial equilibrium to the periodic solutions are computed by either continuation methods and or by time integration and compared with the theoretical and experimental results.

Modeling hydrothermal system: deriving observables and hydrothermal instability in volcanic and non-volcanic setting

Hydrothermal fluids are a fundamental resource for understanding and monitoring volcanic and non-volcanic systems. This thesis is focused on the study of hydrothermal system through numerical modeling with the geothermal simulator TOUGH2. Several simulations are presented, and geophysical and geochemical observables, arising from fluids circulation, are analyzed in detail throughout the thesis. In a volcanic setting, fluids feeding fumaroles and hot spring may play a key role in the hazard evaluation. The evolution of the fluids circulation is caused by a strong interaction between magmatic and hydrothermal systems. A simultaneous analysis of different geophysical and geochemical observables is a sound approach for interpreting monitored data and to infer a consistent conceptual model. Analyzed observables are ground displacement, gravity changes, electrical conductivity, amount, composition and temperature of the emitted gases at surface, and extent of degassing area. Results highlight the different temporal response of the considered observables, as well as the different radial pattern of variation. However, magnitude, temporal response and radial pattern of these signals depend not only on the evolution of fluid circulation, but a main role is played by the considered rock properties. Numerical simulations highlight differences that arise from the assumption of different permeabilities, for both homogeneous and heterogeneous systems. Rock properties affect hydrothermal fluid circulation, controlling both the range of variation and the temporal evolution of the observable signals. Low temperature fumaroles and low discharge rate may be affected by atmospheric conditions. Detailed parametric simulations were performed, aimed to understand the effects of system properties, such as permeability and gas reservoir overpressure, on diffuse degassing when air temperature and barometric pressure changes are applied to the ground surface. Hydrothermal circulation, however, is not only a characteristic of volcanic system. Hot fluids may be involved in several mankind problems, such as studies on geothermal engineering, nuclear waste propagation in porous medium, and Geological Carbon Sequestration (GCS). The current concept for large-scale GCS is the direct injection of supercritical carbon dioxide into deep geological formations which typically contain brine. Upward displacement of such brine from deep reservoirs driven by pressure increases resulting from carbon dioxide injection may occur through abandoned wells, permeable faults or permeable channels. Brine intrusion into aquifers may degrade groundwater resources. Numerical results show that pressure rise drives dense water up to the conduits, and does not necessarily result in continuous flow. Rather, overpressure leads to new hydrostatic equilibrium if fluids are initially density stratified. If warm and salty fluid does not cool passing through the conduit, an oscillatory solution is then possible. Parameter studies delineate steady-state (static) and oscillatory solutions.

Effect of Aspect Ratio on the Three-Dimensional Global Instability Analysis of Incompressible Open Cavity Flows

The stability analysis of open cavity flows is a problem of great interest in the aeronautical industry. This type of flow can appear, for example, in landing gears or auxiliary power unit configurations. Open cavity flows is very sensitive to any change in the configuration, either physical (incoming boundary layer, Reynolds or Mach numbers) or geometrical (length to depth and length to width ratio). In this work, we have focused on the effect of geometry and of the Reynolds number on the stability properties of a threedimensional spanwise periodic cavity flow in the incompressible limit. To that end, BiGlobal analysis is used to investigate the instabilities in this configuration. The basic flow is obtained by the numerical integration of the Navier-Stokes equations with laminar boundary layers imposed upstream. The 3D perturbation, assumed to be periodic in the spanwise direction, is obtained as the solution of the global eigenvalue problem. A parametric study has been performed, analyzing the stability of the flow under variation of the Reynolds number, the L/D ratio of the cavity, and the spanwise wavenumber β. For consistency, multidomain high order numerical schemes have been used in all the computations, either basic flow or eigenvalue problems. The results allow to define the neutral curves in the range of L/D = 1 to L/D = 3. A scaling relating the frequency of the eigenmodes and the length to depth ratio is provided, based on the analysis results.

Instability Analysis of Incompressible Open Cavity Flows

El problema del flujo sobre una cavidad abierta ha sido estudiado en profundidad en la literatura, tanto por el interés académico del problema como por sus aplicaciones prácticas en gran variedad de problemas ingenieriles, como puede ser el alojamiento del tren de aterrizaje de aeronaves, o el depósito de agua de aviones contraincendios. Desde hace muchos a˜nos se estudian los distintos tipos de inestabilidades asociadas a este problema: los modos bidimensionales en la capa de cortadura, y los modos tridimensionales en el torbellino de recirculación principal dentro de la cavidad. En esta tesis se presenta un estudio paramétrico completo del límite incompresible del problema, empleando la herramienta de estabilidad lineal conocida como BiGlobal. Esta aproximación permite contemplar la estabilidad global del flujo, y obtener tanto la forma como las características de los modos propios del problema físico, sean estables o inestables. El estudio realizado permite caracterizar con gran detalle todos los modos relevantes, así como la envolvente de estabilidad en el espacio paramétrico del problema incompresible (Mach nulo, variación de Reynolds, espesor de capa límite incidente, relación altura/profundidad de la cavidad, y longitud característica de la perturbación en la dirección transversal). A la luz de los resultados obtenidos se proponen una serie de relaciones entre los parámetros y características de los modos principales, como por ejemplo entre el Reynolds crítico de un modo, y la longitud característica del mismo. Los resultados numéricos se contrastan con una campaña experimental, siendo la principal conclusión de dicha comparación que los modos lineales están presentes en el flujo real saturado, pero que existen diferencias notables en frecuencia entre las predicciones teóricas y los experimentos. Para intentar determinar la naturaleza de dichas diferencias se realiza una simulación numérica directa tridimensional, y se utiliza un algoritmo de DMD (descomposición dinámica de modos) para describir el proceso de saturación. ABSTRACT The problem of the flow over an open cavity has been studied in depth in the literature, both for being an interesting academical problem and due to the multitude of industrial applications, like the landing gear of aircraft, or the water deposit of firefighter airplanes. The different types of instabilities appearing in this flow studied in the literature are two: the two-dimensional shear layer modes, and the three-dimensional modes that appear in the main recirculating vortex inside the cavity. In this thesis a parametric study in the incompressible limit of the problem is presented, using the linear stability analysis known as BiGlobal. This approximation allows to obtain the global stability behaviour of the flow, and to capture both the morphological features and the characteristics of the eigenmodes of the physical problem, whether they are stable or unstable. The study presented here characterizes with great detail all the relevant eigenmodes, as well as the hypersurface of instability on the parameter space of the incompressible problem (Mach equal to zero, and variation of the Reynolds number, the incoming boundary layer thickness, the length to depth aspect ratio of the cavity and the spanwise length of the perturbation). The results allow to construct parametric relations between the characteristics of the leading eigenmodes and the parameters of the problem, like for example the one existing between the critical Reynolds number and its characteristic length. The numerical results presented here are compared with those of an experimental campaign, with the main conclusion of said comparison being that the linear eigenmode are present in the real saturated flow, albeit with some significant differences in the frequencies of the experiments and those predicted by the theory. To try to determine the nature of those differences a three-dimensional direct numerical simulation, analyzed with Dynamic Mode Decomposition algorithm, was used to describe the process of saturation.

A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets

Numerical simulations are conducted to investigate how a droplet of Newtonian liquid. entrained in a higher viscosity Newtonian liquid, behaves when passing through an axisymmetric microfluidic contraction. Simulations are performed using a transient Volume of Fluid finite volume algorithm, and cover ranges of Reynolds and Weber numbers relevant to microfluidic flows. Results are presented for a droplet to surrounding fluid viscosity ratio of 0.001. In contrast to behaviour at higher viscosity ratios obtained previously by the authors, shear and interfacial tension driven instabilities often develop along the droplet Surface. leading to complex shape development, and in some instances, droplet breakup. (c) 2006 Elsevier Ltd. All rights reserved.

The role of pinning and instability in a class of non-equilibrium growth models

We study the dynamics of a growing crystalline facet where the growth mechanism is controlled by the geometry of the local curvature. A continuum model, in (2+1) dimensions, is developed in analogy with the Kardar-Parisi-Zhang (KPZ) model is considered for the purpose. Following standard coarse graining procedures, it is shown that in the large time, long distance limit, the continuum model predicts a curvature independent KPZ phase, thereby suppressing all explicit effects of curvature and local pinning in the system, in the "perturbative" limit. A direct numerical integration of this growth equation, in 1+1 dimensions, supports this observation below a critical parametric range, above which generic instabilities, in the form of isolated pillared structures lead to deviations from standard scaling behaviour. Possibilities of controlling this instability by introducing statistically "irrelevant" (in the sense of renormalisation groups) higher ordered nonlinearities have also been discussed.