This master thesis treats the modelling of a common-rail direct fuel injection system where pressure generation is decoupled from the injection process. It has been shown that the fuel pressure plays a vital role for the general performance of the engine, affecting both emissions and efficiency, and it is carefully regulated to achieve optimal performance at different operating points. In an attempt to facilitate the development of the responsible control algorithms, a simulation framework has been requested. A model describing the complete work cycle of the high-pressure fuel system is developed and implemented in a Simulink environment. It is to a large extent based on the underlying physics and constructed in a modular manner, which allows for different engine configurations to be simulated. The modelled pressure signal is compared to experimental data at different operating points with promising results in capturing the transient behaviour from a low-level perspective. Additionally, it manages to replicate some of the pressure oscillations which has been observed in the real system and it shows good response to changes in the input signals. However, there are some areas which are subject to improvement since capturing the static pressure levels over longer drive cycles has proved to be a difficult task. Overall, the developed model serves as a starting point for future development and validation of control algorithms.

<div>In order to better understand and quantify the effect of instabilities in systems utilizing flow boiling heat transfer, the present study explores dynamic results for pressure drop, mass velocity, thermodynamic equilibrium quality, and heated wall temperature to ascertain and analyze the dominant modes in which they oscillate. Flow boiling experiments are conducted for a range of mass velocities with both subcooled and saturated inlet conditions in vertical upflow, vertical downflow, and horizontal flow orientations. High frequency pressure measurements are used to investigate the influence of individual flow loop components (flow boiling module, pump, pre-heater, condenser, etc.) on dynamic behavior of the fluid, with fast Fourier transforms of the same used to provide critical frequency domain information. Conclusions from this analysis are used to isolate instabilities present within the system due to physical interplay between thermodynamic and hydrodynamic effects. Parametric analysis is undertaken to better understand the conditions under which these instabilities form and their impact on system performance. Several prior stability maps are presented, with new stability maps provided to better address contextual trends discovered in the present study.</div><div>Further, this study utilizes experimental results for vertical upflow boiling of FC-72 in a rectangular channel with finite inlet quality to investigate Density Wave Oscillations (DWOs) and assess their potential impact on design of two-phase systems for future space missions. High-speed flow visualization image sequences are presented and used to directly relate the cyclical passage of High and Low Density Fronts (HDFs and LDFs) to dominant low-frequency oscillations present in transient pressure signals commonly attributed to DWOs. A methodology is presented to determine frequency and amplitude of DWO induced pressure oscillations, which are then plotted for a wide range of relevant operating conditions. Mass velocity (flow inertia) is seen to be the dominant parameter influencing frequency and amplitude of DWOs. Amplitude of pressure oscillations is at most 7% of the time-averaged pressure level for current operating conditions, meaning there is little risk to space missions. Reconstruction of experimental pressure signals using a waveform defined by frequency and amplitude of DWO induced pressure fluctuations is seen to have only moderate agreement with the original signal due to the oversimplifications of treating DWO induced fluctuations as perfectly sinusoidal in nature, assuming they occur at a constant frequency value, and neglecting other transient flow features. This approach is nonetheless determined to have potential value for use as a boundary condition to introduce DWOs in two-phase flow simulations should a model be capable of accurately predicting frequency and amplitude of oscillation.</div><div>Additionally, this study presents a new mechanistic model for Density Wave Oscillations (DWOs) in vertical upflow boiling using conclusions drawn from analysis of flow visualization images and transient experimental results as a basis from which to begin modeling. Counter to many prior studies attributing DWOs to feedback effects between flow rate, pressure drop, and flow enthalpy causing oscillations in position of the bulk boiling boundary, the present instability mode stems primarily from body force acting on liquid and vapor phases in a separated flow regime leading to liquid accumulation in the near-inlet region of the test section, which eventually departs and moves along the channel, acting to re-wet liquid film along the channel walls and re-establish annular, co-current flow. This process was modeled by dividing the test section into three distinct control volumes and solving transient conservation equations for each, yielding predictions of frequencies at which this process occurs as well as amplitude of associated pressure oscillations. Values for these parameters were validated against an experimental database of 236 FC-72 points and show the model provides good predictive accuracy and capably captures the influence of parametric changes to operating conditions.</div><div>Also, this study shows analysis of pressure signals in condensing systems reveal the presence of relevant oscillatory phenomena during flow condensation as well, which may impact performance in applications concerned with precise system control. Towards this end, the present study presents results for oscillatory behavior observed in pressure measurements during flow condensation of FC-72 in a smooth circular tube in vertical upflow, vertical downflow, and horizontal flow orientations. Dynamic behavior observed within the test section is determined to be independent of other components within the flow loop, allowing it to be isolated and interpreted as resulting from physical aspects of two-phase flow with condensation. The presence of a peak oscillatory mode (one of significantly larger amplitude than any others present) is seen for 72% of</div><div>vertical upflow test cases, 61% of vertical downflow, and 54% of horizontal flow. Relative intensities of this peak oscillatory mode are evaluated through calculation of Q Factor for the corresponding frequency response peak. Frequency and amplitude of peak oscillatory modes are also evaluated. Overall, vertical upflow is seen to exhibit the most significant oscillatory behavior, although in its maximum case amplitude is only seen to be 7.9% of time-averaged module inlet pressure, indicating there is little safety risk posed by oscillations under current operating conditions. Flow visualization image sequences for each orientation are also presented and used to draw parallels between physical characteristics of condensate film behavior under different operating conditions and trends in oscillatory behavior detected in pressure signals</div><div>Further, the present work outlines a new methodology utilizing temperature and pressure measurements to identify condensation flow regimes. For vertical upflow condensation, amplitude of dynamic temperature and pressure oscillations are shown to clearly indicate transition from counter-current flow regimes (i.e., falling film, oscillating film, flooding) to annular, co-current flow (climbing film flow regime). In horizontal flow condensation, standard deviation between multiple thermocouple measurements distributed around the tube circumference was calculated at all axial (stream-wise) measurement locations. High values of standard deviation are present for stratified flow (stratified flow, wavy-stratified, plug flow), while axisymmetric flow regimes (i.e., slug flow, annular flow) yield significantly lower values. Successful development of this technique represents a valuable contribution to literature as it allows condensation flow regime to be identified without the often-costly restriction of designing a test section to allow optical access. Identified flow regimes in both vertical upflow and horizontal flow orientations are compared to regime maps commonly found in the literature in pursuit of optimum performing maps.</div><div>Finally, the present study aims to better analyze the influence of body force on flow condensation heat transfer by conducting tests at multiple orientations in Earth’s gravity. Dielectric FC-72 is condensed in a smooth stainless-steel tube with 7.12 mm diameter and 574.55 mm condensing length by counterflow of cooling water across the outer surface of the tube. Test conditions span FC-72 mass velocities of 50.3 – 360.3 kg/m2s, test section inlet pressures of 127.0 – 132.1 kPa, and test section inlet thermodynamic equilibrium qualities of 0.13 – 1.15. A subset of data gathered corresponding to axisymmetric, annular condensation heat transfer is identified and a detailed methodology for data reduction to calculate heat transfer coefficient presented. Uncertainty analysis is also presented and indicates channel average heat transfer coefficients are calculated within ±3.6% to ±26.7% (depending on operating conditions). Analysis of parametric trends for condensation heat transfer reveals the dominant influence of mass velocity (flow inertia), secondary influence of vapor mass fraction (thermodynamic equilibrium quality), and strong dependence on orientation (body force) at low mass velocities. At higher mass velocities results for all orientations investigated begin to converge, indicating body force independent annular condensation heat transfer is achieved. Separated Flow Model predictions of vertical downflow condensation heat transfer provide reasonable agreement with experimental results, evidence by a Mean Absolute Error (MAE) of 31.2%. Evaluation of condensation heat transfer correlations for horizontal flow reveal most correlations struggle for cases with high liquid content. Specific correlations are identified for superior accuracy in predicting the measured data.</div>

Les moteurs à propergol solide de certains lanceurs spatiaux peuvent présenter des oscillations de poussée provoquées par des oscillations de la pression interne du moteur. Il est désormais connu qu'une instabilité hydrodynamique propre à l'écoulement et mise en évidence par une approche de stabilité linéaire, appelée Vortex Shedding Parietal, est le phénomène à l'origine de ces oscillations. Cependant plusieurs questions subsistent quant à l'apparition de ces oscillations. En particulier, on constate que ces oscillations ne sont observées que dans la seconde moitié du tir. On suppose que la transition laminaire-turbulent de l'écoulement joue un rôle dans ce phénomène. Par ailleurs, l'utilisation possible de la géométrie étoilée pour les blocs de propergol de géométrie étoilée oblige à reprendre l'étude de stabilité linéaire pour l'étude de ces oscillations. Cette thèse porte ainsi à la fois sur l'étude de la stabilité linéaire de l'écoulement induit par une géométrie étoilée et sur la transition laminaire-turbulent d'une configuration circulaire. Une solution analytique de l'écoulement au sein d'un conduit à paroi débitante de géométrie étoilée est d'abord obtenue. Le profil de vitesse associé peut présenter des points d'inflexion en fonction de la déformation radiale ou du nombre de Reynolds de l'écoulement. Ensuite les approches de stabilité linéaire locale et biglobale sont appliquées à cette solution. Par rapport au cas circulaire, cette étude de stabilité linéaire met en évidence la possible existence d'un ou plusieurs modes amplifiés. Enfin, une simulation numérique LES de la transition laminaire-turbulent de la configuration du montage VALDO est effectuée. / Solid rocket motors may exhibit thrust oscillations induced by pressure oscillations inside the engine. In the scientific literature, the so-called vortex shedding parietal, a hydrodynamic instability specific to the internal flow, has been proven to be the triggering mechanism through linear stability analyses. However, some questions still remain. These oscillations are observed during the second half of Ariane 5 launches only, the laminar-turbulent transition of the intern flow being expected to play a significant role in that behavior. Additionally, non-circular grain shapes used in some engines may impact the stability study. The present PhD work extends the base flow linear stability studies to star-shaped geometries and investigate transition mechanisms in a circular configuration. Firstly, an analytical expression of the base flow in a duct with a star-shaped cross-section is determined. Wall-normal injection is assumed, which is representative of solid rocket engines. The velocity profile may exhibit an inflection point depending on both the Reynolds number and wall radial deformation. Secondly, both local and biglobal stability analyses are applied to this flow. In contrast to the circular case, the flow may exhibit one or several unstable modes. Lastly, a large eddy simulation of the VALDO experimental test bench is performed. The influence of the acoustic impedance of the outflow boundary condition on the amplitude of pressure oscillations is highlighted.

Acoustic instabilities in a combustion chamber arise due to the coupling of acoustic pressure with in-phase heat-release, and are characterized by large amplitude oscillations of one or more natural modes of combustor. Even though an array of studies, both theoretical and experimental, has been conducted by a number of authors in this field to extract the flame response, most of these are based on kinematic flame models. In this dissertation, an experimental study of a subsonic flame's intrinsic response to acoustic pressure perturbations is performed for the case of a tube closed at one end and the other end opened to the atmospheric conditions. Pressure fluctuations inside the tube are measured for hot and cold side flows, and their varying trend is explained. The frequencies obtained from Fourier transform analysis exhibit a strong dependence with the distance between the stabilized flame position and open end of the tube. For different values of flame position (xf ), the values of growth constant 's' are calculated from the pressure versus time data readings procured from acoustic pressure transducer and dominant frequencies are analyzed from windowed FFT of the same. The expression for obtaining response function from the measured pressure fluctuations has been derived from the 1-D linearized conservation equations. The undamped response function plot is obtained by adding the decay rates at different frequencies inside the tube to the corresponding growth rates. Finally, the effect of blockage of pre-mixed flow on the growth rates inside the tube and consequently, the flame response values, is studied by repeating the experiment with different types of flame holders. A large number of theoretical flame-response models have been developed in modern literature, and some of these models are compared with the experimentally obtained response. Suggestions are also cited in this study so as to account for the observed deviations in trends. This includes a revisit of the intrinsic flame model by incorporating the effect of flame-area perturbations, with the aid of analyzed steady flame images.

Les moteurs de fusées à ergols solides (SRMs) sont sensibles aux instabilités hydrodynamiques qui peuvent déclencher des oscillations auto-entretenues de pression de grandes amplitudes lorsqu’elles se couplent à l’un des modes acoustiques du système. Le moteur de ces instabilités est la formation de structures tourbillonnaires cohérentes synchronisées par des ondes acoustiques longitudinales. Pour certaines conditions de fonctionnement, les ondes acoustiques générées par l’interaction de ces tourbillons avec la tuyère amorcée du moteur renforcent l’oscillation acoustique. L’objectif des travaux menés dans cette thèse est de déterminer l’amplitude et la fréquence des oscillations de pression au cycle limite des instabilités. Celui-ci est atteint par saturation non linéaire des sources, qui est la conséquence de la formation de grosses structures cohérentes. Dans ce cas l’interaction tourbillon tuyère devient insensible à l’amplitude de l’onde du mode acoustique établi dans le foyer. Dans ces conditions, on peut se concentrer sur l’interaction d’un tourbillon avec la tuyère dans le mécanisme de production sonore. En considérant un écoulement incompressible et l’absence de frottement, un premier modèle analytique est développé permettant de déterminer la production sonore d’un tourbillon ingéré par une tuyère bidimensionnelle plane, lorsque le tourbillon est traité comme une ligne vorticité. Des expériences précédentes indiquent que le volume de la cavité autour de l’entrée d’une tuyère intégrée a une grande influence sur l’amplitude des oscillations de pression dans les grands SRMs. On montre que ceci est dû au champ de vitesse acoustique induit par la compressibilité du gaz dans la cavité qui produit une fluctuation de vitesse transverse à la trajectoire du tourbillon. Une seconde alternative au modèle analytique incompressible est développée en considérant toujours l’absence de frottement, mais un modèle compressible de l’interaction tourbillon-tuyère. Celui-ci repose sur un code aéroacoustique pour les écoulements internes basé sur les équations d’Euler (EIA) qui est utilisé ici pour la simulation de l’interaction tourbillon-tuyère. Une étude systématique de cette interaction a été menée pour une tuyère amorcée. Les résultats ont permis de proposer un modèle de sources localisées pour des ondes planes basé sur une analyse théorique des lois d’échelles de ces phénomènes. Les simulations de ces interactions tourbillons-tuyères ont été réalisées pour différents types de tuyères. En employant un bilan énergétique, un modèle avec un seul paramètre de contrôle est formulé, qui permet de reproduire qualitativement le comportement du cycle limite d’oscillations de pression observées dans des expériences réalisées avec des gaz froids décrites dans la littérature. Finalement le modèle Euler est utilisé pour comparer la production de son par interaction tourbillon-tuyère avec celle due à l’ingestion d’une onde d’entropie, appelée aussi tache d’entropie. Contrairement au cas des tourbillons, le bruit produit par ingestion de taches d’entropie n’est pas sensible au volume de la cavité d’une tuyère intégrée. Ces résultats indiquent que le bruit produit par les tourbillons est dominant dans le cas des SRMs étudiés. L’ensemble de ces travaux permet d’améliorer la compréhension des phénomènes d’interaction entre des non-homogénéités de l’écoulement et la tuyère. Elle permet surtout de déterminer quels sont les facteurs de l’écoulement et les éléments géométriques importants qui pilotent le niveau sonore produit par ces interactions. Les modèles développés dans ces travaux, avec divers degrés d’approximation et de complexité permettent d’enrichir la gamme des outils de conception des SRMs. / Solid Rocket Motors (SRMs) can display self-sustained acoustic oscillations driven by coupling between hydrodynamic instabilities of the internal flow and longitudinal acoustic standing waves. The hydrodynamic instabilities are triggered by the acoustic standing wave and results in the formation of coherent vortical structures. For nominal ranges of flow conditions the sound waves generated by the interaction between these vortices and the choked nozzle at the end of the combustion chamber reinforces the acoustic oscillation. Most available literature on this subject focuses on the threshold of instability using a linear model. The focus of this work is on the prediction of the limit-cycle amplitude. The limit-cycle is reached due to nonlinear saturation of the source, as a consequence of the formation of large coherent vortical structures. In this case the vortex-nozzle interaction becomes insensitive to the amplitude of the acoustic standing wave. Hence, one can focus on the sound generation of a vortex with the nozzle. Sound production can be predicted from an analytical two-dimensional planar incompressible frictionless model using the so-called Vortex Sound Theory. In this model the vorticity is assumed to be concentrated in a line vortex. Experiments indicate that the volume of cavities around so-called “integrated nozzles” have a large influence on the pulsation amplitude for large SRMs. This is due to the acoustical field normal to the vortex trajectory, induced by the compressibility of the gas in this cavity. As an alternative to the incompressible analytical model a compressible frictionless model with an internal Euler Aeroacoustic (EIA) flow solver is used for simulations of vortex-nozzle interaction. A dedicated numerical simulation study focusing on elementary processes such as vortex-nozzle and entropy spot-nozzle interaction allows a systematic variation of relevant parameters and yields insight which would be difficult by means of limit cycle studies of the full engine. A systematic study of the vortex-nozzle interaction in the case of a choked nozzle has been undertaken. The results are summarized by using a lumped element model for plane wave propagation, which is based on theoretical scaling laws. From EIA simulations it appears that sound due to vortex-nozzle interaction is mainly generated during the approach phase and that for the relevant parameter range there is no impingement of the vortex on the nozzle wall as has been suggested in the literature. Using an energy balance approach, a single fit-parameter model is formulated which qualitatively predicts limit-cycle observations in cold gas-scale experiments reported in the literature. Finally the Euler model is used to compare the sound production by vortex-nozzle interaction with that due to the ingestion of an entropy non-uniformity also called entropy spot. In addition to insight, this study provides a systematic procedure to develop a lumped element model for the sound source due to non-homogeneous flow-nozzle interactions in SRMs. Such lumped models based on experimental data or a limited number of flow simulations can be used to ease the design of SRMs.

In the first part of this dissertation, the phenomenon of self-sustained pressure os-cillations due to the flow past a circular, axisymmetric cavity, associated with inline gate valves, was investigated. In many engineering applications, such as flows through open gate valves, there exists potential for coupling between the vortex shedding from the up-stream edge of the cavity and a diametral mode of the acoustic pressure fluctuations. The effects of the internal pipe geometry immediately upstream and downstream of the shal-low cavity on the characteristics of partially trapped diametral acoustic modes were in-vestigated numerically and experimentally on a scaled model of a gate valve mounted in a pipeline that contained convergence-divergence sections in the vicinity of the valve. The resonant response of the system corresponded to the second acoustic diametral mode of the cavity. Excitation of the dominant acoustic mode was accompanied by pressure oscillations, and, in addition to that, as the angle of the converging-diverging section of the main pipeline in the vicinity of the cavity increased, the trapped behavior of the acoustic diametral modes diminished, and additional antinodes of the acoustic pressure wave were observed in the main pipeline. In addition to that, the effect of shallow chamfers, introduced at the upstream and/or downstream cavity edges, was investigated in the experimental system that con-tained a deep, circular, axisymmetric cavity. Through the measurements of unsteady pressure and associated acoustic mode shapes, which were calculated numerically for several representative cases of the internal cavity geometry, it was possible to identify the configuration that corresponded to the most efficient noise suppression. This arrangement also allowed calculation of the azimuthal orientation of the acoustic modes, which were classified as stationary, partially spinning or spinning. Introduction of shallow chamfers at the upstream and the downstream edges of the cavity resulted in changes of azimuthal orientation and spinning behaviour of the acoustic modes. In addition, introduction of splitter plates in the cavity led to pronounced change in the spatial orientation and the spinning behaviour of the acoustic modes. The short splitter plates changed the behaviour of the dominant acoustic modes from partially spinning to stationary, while the long split-ter plates enforced the stationary behaviour across all resonant acoustic modes. Finally, the evolution of fully turbulent, acoustically coupled shear layers that form across deep, axisymmetric cavities and the effects of geometric modifications of the cavity edges on the separated flow structure were investigated using digital particle image velocimetry (PIV). Instantaneous, time- and phase-averaged patterns of vorticity pro-vided insight into the flow physics during flow tone generation and noise suppression by the geometric modifications. In particular, the first mode of the shear layer oscillations was significantly affected by shallow chamfers located at the upstream and, to a lesser degree, the downstream edges of the cavity. In the second part of the dissertation, the performance of aortic heart valve pros-thesis was assessed in geometries of the aortic root associated with certain types of valve diseases, such as aortic valve stenosis and aortic valve insufficiency. The control case that corresponds to the aortic root of a patient without valve disease was used as a reference. By varying the aortic root geometry, it was possible to investigate corresponding changes in the levels of Reynolds shear stress and establish the possibility of platelet activation and, as a result of that, the formation of blood clots. / Graduate / 0541 / 0546 / 0548 / 0986 / alexbn024@gmail.com

In the first part of this dissertation, the phenomenon of self-sustained pressure os-cillations due to the flow past a circular, axisymmetric cavity, associated with inline gate valves, was investigated. In many engineering applications, such as flows through open gate valves, there exists potential for coupling between the vortex shedding from the up-stream edge of the cavity and a diametral mode of the acoustic pressure fluctuations. The effects of the internal pipe geometry immediately upstream and downstream of the shal-low cavity on the characteristics of partially trapped diametral acoustic modes were in-vestigated numerically and experimentally on a scaled model of a gate valve mounted in a pipeline that contained convergence-divergence sections in the vicinity of the valve. The resonant response of the system corresponded to the second acoustic diametral mode of the cavity. Excitation of the dominant acoustic mode was accompanied by pressure oscillations, and, in addition to that, as the angle of the converging-diverging section of the main pipeline in the vicinity of the cavity increased, the trapped behavior of the acoustic diametral modes diminished, and additional antinodes of the acoustic pressure wave were observed in the main pipeline. In addition to that, the effect of shallow chamfers, introduced at the upstream and/or downstream cavity edges, was investigated in the experimental system that con-tained a deep, circular, axisymmetric cavity. Through the measurements of unsteady pressure and associated acoustic mode shapes, which were calculated numerically for several representative cases of the internal cavity geometry, it was possible to identify the configuration that corresponded to the most efficient noise suppression. This arrangement also allowed calculation of the azimuthal orientation of the acoustic modes, which were classified as stationary, partially spinning or spinning. Introduction of shallow chamfers at the upstream and the downstream edges of the cavity resulted in changes of azimuthal orientation and spinning behaviour of the acoustic modes. In addition, introduction of splitter plates in the cavity led to pronounced change in the spatial orientation and the spinning behaviour of the acoustic modes. The short splitter plates changed the behaviour of the dominant acoustic modes from partially spinning to stationary, while the long split-ter plates enforced the stationary behaviour across all resonant acoustic modes. Finally, the evolution of fully turbulent, acoustically coupled shear layers that form across deep, axisymmetric cavities and the effects of geometric modifications of the cavity edges on the separated flow structure were investigated using digital particle image velocimetry (PIV). Instantaneous, time- and phase-averaged patterns of vorticity pro-vided insight into the flow physics during flow tone generation and noise suppression by the geometric modifications. In particular, the first mode of the shear layer oscillations was significantly affected by shallow chamfers located at the upstream and, to a lesser degree, the downstream edges of the cavity. In the second part of the dissertation, the performance of aortic heart valve pros-thesis was assessed in geometries of the aortic root associated with certain types of valve diseases, such as aortic valve stenosis and aortic valve insufficiency. The control case that corresponds to the aortic root of a patient without valve disease was used as a reference. By varying the aortic root geometry, it was possible to investigate corresponding changes in the levels of Reynolds shear stress and establish the possibility of platelet activation and, as a result of that, the formation of blood clots. / Graduate / 0541 / 0546 / 0548 / 0986 / alexbn024@gmail.com

Cette thèse vise à modéliser les instabilités hydrodynamiques générant des détachements tourbillonnaires pariétaux (ou VSP) responsables des Oscillations De Pression dans les moteurs à propergol solide longs et segmentés par interaction avec l’acoustique du moteur. Ces instabilités sont modélisées en tant que modes de stabilité linéaire globaux de l’écoulement d’un conduit à parois débitantes. En supposant que les structures pariétales émergent d’une perturbation de l’écoulement de base, des modes discrets et indépendants du maillage utilisé sont calculés. Dans ce but, une discrétisation par collocation spectrale multi-domaine est implémentée dans un solveur parallèle afin de s’affranchir de la croissance polynomiale des fonctions propres et de la présence de couches limites. Les valeurs propres ainsi calculées dépendent explicitement des frontières du domaine, à savoir la position de la perturbation et celle de la sortie, et sont ensuite validées par simulation numérique directe. On montre alors qu’elles permettent bien de décrire la réponse à une perturbation initiale de l’écoulement modifié par une rupture de débit pariétale. Ensuite, la simulation d’une réponse forcée par l’acoustique se fait sous forme de structures tourbillonnaires dont les fréquences discrètes sont en accord avec celles des modes de stabilité. Ces structures sont réfléchies en ondes de pression de même fréquences remontant l’écoulement. Finalement, la simulation numérique et la théorie de la stabilité permettent de montrer que le VSP, dont la réponse est linéaire vis-à-vis d’un forçage compressible comme l’acoustique, est le phénomène moteur des Oscillations De Pression. / The current work deals with the modeling of the hydrodynamic instabilities that play a major role in the triggering of the Pressure Oscillations occurring in large segmented solid rocket motors. These instabilities are responsible for the emergence of Parietal Vortex Shedding (PVS) and they interact with the boosters acoustics. They are first modeled as eigenmodes of the internal steady flowfield of a cylindrical duct with sidewall injection within the global linear stability theory framework. Assuming that the related parietal structures emerge from a baseflow disturbance, discrete meshindependant eigenmodes are computed. In this purpose, a multi-domain spectral collocation technique is implemented in a parallel solver to tackle numerical issues such as the eigenfunctions polynomial axial amplification and the existence of boundary layers. The resulting eigenvalues explicitly depend on the location of the boundaries, namely those of the baseflow disturbance and the duct exit, and are then validated by performing Direct Numerical Simulations. First, they successfully describe flow response to an initial disturbance with sidewall velocity injection break. Then, the simulated forced response to acoustics consists in vortical structures wihich discrete frequencies that are in good agreement with those of the eigenmodes. These structures are reflected into upstream pressure waves with identical frequencies. Finally, the PVS, which response to a compressible forcing such as the acoustic one is linear, is understood as the driving phenomenon of the Pressure Oscillations thanks to both numerical simulation and stability theory.