Les lanceurs ont un impact sur la composition de l'atmosphere, et en particulier sur l'ozone stratospherique. Parmi tous les types de propulsion, les moteurs à propergol solide ont fait l'objet d'une attention particulière car leurs émissions sont responsables d'un appauvrissement significatif d'ozone dans le panache des lanceurs lors des premières heures suivant le lancement. Ce phénomène est principalement dû à la conversion de l'acide chlorhydrique, un composé chimique présent en grandes quantités dans les émissions de ce type de moteur, en chlore actif qui réagit par la suite avec l'ozone dans un cycle catalytique similaire à celui responsable du "trou de la couche d'ozone", cette diminution périodique de l'ozone en Antarctique. Cette conversion se produit dans le panache supersonique, où les hautes températures favorisent une seconde combustion entre certaines espèces chimiques du panache et l'air ambiant. L'objectif de cette étude est d'évaluer la concentration de chlore actif dans le panache d'un moteur à propergol solide en utilisant la technique des Simulations aux Grandes Echelles (SGE). Le gaz est injecté à travers la tuyère d'un moteur et une méthode de couplage entre deux instances du solveur de mécanique des fluides est utilisée pour étendre autant que possible le domaine de calcul derrière la tuyère (jusqu'à l'équivalent de 400 diamètres de sortie de la tuyère). Cette méthodologie est validée par une première SGE sans chimie, en analysant les caractéristiques de l'écoulement supersonique avec co-écoulement obtenu par ce calcul. Ensuite, le chimie mettant en jeu la conversion des espèces chlorées a été étudiée au moyen d'un modèle "hors-ligne" permettant de résoudre une chimie complexe le long de lignes de courant extraites d'un écoulement moyenné dans le temps résultant du calcul précédent (non réactif). Enfin, une SGE multi-espèces est réalisée, incluant un schéma chimique auparavant réduit afin de limiter le coût de calcul. Cette simulation représente une des toutes premières SGE d'un jet supersonique réactif, incluant la tuyère, effectuée sur un domaine de calcul aussi long. En capturant avec précision le mélange du panache avec l'air ambiant ainsi que les interactions entre turbulence et combustion, la technique des simulations aux grandes échelles offre une évaluation des concentrations des espèces chimiques dans le jet d'une precision inédite. Ces résultats peuvent être utilisés pour initialiser des calculs atmosphériques sur de plus larges domaines, afin de modéliser les réactions entre chlore actif et ozone et de quantifier l'appauvrissement en ozone dans le panache. / Rockets have an impact on the chemical composition of the atmosphere, and particularly on stratospheric ozone. Among all types of propulsion, Solid-Rocket Motors (SRMs) have given rise to concerns since their emissions are responsible for a severe decrease in ozone concentration in the rocket plume during the first hours after a launch. The main source of ozone depletion is due to the conversion of hydrogen chloride, a chemical compound emitted in large quantities by ammonium perchlorate based propellants, into active chlorine compounds, which then react with ozone in a destructive catalytic cycle, similar to those responsible for the Antartic "Ozone hole". This conversion occurs in the hot, supersonic exhaust plume, as part of a strong second combustion between chemical species of the plume and air. The objective of this study is to evaluate the active chlorine concentration in the far-field plume of a solid-rocket motor using large-eddy simulations (LES). The gas is injected through the entire nozzle of the SRM and a local time-stepping method based on coupling multi-instances of the fluid solver is used to extend the computational domain up to 400 nozzle exit diameters downstream of the nozzle exit. The methodology is validated for a non-reactive case by analyzing the flow characteristics of the resulting supersonic co-flowing under-expanded jet. Then the chemistry of chlorine is studied off-line using a complex chemistry solver applied on trajectories extracted from the LES time-averaged flow-field. Finally, the online chemistry is analyzed by means of the multi-species version of the LES solver using a reduced chemical scheme. To the best of our knowledge, this represents one of the first LES of a reactive supersonic jet, including nozzle geometry, performed over such a long computational domain. By capturing the effect of mixing of the exhaust plume with ambient air and the interactions between turbulence and combustion, LES offers an evaluation of chemical species distribution in the SRM plume with an unprecedented accuracy. These results can be used to initialize atmospheric simulations on larger domains, in order to model the chemical reactions between active chlorine and ozone and to quantify the ozone loss in SRM plumes.

La modélisation par simulation des grandes échelles (Large-Eddy Simulation - LES) des processus physiques régissant la couche limite atmosphérique (CLA) demeure complexe de part la difficulté des modèles à capter l'évolution de la turbulence entre différentes conditions de stratification. De ce fait, l'étude LES du cycle diurne complet de la CLA comprenant des situations convectives la journée et des conditions stables la nuit est très peu documenté. La simulation de la couche limite stable où la turbulence est faible, intermittente et qui est caractérisée par des structures turbulentes de petite taille est tout particulièrement compliquée. En conséquence, la capacité de la LES à bien reproduire les conditions météorologiques de la CLA, notamment en situation stable, est étudiée à l'aide du code de mécanique des fluides développé par EDF R&D, Code_Saturne. Dans une première étude, le modèle LES est validé sur un cas de couche limite convective quasi stationnaire sur terrain homogène. L'influence des modèles sous-maille de Smagorinsky, Germano-Lilly, Wong-Lilly et WALE (Wall-Adapting Local Eddy-viscosity) ainsi que la sensibilité aux méthodes de paramétrisation sur les champs moyens, les flux et les variances est discutées. Dans une seconde étude le cycle diurne complet de la CLA pendant la campagne de mesure Wangara est modélisé. L'écart aux mesures étant faible le jour, ce travail se concentre sur les difficultés rencontrées la nuit à bien modéliser la couche limite stable. L'impact de différents modèles sous-maille ainsi que la sensibilité au coefficient de Smagorinsky ont été analysés. Par l'intermédiaire d'un couplage radiatif réalisé en LES, les répercussions du rayonnement infrarouge et solaire sur le jet de basse couche nocturne et le gradient thermique près de la surface sont exposées. De plus l'adaptation de la résolution du domaine à l'intensité de la turbulence et la forte stabilité atmosphérique durant l'expérience Wangara sont commentées. Enfin un examen des oscillations numériques inhérentes à Code_Saturne est réalisé afin d'en limiter les effets

O presente trabalho tem como objetivo estudar escoamentos de tornados e sua ação sobre corpos imersos empregando ferramentas numéricas da Engenharia do Vento Computacional (EVC). Os tornados constituem-se atualmente em uma das causas de desastres naturais no Brasil, especialmente nas regiões sul e sudeste do país, como também em alguns países vizinhos. Os efeitos gerados são geralmente localizados e de curta duração, podendo ser devastadores dependendo da escala do tornado. Tais características dificultam a realização de estudos detalhados a partir de eventos reais, o que levou ao desenvolvimento de modelos experimentais e numéricos. A abordagem numérica é utilizada neste trabalho para a simulação de tornados, a qual se baseia nas equações de Navier-Stokes e na equação de conservação de massa, considerando a hipótese de pseudo-compressibilidade e condições isotérmicas. Para escoamentos com turbulência utiliza-se a Simulação Direta de Grandes Escalas com o modelo clássico de Smagorinsky para as escalas inferiores à resolução da malha (Large Eddy Simulation ou LES em inglês). A discretização das equações fundamentais do escoamento se realiza com um esquema explícito de dois passos de Taylor-Galerkin, onde o Método dos Elementos Finitos é empregado na discretização espacial utilizando-se o elemento hexaédrico trilinear isoparamétrico com um ponto de integração e controle de modos espúrios Na presença de corpos imersos que se movem para simular os deslocamentos dos tornados, o escoamento é descrito cinematicamente através de uma formulação Arbitrária Lagrangeana-Euleriana (ALE) que inclui um esquema de movimento de malha. Tornados são reproduzidos através da simulação numérica de dispositivos experimentais e do Modelo de Vórtice Combinado de Rankine (RCVM). Exemplos clássicos da Dinâmica dos Fluidos Computacional são apresentados inicialmente para a verificação das ferramentas numéricas implementadas. Finalmente, problemas envolvendo tornados móveis e estacionários são analisados, incluindo sua ação sobre corpos imersos. Nos modelos baseados em experimentos, a variação da relação de redemoinho determinou os diferentes padrões de escoamento observados no laboratório. Nos exemplos de modelo de vórtice, quando o tornado impactou o corpo imerso gerou picos de forças em todas as direções e, após a passar pelo mesmo, produziu uma alteração significativa na estrutura do vórtice. / Analyses of tornado flows and its action on immersed bodies using numerical tools of Computational Wind Engineering (CWE) are the main aims of the present work. Tornadoes are currently one of the causes of natural disasters in Brazil, occurring more frequently in the southern and southeastern regions of the country, as well as in some neighboring countries. Effects are usually localized, presenting a short time interval, which can be devastating depending on the scale of the tornado. These characteristics difficult to carry out detailed studies based on real events, leading to the development of experimental and numerical models. The numerical approach is used in this work for the simulation of tornadoes, which is based on the Navier-Stokes equations and the mass conservation equation, considering the hypothesis of pseudo-compressibility and isothermal conditions. For turbulent flows, Large Eddy Simulation (LES) is used with the classical Smagorinsky model for sub-grid scales Discretization is performed the explicit two-step Taylor-Galerkin scheme, where the Finite Element Method is used in spatial discretization using isoparametric trilinear hexahedral elements with one-point quadrature and hourglass control. In the presence of immersed bodies that are moving in order to simulate translating tornadoes, the flow is kinematically described through a Lagrangian-Eulerian Arbitrary (ALE) formulation, which includes a mesh motion scheme. Tornadoes are reproduced using numerical simulation of experimental devices and the Rankine Combined Vortex Model (RCVM). Classical examples of Computational Fluid Dynamics are presented initially for the verification of the numerical tools implemented here. Finally, problems involving moving and stationary tornadoes are analyzed, including their actions on immersed bodies. For models based on experiments, the variation of the swirl ratio determined the different flow patterns observed in the laboratory. In the vortex model examples, when the tornado impacted on the immersed body, peaks of forces were generated in all directions and, after passing over it, a significant change in the structure of the vortex was produced.

QC 20110926

Cette thèse se focalise sur la compréhension et la prédiction des instabilités de combustion dans les systèmes à haute pression. Elle s'oriente autour de la simulation numérique d’un banc d'essai, opéré à l'université de Purdue, comprenant un injecteur caractéristique des moteurs-fusées et dont les propriétés acoustiques peuvent varier à l'aide d’un tube d'injection mobile. Une méthode d'initialisation et d'allumage pour les calculs LES de chambres de combustions terminées par une tuyère est présentée. Un point de fonctionnement instable est choisi pour étudier le mécanisme de l'instabilité. Les simulations sont comparées aux résultats expérimentaux en terme de fréquence et structure du mode instable. La fonction de transfert de flamme est calculée à l'aide du modèle n − τ puis implémentée dans un solveur acoustique (ne résolvant que les perturbations acoustiques à partir de l'équation de Helmholtz en écoulement réactif). Différents modèles d'impédance de tuyère, extraits de la littérature, sont comparés et leurs impacts sur les résultats de stabilité sont analysés. Le théorème d’impédance translatée est implémenté dans le solveur acoustique pour analyser, à faible coût de calcul, l’influence de la variation de la longueur du tube d'injection. Des écarts entre les fréquences prédites et celles trouvées expérimentalement subsistent mais la carte de stabilité de l’expérience est bien reproduite. / This thesis focuses on the understanding and the prediction of combustion instability in high-pressure devices. A model rocket combustor, tested experimentally at Purdue University, with continuously variable acoustic properties, thanks to a variable-length injector tube, is simulated. A method to initialize and ignite Large-Eddy-Simulation (LES) calculation of combustion chamber surrounded by nozzle is proposed. An unstable operating point is then chosen to investigate the mechanism of the instability. The simulations are compared to experimental results in terms of frequency and mode structure. The flame transfer function is calculated using the n − τ model to feed an acoustic solver which solves only the acoustic perturbation using a Helmholtz equation in reacting flows. The importance of the modeling of the nozzles impedance is studied through the main theories in the literature. The impedance translation theorem is implemented in the acoustic solver to analyze at low cost the influence of the variation of the injector tube. Despite differences in frequency of the instability, the stability map of the experiment is well reproduced.

O presente trabalho tem como objetivo estudar escoamentos de tornados e sua ação sobre corpos imersos empregando ferramentas numéricas da Engenharia do Vento Computacional (EVC). Os tornados constituem-se atualmente em uma das causas de desastres naturais no Brasil, especialmente nas regiões sul e sudeste do país, como também em alguns países vizinhos. Os efeitos gerados são geralmente localizados e de curta duração, podendo ser devastadores dependendo da escala do tornado. Tais características dificultam a realização de estudos detalhados a partir de eventos reais, o que levou ao desenvolvimento de modelos experimentais e numéricos. A abordagem numérica é utilizada neste trabalho para a simulação de tornados, a qual se baseia nas equações de Navier-Stokes e na equação de conservação de massa, considerando a hipótese de pseudo-compressibilidade e condições isotérmicas. Para escoamentos com turbulência utiliza-se a Simulação Direta de Grandes Escalas com o modelo clássico de Smagorinsky para as escalas inferiores à resolução da malha (Large Eddy Simulation ou LES em inglês). A discretização das equações fundamentais do escoamento se realiza com um esquema explícito de dois passos de Taylor-Galerkin, onde o Método dos Elementos Finitos é empregado na discretização espacial utilizando-se o elemento hexaédrico trilinear isoparamétrico com um ponto de integração e controle de modos espúrios Na presença de corpos imersos que se movem para simular os deslocamentos dos tornados, o escoamento é descrito cinematicamente através de uma formulação Arbitrária Lagrangeana-Euleriana (ALE) que inclui um esquema de movimento de malha. Tornados são reproduzidos através da simulação numérica de dispositivos experimentais e do Modelo de Vórtice Combinado de Rankine (RCVM). Exemplos clássicos da Dinâmica dos Fluidos Computacional são apresentados inicialmente para a verificação das ferramentas numéricas implementadas. Finalmente, problemas envolvendo tornados móveis e estacionários são analisados, incluindo sua ação sobre corpos imersos. Nos modelos baseados em experimentos, a variação da relação de redemoinho determinou os diferentes padrões de escoamento observados no laboratório. Nos exemplos de modelo de vórtice, quando o tornado impactou o corpo imerso gerou picos de forças em todas as direções e, após a passar pelo mesmo, produziu uma alteração significativa na estrutura do vórtice. / Analyses of tornado flows and its action on immersed bodies using numerical tools of Computational Wind Engineering (CWE) are the main aims of the present work. Tornadoes are currently one of the causes of natural disasters in Brazil, occurring more frequently in the southern and southeastern regions of the country, as well as in some neighboring countries. Effects are usually localized, presenting a short time interval, which can be devastating depending on the scale of the tornado. These characteristics difficult to carry out detailed studies based on real events, leading to the development of experimental and numerical models. The numerical approach is used in this work for the simulation of tornadoes, which is based on the Navier-Stokes equations and the mass conservation equation, considering the hypothesis of pseudo-compressibility and isothermal conditions. For turbulent flows, Large Eddy Simulation (LES) is used with the classical Smagorinsky model for sub-grid scales Discretization is performed the explicit two-step Taylor-Galerkin scheme, where the Finite Element Method is used in spatial discretization using isoparametric trilinear hexahedral elements with one-point quadrature and hourglass control. In the presence of immersed bodies that are moving in order to simulate translating tornadoes, the flow is kinematically described through a Lagrangian-Eulerian Arbitrary (ALE) formulation, which includes a mesh motion scheme. Tornadoes are reproduced using numerical simulation of experimental devices and the Rankine Combined Vortex Model (RCVM). Classical examples of Computational Fluid Dynamics are presented initially for the verification of the numerical tools implemented here. Finally, problems involving moving and stationary tornadoes are analyzed, including their actions on immersed bodies. For models based on experiments, the variation of the swirl ratio determined the different flow patterns observed in the laboratory. In the vortex model examples, when the tornado impacted on the immersed body, peaks of forces were generated in all directions and, after passing over it, a significant change in the structure of the vortex was produced.

O presente trabalho tem como objetivo estudar escoamentos de tornados e sua ação sobre corpos imersos empregando ferramentas numéricas da Engenharia do Vento Computacional (EVC). Os tornados constituem-se atualmente em uma das causas de desastres naturais no Brasil, especialmente nas regiões sul e sudeste do país, como também em alguns países vizinhos. Os efeitos gerados são geralmente localizados e de curta duração, podendo ser devastadores dependendo da escala do tornado. Tais características dificultam a realização de estudos detalhados a partir de eventos reais, o que levou ao desenvolvimento de modelos experimentais e numéricos. A abordagem numérica é utilizada neste trabalho para a simulação de tornados, a qual se baseia nas equações de Navier-Stokes e na equação de conservação de massa, considerando a hipótese de pseudo-compressibilidade e condições isotérmicas. Para escoamentos com turbulência utiliza-se a Simulação Direta de Grandes Escalas com o modelo clássico de Smagorinsky para as escalas inferiores à resolução da malha (Large Eddy Simulation ou LES em inglês). A discretização das equações fundamentais do escoamento se realiza com um esquema explícito de dois passos de Taylor-Galerkin, onde o Método dos Elementos Finitos é empregado na discretização espacial utilizando-se o elemento hexaédrico trilinear isoparamétrico com um ponto de integração e controle de modos espúrios Na presença de corpos imersos que se movem para simular os deslocamentos dos tornados, o escoamento é descrito cinematicamente através de uma formulação Arbitrária Lagrangeana-Euleriana (ALE) que inclui um esquema de movimento de malha. Tornados são reproduzidos através da simulação numérica de dispositivos experimentais e do Modelo de Vórtice Combinado de Rankine (RCVM). Exemplos clássicos da Dinâmica dos Fluidos Computacional são apresentados inicialmente para a verificação das ferramentas numéricas implementadas. Finalmente, problemas envolvendo tornados móveis e estacionários são analisados, incluindo sua ação sobre corpos imersos. Nos modelos baseados em experimentos, a variação da relação de redemoinho determinou os diferentes padrões de escoamento observados no laboratório. Nos exemplos de modelo de vórtice, quando o tornado impactou o corpo imerso gerou picos de forças em todas as direções e, após a passar pelo mesmo, produziu uma alteração significativa na estrutura do vórtice. / Analyses of tornado flows and its action on immersed bodies using numerical tools of Computational Wind Engineering (CWE) are the main aims of the present work. Tornadoes are currently one of the causes of natural disasters in Brazil, occurring more frequently in the southern and southeastern regions of the country, as well as in some neighboring countries. Effects are usually localized, presenting a short time interval, which can be devastating depending on the scale of the tornado. These characteristics difficult to carry out detailed studies based on real events, leading to the development of experimental and numerical models. The numerical approach is used in this work for the simulation of tornadoes, which is based on the Navier-Stokes equations and the mass conservation equation, considering the hypothesis of pseudo-compressibility and isothermal conditions. For turbulent flows, Large Eddy Simulation (LES) is used with the classical Smagorinsky model for sub-grid scales Discretization is performed the explicit two-step Taylor-Galerkin scheme, where the Finite Element Method is used in spatial discretization using isoparametric trilinear hexahedral elements with one-point quadrature and hourglass control. In the presence of immersed bodies that are moving in order to simulate translating tornadoes, the flow is kinematically described through a Lagrangian-Eulerian Arbitrary (ALE) formulation, which includes a mesh motion scheme. Tornadoes are reproduced using numerical simulation of experimental devices and the Rankine Combined Vortex Model (RCVM). Classical examples of Computational Fluid Dynamics are presented initially for the verification of the numerical tools implemented here. Finally, problems involving moving and stationary tornadoes are analyzed, including their actions on immersed bodies. For models based on experiments, the variation of the swirl ratio determined the different flow patterns observed in the laboratory. In the vortex model examples, when the tornado impacted on the immersed body, peaks of forces were generated in all directions and, after passing over it, a significant change in the structure of the vortex was produced.

This thesis is concerned with studies on spatially developing high speed mixing layers with twin objectives: (a) to provide enhanced and detailed understanding of spatial development of two-dimensional mixing layer emanating from splitter plate through large eddy simulation (LES, from now on) technique and (b) to evolve a consistent strategy for Unsteady Reynolds Averaged Navier-Stokes (URANS) approach to mixing layer calculations. The inspiration for this work arose out of the explanations that were being developed for the reduction in the mixing layer thickness with compressibility (measured by a parameter called convective Mach number, Mc). The reasons centered around increased stability, increase in compressible dissipation that was later discounted in favor of reduction in production and pressure-strain terms (with Mc, of course). These were obtained with direct numerical simulations (DNS) or LES techniques with homogeneous shear flow or temporal mixing layer. As apart, there was also a wide held view that using RANS (steady) techniques did not capture the compressibility effects when used in a way described above and so classical industrial codes for computing mixing- layer-embedded flows are unsuitable for such applications. Other important aspects that come out of the examination of literature are: the mixing layer growth is controlled in the initial stages by the double- boundary layer profile over the splitter plate and results in the mixing layer growth that is somewhat irregular due to doubling and merging of vertical structures. The view point of a smooth growth of the mixing layer is a theo- retical approximation arising out of the use of a smooth tan-hyperbolic profile that results at larger distances from the splitter plate. For all practical applications, it is inferred that the initial development is what is important because the processes of ignition and stable combustion occur close to the splitter plate. For these reasons, it was thought that understanding the development of the mixing layer is best dealt with using accurate spatial simulation with the appropriate initial profile. The LES technique used here is drawn from an OpenFOAM approach for dissimilar gases and uses one-equation Eddy Model for SGS stresses. The temporal discretization is second order accurate backward Euler and spatial discretization is fourth order least squares; the algorithm used for solving the equations is PISO and the parallelized code uses domain decomposition approach to cover large spatial domain. The calculations are performed with boundary layer profiles over the splitter plate and an initial velocity field with white noise-like fluctuations to simulate the turbulence as in the experiments. Grid independence studies are performed and several experimental cases are considered for comparison with measured data on the velocity and temperature fields as well as turbulent statistics. These comparisons are excellent for the mean field behavior and moderately acceptable for turbulent kinetic energy and shear stress. To further benefit from the LES approach, the details of the mixing layer are calculated as a function of four independent parameters on which the growth depends: convective Mach number (Mc = (U1 -U2)/ (a1 +a2)), stream speed ratio (r = U2=U1), stream density ratio (s = p2/p1) and the average velocity of the two streams ((U1+U2)=2) and examine the various terms in the equations to enable answering the questions discussed earlier. It is uncovered that r has significant influence on the attainment of self similarity (which also implies on the rate of removal of velocity defect in the double-boundary layer profile) and other parameters have a very weak influence. The minimum velocity variation with distance from the splitter plate has the 1/paxial distance behavior like in wakes; however, after a distance, departure to linear rise occurs and the distance it takes for this to appear is delayed with Mc. Other features such as the coherent structures, their merger or break up, the area of the structures, convective velocity information extraction from the coherent structures, the behavior of the pressure field in the mixing layer through the field are elucidated in detail; the behavior of the correlations between parameters (like pressure, velocity etc) at different points is used to elucidate the coherence of their fluctuating field. The effects of the parameters on the energy spectra have expected trends. An examination of the kinetic energy budget terms reveals that • the production term is the main source of the xx turbulence stress, whereas it is not significant in the yy component. • A substantial portion of this is carried by the pressure-velocity coupling from the xx direction to the yy direction, which becomes the main source term in the yy component. • Both, the production term as well as the pressure-velocity term show a clear decrease with increase in Mc. The high point of the thesis is related to using the understanding derived from an analysis of various source terms in the kinetic energy balance to evolve an unsteady Reynolds Averaged Navier Stokes (URANS) model for calculating high speed mixing layers, a subject that has eluded international research till now. It recognizes that the key feature affected by ompressibility is related to the anisotropy of the stress tensor. The relationship between stress component (_Txy) and the velocity gradient (Sxy) as obtained from LES is set out in the form of a simple relationship accounting for the effects of other parameters obtained earlier in this thesis. A minor influence due to _Tyy is extracted by describing its dependence on Sxy again as gleaned from LES studies. The needed variation of Prandtl and Schmidt numbers through the field is extracted. While the detailed variations can in fact be taken into account in URANS simulations, a simple assumption of these values being around 0.3 is chosen for the present simulations of URANS. Introduction of these features into the momentum equation gives the much expected variation of the reduction in the growth rate of the mixing layer with convective Mach number as in experiments. The relationships that can be used in high speed mixing layers are Introduction of these features into the momentum equation gives the much expected variation of the reduction in the growth rate of the mixing layer with convective Mach number as in experiments. This is then a suggested new approach to solve high speed mixing layers. While it can be thought that the principal contributions of the thesis are complete here, an additional segment is presented related to entropy view of the mixing layer. This study that considers the mixing layer with two different species expresses various terms involved in the entropy conservation equation and obtains the contribution of various terms on the entropy change for various Mc. It is first verified that the entropy derived from the conservation equation matches with those calculated from fluid properties, entropy being a state variable. It is shown that irreversible diffusion comes down the most with convective Mach number. Left: This image shows pictorially the flow of source of turbulent stress from the axial to the cross wise turbulent stress. Production (Σ) of turbulence happens mainly in the xx direction, a part of it is carried by the pressure-velocity correlation to the yy direction, which itself has a low production. With increasing Mc, both the production as well as the pressure-velocity correlation decrease. Right: This image shows the growth rate obtained from simulations scaled with the incompressible growth rate, of LES and RANS in the background of experiments (others). As is clear, the growth rate obtained is well within the band of experimental results.

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第21582号 / 理博第4489号 / 新制||理||1644(附属図書館) / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)准教授 竹見 哲也, 准教授 石岡 圭一, 教授 秋友 和典 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

<div>The growing need to improve the power density of compact thermal systems necessitates developing new techniques to modulate the convective heat transfer efficiently. In the present research, acoustic streaming is evaluated as a potential technology to achieve this objective. Numerical simulations using the linearized and fully non-linear Navier-Stokes equations are employed to characterize the physics underlying this process. The linearized Navier-Stokes equations accurately replicate the low-frequency flow unsteadiness, which is used to find the optimal control parameters. Local and global stability analysis tools were developed to identify the modes with a global and positive heat transfer effect.</div><div><br></div><div>High-fidelity numerical simulations are performed to evaluate the effect of the excitation at selected frequencies, directed by the linear stability analysis, on the heat and momentum transport in the flow. Results indicate that, under favorable conditions, superimposing an acoustic wave, traveling along with the flow, can <i>resonate</i> within the domain and lead to a significant heat transfer enhancement with minimal skin friction losses. Two main flow configurations are considered; at the fixed Reynolds number Re_b=3000, in the supersonic case, 10.1% heat transfer enhancement is achieved by an 8.4% skin friction increase; however, in the subsonic case, 10% enhancement in heat transfer only caused a 5.3% increase to the skin friction. The deviation between these two quantities suggests a violation of the Reynolds analogy. This study is extended to include a larger Reynolds number, namely Re_b=6000 at M_b=0.75 and a similar response is observed. The effect of excitation amplitude and frequency on the resonance, limit-cycle oscillations, heat transfer, and skin friction are also investigated here.</div><div><br></div><div>Applying acoustic waves normal to the flow in the spanwise direction disrupts the near-wall turbulent structures that are primarily responsible for heat and momentum transport near the solid boundary. Direct numerical simulations were employed to investigate this technique in a supersonic channel flow at M_b=1.5 and Re_b=3000. The external excitation is applied through a periodic body force in the spanwise direction, mimicking loudspeakers placed on both walls that are operating with a 180^o phase shift. By keeping the product of forcing amplitude A_f and pulsation period (<i>T</i>) constant, spanwise velocity perturbations are generated with a similar amplitude at different frequencies. Under this condition, spanwise pulsations at <i>T</i>=20 and <i>T</i>=10 show up to 8% reduction in Nusselt number as well as the skin friction coefficient. Excitation at higher or lower frequencies fails to achieve such high level of modulations in heat and momentum transport processes near the walls.<br> <br>In configurations involving a spatially-developing boundary layer, a computational setup that includes laminar, transitional, and turbulent regions inside the domain is considered and the impact of acoustic excitation on this flow configuration has been characterized. Large-eddy simulations with dynamic Smagorinsky sub-grid scale modeling has been implemented, due to the excessive computational cost of DNS calculations at high-Reynolds numbers. The optimal excitation frequency that resembles the mode chosen for the fully-developed case has been identified via global stability analysis. Fully non-linear simulations of the spatially-developing boundary layer subjected to the excitation at this frequency reveal an interaction between the <i>pulsations</i> and the perturbations originated from the tripping which creates a re-laminarization zone traveling downstream. Such technique can locally enhance or reduce the heat transfer along the walls.<br></div>