京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第22055号 / 工博第4636号 / 新制||工||1723(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 細田 尚, 教授 戸田 圭一, 准教授 音田 慎一郎 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Slush is a two phase flow of solid particles (crystals) and liquid at the triple point temperature, and constitutes an appealing alternative to liquid propellants for space launchers. The crystals give to the mixture higher density and lower specific enthalpy than liquid, enabling reduced tank volume storage and larger fuel holding time. However, the presence of solid crystals significantly modifies the thermo-hydraulics of the fuel transport, and requires novel predictive tools and diagnostic techniques for efficiently exploiting slush propellants. This thesis contributes to both aspects. In particular, this work studied the flow pressure losses and the heat transfer of solid-liquid mixtures in pipelines, combining experimental and numerical methods. Hydraulic and thermal flow features were analyzed separately with substitute mixtures chosen to mimic the behavior of slush flows in engine fuel feed systems. A dedicated facility was designed and built. The pipeline mounted conventional probes for pressure, temperature and mass flow rate measurements. Moreover, a capacitance-based density meter was developed and validated to measure the mixture's solid content. Optical flow visualization and image processing routines were combined to retrieve particulate phase distribution and velocity fields. The experimental work was complemented with 3D Unsteady Reynolds Averaged Navier Stokes simulations in OpenFOAM. The simulations coupled the Euler-Euler approach with the granular kinetic theory for the treatment of the solid dispersed phase. The model was validated with the experimental results on the pressure drop, heat transfer and solid volume fraction.The resulting physical insights and the proposed empirical correlations on the pressure drop and heat transfer in solid-liquid flows contribute to move a step forward towards slush propelled space launchers. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

<div>High speed flow problems of practical interest require a solution of nonequilibrium aerothermochemistry to accurately predict important flow phenomena including surface heat transfer and stresses. As a majority of these flow problems are in the continuum regime, Computational Fluid Dynamics (CFD) is a useful tool for flow modeling. This work presents the development of a nonequilibrium add-on solver to ANSYS Fluent utilizing user-defined-functions to model salient aspects of nonequilibrium flow in air. The developed solver was verified for several benchmark nonequilibrium flow problems and compared with the available experimental data and other nonequilibrium flow simulations. <br></div><div><br></div><div>The rate of dissociation behind a strong shock in thermochemical nonequilibrium depends on the vibrational excitation of molecules. The Macheret-Fridman (MF) classical impulsive model provides analytical expressions for nonequilibrium dissociation rates. The original form of the model was limited to the dissociation of homonuclear molecules. In this work, a general form of the MF model has been derived and present macroscopic rates applicable for modeling dissociation in CFD. Additionally, some improvements to the prediction of mean energy removed in dissociation in the MF-CFD model has been proposed based on the comparisons with available QCT data. In general, the results from the MF-CFD model upon investigating numerous nonequilibrium flows are promising and the model shows a possibility of becoming the standard tool for investigating nonequilibrium flows in CFD.</div><div><br></div><div>The aerodynamic deorbit experiment (ADE) CubeSat has dragsail to accompany accelerated deorbiting of a CubeSat post-mission. A good estimation of the aerothermal load on a reentry CubeSat is paramount to ensure a predictable reentry. This study investigates the aerothermal load on an ADE CubeSat using the direct simulation Monte Carlo (DSMC) methods and Navier-Stokes-Fourier continuum based methods with slip boundary conditions. The aerothermal load on an ADE CubeSat at 90 km altitude from the DSMC and continuum methods were consistent with each other. The continuum breakdown at a higher altitude of 95 km resulted in a strong disagreement between the continuum and DSMC solutions. Overall, the continuum methods could offer a considerable computational cost saving to the DSMC methods in predicting aerothermal load on an ADE CubeSat at low altitudes.<br> </div>

Les travaux numériques et expérimentaux de cette thèse contribuent à une meilleure compréhension de la dynamique de grosses particules dans un écoulement turbulent. Un premier volet m’a permis de quantifier leur mouvement relatif au fluide, ainsi que leur influence locale sur l’écoulement turbulent. Dans un second volet, j'ai trouvé que l'effet collectif des particules est d'atténuer l’amplitude des fluctuations turbulentes. En revanche, celles-ci n’ont pas d’influence sur les propriétés statistiques fines de l’écoulement. Aussi, ces mesures suggèrent l’existence d’une transition de phase dans les grandes échelles de l’écoulement au-delà d’un seuil critique du nombre de particules. / The numerical and experimental work of this thesis contribute to a better understanding of the dynamics of large particles in a turbulent flow. The first part allowed me to quantify their relative motion to the flow and their local influence on the surrounding flow. In a second part, I found that the collective effect of particles is to reduce the amplitude of turbulent fluctuations. In revanche, they have no influence on the fine statistical properties of the flow. Also, these measurements suggest the existence of a phase transition in the larger scales of the flow beyond a critical threshold of the number of particles.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第19288号 / 工博第4085号 / 新制||工||1630(附属図書館) / 32290 / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 細田 尚, 教授 戸田 圭一, 教授 後藤 仁志 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

No description available.

Modeling of complex flows involving the combined effects of flow transition and streamline curvature using two advanced turbulence models, one in the Reynolds-averaged Navier-Stokes (RANS) category and the other in the hybrid RANS-Large eddy simulation (LES) category is considered in this research effort. In the first part of the research, a new scalar eddy-viscosity model (EVM) is proposed, designed to exhibit physically correct responses to flow transition, streamline curvature, and system rotation effects. The four equation model developed herein is a curvature-sensitized version of a commercially available three-equation transition-sensitive model. The physical effects of rotation and curvature (RC) enter the model through the added transport equation, analogous to a transverse turbulent velocity scale. The eddy-viscosity has been redefined such that the proposed model is constrained to reduce to the original transition-sensitive model definition in nonrotating flows or in regions with negligible RC effects. In the second part of the research, the developed four-equation model is combined with a LES technique using a new hybrid modeling framework, dynamic hybrid RANS-LES. The new framework is highly generalized, allowing coupling of any desired LES model with any given RANS model and addresses several deficiencies inherent in most current hybrid models. In the present research effort, the DHRL model comprises of the proposed four-equation model for RANS component and the MILES scheme for LES component. Both the models were implemented into a commercial computational fluid dynamics (CFD) solver and tested on a number of engineering and generic flow problems. Results from both the RANS and hybrid models show successful resolution of the combined effects of transition and curvature with reasonable engineering accuracy, and for only a small increase in computational cost. In addition, results from the hybrid model indicate significant levels of turbulent fluctuations in the flowfield, improved accuracy compared to RANS models predictions, and are obtained at a significant reduction of computational cost compared to full LES models. The results suggest that the advanced turbulence modeling techniques presented in this research effort have potential as practical tools for solving low/high Re flows over blunt/curved bodies for the prediction of transition and RC effects.

We develop a viscoelastic version of the volume of fluid algorithm for tracking deformable interfaces. The code uses the piecewise linear interface calculation method to reconstruct the interface, the continuous surface force formulation to model interfacial tension forces and utilizes the semi-implicit Stokes solver (enabling computations at low Reynolds numbers). The algorithm is primarily designed to simulate the flow of superposed fluids and the drop in a flow problem in 2D shear flows of viscoelastic and/or Newtonian fluids. The code is validated against linear stability theory for the two-layer flow case and against experimental and other documented numerical investigations for the droplet-matrix case. / Ph. D.

Viscous drag reduction plays a vital role in increasing the performance of vehicles. However, there are only so many measurement techniques that can quickly and accurately measure this when compared to pressure drag measurement techniques. The current study makes use of one of the direct and robust measurement techniques that exist, called the Oil Film Interferometry (OFI) to estimate skin friction on the NASA/Virginia Tech BeVERLI (Benchmark Validation Experiment for RANS and LES Investigations) hill. This project aims to develop a detailed database of non-equilibrium, separated turbulent boundary layer flows obtained through wind tunnel experiments for CFD validation. Skin friction measurements are obtained at specific critical locations on the hill and in its close proximity. The challenges involved in obtaining skin friction data from these locations are discussed in detail. Detailed discussions on the experimental setup and data processing methodology are presented. Qualitative and quantitative results from each measurement location are discussed along with uncertainties to explain certain key flow physics. Additionally, skin friction coefficients from selected overlapping measurement locations from another experimental flow measurement technique called Laser Doppler Velocimetry (LDV) are compared with OFI, and a cross-instrument study is performed. Finally, results from well-refined RANS CFD simulations are assessed with the experimental results, and critical improvement areas are identified. / Master of Science / Drag force is a parameter that significantly contributes to the performance efficiency of any vehicle moving in a fluid. This force is categorised into two types - pressure and viscous drag- both of which need to be minimised as much as possible to contribute towards higher vehicle performance. While there are numerous measurement techniques and documentation currently available to measure pressure drag, this is not the case with the measurement of viscous drag. Skin friction measurement directly relates to the estimation of viscous drag, but accurate and quick measurement of this quantity highly challenging with countable measurement techniques currently available. Through this project, BeVERLI (Benchmark Validation Experiment for RANS and LES Investigations), a detailed documentation is developed for accurate measurement of skin friction through Oil Film Interferometry (OFI). The results obtained through this measurement is explained with a detailed experimental procedure as well as using a data processing code. The accuracy of these results are then discussed with the results from another flow measurement technique called Laser Doppler Velocimetry (LDV) and from Computational Fluid Dynamics (CFD).

Yes / A set of shallow-water equations (SWEs) based on a k-epsilon Reynold stress model is established to simulate the turbulent flows over a complex roughness bed. The fundamental equations are discretized by the second-order finite-difference method (FDM), in which spatial and temporal discretization are conducted by staggered-grid and leap-frog schemes, respectively. The turbulent model in this study stems from the standard k-epsilon model, but is enhanced by replacing the conventional vertical production with a more rigorous and precise generation derived from the energy spectrum and turbulence scales. To verify its effectiveness, the model is applied to compute the turbulence in complex flow surroundings (including a rough bed) in an abrupt bend and in a natural waterway. The comparison of the model results against experimental data and other numerical results shows the robustness and accuracy of the present model in describing hydrodynamic characteristics, especially turbulence features on the complex roughness bottom. / National Key Research and Development Program of China (Grant No: 2016YFE0122500, 2013CB036401 and 2013CB036402), China Postdoctoral Science Foundation (Grant No: 2016M591184) and Programme of Introducing Talents of Discipline to Universities (Grant No: BC2018038) / Research Development Fund Publication Prize Award winner, June 2019.