Numerical simulation of fresh SCC flow in wall and beam elements using flow dynamics models

Abstract : Recently, there is a great interest to study the flow characteristics of suspensions in different environmental and industrial applications, such as snow avalanches, debris flows, hydrotransport systems, and material casting processes. Regarding rheological aspects, the majority of these suspensions, such as fresh concrete, behave mostly as non-Newtonian fluids. Concrete is the most widely used construction material in the world. Due to the limitations that exist in terms of workability and formwork filling abilities of normal concrete, a new class of concrete that is able to flow under its own weight, especially through narrow gaps in the congested areas of the formwork was developed. Accordingly, self-consolidating concrete (SCC) is a novel construction material that is gaining market acceptance in various applications. Higher fluidity characteristics of SCC enable it to be used in a number of special applications, such as densely reinforced sections. However, higher flowability of SCC makes it more sensitive to segregation of coarse particles during flow (i.e., dynamic segregation) and thereafter at rest (i.e., static segregation). Dynamic segregation can increase when SCC flows over a long distance or in the presence of obstacles. Therefore, there is always a need to establish a trade-off between the flowability, passing ability, and stability properties of SCC suspensions. This should be taken into consideration to design the casting process and the mixture proportioning of SCC. This is called “workability design” of SCC. An efficient and non-expensive workability design approach consists of the prediction and optimization of the workability of the concrete mixtures for the selected construction processes, such as transportation, pumping, casting, compaction, and finishing. Indeed, the mixture proportioning of SCC should ensure the construction quality demands, such as demanded levels of flowability, passing ability, filling ability, and stability (dynamic and static). This is necessary to develop some theoretical tools to assess under what conditions the construction quality demands are satisfied. Accordingly, this thesis is dedicated to carry out analytical and numerical simulations to predict flow performance of SCC under different casting processes, such as pumping and tremie applications, or casting using buckets. The L-Box and T-Box set-ups can evaluate flow performance properties of SCC (e.g., flowability, passing ability, filling ability, shear-induced and gravitational dynamic segregation) in casting process of wall and beam elements. The specific objective of the study consists of relating numerical results of flow simulation of SCC in L-Box and T-Box test set-ups, reported in this thesis, to the flow performance properties of SCC during casting. Accordingly, the SCC is modeled as a heterogeneous material. Furthermore, an analytical model is proposed to predict flow performance of SCC in L-Box set-up using the Dam Break Theory. On the other hand, results of the numerical simulation of SCC casting in a reinforced beam are verified by experimental free surface profiles. The results of numerical simulations of SCC casting (modeled as a single homogeneous fluid), are used to determine the critical zones corresponding to the higher risks of segregation and blocking. The effects of rheological parameters, density, particle contents, distribution of reinforcing bars, and particle-bar interactions on flow performance of SCC are evaluated using CFD simulations of SCC flow in L-Box and T-box test set-ups (modeled as a heterogeneous material). Two new approaches are proposed to classify the SCC mixtures based on filling ability and performability properties, as a contribution of flowability, passing ability, and dynamic stability of SCC.

Natural convection from a vertical cone in a porous medium due to the combined effects of heat and mass diffusion with non-uniform wall temperature/concentration or heat/mass flux and suction/injection

An analysis has been carried out to study the non-Darcy natural convention flow of Newtonian fluids on a vertical cone embedded in a saturated porous medium with power-law variation of the wall temperature/concentration or heat/mass flux and suction/injection with the streamwise distance x. Both non-similar and self-similar solutions have been obtained. The effects of non-Darcy parameter, ratio of the buoyancy forces due to mass and heat diffusion, variation of wall temperature/concentration or heat/mass flux and suction/injection on the Nusselt and Sherwood numbers have been studied.

Combined convection in power-law fluids along a nonisothermal vertical plate in a porous medium

The problem of combined convection from vertical surfaces in a porous medium saturated with a power-law type non-Newtonian fluid is investigated. The transformed conservation laws are solved numerically for the case of variable surface heat flux conditions. Results for the details of the velocity and temperature fields as well as the Nusselt number have been presented. The viscosity index ranged from 0.5 to 2.0.

Free convection in power-law fluids near a three-dimensional stagnation point

Numerical results are presented for the free-convection boundary-layer equations of the Ostwald de-Waele non-Newtonian power-law type fluids near a three-dimensional (3-D) stagnation point of attachment on an isothermal surface. The existence of dual solutions that are three-dimensional in nature have been verified by means of a numerical procedure. An asymptotic solution for very large Prandtl numbers has also been derived. Solutions are presented for a range of values of the geometric curvature parameter c, the power-law index n, and the Prandtl number Pr.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Run off Conditions on Rimming Power Law Fluids

The rimming ?ow of a power-law ?uid in the inner surface of a horizontal rotating cylinder is investigated. Exploiting the fact that the liquid layer is thin, the simplest lubrication theory is applied. The generalized run-off condition for the steady-state ?ow of the power-law liquid is derived. In the bounds implied by this condition, ?lm thickness admits a continuous solution. In the supercritical case when the mass of non-Newtonian liquid exceeds a certain value or the speed of rotation is less than an indicated limit, a discontinuous solution is possible and a hydraulic jump may occur in the steady-state regime. The location and height of the hydraulic jump for the power-law liquid is determined.

The Inertio-Elastic Planar Entry Flow of Low-Viscosity Elastic Fluids in Micro-fabricated Geometries

The non-Newtonian flow of dilute aqueous polyethylene oxide (PEO) solutions through microfabricated planar abrupt contraction-expansions is investigated. The contraction geometries are fabricated from a high-resolution chrome mask and cross-linked PDMS gels using the tools of soft-lithography. The small length scales and high deformation rates in the contraction throat lead to significant extensional flow effects even with dilute polymer solutions having time constants on the order of milliseconds. The dimensionless extra pressure drop across the contraction increases by more than 200% and is accompanied by significant upstream vortex growth. Streak photography and videomicroscopy using epifluorescent particles shows that the flow ultimately becomes unstable and three-dimensional. The moderate Reynolds numbers (0.03 ≤ Re ≤ 44) associated with these high Deborah number (0 ≤ De ≤ 600) microfluidic flows results in the exploration of new regions of the Re-De parameter space in which the effects of both elasticity and inertia can be observed. Understanding such interactions will be increasingly important in microfluidic applications involving complex fluids and can best be interpreted in terms of the elasticity number, El = De/Re, which is independent of the flow kinematics and depends only on the fluid rheology and the characteristic size of the device.

Effect of cosurfactant on the supramolecular structure and physicochemical properties of non-ionic biocompatible microemulsions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The carreau-yasuda fluids: a skin friction equation for turbulent flow in pipes and kolmogorov dissipative scales

In this work the turbulent flow of the Non-Newtonian Carreau-Yasuda fluid will be studied. A skin friction equation for the turbulent flow of Carreau-Yasuda fluids will be derived assuming a logarithmic behavior of the turbulent mean velocity for the near wall flow out of the viscous sub layer. An alternative near wall characteristic length scale which takes into account the effects of the relaxation time will be introduced. The characteristic length will be obtained through the analysis of viscous region near the wall. The results compared with experimental data obtained with Tylose (methyl hydroxil cellulose) solutions showing good agreement. The relations between scales integral and dissipative obtained for length, time, velocity, kinetic energy, and vorticity will be derived for this type of fluid. When the power law index approach to unity the relations reduces to Newtonian case.

Numerical study of ultrasound induced non-linear shape and size bubble oscillations in viscoelastic media

The theoretical study of forced bubble oscillations is motivated by the importance of cavitation bubbles and oscillating encapsulated microbubbles (i.e. contrast agents) in medical sciences. In more details,theoretical studies on bubble dynamics addressing the sound-bubble interaction phenomenon provide the basis for understanding the dynamics of contrast agent microbubbles used in medical diagnosis and of non-linearly oscillating cavitation bubbles in the case of high-intensity ultrasound therapy. Moreover, the inclusion of viscoelasticity is of vital importance for an accurate theoretical analysis since most biological tissues and fluids exhibit non-Newtonian behavior.

One-dimensional Model for Fluids of Third-grade in Tubes with Constant Radius

Our goal in this paper is to extend previous results obtained for Newtonian and secondgrade fluids to third-grade fluids in the case of an axisymmetric, straight, rigid and impermeable tube with constant cross-section using a one-dimensional hierarchical model based on the Cosserat theory related to fluid dynamics. In this way we can reduce the full threedimensional system of equations for the axisymmetric unsteady motion of a non-Newtonian incompressible third-grade fluid to a system of equations depending on time and on a single spatial variable. Some numerical simulations for the volume flow rate and the the wall shear stress are presented.

Nonsimilar incompressible laminar boundary-layer flows in micropolar fluids

The solution of the steady laminar incompressible nonsimilar boundary-layer problem for micropolar fluids over two-dimensional and axisymmetric bodies has been presented. The partial differential equations governing the flow have been transformed into new co-ordinates having finite range. The resulting equations have been solved numerically using implicit finite-difference scheme. The computations have been carried out for a cylinder and a sphere. The results indicate that the separation in micropolar fluids occurs at earlier streamwise locations as compared to Newtonian fluids. The skin friction and velocity profiles depend on the shape of the body and are almost insensitive to microrotation or coupling parameter, provided the coupling parameter is small. On the other hand, the microrotation profiles and microrotation gradient depend on the microrotation parameter and they are insensitive to the coupling parameter.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.