A mathematical model for a consolidation process of a highlyconcentrated, flocculated suspension is developed.Thesuspension is treated as a mixture of a fluid and solidparticles by an Eulerian two-phase fluid model.W e characterizethe suspension by constitutive relations correlating thestresses, interaction forces, and inter-particle forces toconcentration and velocity gradients.This results in threeempirically determined material functions: a hystereticpermeability, a non-Newtonian viscosity and a non-reversibleparticle interaction pressure.P arameters in the models arefitted to experimental data. A simulation program using finite difference methods both intime and space is applied to one and two dimensional testcases.Numer ical experiments are performed to study the effectof different viscosity and permeability models. The effect ofshear on consolidation rate is studied and it is significantwhen the permeability hysteresis model is employed.

The prediction of two-phase fluid-solid (gas-solid and liquid-solid) flow remains a major challenge in many engineering and industrial applications. Numerical modeling of these flows is complicated and various studies have been conducted to improve the model performance. In the present work, the two-fluid model of Bolio et al. (1995), developed for dilute turbulent gas-solid flows, is employed to investigate turbulent two-phase liquid-solid flows in both a vertical pipe and a horizontal channel. <p>Fully developed turbulent gas-solid and liquid-solid flows in a vertical pipe and liquid-solid (slurry) flow in a horizontal channel are numerically simulated. The momentum equations for the fluid and solid phases were solved using the finite volume technique developed by Patankar (1980). Mean and fluctuating velocities for both phases, solids concentration, and pressure drop were predicted and compared with the available experimental data. In general, the mean velocity predictions for both phases were in good agreement with the experimental data for vertical flow cases, considered in this work. <p>For dilute gas-solid vertical flows, the predictions were compared with the experimental data of Tsuji et al. (1984). The gas-phase fluctuating velocity in the axial direction was significantly under-predicted while the results for the solids fluctuating velocity were mixed. There was no data to compare the solids concentration but the profiles looked realistic. The pressure drop was observed to increase with increasing Reynolds number and mass loading when compared with the data of Henthorn et al. (2005). The pressure drop first decreased as particle size increased and then started increasing. This behaviour was shown by both experimental data and model predictions. <p>For the liquid-solid flow simulations the mean velocity profiles for both phases, and the liquid-phase turbulence kinetic energy predictions (for dilute flow case), were in excellent agreement with the experimental data of Alejbegovic et al. (1995) and Sumner et al. (1990). The solids concentration profiles were poorly predicted, especially for the lighter particles. The granular temperature profiles, accounting for the solids velocity fluctuations, for the dilute flow case failed to agree with the data, although they captured the overall trend. The liquid-solid pressure drop predictions, using the present model, were only successful for some particles. <p>The solids concentration predictions for the horizontal flow case were similar to the experimental measurements of Salomon (1965), except for a sharp peak at the bottom wall and the opposite curvature. The mixture velocity profiles were asymmetric, due to the addition of particles, and were similar to the experimental data, though only a partial agreement was observed between the predictions and the data.<p>A conclusion from this work is that the present model, which was developed for dilute gas-solid flows, is inadequate when liquid-solid flows are considered. Further improvements, such as including the interstitial fluid effects while computing the liquid-phase stress, are needed to improve the predictive capability of this two-fluid model.

The prediction of two-phase fluid-solid (gas-solid and liquid-solid) flow remains a major challenge in many engineering and industrial applications. Numerical modeling of these flows is complicated and various studies have been conducted to improve the model performance. In the present work, the two-fluid model of Bolio et al. (1995), developed for dilute turbulent gas-solid flows, is employed to investigate turbulent two-phase liquid-solid flows in both a vertical pipe and a horizontal channel. <p>Fully developed turbulent gas-solid and liquid-solid flows in a vertical pipe and liquid-solid (slurry) flow in a horizontal channel are numerically simulated. The momentum equations for the fluid and solid phases were solved using the finite volume technique developed by Patankar (1980). Mean and fluctuating velocities for both phases, solids concentration, and pressure drop were predicted and compared with the available experimental data. In general, the mean velocity predictions for both phases were in good agreement with the experimental data for vertical flow cases, considered in this work. <p>For dilute gas-solid vertical flows, the predictions were compared with the experimental data of Tsuji et al. (1984). The gas-phase fluctuating velocity in the axial direction was significantly under-predicted while the results for the solids fluctuating velocity were mixed. There was no data to compare the solids concentration but the profiles looked realistic. The pressure drop was observed to increase with increasing Reynolds number and mass loading when compared with the data of Henthorn et al. (2005). The pressure drop first decreased as particle size increased and then started increasing. This behaviour was shown by both experimental data and model predictions. <p>For the liquid-solid flow simulations the mean velocity profiles for both phases, and the liquid-phase turbulence kinetic energy predictions (for dilute flow case), were in excellent agreement with the experimental data of Alejbegovic et al. (1995) and Sumner et al. (1990). The solids concentration profiles were poorly predicted, especially for the lighter particles. The granular temperature profiles, accounting for the solids velocity fluctuations, for the dilute flow case failed to agree with the data, although they captured the overall trend. The liquid-solid pressure drop predictions, using the present model, were only successful for some particles. <p>The solids concentration predictions for the horizontal flow case were similar to the experimental measurements of Salomon (1965), except for a sharp peak at the bottom wall and the opposite curvature. The mixture velocity profiles were asymmetric, due to the addition of particles, and were similar to the experimental data, though only a partial agreement was observed between the predictions and the data.<p>A conclusion from this work is that the present model, which was developed for dilute gas-solid flows, is inadequate when liquid-solid flows are considered. Further improvements, such as including the interstitial fluid effects while computing the liquid-phase stress, are needed to improve the predictive capability of this two-fluid model.

<p>A mathematical model for a consolidation process of a highlyconcentrated, flocculated suspension is developed.Thesuspension is treated as a mixture of a fluid and solidparticles by an Eulerian two-phase fluid model.W e characterizethe suspension by constitutive relations correlating thestresses, interaction forces, and inter-particle forces toconcentration and velocity gradients.This results in threeempirically determined material functions: a hystereticpermeability, a non-Newtonian viscosity and a non-reversibleparticle interaction pressure.P arameters in the models arefitted to experimental data.</p><p>A simulation program using finite difference methods both intime and space is applied to one and two dimensional testcases.Numer ical experiments are performed to study the effectof different viscosity and permeability models. The effect ofshear on consolidation rate is studied and it is significantwhen the permeability hysteresis model is employed.</p>

Cette travail de thèse est lié au besoin de modélisation des écoulements diphasiques en générateurs de vapeur (entrée liquide et sortie vapeur). La démarche proposée consiste à faire le choix d'une modélisation hybride de l'écoulement, en scindant la phase gaz en deux champs, modélisés de manières différentes. Ainsi, les petites bulles sphériques sont modélisées avec une approche dispersée classique avec le modèle eulérien à deux fluides, et les bulles déformées sont simulées à l'aide d'une méthode de localisation d'interface.Le travail effectué porte sur la mise en place, la vérification et la validation du modèle dédié aux larges bulles déformées, ainsi que le couplage entre les deux approches pour le gaz gaz, permettant des premiers calculs de démonstration utilisant l'approche hybride complète. / The present study is related to the need of modeling the two-phase flows occuring in a steam generator (liquid at inlet and vapour at outlet). The choice is made to investigate a hybrid modeling of the flow, considering the gas phase as two separated fields, each one being modeled with different closure laws. In so doing, the small and spherical bubbles are modeled through a dispersed approach within the two-fluid model, and the distorted bubbles are simulated with an interface locating method.The main outcome is about the implementation, the verification and the validation of the model dedicated to the large and distorted bubbles, as well as the coupling of the two approaches for the gas, allowing the presentation of demonstration calculations using the so-called hybrid approach.

Ce travail a été consacré à la simulation numérique des équations de la mécanique des fluides par des méthodes de volumes finis implicites. Tout d’abord, nous avons étudié et mis en place une version implicite du schéma de Roe pour les écoulements monophasiques et diphasiques compressibles. Grâce à la méthode de Newton utilisée pour résoudre les systèmes nonlinéaires, nos schémas sont conservatifs. Malheureusement, la résolution de ces systèmes est très coûteuse. Il est donc impératif d’utiliser des algorithmes de résolution performants. Pour des matrices de grande taille, on utilise souvent des méthodes itératives dont la convergence dépend de leur spectre. Nous avons donc étudié le spectre du système linéaire et proposé une stratégie de Scaling pour améliorer le conditionnement de la matrice. Combinée avec le préconditionneur classique ILU, notre stratégie de Scaling a réduit de façon significative le nombre d’itérations GMRES du système local et le temps de calcul. Nous avons également montré l’intérêt du schéma centré pour la simulation de certains écoulements à faible nombre de Mach. Nous avons ensuite étudié et implémenté la méthode de décomposition de domaine pour les écoulements compressibles. Nous avons proposé une nouvelle variable interface qui rend la méthode du complément de Schur plus facile à construire et nous permet de traiter les termes de diffusion. L’utilisation du solveur itératif GMRES plutôt que Richardson pour le système interface apporte aussi une amélioration des performances par rapport aux autres méthodes. Nous pouvons également découper notre domaine de calcul en un nombre quelconque de sous-domaines. En utilisant la stratégie de Scaling pour le système interface, nous avons amélioré le conditionnement de la matrice et réduit le nombre d’itérations GMRES de ce système. En comparaison avec le calcul distribué classique, nous avons montré que notre méthode est robuste et efficace. / This thesis deals with numerical simulations of compressible fluid flows by implicit finite volume methods. Firstly, we studied and implemented an implicit version of the Roe scheme for compressible single-phase and two-phase flows. Thanks to Newton method for solving nonlinear systems, our schemes are conservative. Unfortunately, the resolution of nonlinear systems is very expensive. It is therefore essential to use an efficient algorithm to solve these systems. For large size matrices, we often use iterative methods whose convergence depends on the spectrum. We have studied the spectrum of the linear system and proposed a strategy, called Scaling, to improve the condition number of the matrix. Combined with the classical ILU preconditioner, our strategy has reduced significantly the GMRES iterations for local systems and the computation time. We also show some satisfactory results for low Mach-number flows using the implicit centered scheme. We then studied and implemented a domain decomposition method for compressible fluid flows. We have proposed a new interface variable which makes the Schur complement method easy to build and allows us to treat diffusion terms. Using GMRES iterative solver rather than Richardson for the interface system also provides a better performance compared to other methods. We can also decompose the computational domain into any number of subdomains. Moreover, the Scaling strategy for the interface system has improved the condition number of the matrix and reduced the number of GMRES iterations. In comparison with the classical distributed computing, we have shown that our method is more robust and efficient.

We develop a two-fluid model (TFM) for heat transfer in dense non-Brownian suspensions. Specifically, we propose closure relations for the inter-phase heat transfer coefficient and the thermal diffusivity of the particle phase based on calibration against experimental data. The model is then employed to simulate non-isothermal flow in an annular Couette cell. We find that, when the shear rate is controlled by the rotation of the inner cylinder, both the shear and thermal gradients are responsible for particle migration. Within the TFM framework, we identify the origin and functional form of a "thermo-rheological" migration force that rationalizes our observations. Furthermore, we apply our model to flow in eccentric Couette cells. Our simulations reveal that the system's heat transfer coefficient is affected by both the classic shear-induced migration of particles and the newly identified thermo-rheological migration effect. Finally, we employed the proposed computational TFM framework to analyze electronics cooling by forced convection for microchannel cooling. We used a suspensions of high thermal conductivity (Boron Nitride) particles in a 3M Fluorinert FC-43 cooling fluid. Three-dimensional simulations were run to quantify the temperature distributions under uniform heating (5 W) and under hot-spot heating (2 W/cm^2) conditions. A 100 K junction level temperature improvement (enhanced thermal spreading) was seen for hot-spot heating and 15 K was observed for uniform heating, demonstrating the enhanced cooling capabilities of dense particulate suspensions of high-conductivity particles, over a clear FC-43 fluid.

abstract: The applications utilizing nanoparticles have grown in both industrial and academic areas because of the very large surface area to volume ratios of these particles. One of the best ways to process and control these nanoparticles is fluidization. In this work, a new microjet and vibration assisted (MVA) fluidized bed system was developed in order to fluidize nanoparticles. The system was tested and the parameters optimized using two commercially available TiO2 nanoparticles: P25 and P90. The fluidization quality was assessed by determining the non-dimensional bed height as well as the non-dimensional pressure drop. The non-dimensional bed height for the nanosized TiO2 in the MVA system optimized at about 5 and 7 for P25 and P90 TiO2, respectively, at a resonance frequency of 50 Hz. The non-dimensional pressure drop was also determined and showed that the MVA system exhibited a lower minimum fluidization velocity for both of the TiO2 types as compared to fluidization that employed only vibration assistance. Additional experiments were performed with the MVA to characterize the synergistic effects of vibrational intensity and gas velocity on the TiO2 P25 and P90 fluidized bed heights. Mathematical relationships were developed to correlate vibrational intensity, gas velocity, and fluidized bed height in the MVA. The non-dimensional bed height in the MVA system is comparable to previously published P25 TiO2 fluidization work that employed an alcohol in order to minimize the electrostatic attractions within the bed. However, the MVA system achieved similar results without the addition of a chemical, thereby expanding the potential chemical reaction engineering and environmental remediation opportunities for fluidized nanoparticle systems. In order to aid future scaling up of the MVA process, the agglomerate size distribution in the MVA system was predicted by utilizing a force balance model coupled with a two-fluid model (TFM) simulation. The particle agglomerate size that was predicted using the computer simulation was validated with experimental data and found to be in good agreement. Lastly, in order to demonstrate the utility of the MVA system in an air revitalization application, the capture of CO2 was examined. CO2 breakthrough time and adsorption capacities were tested in the MVA system and compared to a vibrating fluidized bed (VFB) system. Experimental results showed that the improved fluidity in the MVA system enhanced CO2 adsorption capacity. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2019

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第18590号 / 工博第3951号 / 新制||工||1607(附属図書館) / 31490 / 京都大学大学院工学研究科原子核工学専攻 / (主査)教授 功刀 資彰, 教授 中部 主敬, 准教授 横峯 健彦 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

No description available.