Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

This thesis studies advection modeling for heat transfer and fluid flow problems using a new Non--Inverted Skew Upwind Scheme (called NISUS). Variants of the new scheme are formulated and developed with 8-noded hexahedral elements using the Finite Element Method (FEM)and rectangular elements based on a Finite Volume Method (FVM). A new method of mass weighting to predict convective fluxes of each scalar from the nodal point values is developed. Due to an explicit representation in terms of nodal variables, local inversion of the upwind coefficient matrix is not needed. Also, this thesis evaluates two variants of the new scheme (i.e., 3-node / 3-point and 4-node / 8-point formulations) within a 3--D FEM and a third variant within a 2--D FVM. The 3--D FEM variants are applied to a variety of test problems involving the transport of a scalar variable, while the 2--D FVM variant is applied to fluid flow problems including natural convection in an enclosure and micro--channel flow simulations. The promising performance of NISUS, as compared with exact and previous solutions, is demonstrated both in terms of accuracy and stability. Furthermore, a new data storage format called Compressed Banded Data (CBD) is developed for sparse banded matrices generated by the control volume finite element method (CVFEM). The platform of the new CBD structure permits dynamic switching between various solvers, without any procedural change in the implementation of existing simulation software. The performance of different Krylov techniques with an ILU(0) preconditioner is observed and compared in three test problems with a direct solver. / October 2006

Electric breakdown in power transformers is preceded by pre-breakdown events such as streamers. The understanding of these phenomena is important in order to optimize liquid insulation systems. Earlier works have derived a model that describes streamers in transformer oil and utilized a finite element method to produce numerical solutions. This research investigates the consequences of changing the numerical method to a finite volume-based solver implemented in OpenFOAM. Using a standardized needle-sphere geometry, a number of oil and voltage combinations were simulated, and the results are for the most part similar to those produced by the previous method. In cases with differing results the change is attributed to the more stable numerical performance of the OpenFOAM solver. A proof of concept for the extension of the simulation from a two-dimensional axial symmetry to three dimensions is also presented. / Elektriska genomslag i högspänningstransformatorer föregås av bildandet av elektriskt ledande kanaler som kallas streamers. En god förståelse av detta fenomen är viktigt vid konstruktionen av oljebaserad elektrisk isolation. Tidigare forskning i ämnet har tagit fram en modell för fortplantningen av streamers. Denna modell har sedan lösts numeriskt av ett beräkningsverktyg baserat på finita elementmetoden. I denna uppsats undersöks konsekvenserna av att byta metod till finita volymsmetoden genom att implementera en lösare i OpenFOAM. En standardiserad nål-sfär-geometri har ställts upp och ett flertal kombinationer av oljor och spänningar har simulerats. De flesta resultaten visar god överensstämmande med tidigare forskning medan resultat som avviker har tillskrivits de goda numeriska egenskaperna hos OpenFOAM-lösaren. En ny typ av simulering har även genomförts där simulationen utökas från en tvådimensionell axisymmetrisk geometri til tre dimensioner.

This master thesis was provided by Scania AB. The objective of this thesis was to modify an application in the free Computational Fluid Dynamics software OpenFOAM to be able to handle spray and wall film modeling of a Urea Water Solution together with Conjugate Heat Transfer. The basic purpose is to widen the knowledge of the vaporization process of a Urea Water Solution in the exhaust gas after treatment system for a diesel engine by using OpenFOAM. First, urea has been modeled as a very viscous liquid at low temperature to mimic the solidication process of urea. Second, the development of the new application has been done. At last, test simulations of a simple test case are performed with the new application. The results are then compared with simplied hand calculations to verify a correct behavior of certain exposed source terms. The new application is working properly for the test case but to ensure the reliability, the results need to be compared with another Computational Fluid Dynamics software or more preferable, real experiments. For more advanced geometries, the continued development presented last in this thesis is highly recommended to follow.

Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

An axisymmetric, Navier-Stokes finite volume flow solver, which uses Harten, Lax and van Leer (HLL) and Harten, Lax and van Leer&ndash / Contact (HLLC) upwind flux differencing scheme for spatial and uses Runge-Kutta explicit multi-stage time stepping scheme for temporal discretization on unstructured meshe is developed. Developed solver can solve the compressible axisymmetric flow. The spatial accuracy of the solver can be first or second order accurate. Second order accuracy is achieved by piecewise linear reconstruction. Gradients of flow variables required for piecewise linear reconstruction are calculated by Green-Gauss theorem. Baldwin-Lomax turbulent model is used to compute the turbulent viscosity. Approximate Riemann solver of HLL and HLLC implemented in solver are validated by solving a cylindrical explosion case. Also the solver&rsquo / s capability of solving unstructured, multi-zone domain is investigated by this problem. First and second order results of solver are compared by solving the flow over a circular bump. Axisymmetric flow in solid propellant rocket motor is solved in order to validate the axisymmetric feature of solver. Laminar flow over flat plate is solved for viscous terms validation. Turbulent model is studied in the flow over flat plate and flow with mass injection test cases.

We connect the Mimetic Finite Difference method (MFD) with the finite-volume two-point flux scheme (TPFA) for Voronoi meshes. The main effect is reducing the saddle-point system to a much smaller symmetric-positive definite matrix. In addition, the generalization allows MFD to seamlessly integrate with existing porous media modeling technology. The generalization also imparts the monotonicity property of the TPFA method on MFD. The connection is achieved by altering the consistency condition of the velocity bilinear operator. First-order convergence theory is presented as well as numerical results that support the claims. We demonstrate a methodology for using MFD in modeling fluid flow in fractures coupled with a reservoir. The method can be used for nonplanar fractures. We use the method to demonstrate the effects of fracture curvature on single-phase and multi-phase flows. Standard benchmarks are used to demonstrate the accuracy of the method. The approach is coupled with existing reservoir simulation technology. / text

Gradient-based aerodynamic shape optimization, based on Computational Fluid Dynamics analysis of the flow, is a method that can automatically improve designs of aircraft components. The prospect is to reduce a cost function that reflects aerodynamic performances. When the shape is described by a large number of parameters, the calculation of one gradient of the cost function is only feasible by recourse to techniques that are derived from the theory of optimal control. In order to obtain the best computational efficiency, the so called adjoint method is applied here on the complete mapping, from the parameters of design to the values of the cost function. The mapping considered here includes the Euler equations for compressible flow discretized on unstructured meshes by a median-dual finite-volume scheme, the primal-to-dual mesh transformation, the mesh deformation, and the parameterization. The results of the present research concern the detailed derivations of expressions, equations, and algorithms that are necessary to calculate the gradient of the cost function. The discrete adjoint of the Euler equations and the exact dual-to-primal transformation of the gradient have been implemented for 2D and 3D applications in the code Edge, a program of Computational Fluid Dynamics used by Swedish industries. Moreover, techniques are proposed here in the aim to further reduce the computational cost of aerodynamic shape optimization. For instance, an interpolation scheme is derived based on Radial Basis Functions that can execute the deformation of unstructured meshes faster than methods based on an elliptic equation. In order to improve the accuracy of the shape, obtained by numerical optimization, a moving mesh adaptation scheme is realized based on a variable diffusivity equation of Winslow type. This adaptation has been successfully applied on a simple case of shape optimization involving a supersonic flow. An interpolation technique has been derived based on a mollifier in order to improve the convergence of the coupled mesh-flow equations entering the adaptive scheme. The method of adjoint derived here has also been applied successfully when coupling the Euler equations with the boundary-layer and parabolized stability equations, with the aim to delay the laminar-to-turbulent transition of the flow. The delay of transition is an efficient way to reduce the drag due to viscosity at high Reynolds numbers.

Lorsqu'une suspension concentrée est en écoulement, il est fréquent d'observer que la concentration en particules ne reste pas homogène mais que les particules migrent vers des régions préférentielles de l'écoulement. Globalement, il existe deux types de modèles pour décrire cette interaction entre l'écoulement et la structure qui apparaît dans la suspension. Les premiers sont assez phénoménologiques et reposent sur l'étude des collisions qui surviennent entre les particules en écoulement. Ils donnent lieu à une description de la migration en terme de diffusion des particules. Le second modèle, appelé "Suspension Balance Model" (SBM), fait appel à l'action des contraintes normales d'origine particulaire engendrées par l'écoulement. Ce modèle semble très pertinent mais son utilisation souffre du manque de données sur les contraintes normales dans les suspensions. Ce manuscrit rapporte une étude expérimentale et numérique de la rhéologie et de la migration induite par un écoulement dans des suspensions non-browniennes concentres. La partie expérimentale consiste à mesurer, en géométrie torsionnelle plan-plan, la viscosité, les deux différences de contraintes normales et le tenseur des contraintes particulaires. Les contraintes particulaires déterminées expérimentalement sont alors injectées dans le "Suspension Balance Model" qui relie le flux de particules à la divergence du tenseur des contraintes particulaires, dans le cas d'un écoulement de cisaillement simple. Les équations couplées de la conservation de la masse, des particules et du moment d'inertie sont implémentées dans OpenFOAM et résolues par la méthode des volumes finis. Les résultats numériques sont comparés à des résultats numériques et expérimentaux de la littérature. Enfin, le SBM est généralisé pour être utilisé dans tout type d'écoulement à 2 dimensions ; les cas du cisaillement d'un nuage de particules et de l'effet de la gravité dans un écoulement de Couette horizontal sont traités.

The aim of this research dissertation is at studying natural and mixed convections of fluid flows, and to develop and validate numerical schemes for interface tracking in order to treat incompressible and immiscible fluid flows, later. In a first step, an original numerical method, based on Finite Volume discretizations, is developed for modeling low Mach number flows with large temperature gaps. Three physical applications on air flowing through vertical heated parallel plates were investigated. We showed that the optimum spacing corresponding to the peak heat flux transferred from an array of isothermal parallel plates cooled by mixed convection is smaller than those for natural or forced convections when the pressure drop at the outlet keeps constant. We also proved that mixed convection flows resulting from an imposed flow rate may exhibit unexpected physical solutions; alternative model based on prescribed total pressure at inlet and fixed pressure at outlet sections gives more realistic results. For channels heated by heat flux on one wall only, surface radiation tends to suppress the onset of recirculations at the outlet and to unify the walls temperature. In a second step, the mathematical model coupling the incompressible Navier-Stokes equations and the Level-Set method for interface tracking is derived. Improvements in fluid volume conservation by using high order discretization (ENO-WENO) schemes for the transport equation and variants of the signed distance equation are discussed