The work presented in this thesis studies acoustic perturbations in slowly varying pipes. The slow variation is introduced in the form of a small parameter ${\epsilon}$ and through this in turn gives rise to a slow axial scale $X$ such that $X = {\epsilon}x$ where $x$ is the normal axial coordinate. This allows an asymptotic approach and the WKB method is used to solve the subsequent mathematical problems. The first deals with the existence of a trapped mode in a hard-walled pipe of varying radius conveying fluid. For the derived leading order propagating mode solution, its amplitude becomes singular at transition points $X_{t}$ and $X_{t'}$ where $X_{t} > 0$ and $X_{t'} < 0$ and thus is unable to propagate past these points. Because of the break down in the solution, this leads to the theory that in the neighbourhood of these points there exists a boundary layer in which the original assumption about having slow variation does not hold. By first seeking the thickness of the layer, valid solutions can then be derived and then matched to the outer solutions in order to produce a uniform solution which holds for the entire axial domain. Once this is achieved, it is then used to derive trapped mode solutions. In this case, the theory used is that of two single turning points which are then combined to obtain the full solution. It is illustrated through consideration of examples and the dependence on ${\epsilon}$ is also shown through various plots. This problem will be considered for a symmetric and asymmetric duct and for differing duct parameters. Problems may arise when the two turning points lie close together and so we seek to improve on the method used by deriving a solution to trapped modes encompassing both turning points, which will be proposed together with some illustrations in order to justify its use and reliability. Next, the case of mode propagations on a thin elastic shell of varying radius conveying fluid is studied. The acoustic solutions of a straight shell in vacuo are first briefly reviewed and then built up by the addition of radius variation and the presence of a stationary fluid. The work presented first outlines the analysis for wave propagation in a slowly-varying thin elastic shell in vacuo. It is found that the shell and the fluid terms are coupled through the fluid pressure term, which is added to the equation governing the radial shell displacements since the pressure is assumed to affect radial motion only. Once the newly corrected equation for the radial shell displacements has been obtained, together with the axial and azimuthal displacements equations, this new system of governing equations is then separated into leading order ${\epsilon}^{0}$ and first order ${\epsilon}^{1}$ systems. In order to simplify the calculations, only the zeroth azimuthal order $m = 0$ will be studied here. With this simplification, a notable result is that the solutions of the torsional motion is decoupled from the axial and radial solutions. Once the dispersion equation is extracted from the leading order system, it can be seen that the axial and radial solutions are in fact coupled. The solution to the in vacuo with varying radius problem is first briefly presented and it is then followed by the solution to the fluid inclusion problem with varying radius, which makes up the main part of this section. The solution is studied for various frequencies and at various points along the shell. In addition, the axial and radial components of the first three modes are examined along with their amplitudes and energy distributions. Finally, mean flow is added and the same analysis is carried out, paying particular attention to the differences which arise in comparison to the stationary flow case.

Dans cette thèse, nous étudions les écoulements de fluides compressibles décrits par les équations de Navier-Stokes-Fourier dans les cas stationnaire et instationnaire et avec des conditions de bord assurant l’isolation thermique et mécanique du fluide. On commence par le cas stationnaire barotrope et des conditions de Navier à la frontière du domaine. La pression est donc de la forme p(%) = % où est appelé coefficient adiabatique et nous arrivons à montrer l’existence de solutions faibles pour > 1.On généralise ensuite ce résultat aux équations de Navier-Stokes-Fourier avec conduction de la chaleur et glissement (partiel ou total) au bord, toujours dans le cas stationnaire. On montre cette fois-ci l’existence de solutions faibles particulières appelées solutions entropiques variationnelles respectant l’inégalité d’entropie pour > 1 et l’existence de solutions faibles respectant le bilan de l’énergie totale au sens faible pour > 5/4. On travaille ensuite sur les écoulements instationnaires décrits par les équations de Navier-Stokes-Fourier sur une large variété de domaines non bornés, tout d’abord pour des conditions de bord d’adhérence puis pour des conditions de Navier à la frontière (ce qui restreintquelque peu la diversité des domaines non bornés admissibles). On arrive à montrer l’existence de solutions faibles particulières respectant l’inégalité d’entropie et une inégalité de dissipation remplaçant l’égalité de conservation d’énergie totale dans le volume qui n’a plus de sens dans les domaines non bornés. Par après, on met en place une inégalité dite d’entropie relative dont on montre qu’elle est respectée par certaines des solutions faibles exhibées auparavant. Ces solutions sont appelées solutions dissipatives. On parvient à prouver que pour chaque donnée initiale, il existe au moins une solution dissipative. Cette inégalité d’entropie relative nous permet de démontrer le principe d’unicité forte-faiblepour nos solutions dissipatives. Précisément, cela signifie qu’une solution dissipative et une solution forte issues des mêmes données initiales coïncident sur le temps maximal d’existence de la solution forte. La propriété d’unicité forte-faible donne un fondement à la notion de solution dissipative pour les domaines non bornés. / In this thesis, we study the Navier-Stokes-Fourier system describing the flow of compressible fluids both in the steady and unsteady case and we suppose that the fluid is thermally and mechanically isolated. We start with the case of a steady barotropic fluid and Navier boundary conditions. In this situation, the pressure law considered is of the form p(%) = %, where is called the adiabatic constant. We show the existence of weak solutions for > 1. We then extend this result to the complete Navier-Stokes-Fourier system with heat conductivity and slip or partially slip boundary conditions, once again in thesteady case. In this setup, we prove the existence of a specific type of weak solutions, called variationnal entropy solutions, which satisfy the entropy inequality for > 1 and the existence of weak solutions satisfying the conservation of total energy in its weak formulation for > 5/4. We then treat the unsteady flows described by the complete Navier-Stokes-Fourier system on a large class of unbouded domains, first with no-slip boundary conditions and then with the Navier boundary conditions which reduce the class of the admissible unbounded domains. We manage to prove the existence of a specific type of weak solutions verifying the entropy inequality and a dissipation inequality instead of the global conservation of total energy which is no more relevant in the unbounded domains. Afterwards, we establish a new inequality called relative entropy inequality and we show that it is satisfied by some of the weak solutions presented previously. These are called dissipative solutions. Next we show that for any given initial data there exists at least one dissipative solution. This observation allows us toperform the proof of the weak-strong uniqueness principle in the class of dissipative solutions. Precisely, it means that a dissipative solution and a classical one emanating from the same initial data coincide as long as the latter exists. The weak-strong uniqueness property justifies the concept of dissipative solutions in the situation of unbounded domains.

abstract: Rapid expansion of dense beds of fine, spherical particles subjected to rapid depressurization is studied in a vertical shock tube. As the particle bed is unloaded, a high-speed video camera captures the dramatic evolution of the particle bed structure. Pressure transducers are used to measure the dynamic pressure changes during the particle bed expansion process. Image processing, signal processing, and Particle Image Velocimetry techniques, are used to examine the relationships between particle size, initial bed height, bed expansion rate, and gas velocities. The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height. Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds. To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow. / Dissertation/Thesis / Rapid expansion of bed composed of [212, 297]micron particles. / Rapid expansion of bed composed of [44, 90]micron particles. / Rapid expansion of bed composed of [150, 212]micron particles. / Doctoral Dissertation Engineering 2019

Dynamic response of materials to high-speed and high-intensity loading conditions is important in several applications including high-speed flows with droplets, bubbles and particles, and hyper-velocity impact and penetration processes. In such high-pressure physics problems, simulations encounter challenges associated with the treatment of material interfaces, particularly when strong nonlinear waves like shock and detonation waves impinge upon them. To simulate such complicated interfacial dynamics problems, a fixed Cartesian grid approach in conjunction with levelset interface tracking is attractive. In this regard, a sharp interface Cartesian grid-based, Ghost Fluid Method (GFM) is developed for resolving embedded fluid, elasto-plastic solid and rigid (solid) objects in hyper-velocity impact and high-intensity shock loaded environment. The embedded boundaries are tracked and represented by virtue of the level set interface tracking technique. The evolving multi-material interface and the flow are coupled by meticulously enforcing the boundary conditions and jump relations at the interface. In addition, a tree-based Local Mesh Refinement scheme is employed to efficiently resolve the desired physics. The framework developed is generic and is applicable to interfaces separating a wide range of materials and for a broad spectrum of speeds of interaction (O(km/s)). The wide repertoire of problems solved in this work demonstrates the flexibility, stability and robustness of the method in accurately capturing the dynamics of the embedded interface. Shocks interacting with large ensembles of particles are also computed.

Cette thèse porte sur la simulation numérique d'un écoulement interne compressible, ins- tationnaire et turbulent à l'aide de la Simulation des Grandes Echelles (SGE). Dans cette approche, les grandes échelles énergétiques et instationnaires de l'écoulement sont calcu- lées, tandis que les petites échelles sont modélisées. L'objectif de ce travail, réalisé en collaboration avec le Consortium Industrie-Recherche dans les Turbomachines (CIRT) et le LEMFI, est d'analyser les potentialités de la SGE pour le calcul d'écoulements confinés en géométrie 3D en vue d'applications aux turboma- chines. Le cas d'étude retenu pour la validation de la méthodologie est l'écoulement dans une tuyère 3D transsonique pour lequel il existe de nombreux résultats statistiques (LEMFI) et expérimentaux (ONERA). La tuyère est caractérisée par une bosse en flèche sur la paroi basse génératrice d'effets 3D. L'écoulement, sans rotation, présente cependant les phénomènes physiques complexes d'interaction onde de choc/couche limite et de large dé- collement, comme ceux rencontrés au sein des turbomachines. La SGE est obtenue à partir d'un modèle de sous maille d'échelles mixtes développé au LIMSI-CNRS. La discrétisa- tion temporelle des équations est réalisée par un schéma de Runge-Kutta explicite d'ordre deux. Les flux Euler sont discrétisés par un schéma TVD d'Harten-Yee d'ordre deux tandis que les flux visqueux le sont par un schéma centré. La SGE a permis d'obtenir des informations instationnaires sur l'écoulement, et de mettre en évidence la formation et le lâcher de tourbillons dans la tuyère. La solution insta- tionnaire est différente de la solution stationnaire RANS obtenue avec une modélisation statistique classique et montre l'oscillation du choc et la déstabilisation du décollement au cours du temps. Les résultats de la SGE obtenus sur le nombre de Mach isentropique, les profils de la vitesse moyenne et les tensions de Reynolds sont discutés et comparés aux résultats expérimentaux (ONERA) et statistiques, obtenus avec le modèle de Launder- Shima (LEMFI). Le bilan de l'énergie cinétique turbulente (k) est analysé et comparé à celui donné par la modélisation statistique, dans le décollement. Les résultats de la SGE montrent des différences notables avec les résultats statistiques dans le c÷ur de la tuyère : les termes de fluctuations de pression importants mesurés dans le choc sont pris en compte par la SGE.

Accurate measurement of turbine engine compressor inlet total temperature is paramount for controlling engine speed and pressure ratio. Various methods exist for measuring compressor inlet total temperature on turbojet engines with hydromechanical control. One method involves the use of an ejector-diffuser system (eductor) to pull air from the engine inlet in order to measure the incoming total temperature. Analysis of historical test data has revealed that the inlet temperature measurement can be biased at certain flight conditions causing engine mis-scheduling and off-nominal engine operation. This bias is characterized primarily by adverse heat transfer effects and secondly by poor flow quality in the eductor tubing. Alternate eductor system configurations have been proposed to mitigate temperature bias. A one-dimensional engineering model of the eductor system was developed to facilitate the analysis of baseline and alternate eductor configurations. The model is calibrated with results from Computational Fluid Dynamics and validated with ground test data. The validated model is used to quantify the performance of several eductor configurations throughout the range of expected operating conditions and to quantify the amount of compressor inlet temperature measurement bias mitigation each configuration provides.

We consider the application of least-squares variational principles and the finite element method to the numerical solution of boundary value problems arising in the fields of solidand fluidmechanics.For manyof these problems least-squares principles offer many theoretical and computational advantages in the implementation of the corresponding finite element model that are not present in the traditional weak form Galerkin finite element model.Most notably, the use of least-squares principles leads to a variational unconstrained minimization problem where stability conditions such as inf-sup conditions (typically arising in mixed methods using weak form Galerkin finite element formulations) never arise. In addition, the least-squares based finite elementmodelalways yields a discrete system ofequations witha symmetric positive definite coeffcientmatrix.These attributes, amongst manyothers highlightedand detailed in this work, allow the developmentofrobust andeffcient finite elementmodels for problems of practical importance. The research documented herein encompasses least-squares based formulations for incompressible and compressible viscous fluid flow, the bending of thin and thick plates, and for the analysis of shear-deformable shell structures.

Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements compressibles comportant des interfaces. Ces interfaces peuvent séparer un fluide et un solide rigide, deux fluides de lois d'état différentes ou encore un fluide et un solide élastique. Dans un premier temps, nous avons élaboré une méthode de type frontières immergées afin d'imposer une condition de glissement au bord d'un obstacle rigide de manière précise. Nous avons ensuite étudié et validé un schéma de type interface non diffuse pour les écoulements multi-matériaux en vue d'appliquer la méthode de frontières immergées aux solides déformables.

Aquesta tesi doctoral proposa formulacions de Galerkin Discontinu (DG) d’alt ordre per la captura de shocks, obtenint alhora solucions altament precises per problemes de flux compressible. En les últimes dècades, la investigació en els mètodes de DG ha estat en constant creixement. L'èxit dels mètodes DG en problemes hiperbòlics ha conduit el seu desenvolupament en lleis de conservació no lineals i problemes de convecció dominant. Entre els avantatges dels mètodes DG, destaquen la seva estabilitat inherent i les propietats locals de conservació. D'altra banda, els mètodes DG estan especialment dissenyats per l’ús aproximacions d'ordre superior. De fet, en els últims anys s'ha demostrat que la resolució de problemes de convecció dominant ja no es restringeix només a elements d'ordre inferior. De fet, es necessiten models numèrics d'alta precisió per aconseguir prediccions altament fiables dins la dinàmica de fluids computacional (CFD). En aquest context es presenten i discuteixen dos tècniques de captura de shocks. En primer lloc, es presenta una tècnica novedosa i senzilla basada en la introducció d'una nova base de funcions de forma. Aquesta base té la capacitat de canviar a nivell local entre una interpolació contínua o discontínua, depenent de la suavitat de la funció que es vol aproximar. En presència de xocs, les discontinuïtats introduïdes dins l’element permeten incloure l'estabilització necessària gràcies a l’ús dels fluxos numèrics, i alhora exploten les propietats intrínsiques del mètodes DG. En conseqüència, es poden utilitzar malles grolleres amb elements d’ordre superior. Amb aquestes discretitzacions i, utilitzant el mètode proposats, els xocs queden continguts a l’interior de l’element i per tant, és possible evitar l’ús de tècniques de refinament adaptatiu de la malla, alhora que es manté la localitat i compacitat dels esquemes DG. En segon lloc, es proposa una tècnica clàssica i, aparentment simple: la introducció de la viscositat artificial. Primerament es realitza un estudi detallat per al cas unidimensional. S’obté una viscositat d’alta precisió que escala segons el valor hk amb 1 ≤ k ≤ p i essent h la mida de l’element. En conseqüència, s’obté un xoc amb amplitud del mateix ordre. Seguidament, l'estudi de la viscositat unidimensional obtenida s'extén al cas multidimensional per a malles triangulars. L'extensió es basa en la projecció de la viscositat unidimensional en unes determinades direccions espacials dins l’element. Es demostra de manera consistent que la viscositat introduïda és, com a molt, del mateix ordre que la resolució donada per la discretització espacial, és a dir, h/p. El mètode és especialment eficient per aproximacions de Galerkin discontinu d’alt ordre, per exemple p≥ 3. Les dues metodologies es validen mitjançant una àmplia selecció d’exemples numèrics. En alguns exemples, els mètodes proposats permeten una reducció en el nombre de graus de llibertat necessaris per capturar xocs acuradament de fins i tot un ordre de magnitud, en comparació amb mètodes estàndar de refinament adaptatiu amb aproximacions de baix ordre. / This thesis proposes shock-capturing methods for high-order Discontinuous Galerkin (DG) formulations providing highly accurate solutions for compressible flows. In the last decades, research in DG methods has been very active. The success of DG in hyperbolic problems has driven many studies for nonlinear conservation laws and convection-dominated problems. Among all the advantages of DG, their inherent stability and local conservation properties are relevant. Moreover, DG methods are naturally suited for high-order approximations. Actually, in recent years it has been shown that convection-dominated problems are no longer restricted to low-order elements. In fact, highly accurate numerical models for High-Fidelity predictions in CFD are necessary. Under this rationale, two shock-capturing techniques are presented and discussed. First, a novel and simple technique based on on the introduction of a new basis of shape functions is presented. It has the ability to change locally between a continuous or discontinuous interpolation depending on the smoothness of the approximated function. In the presence of shocks, the new discontinuities inside an element introduce the required stabilization thanks to the numerical fluxes, thus exploiting DG inherent properties. Large high-order elements can therefore be used and shocks are captured within a single element, avoiding adaptive mesh refinement and preserving the locality and compactness of the DG scheme. Second, a classical and, apparently simple, technique is advocated: the introduction of artificial viscosity. First, a one-dimensional study is perfomed. Viscosity of the order O(hk) with 1≤ k≤ p is obtained, hence inducing a shock width of the same order. Second, the study extends the accurate one-dimensional viscosity to triangular multidimensional meshes. The extension is based on the projection of the one-dimensional viscosity into some characteristic spatial directions within the elements. It is consistently shown that the introduced viscosity scales, at most, withthe DG resolutions length scales, h/p. The method is especially reliable for highorder DG approximations, say p≥3. A wide range of different numerical tests validate both methodologies. In some examples the proposed methods allow to reduce by an order of magnitude the number of degrees of freedom necessary to accurately capture the shocks, compared to standard low order h-adaptive approaches.

Simulating transient compressible flows involving shock waves presents challenges to the CFD practitioner in terms of the mesh quality required to resolve discontinuities and prevent smearing. This document discusses a novel two-dimensional Cartesian anisotropic mesh adaptation technique implemented for transient compressible flow. This technique, originally developed for laminar incompressible flow, is efficient because it refines and coarsens cells using criteria that consider the solution in each of the cardinal directions separately. In this document the method will be applied to compressible flow. The procedure shows promise in its ability to deliver good quality solutions while achieving computational savings. Transient shock wave diffraction over a backward step and shock reflection over a forward step are considered as test cases because they demonstrate that the quality of the solution can be maintained as the mesh is refined and coarsened in time. The data structure is explained in relation to the computational mesh, and the object-oriented design and implementation of the code is presented. Refinement and coarsening algorithms are outlined. Computational savings over uniform and isotropic mesh approaches are shown to be significant.