No description available.

We consider the flow of two immiscible fluids in a thin inclined channel subject to gravity and a change in pressure. In particular, we focus on the effects of Navier slip along the channel walls on the long-wave linear stability. Of interest are two different physical scenarios. The first corresponds to two incompressible fluid layers separated by a sharp interface, while the second focuses on a more dense fluid below a compressible gas. From a lubrication analysis, we find in the first scenario that the system is stable in the zero-Reynolds number limit with the slip effects modifying the decay rate of the stable perturbation. In the case of the Rayeligh-Taylor problem, slip along the less dense fluid wall has a destabilizing effect. In the second scenario, fluid inertia is pertinent, and we find neutral stability criteria are not significantly affected with the presence of slip.

This thesis is dedicated to the derivation and the validation of a new collision model as a stabilization technique for high-order lattice Boltzmann methods (LBM). More specifically, it intends to stabilize simulations of: (1) isothermal and weakly compressible flows at high Reynolds numbers, and (2) fully compressible flows including discontinuities such as shock waves. The new collision model relies on an enhanced regularization step. The latter includes a recursive computation of nonequilibrium Hermite polynomial coefficients. These recursive formulas directly derive from the Chapman-Enskog expansion, and allow to properly filter out second- (and higher-) order nonhydrodynamic contributions in underresolved conditions. This approach is even more interesting since it is compatible with a very large number of velocity sets. This high-order LBM is first validated in the isothermal case, and for high-Reynolds number flows. The coupling with a shock-capturing technique allows to further extend its validity domain to the simulation of fully compressible flows including shockwaves. The present work ends with the linear stability analysis(LSA) of the new approach, in the isothermal case. This leads to a proper quantification of the impact induced by each discretization (velocity and numerical) on the spectral properties of the related set of equations. The LSA of the recursive regularized LBM finally confirms the drastic stability gain obtained with this new approach.

Numerical methods for solving linearly damped Hamiltonian ordinary differential equations are analyzed and compared. The methods are constructed from the well-known Störmer-Verlet and implicit midpoint methods. The structure preservation properties of each method are shown analytically and numerically. Each method is shown to preserve a symplectic form up to a constant and are therefore conformal symplectic integrators, with each method shown to accurately preserve the rate of momentum dissipation. An analytical linear stability analysis is completed for each method, establishing thresholds between the value of the damping coefficient and the step-size that ensure stability. The methods are all second order and the preservation of the rate of energy dissipation is compared to that of a third order Runge-Kutta method that does not preserve conformal properties. Numerical experiments will include the damped harmonic oscillator and the damped nonlinear pendulum.

The linear receptivity and stability of plane idealized detonation with one-step Arrhenius type reaction kinetics is explored in the case of three-dimensional perturbations to a Zel'dovich-von Neumann-Doering base flow. This is explored in both overdriven and explicitly Chapman-Jouguet detonation. Additionally, the use of a multi-domain spectral collocation method for solving the conventional stability problem is explored within the context of normal-mode detonation. An extension of the stability analysis to confined detonations in a slightly porous walled tube is also carried out. Finally, an asymptotic analysis of a detonation with two-step reaction kinetics in the limit of large activation energy and for general overdrive and reaction order is performed yielding a nonlinear evolution equation for perturbations that produce stable limit cycle solutions.

This thesis provides a numerical and theoretical investigation of transitional and turbulent enclosed rotating flows, with a focus on the formation of macroscopic coherent flow structures. The underlying processes are strongly threedimensional due to the presence of boundary layers on the discs and on the walls of the outer (resp. inner) cylindrical shroud (resp. shaft). The complexity of these flows poses a great challenge in fundamental research however the present work is also of importance for industrial rotating machinery, from hard-drives to space engines turbopumps - the design issues of the latter being behind the motivation for this thesis. The present work consists of two major investigations. First, industrial cavities are modeled by smooth rotor/stator cavities and therein the dominant flow dynamics is investigated. For the experimental campaigns on industrial machinery revealed dangerous unsteady phenomena within the cavities, the emphasis is put on the reproduction and monitoring of unsteady pressure fluctuations within the smooth cavities. Then, the LES of three configurations of real industrial turbines are conducted to study in situ the pressure fluctuations and apply the diagnostics already vetted on academic problems.

At the microscopic level, many materials are made of smaller and randomly oriented grains. These grains are separated by grain boundaries which tend to decrease the electrical and thermal conductivity of the material. The motion of grain boundaries is an important phenomenon controlling the grain growth in materials processing and synthesis. Mathematical modeling and simulation is a powerful tool for studying the motion of grain boundaries. The research reported in this thesis is focused on the numerical simulation and analysis of a coupled surface and grain boundary motion which models the evolution of grain boundary and the diffusion of the free surface during the process of grain growth. The “quarter loop” geometry provides a convenient model for the study of this coupled motion. Two types of normal curve velocities are involved in this model: motion by mean curvature and motion by surface diffusion. They are coupled together at a triple junction. A front tracking method is used to simulate the migration. To describe the problem, different formulations are presented and discussed. A new formulation that comprises partial differential equations and algebraic equations is proposed. It preserves arc length parametrization up to scaling and exhibits good numerical performance. This formulation is shown to be well-posed in a reduced, linear setting. Numerical simulations are implemented and compared for all formulations. The new formulation is also applied to some other related problems. We investigate numerically the linear stability of the travelling wave solutions for the quarter loop problem and a simple grain boundary motion problem for both curves in two dimensions and surfaces in three dimensions. The numerical results give evidence that they are convectively stable. A class of high order three-phase boundary motion problems are also studied. We consider a region where three phase boundaries meet at a triple junction and evolve with specified normal velocities. A system of partial differential algebraic equations (PDAE) is proposed to describe this class of problems by extending the discussion for the coupled surface and grain boundary motion. The linear well-posedness of the system is analyzed and numerical simulations are performed.

At the microscopic level, many materials are made of smaller and randomly oriented grains. These grains are separated by grain boundaries which tend to decrease the electrical and thermal conductivity of the material. The motion of grain boundaries is an important phenomenon controlling the grain growth in materials processing and synthesis. Mathematical modeling and simulation is a powerful tool for studying the motion of grain boundaries. The research reported in this thesis is focused on the numerical simulation and analysis of a coupled surface and grain boundary motion which models the evolution of grain boundary and the diffusion of the free surface during the process of grain growth. The “quarter loop” geometry provides a convenient model for the study of this coupled motion. Two types of normal curve velocities are involved in this model: motion by mean curvature and motion by surface diffusion. They are coupled together at a triple junction. A front tracking method is used to simulate the migration. To describe the problem, different formulations are presented and discussed. A new formulation that comprises partial differential equations and algebraic equations is proposed. It preserves arc length parametrization up to scaling and exhibits good numerical performance. This formulation is shown to be well-posed in a reduced, linear setting. Numerical simulations are implemented and compared for all formulations. The new formulation is also applied to some other related problems. We investigate numerically the linear stability of the travelling wave solutions for the quarter loop problem and a simple grain boundary motion problem for both curves in two dimensions and surfaces in three dimensions. The numerical results give evidence that they are convectively stable. A class of high order three-phase boundary motion problems are also studied. We consider a region where three phase boundaries meet at a triple junction and evolve with specified normal velocities. A system of partial differential algebraic equations (PDAE) is proposed to describe this class of problems by extending the discussion for the coupled surface and grain boundary motion. The linear well-posedness of the system is analyzed and numerical simulations are performed.

At the microscopic level, many materials are made of smaller and randomly oriented grains. These grains are separated by grain boundaries which tend to decrease the electrical and thermal conductivity of the material. The motion of grain boundaries is an important phenomenon controlling the grain growth in materials processing and synthesis. Mathematical modeling and simulation is a powerful tool for studying the motion of grain boundaries. The research reported in this thesis is focused on the numerical simulation and analysis of a coupled surface and grain boundary motion which models the evolution of grain boundary and the diffusion of the free surface during the process of grain growth. The “quarter loop” geometry provides a convenient model for the study of this coupled motion. Two types of normal curve velocities are involved in this model: motion by mean curvature and motion by surface diffusion. They are coupled together at a triple junction. A front tracking method is used to simulate the migration. To describe the problem, different formulations are presented and discussed. A new formulation that comprises partial differential equations and algebraic equations is proposed. It preserves arc length parametrization up to scaling and exhibits good numerical performance. This formulation is shown to be well-posed in a reduced, linear setting. Numerical simulations are implemented and compared for all formulations. The new formulation is also applied to some other related problems. We investigate numerically the linear stability of the travelling wave solutions for the quarter loop problem and a simple grain boundary motion problem for both curves in two dimensions and surfaces in three dimensions. The numerical results give evidence that they are convectively stable. A class of high order three-phase boundary motion problems are also studied. We consider a region where three phase boundaries meet at a triple junction and evolve with specified normal velocities. A system of partial differential algebraic equations (PDAE) is proposed to describe this class of problems by extending the discussion for the coupled surface and grain boundary motion. The linear well-posedness of the system is analyzed and numerical simulations are performed. / Science, Faculty of / Mathematics, Department of / Graduate

Previous experimental and numerical studies showed that two-dimensional roughness elements can stabilize disturbances inside a hypersonic boundary layer, and eventually delay the transition onset. The objective of the thesis is to evaluate the response of disturbances propagating inside a hypersonic boundary layer to various two-dimensional surface deformations of different shapes. The proposed deformations consist of a backward step, forward step, a combination of backward and forward steps, two types of wavy surfaces, surface dips or surface humps. Disturbances inside of a Mach 5.92 flat-plate boundary layer are excited by pulsed or periodic wall blowing and suction at an upstream location. The numerical tools consist of the Navier-Stokes equations in curvilinear coordinates and a linear stability analysis tool. Results show that all types of surface deformations are able to reduce the amplitude of boundary layer disturbances to a certain degree. The amount of reduction in the disturbance energy is related to the type of pressure gradient created by the deformation, adverse or favorable.