京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21100号 / 工博第4464号 / 新制||工||1694(附属図書館) / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 黒瀬 良一, 教授 中部 主敬, 教授 吉田 英生 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Temporally developing direct numerical simulations (T-DNS) are performed and validated for bypass transition of a zero pressure gradient flat plate boundary layer to understand the interplay between pressure-strain terms and flow instability mechanisms, and to propose and validate a phenomenological hypothesis for the identification of a robust transition onset marker for use in transition-sensitive Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations. Results show that transition initiates at a location where the slow pressure-strain term becomes more dominant than the rapid term in the pre-transitional boundary layer region. The slow pressure strain term is responsible for the transfer of turbulence energy from the streamwise component to other components while the rapid pressure strain term counteracts with the slow term in the pre-transitional regime before transition onset akin to a shear sheltering like effect. The relative magnitudes of the slow and rapid terms thus provide a basis for the development of physically meaningful large-scale parameters that can be used as a transition onset marker for Reynolds averaged Navier-Stokes RANS simulations.

The present work seeks to fully investigate, describe and characterize the distinct flow regimes existing within a mixing vessel at various rotational speeds. This investigation is computational in nature and simulates the flow within a baffled tank containing a Rushton turbine of the standard configuration. For a Re based on impeller diameter and blade rotational speed (Re â ¡ Ï ND2/μ) the following flow regimes were identified and investigated in detail: Reverse/reciprocating flows at very low Re (<10); stalled flows at low Re (â 10); laminar pumping flow for higher Re and transitional pumping flow (10 squared < Re <10 to the 4th). For the three Re numbers 1, 10 and 28, it was found that for the higher Re number (28), the flow exhibited the familiar outward pumping action associated with radial impellers under turbulent flow conditions. However, as the Re number decreases, the net radial flow during one impeller revolution was reduced and for the lowest Re number a reciprocating motion with negligible net pumping was observed. In order to elucidate the physical mechanism responsible for the observed flow pattern at low Re, the forces acting on a fluid element in the radial direction were analyzed. Based on this analysis, a simplified quasi-analytic model of the flow was developed that gives a satisfactory qualitative, as well as quantitative representation of the flow at very low Re. Investigation of the transitional flow regime (Re â 3000) includes a compilation and characterization of ensemble and turbulent quantities such as the Reynolds stress components, dissipation length η and time scales Ï , as well a detailed investigation of the near-impeller flow and trailing vortex. Calculation and compilation of all terms in the turbulent kinetic energy transport equation was performed (including generation and the illusive turbulent pressure work). Specifically, the most important transport mechanism was turbulent convection/diffusion from the impeller disk-plane/trailing vortex region. Mean flow transport of turbulent kinetic energy was primarily towards the impeller disk-plane and radially outward from the trailing vortex region. The turbulent pressure work was found to partially counteract turbulent convection. Turbulent dissipation followed by turbulent viscous work were found to be the least important mechanism responsible for turbulent transport with both terms being maximized within the vortex region and at the disk-plane down-stream from the vortices. / Ph. D.

The collisions frequency of dispersed phases (particles, droplets, bubbles) in a turbulent carrier phase is a fundamental quantity that is needed for modeling multiphase flows with applications to chemical processes, minerals flotation, food science, and many other industries. In this dissertation, numerical simulations are performed to determine collisions frequency of bi-dispersed particles (solid particles and bubbles) in homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and Large Eddy simulations (LES) are conducted to determine velocity fluctuations of the carrier phase. The DNS results are used to validate existing theoretical models as well as the LES results. The dissertation also presents a CFD-based flotation model for predicting the pulp recovery rate in froth flotation machines. In the direct numerical simulations work, particles and bubbles suspended in homogeneous isotropic turbulence are tracked and their collisions frequency is determined as a function of particle Stokes number. The effects of the dispersed phases on the carrier phase are neglected. Particles and bubbles of sizes on the order of Kolmogorov length scale are treated as point masses. Equations of motion of dispersed phases are integrated simultaneously with the equations of the carrier phase using the same time stepping scheme. In addition to Stokes drag, the pressure gradient in the carrier phase and added-mass forces are also included. The collision model used here allows overlap of particles and bubbles. Collisions kernel, radial relative velocity, and radial distribution function found by DNS are compared to theoretical models over a range of particle Stokes number. In general, good agreement between DNS and recent theoretical models is obtained for radial relative velocity for both particle-particle and particle-bubble collisions. The DNS results show that around Stokes number of unity particles of the same group undergo expected preferential concentration while particles and bubbles are segregated. The segregation behavior of particles and bubbles leads to a radial distribution function that is less than one. Existing theoretical models do not account for effects of this segregation behavior of particles and bubbles on the radial distribution function. In the large-eddy simulations efforts, the dissertation addresses the importance of the subgrid fluctuations on the collisions frequency and investigates techniques for predicting those fluctuations. The cases studied are of particles-particles and particles-bubbles collisions at Reynolds number Re_λ = 96. A study is conducted first by neglecting the effects of subgrid velocity fluctuations on particles and bubbles motions. It is found that around Stokes number of unity solid particles of the same group undergo the well known preferential concentration as observed in the DNS. Effects of pressure gradient on the particles are negligible due to their small sizes. Bubbles as a low inertia particles are very sensitive to subgrid velocity and acceleration fields where the effects of pressure gradient in the carrier phase are dominant. However, particle-bubble radial distribution functions from LES are not as low as that from DNS. To account for the effects of subgrid field on the dispersion of particles and bubbles, a new multifractal methodology has been developed to construct a subgrid vorticity field from the resolved vorticity field in frame work of LES. A Poisson's solver is used to obtain the subgrid velocity field from the subgrid vorticity field. Accounting for the subgrid velocity fluctuations (but neglecting pressure gradient) produced minor changes in the radial distribution function for particle-particle and particle-bubble collisions. We conclude from this study that for accurate particle tracking in LES the subgrid velocity fluctuations must be dynamically realizable field (temporally and spatially correlated with the large scale motion). Adding random SGS velocity fluctuations is not enough to capture the correct radial distribution functions of dispersed phases especially for bubbles-particles collisions where the pressure gradient term ( or acceleration Du_f′/Dt) is responsible for particle-bubble segregation around particle Stokes number near one. A CFD-based model for minerals flotation machines has been developed in this dissertation. The objective of flotation models is to predict the recovery rate of minerals from a flotation cell. The developed model advances the state-of-the-art of pulp recovery rate prediction by incorporating validated theoretical collisions frequency models and detailed hydrodynamics from two-phase flow simulations. Spatial distributions of dissipation rate and air volume fraction are determined by the two-phase hydrodynamic simulations. Knowing these parameters throughout the machine is essential in understanding the effectiveness of different components of flotation machine (rotor, stator or disperser, jets) on the flotation efficiency. The developed model not only predicts the average pulp recovery rate but also it indicates regions of high/low recovery rates. The CFD-based flotation model presented here can be used to determine the dependence of recovery rate constant at any locality within the pulp based on particle diameter, particle specfic gravity, contact angle, and surface tension. / Ph. D.

The focus of this thesis is the numerical study of subcritical transition to turbulence in boundary-layer flows. For the most part, boundary layers with uniform suction are considered. Constant homogeneous suction counteracts the spatial growth of the boundary layer, rendering the flow parallel. This enables research approaches which are not feasible in the context of spatially developing flows. In the first part, the laminar–turbulent separatrix of the asymptotic suction boundary layer (ASBL) is investigated numerically by means of an edge-tracking algorithm. The obtained edge states experience recurrent dynamics, going through calm and bursting phases. The self-sustaining mechanism bears many similarities with the classical regeneration cycle of near-wall turbulence. The recurrent simple structure active during calm phases is compared to the nucleation of turbulence events in bypass transition originating from delocalised initial conditions. The implications on the understanding of the bypass-transition process and the edge state's role are discussed. Based on this understanding, a model is constructed which predicts the position of the nucleation of turbulent spots during free-stream turbulence induced transition in spatially developing boundary-layer flow. This model is used together with a probabilistic cellular automaton (PCA), which captures the spatial spreading of the spots, correctly reproducing the main statistical characteristics of the transition process. The last part of the thesis is concerned with the spatio-temporal aspects of turbulent ASBL in extended numerical domains near the onset of sustained turbulence. The different behaviour observed in ASBL, i.e. absence of sustained laminar–turbulent patterns, which have been reported in other wall-bounded flows, is associated with different character of the large-scale flow. In addition, an accurate quantitative estimate for the lowest Reynolds number with sustained turbulence is obtained / <p>QC 20160429</p>

<p>Transition from laminar to turbulent flow has beentraditionally studied in terms of exponentially growingeigensolutions to the linearized disturbance equations.However, experimental findings show that transition may occuralso for parameters combinations such that these eigensolutionsare damped. An alternative non-modal growth mechanism has beenrecently identified, also based on the linear approximation.This consists of the transient growth of streamwise elongateddisturbances, mainly in the streamwise velocity component,called streaks. If the streak amplitude reaches a thresholdvalue, secondary instabilities can take place and provoketransition. This scenario is most likely to occur in boundarylayer flows subject to high levels of free-stream turbulenceand is the object of this thesis. Different stages of theprocess are isolated and studied with different approaches,considering the boundary layer flow over a flat plate. Thereceptivity to free-stream disturbances has been studiedthrough a weakly non-linear model which allows to disentanglethe features involved in the generation of streaks. It is shownthat the non-linear interaction of oblique waves in thefree-stream is able to induce strong streamwise vortices insidethe boundary layer, which, in turn, generate streaks by thelift-up effect. The growth of steady streaks is followed bymeans of Direct Numerical Simulation. After the streaks havereached a finite amplitude, they saturate and a new laminarflow, characterized by a strong spanwise modulation isestablished. Using Floquet theory, the instability of thesestreaks is studied to determine the features of theirbreakdown. The streak critical amplitude, beyond which unstablewaves are excited, is 26% of the free-stream velocity. Theinstability appears as spanwise (sinuous-type) oscillations ofthe streak. The late stages of the transition, originating fromthis type of secondary instability, are also studied. We foundthat the main structures observed during the transition processconsist of elongated quasi-streamwise vortices located on theflanks of the low speed streak. Vortices of alternating signare overlapping in the streamwise direction in a staggeredpattern.</p><p><strong>Descriptors:</strong>Fluid mechanics, laminar-turbulenttransition, boundary layer flow, transient growth, streamwisestreaks, lift-up effect, receptivity, free-stream turbulence,nonlinear mechanism, streak instability, secondary instability,Direct Numerical Simulation.</p> / QC 20100518

Le développement de nouvelles technologies pour le transport aérien moins polluant est de plus en plus basé sur la simulation numérique, qui nécessite alors une description fiable de la chimie. Pour la plupart des carburants, la description de la combustion nécessite des mécanismes détaillés mais leur utilisation dans une simulation numérique de combustion turbulente est limitée par le coût calcul. Des mécanismes cinétiques réduits et des méthodes de tabulation ont été proposés pour surmonter ce problème. Ces descriptions chimiques simplifiées ayant été développées dans le cadre de configurations laminaires, cette thèse propose de les évaluer dans des configurations turbulentes: une DNS de flamme prémélangée méthane/air de type Bunsen et une LES d’un brûleur expérimental. Les mécanismes sont analysés en termes de structure de flamme, paramètres de flamme globaux, longuer de flamme, prediction des concentrations en espèces majoritaires et des émissions polluantes. Une méthodologie pour évaluer a priori la capacité d’un mécanisme à prédire correctement des phénomènes chimiques tridimensionnels est proposée en se basant sur les résultats de flammes laminaires monodimensionnelles non étirées et étirées. Il ressort que, d’une part, pour construire un mécanisme réduit, il est nécessaire de faire un compromis entre coût calcul, robustesse et qualité des résultats. D’autre part, la qualité des résultats de DNS et LES de configurations tridimensionnelles turbulentes peut être anticipée par une analyse du comportement des schémas réduits dans des configurations simplifiées de flammes monodimensionnelles laminaires non étirées et étirées. / A growing need for numerical simulations based on reliable chemistries has been observed in the last years in order to develop new technologies which could guarantee the reduction of the enviromental impact on air transport. The description of combustion requires the use of detailed kinetic mechanisms for most hydro-carbons. Their use in turbulent combustion simulation is still prohibitive because of their high computational cost. Reduced chemistries and tabulation methods have been proposed to over-come this problem. Since all these reductions have been developed for laminar configurations, this thesis proposes to evaluate their performances in simulations of turbulent configurations such as a DNS of a premixed Bunsen methane/air flame and a LES of an experimental PREC-CINSTA burner. The mechanisms are analysed in terms of flame structure, global burning parameters, flame length, prediction of major species concentrations and pollutant emissions. An a priori methodology based on one-dimensional unstrained and strained laminar flames to evaluate the mechanism capability to predict three-dimensional turbulent flame features is therefore proposed. On the one hand when building a new reduced scheme, its requirements should be fixed compromising the computational cost, the robustness of the chemical description and the desired quality of results. On the other hand, the quality of DNS or LES results in three-dimensional configurations could be anticipated testing the reduced mechanism on laminar one-dimensional premixed unstrained and strained flames.

Ces cinquantes dernières années, la majorité des paramètres de conception des moteurs cryotechniques ont été ajustés en l'absence d'une compréhension détaillée des phénomènes de combustion, en raison des limites des diagnostiques expérimentaux et des capacités de calcul. L'objectif de cette thèse est de réaliser des simulations numériques instationnaires d'écoulements réactifs transcritiques de haute fidélité, pour permettre une meilleure compréhension de la dynamique de flamme dans les moteurs cryotechniques et finalement guider leur amélioration. Dans un premier temps, la thermodynamique gaz-réel et son impact sur les schémas numériques sont présentés. Comme la Simulation aux Grandes Echelles (SGE) comporte des équations filtrées, les effets de filtrages induits par la thermodynamique gaz-réel sont ensuite mis en évidence dans une configuration transcritique type et un opérateur de diffusion artificiel, spécifique au gaz réel, est proposé pour lisser les gradients transcritiques en SGE. Dans un deuxième temps, une étude fondamentale du mélange turbulent et de la combustion dans la zone proche-injecteur des moteurs cryotechniques est menée grâce à la Simulation Numérique Directe (SND). Dans le cas non-réactif, les lâchers tourbillonnaires dans le sillage de la lèvre de l’injecteur jouent un rôle majeur dans le mélange turbulent et provoquent la formation de structures en peigne déjà observées expérimentalement dans des conditions similaires. Dans le cas réactif, la flamme reste attachée à la lèvre de l'injecteur, sans extinction locale, et les structures en peigne disparaissent. La structure de flamme est analysée et différents modes de combustion sont identifiés. Enfin, une étude de flamme-jet transcritique H2/O2, accrochée à un injecteur coaxial avec et sans retrait interne, est menée. Les résultats numériques sont d'abord validés par des données expérimentales pour l'injecteur sans retrait. Ensuite, la configuration avec retrait est comparée à la solution de référence sans retrait et à des données experimentales pour observer les effets de ce paramètre de conception sur l'efficacité de combustion. / In the past fifty years, most design parameters of the combustion chamber of Liquid Rocket Engines (LREs) have been adjusted without a detailed understanding of combustion phenomena, because of both limited experimental diagnostics and numerical capabilities. The objective of the present thesis work is to conduct high-fidelity unsteady numerical simulations of transcritical reacting flows, in order to improve the understanding of flame dynamics in LRE, and eventually provide guidelines for their improvement. First real-gas thermodynamics and its impact on numerical schemes are presented. As Large-Eddy Simulation (LES) involves filtered equations, the filtering effects induced by real-gas thermodynamics are then highlighted in a typical 1D transcritical configuration and a specific real-gas artificial dissipation is proposed to smooth transcritical density gradients in LES. Then, a Direct Numerical Simulation (DNS) study of turbulent mixing and combustion in the near-injector region of LREs is conducted. In the non-reacting case, vortex shedding in the wake of the lip of the injector is shown to play a major role in turbulent mixing, and induces the formation of finger-like structures as observed experimentally in similar operating conditions. In the reacting case, the flame is attached to the injector rim without local extinction and the finger-like structures disappear. The flame structure is analyzed and various combustion modes are identified. Finally, a LES study of a transcritical H2/O2 jet flame, issuing from a coaxial injector with and without inner recess, is conducted. Numerical results are first validated against experimental data for the injector without recess. Then, the recessed configuration is compared to the reference solution and to experimental results, to scrutinize the effects of this design parameter on combustion efficiency.

La perspective lagrangienne, décrivant un écoulement selon les trajectoires de traceurs fluides, est une approche naturelle pour étudier les phénomènes de dispersion dans les écoulements turbulents. En turbulence de paroi, le mouvement des traceurs est influencé par le cisaillement moyen et par une forte inhomogénéité et anisotropie en proche paroi. On étudie les propriétés lagrangiennes d’un écoulement de canal turbulent par simulation numérique directe à un nombre de Reynolds modéré. Les statistiques d’accélération lagrangienne sont comparées aux expériences de suivi de particules réalisées en parallèle à ce travail. Comme en turbulence homogène isotrope (THI), les composantes d’accélération le long des trajectoires lagrangiennes se décorrèlent sur des temps comparables aux plus petites échelles de l’écoulement, tandis que la norme de l’accélération reste corrélée plus longtemps. La persistance d’anisotropie à petite échelle loin de la paroi est constatée par l’existence d’une corrélation croisée non nulle entredeux composantes de l’accélération. On montre que, en conséquence des flux moyens d’énergie cinétique en turbulence de paroi, près des parois les traceurs se déplacent et s’étalent sur des plus grandes distances quand ils sont suivis en arrière dans le temps qu’en avant. La dispersion relative de paires de traceurs est aussi étudiée. Aux temps courts, la séparation des paires est balistique pour toutes les distances à la paroi. Comme en THI, les traceurs se séparent plus rapidement lorsqu’ils sont suivis en arrière dans le temps. Aux temps plus longs, le cisaillement moyen accélère la séparation dans la direction de l’écoulement moyen. Un modèle de cascade balistique initialement proposé pour la THI est adapté aux écoulements inhomogènes / The Lagrangian perspective, describing a flow from the trajectories of fluid tracers, isa natural framework for studying dispersion phenomena in turbulent flows. In wall-boundedturbulence, the motion of fluid tracers is affected by mean shear and by strong inhomogeneityand anisotropy near walls. We investigate the Lagrangian properties of a turbulent channel flowusing direct numerical simulations at a moderate Reynolds number. Lagrangian accelerationstatistics are compared to particle tracking experiments performed in parallel to this work. Asin homogeneous isotropic turbulence (HIT), the acceleration components along Lagrangianpaths decorrelate over time scales representative of the smallest scales of the flow, while theacceleration norm stays correlated for much longer. The persistence of small-scale anisotropy farfrom the wall is demonstrated in the form of a non-zero cross-correlation between accelerationcomponents. As a result of the average fluxes of kinetic energy in wall turbulence, tracers initiallylocated close to the wall travel and spread over longer distances when tracked backwardsin time than forwards. The relative dispersion of tracer pairs is finally investigated. At shorttimes, pair separation is ballistic for all wall distances. As in HIT, relative dispersion is timeasymmetric, with tracers separating faster when tracked backwards in time. At longer times,mean shear dominates leading to rapid separation in the mean flow direction. A ballisticcascade model previously proposed for HIT is adapted to inhomogeneous flows

Escoamentos sobre superfícies côncovas estão sujeitos à instabilidade centrífuga, dando origem a vórtices longitudinais, conhecidos como vórtices de Görtler. Esses vórtices são responsáveis por gerar distorções fortes nos perfis de velocidade. Como os vórtices são contra-rotativos, duas regiões surgem entre os mesmos: uma região de upwash e uma região de downwash. Na região de upwash o fluido próximo à parede é jogado para longe da mesma. Na região de downwash acontece o contrário, o fluido que se desloca a uma velocidade maior é jogado em direção à parede. Os vórtices se amplificam inicialmente de forma linear. À jusante na região não linear de desenvolvimento dos vórtices, a amplitude dos mesmos já é elevada, e há a formação de uma estrutura do tipo cogumelo com a distribuição da componente de velocidade na direção principal do escoamento . Essa nova distribuição de velocidade é tridimensional e difere em muito da camada limite obtida com a solução das equações de Blasius. Levando-se em consideração a camada limite térmica, já foi observado que, na média, há um aumento de transferência de calor na direção da parede. No presente trabalho, é verificado numericamente a transferência de calor na presença de vórtices de Görtler. Para tal, foi desenvolvido e implementado um código de simulação numérica direta espacial (DNS - do inglês Direct Numerical Simulation). Os resultados deste trabalho mostram a intensificação da transferência de calor através dos vórtices de Görtler, tanto no regime não-linear como na instabilidade secundária / Flows over concave surfaces are subject to centrifugal instability. It gives rise to stramwise vortices known as Görtler vortices. These vortices are responsible for generating strong distortions in the velocity profiles. As the vortices are counterrotating, two regions arise between them: a region of uowash and a region of downwash. In the upwash region, the fluid near the wall is convected away from it. In the downwash region the opposite happens, the fluid moving at a faster speed is moved towards the wall. The vortices initially amplify linearly in the downstream. When their amplitude is already high, in the non-linear development region, a mushroom-type structure, with the velocity distribution in the main flow direction, is formed. This new three-dimensional velocity distribution is different from the boundary layer obtained with the solution of Blasius equations. Taking into account a thermal boundary layer, on average, an increase in the heat transfer in the wall direction has been observed. In the present work, it is verified numerically the heat transfer in the presence of Görtler vortices. A simulation code was developed and implemanted usin Direct Numerical Simulation (DNS). The results of this work show the intensification of heat transfer through the Görtler vortices both in the non-linear regime and in the secondary instability