Experimental observation of supersonic non-ideal compressible-fluid flows

The branch of fluid mechanics devoted to the study of compressible flows whose behaviour deviates from the one predicted by the ideal-gas model is termed Non-Ideal Compressible-Fluid Dynamics (NICFD). Departure from ideal behaviour is typically observed in fluids made of complex molecules, operating in thermodynamic conditions close to the liquid-vapour saturation curve and critical point. Differently from the ideal-gas case, non-ideal flow dynamics is strongly influenced by the process conditions, and peculiar flow patterns such as rarefaction shock waves, shock waves with either upstream or downstream sonic states, and split shocks are physically admissible. Complex thermodynamic models were devised and applied to the theoretical analysis and numerical simulation of non-ideal flows. However, only few experimental data in the NICFD regime were available up to date. In the present thesis, experimental benchmark tests are devised and measurements are performed to provide assessment to NICFD theory, thermodynamic models and numerical simulation tools. Experiments are carried out in the Test Rig for Organic VApors (TROVA), the blow-down wind tunnel of the Laboratory of Compressible-fluid dynamics for Renewable Energy Applications (CREALab) of Politecnico di Milano, while numerical simulations are performed using the open-source NICFD solver SU2. Two main flow configurations are considered, namely adapted flows in planar converging-diverging nozzles and uniform supersonic flows with oblique shock waves and expansion fans. An experimental campaign is carried out to study the NICFD flow of siloxane vapour MDM (octamethyltrisiloxane, C8H24O2Si3) in a planar converging-diverging nozzle. The nozzle is designed using the method of characteristics, complemented with state-of-the-art Equations of State (EoS) dealing with non-ideal flow behaviour. Operation of the nozzle is characterised through pressure and temperature measurements, schlieren visualizations, and numerical simulations. Numerical simulations and theoretical predictions based on the isentropic hypothesis for nozzle expansions agree well with experimental data. Consistently with NICFD theory, reservoir conditions are found to significantly influence both the pressure ratio and the Mach number distribution along the nozzle axis. Numerical simulations and the NICFD theory are then extensively applied to the study of supersonic expansions in converging-diverging nozzles, both in the ideal and NICF regimes. Two different topics related to the design and operation of converging-diverging nozzles are addressed. First, a design solution aimed at fixing the location of the minimum-area section of the nozzle is discussed. The geometry of planar converging-diverging nozzles operating with air and MDM is modified by the introduction of a small recessed step at the geometrical throat. The recessed step triggers boundary layer separation and fixes the location of the minimum-area section of the nozzle. Experiments and numerical simulations are applied to the analysis of the adapted flow. A complex perturbation wave pattern originates at the step location, and propagates up to the exhaust section of the nozzle through multiple reflections at the nozzle walls and interactions at the nozzle symmetry axis. Pressure measurements performed along the nozzle axis show that the perturbation introduced by the step influences the flow only locally. For the nozzle operating with MDM in the NICFD regime, the configuration of the perturbation wave pattern in the throat region is found to depend on reservoir conditions. Second, a key parameter for the design and performance analysis of supersonic nozzles is investigated, namely the discharged massflow rate. A closed-form expression for this parameter is derived analytically. The massflow rate is found to be dependent on the radius of curvature of the nozzle profile at the throat section, on the molecular complexity of the working fluid, and on process conditions. A numerical verification of the dependence of the discharged massflow rate on relevant design parameters is performed, both in the ideal and non-ideal flow regimes. The dependence of the massflow rate on the fluid properties and on process conditions is confirmed. Moreover, the discharged massflow rate is found to be affected by the overall shape of the converging portion, and not only by the local curvature of the nozzle profile at the throat section. This is true even in the case of smooth convergents. A second experimental campaign is carried out to characterize oblique shock waves and expansion fans occurring in the NICFD regime. A diamond-shaped airfoil with semi-aperture 7.5° at leading edge and 10° at the trailing edge is placed in a uniform supersonic stream of MDM, which in test conditions is a single-phase vapour lying in the NICFD region. Two oblique shocks at the airfoil leading-edge and centered expansion fans at the shoulder are observed. Oblique shock waves and expansion fans are studied at varying upstream stagnation conditions, for six deviation angles in the range 6.5°-10° obtained by changing the attitude of the model with respect to the wind tunnel axis. Experimental results are assessed against the inviscid shock-expansion theory for two-dimensional steady flows, complemented with state-of-the-art EoS. All the geometrical configurations and operating conditions explored in wind-tunnel tests are reproduced numerically. In the inviscid core of the flow, where measurements are performed, numerical results agree fairly well with experimental data. A significant dependence of the pre-shock, post-shock and post-expansion states on the corresponding upstream flow conditions is found, and the speed of sound is seen to decrease across shock waves and increase upon expansion fans, consistently with the NICFD theory applied to the tested thermodynamic conditions.

La branca della fluidodinamica che si occupa dello studio di correnti comprimibili il cui comportamento si discosta da quanto predetto dal modello di gas ideale è denominata Non-Ideal Compressible-Fluid Dynamics (NICFD). Deviazioni significative dal comportamento ideale si osservano tipicamente in flussi di sostanze costituite da molecole complesse, in condizioni termodinamiche vicine alla curva di saturazione liquido-vapore e al punto critico. A differenza del caso ideale, la dinamica di un flusso non ideale è influenzata sensibilmente dalle condizioni di processo. Complessi modelli termodinamici sono stati applicati all'analisi teorica e alla simulazione numerica di tali correnti. Tuttavia, solo pochi dati sperimentali relativi a questo tipo di correnti sono attualmente disponibili. In questo lavoro di tesi vengono presentati esperimenti su correnti non-ideali di riferimento, effettuati per verificare gli strumenti teorici e numerici e i modelli termodinamici attualmente disponibili. Gli esperimenti sono stati effettuati presso la galleria del vento supersonica TROVA (Test Rig for Organic VApors) del CREALab (Laboratory of Compressible-fluid dynamics for Renewable Energy Applications), Politecnico di Milano. Simulazioni numeriche di correnti non-ideali a supporto degli esperimenti sono state eseguite con il solutore open-source SU2. Nella presente tesi vengono analizzate due principali tipologie di correnti, ossia il flusso adattato in ugelli piani convergenti-divergenti e flussi supersonici uniformi con onde d'urto oblique e ventagli di espansione. Una prima campagna sperimentale viene eseguita per studiare il flusso non-ideale del fluido silossano MDM (octamethyltrisiloxane, C8H24O2Si3) in un ugello piano convergente-divergente. L'ugello è progettato utilizzando il metodo delle caratteristiche e complesse equazioni di stato in grado di trattare il comportamento non-ideale della corrente. Il funzionamento dell'ugello viene caratterizzato tramite misure di pressione e temperatura, visualizzazioni schlieren e simulazioni numeriche. Le simulazioni numeriche e il modello di espansione isentropica per il flusso nell'ugello mostrano un buon accordo con i risultati sperimentali. Coerentemente con la teoria NICFD, le condizioni di serbatoio influenzano significativamente l'andamento in asse ugello sia del rapporto di pressione statica/totale, sia del numero di Mach. Simulazioni numeriche e la teoria NICFD vengono in seguito applicate allo studio delle espansioni supersoniche in ugelli convergenti-divergenti, sia in regime ideale, sia in regime non-ideale. In particolare, vengono approfonditi due aspetti di particolare interesse per il progetto e il funzionamento di ugelli convergenti-divergenti. Per prima cosa viene discussa una soluzione di progetto volta a fissare la posizione della sezione minima dell'ugello. Le geometrie di ugelli piani convergenti-divergenti per aria e MDM vengono modificate con l'introduzione di un piccolo scalino in corrispondenza della gola geometrica. Lo scalino induce la separazione dello strato limite e fissa la posizione della sezione di area minima dell'ugello. La corrente adattata all'interno di tali ugelli viene studiata mediante esperimenti e simulazioni numeriche. In corrispondenza dello scalino sulla sezione di gola ha origine un complesso sistema di onde. Tale sistema di onde si propaga fino alla sezione di efflusso dell'ugello mediante una serie di riflessioni sulle pareti dell'ugello e interazioni lungo l'asse di simmetria. Le misure di pressione effettuate lungo l'asse dell'ugello mostrano che la perturbazione introdotta dallo scalino influenza la corrente solo localmente. Per l'ugello funzionante con MDM in regime non ideale, la configurazione del sistema di onde nella regione di gola dipende dalle condizioni di serbatoio. In seguito viene studiato un parametro chiave per il progetto e l'analisi delle prestazioni di un ugello supersonico, ossia la portata massica scaricata, di cui viene ricavata un'espressione analitica in forma chiusa. La portata massica dipende dal raggio di curvatura dell'ugello in corrispondenza della sezione di gola, dalla complessità molecolare del fluido di lavoro e dalle condizioni di processo. La dipendenza della portata massica scaricata da vari parametri di progetto viene verificata numericamente, sia in regime ideale, sia in regime non ideale, confermando la dipendenza dalle proprietà del fluido e dalle condizioni di processo. Inoltre, viene mostrato come la portata massica scaricata dipenda dalla forma del convergente, e non soltanto dalla curvatura locale del profilo dell'ugello in corrispondenza della sezione di gola. Una seconda campagna sperimentale viene eseguita per caratterizzare onde d'urto oblique e ventagli di espansione in regime non ideale. Un profilo a losanga di semiapertura 7.5° al bordo d'attacco e 10° al bordo d'uscita viene posto in una corrente supersonica uniforme di MDM, che nelle condizioni di prova è in fase vapore, in prossimità della curva di saturazione liquido-vapore e del punto critico. Negli esperimenti vengono caratterizzate le due onde d'urto oblique e i ventagli di espansione sul profilo alare. Le onde d'urto oblique vengono studiate al variare dell'angolo di deflessione della corrente (da 6.5° a 10°) tramite regolazione dell'incidenza del profilo. I risultati sperimentali vengono confrontati con teoria e simulazioni numeriche, utilizzando in entrambi i casi equazioni di stato non ideali. Tutte le configurazioni geometriche e le condizioni operative studiate sperimentalmente vengnono riprodotte numericamente. I risultati evidenziano una forte dipendenza delle condizioni post-urto e post-espansione dai rispettivi stati di monte, e la velocità del suono diminuisce attraverso le onde d'urto e aumenta attraverso i ventagli di espansione, coerentemente con la teoria NICFD applicata alle condizioni operative esaminate.