Rotating-surface-driven non-newtonian flow in a cylindrical enclosure
|
|
- Tiffany Loraine Simmons
- 5 years ago
- Views:
Transcription
1 Korea-Australia Rheology Journal Vol. 22, No. 4, December 2010 pp Rotating-surface-driven non-newtonian flow in a cylindrical enclosure Yalin Kaptan 1, *, Ali Ecder 2 and Kunt Atalik 2 1 Hansung University, Department of Mechanical Systems Engineering, Seoul, Korea 2 Bogazici University, Department of Mechanical Engineering, Istanbul, Turkey (Received April 30, 2010; final revision received July 26, 2010; accepted July 29, 2010) Abstract The numerical simulations of non-newtonian flows become difficult as the Weissenberg number increases. The main objective of this study is to generate a robust and efficient matrix-free Inexact Newton-Krylov solver (IN-GMRES) which can deal with the high nonlinearities arising from the increased Weissenberg numbers. In order to achieve this aim, non-newtonian flows of rotating surface driven cylindrical enclosure problem is numerically investigated by using three differential viscoelastic constitutive relations namely Upper Convected Maxwell (UCM), Oldroyd B and Giesekus. The results obtained by IN-GMRES solver are validated and compared with the POLYFLOW simulations. Additionally, the selection of the constitutive relation, effects of the Weissenberg number and effects of the Reynolds number are studied. The simulations indicate that the generated algorithm is capable of solving higher Weissenberg number problems (up to the Elasticity number limit of 130) when compared to the previous studies. Furthermore, it is shown that with the increasing Weissenberg number, the reversed flow can be observed in the flow domain and in some cases, depending on the Reynolds number, re-formation of the Newtonian like flow is possible at high Weissenberg numbers. Keywords : inexact Newton, GMRES, Giesekus, Oldroyd B, Upper Convected Maxwell 1. Introduction The investigation of the vortex flow in moving edge problems has been popular because of their physical applications such as lubrication, ocean circulation, turbo machinery, electronic cooling, rheology and coating. One of the most investigated types of such flows is rotating disc in a cylindrical enclosure problem, which continues to be an immensely-studied subject by many scientists. This problem is one of the main geometries of rheometry applications and it is modeled by Pao (1970, 1972) for Newtonian fluids. In the cited study, one of the plates (bottom or top) is held stationary while the other plate and the casing are rotated. The solutions for small Reynolds numbers are obtained for fixed aspect ratio (δ=1). Likewise, similar problems are considered by Lopez (1996) and Saci et al. (1998) and results for high Reynolds numbers and various aspect ratios for Newtonian fluids are achieved. The non-newtonian case of this problem is investigated by Moroi et al. (2001) both numerically and experimentally. Upper Convected Maxwell, Giesekus, Phan Thien-Tanner and power law models are used for the numerical part *Corresponding author: yalinkaptan@gmail.com 2010 by The Korean Society of Rheology where Particle Tracking Velocimetry (PTV) technique is adapted for the experimental part. The elasticity number limit for the numerical part of the cited study is Recently, Itoh et al. (2006) worked on the same problem by adopting a similar approach to that of Moroi et al. (2001). In the cited study, the solutions for variable aspect ratios and Weissenberg numbers are obtained by using Giesekus and Power Law constitutive models whereas the results are verified with the experimental data collected by laser Doppler anemometer. The dependence of flow structures on Reynolds number and aspect ratio is investigated up to the Elasticity number limit of 0.1 for the numerical part. Another numerical study is realized by Xue et al. (1999), in which UCM flows in a 3-D geometry for different aspect ratios are simulated. In the mentioned study, the Finite Volume Method is used for discretization while Semi-Implicit Method for Pressure Linked Equations (SIMPLE) method is utilized for the solution. The results point out that the flow patterns are altered from Newtonian like to reverse flow with the increasing Elasticity numbers where the limit of the elasticity number is for the numerical part of the study. The main objective of this study is to generate a robust and efficient matrix-free Inexact Newton-Krylov solver which can deal with the high Weissenberg number cases of the non-newtonian flow computations, specifically rotat- Korea-Australia Rheology Journal December 2010 Vol. 22, No
2 Yalin Kaptan, Ali Ecder and Kunt Atalik Fig. 1. Geometry of the problem. ing disc in a cylindrical enclosure problem. What is more, the effects of the Weissenberg number, the Reynolds number and the selection of the constitutive relation for the flow of the rotating disc in cylindrical enclosure geometry will be studied in the following sections. The geometry of the problem can be seen in Fig. 1 which consists of a stationary cylindrical casing and two circular discs at the top and the bottom of the casing (stationary top disc and the rotating bottom disc). By utilizing the axial symmetry, the computation area is reduced to a rectangle which can be observed in the same figure with a hatched region. The aspect ratio for parallel plate problem is defined as δ=h/r and the bottom disc is rotating with an angular velocity Ω. 2. Governing Equations The vector form of the continuity, momentum and constitutive equations for the steady and incompressible flow can be given as follows: u = 0 Re( u )u = p + β 2 u + ( 1 β) τ τ+ We[ u τ ( u) T τ τ u + ατ τ] = u + ( u) T (3) Here τ is the viscoelastic stress tensor, u is the velocity vector, p is the pressure, β is the viscosity ratio (given that µ s is the solvent viscosity and µ p is the polymer viscosity, β = µ s /(µ s + µ p )) while α is the mobility factor. The realistic range of the mobility factor is between 0 and 0.5 (Bird et al. (1987)); however, the values between 0.05 and 0.15 are generally used to model the non-newtonian behavior in the literature (Bird et al. (1987), Itoh et al. (2006)). The definitions for the Reynolds number and the Weissenberg number are; Re = ΩR 2 /ν, We = λu/r where ν is the kinematic viscosity, λ is the relaxation time, U is the characteristic velocity (U = RΩ) and R is the radius. Increased Weissenberg number can be interpreted as the increasing (1) (2) elastic effects in the fluid. The elastic effects are compared by means of a dimensionless number, namely Elasticity number, which is the ratio of Weissenberg number to the Reynolds number (El = We/Re). In the equations above if β and α are zero, UCM constitutive relation can be obtained. If β is non-zero and α is zero Oldroyd B is reached. Lastly, if both β and α are nonzero the constitutive relation turns out to be Giesekus. With the aid of the axisymmetry, the problem is considered as 2.5 dimensional which means that the velocities in r, θ and z directions are solved for the flow domain shown in Fig. 1. The same approach is adopted by Pao (1970, 1972), Lopez (1996) and Saci et al. (1998) for the Newtonian problems. To ease up the solution stream function ( ψ ) vorticity ( ξ) circulation ( Γ) formulation is used. In addition, the viscoelastic stress tensor is splitted into Newtonian and non-newtonian parts as τ = S+2D. When this approach is used, the type of the vorticity equation is transformed to an elliptic one by which the convergence and the solution limit of the computations are enhanced (Crochet et al. (1991)). In the equation for splitting, S is the non-newtonian part of the stress tensor where D is the rate of deformation tensor. After some rearrangements, we need to solve nine equations for the rotating disc in a cylindrical enclosure flow (equations for ψ, ξ, Γ, S rr, S rθ, S rz, S θθ, S θz and S zz ). The full versions of the equations can be found in Appendix Boundary Conditions Zero stream function value is used for all of the edges. Vorticity boundary conditions are taken as; ξ = 2 ( ψ wall+1 ψ wall) / (dz 2 r) for the top and bottom discs, ξ= 2 ( ψ wall+1 ψ wall) / (dr 2 ) for the right wall and zero for the axis. In the vorticity boundary conditions, dr is the distance between the nodes in r-direction whereas dz is the distance in the z-direction in the similar way. The boundary condition for the circulation on the rotating disc is Γ = r 2 where all of the other edges have zero circulation value. Stress boundary conditions at the walls are applied in the way they are used in the study of Kawabata et al. (1990). It can be explained as equating the second derivatives of the non- Newtonian parts of the stresses to zero except the axis. At the symmetry axis, axially symmetric boundary condition is used for the stresses. 3. Solution Method The differential equations and the boundary conditions obtained at the previous section are discretized by using the finite difference method and as a result non-linear, sparse, coupled, large scale system (Equation 7) is obtained. Fx ( ) = 0 (7) 266 Korea-Australia Rheology Journal
3 Rotating-surface-driven non-newtonian flow in a cylindrical enclosure M 1 k Js k ( ) M 1 Fx ( k) k = x ( k+ 1) = x ( k+ 1) s ( k) ( ) (8) (9) Fx ( k) ( ) Js ( k) 2 η k Fx ( ( k) ) 2 (10) Here J is the Jacobian matrix, M k is the diagonal scaling matrix and s (k) is the linear step solution in the k th step and η is the inexact Newton parameter. For the solution of this system, Inexact Newton- GMRES(m) (Generalized Minimal Residual method with restart m by Saad et al. (1986)) algorithm is used (IN- GMRES(m)). In each step, Equation (7) is linearized by using Newton s method and the obtained system is scaled with respect to the diagonal term to enhance the convergence behavior. The result produced by this procedure is Equation (8). The system in Equation (8) is solved via GMRES(m) and the solution is updated as in Equation (9). This operation is carried out until the inequality in the Equation (10) is satisfied. Since the convergence of the Newton s method is strictly dependent on the initial guess, the continuation method is preferred. The continuation method can be described as using the previous solution as the initial guess. The reason for using GMRES(m) as the linear solver is its superiority over the other methods such as Gauss-Seidel and Successive Over Relaxation (SOR) in terms of its convergence rate, non-increasing residual property and easy application of the matrix-free versions. Since the computation of the Jacobian matrix is costly in terms of computation time and the matrix itself is sparse, the code is written as matrix-free. In the solution of the linear system, GMRES requires merely the products of Jacobian matrix and a vector instead of the Jacobian matrix itself. This operation is carried out by using the formula which is called Directional Differencing Discrete Newton. Fx ( + εv) Fx ( ) Jx ( )v = (11) ε In the Equation 11, ε is the perturbation parameter, the selection of which has an effect on the robustness of the nonlinear solver. This parameter can be selected either as a variable like in the studies of Brown et al. (1990) and Qin et al. (2000) or as a constant. In this study, it is chosen as a constant which is the square root of the machine epsilon. Upwinding is also used in the discretization processes. Additionally, in order to enhance the accuracy of the results, the fourth order difference stencils are used at the interior points. The convergence parameters for the linear and non-linear systems are taken as 1E Results Fig. 2. Streamline contours of the UCM flow with Re=0.32 and δ= Validation In this part of the study validation of the code is performed by comparing the results with the data available in the literature and the data obtained by using POLYFLOW which is commercial software designed for the flow of non-newtonian fluids. The study by Xue et al. (1999), in which UCM flows of rotating bottom disc in a cylindrical enclosure problem is selected as the first comparison case. In this problem, Reynolds number is fixed as 0.32 where aspect ratio is 1 and the constitutive relation is selected as UCM. Fig. 2 shows the streamline contours for different Weissenberg numbers. The streamline contours are in good agreement, therefore the contours of the study of Xue et al. (1999) are not given here. In Fig. 2, the positive values of the streamlines indicate the Newtonian like flow, in which the secondary flow in the r-z plane rotates radially outward from the disk (counter clockwise direction in the figure) whereas the negative values indicate the reversed flow pattern (clock- Korea-Australia Rheology Journal December 2010 Vol. 22, No
4 Yalin Kaptan, Ali Ecder and Kunt Atalik Table 1. Maximum weissenberg numbers computed for different grids Grid Max. We stresses counterbalance the centrifugal force and generate a reversed flow field. In the Fig. 2(b), this reversed flow appears near the right bottom corner of the flow domain. With the further increase in Weissenberg number, elastic forces become dominant and the reversed flow field grows into the domain where it finally covers the whole domain as illustrated in Fig. 2(e). Furthermore, Fig. 3 reveals the comparison of the velocity plots with the POLYFLOW for the same problem. The reversing of the flow can also be observed in these figures. Moreover the solutions for UCM, Oldroyd B and Giesekus models are validated with POLYFLOW up to its limits and the results are turned out to be comparable. Grid dependence is also checked with these analysis and grid is found out to be useful for δ=0.5 cases where grid turns out to be useful for δ=1 cases. Table 2 reveals the grid dependence of the maximum Weissenberg numbers for the Giesekus flow with the parameters of; Re=0, δ=0.5, β = 0.5 and α = 0.1. The maximum Weissenberg number limits are decreasing with the refinement of the grid as expected. Fig. 3. Velocity plots of the UCM flow with Re=0.32 and δ=1. wise direction in the figure). The secondary rotation in the r-z plane for Newtonian fluids is formed because of the centrifugal force applied by the disk. Increasing Weissenberg numbers generate growing normal stresses in the opposite direction of the centrifugal force. After a certain point, approximately We = for this problem, these 4.2. Effects of the constitutive relation, weissenberg number and reynolds number The constitutive relation is the most prominent part of the modeling of the non-newtonian phenomenon. In this study, three differential type constitutive relations, namely UCM, Oldroyd B and Giesekus relations are used. Table 2 points out the comparison of maximum Weissenberg number limits of the matrix-free IN-GMRES code generated in this study and POLYFLOW for different constitutive relations when δ = 0.5. In this problem, the identical boundary conditions and continuation parameters are selected both for IN-GMRES(m) and POLYFLOW codes. This table demonstrates that the solver generated in this study is more efficient than the steady solver in POLYFLOW and they both exhibit similar behavior of maximum Weissenberg number limits in terms of increase or decrease depending on the Reynolds number. The maximum Elasticity number obtained by IN- GMRES code developed in this study is for a nonzero Reynolds number and it is possible to achieve higher Elasticity numbers for different aspect ratios and Reynolds numbers as well. This value is proven to be high when compared with the previous studies of Moroi et al. (2001) (Elasticity number limit is 0.11), Itoh et al. (2006) (Elasticity number limit is 0.1) and Xue et al. (1999) (Elasticity 268 Korea-Australia Rheology Journal
5 Rotating-surface-driven non-newtonian flow in a cylindrical enclosure Table 2. Maximum weissenberg numbers computed for different constitutive models UCM Oldroyd B β = 0.5 Giesekus β = 0.5, α = 0.1 This Study Polyflow This Study Polyflow This Study Polyflow Re = Re = Re = Re = number limit is 0.072). What is more, it exceeds the Elasticity number limits of the steady solver of POLYFLOW. Although the limits of IN-GMRES code are beyond the experimental limits, it is significant to generate a solver with this capability since the application of another problem is always possible. Although UCM is one of the earliest and simplest constitutive relations, its numerical solution is the most challenging one of all. Thus, the maximum Weissenberg number limit is very low compared to the others. For Oldroyd B model, β = 0 value coincides with the UCM relation where the β = 1 value coincides with the Newtonian flow. Numerical solutions of Oldroyd B model become easier when viscosity ratio approaches 1 because the elliptic property of the Oldroyd B model increases. The Giesekus constitutive relation is one of the most realistic constitutive relations in the area of the non-newtonian flow modeling and it can model the shear thinning behavior of non-newtonian fluids which is useful in the modeling of rotating disk in a cylindrical enclosure problem. Beneficially, the maximum Weissenberg number limits are higher for Giesekus constitutive relation compared to the UCM and Oldroyd B models. Since the Giesekus model is the most advantageous one compared with the other two, the simulations carried out in the following sections are performed by using the Giesekus constitutive relation. It is not advisable to compare the velocity fields of three constitutive relations mentioned above for a fixed Weissenberg number because the different viscosity ratios and mobility factors will result in dissimilar elastic effects and consequently the flow fields will not be similar. Fig. 4 reveals the streamline contours of non-newtonian flow computed by using Giesekus constitutive relation with the parameters β = 0.5, α = 0.1, δ = 0.5 and Re = 1. It can be observed from this figure that the elastic forces counterbalance the centrifugal force and the flow transformed from Newtonian like to the fully reversed flow with We = 0.5 case. After this point, the Weissenberg number is further increased and at nearly We = 25, a new flow spot which rotates in the direction of the Newtonian flow Fig. 4. Streamline contours of the Giesekus flow with Re = 1, δ = 0.5, α = 0.1 and δ = 0.5. (counterclockwise) appears. This re-formation means that the elastic forces now can not overcome the centrifugal force arising from the rotation of the disk. Continuous rise in Weissenberg number enlarges this new spot where it covers nearly 75% of the flow domain at We = 130. This phenomenon can be considered as a high Weissenberg number occurrence. The results are validated with POLY- FLOW up to We = 90 and they are in good agreement. Since the re-formation of the Newtonian like spot occurs nearly at We = 25 (El = 25) and the Elasticity limit in the Korea-Australia Rheology Journal December 2010 Vol. 22, No
6 Yalin Kaptan, Ali Ecder and Kunt Atalik Fig. 6. Streamline contours of the Giesekus flow with Re = 5, β = 0.5, α = 0.1 and δ = 0.5. Fig. 5. Velocity plots of the Giesekus flow with Re = 1, β = 0.5, α = 0.1 and δ = 0.5. literature is nearly 0.1, there is no result for the re-formation of the Newtonian flow spot in the literature for the authors knowledge. Although this limit is beyond the experimental studies, it is important to solve this type of problems by use of robust numerical methods such as matrix-free Inexact Newton-Krylov methods used in this study. Fig. 5 reveals the comparison of the velocity plots with the POLYFLOW for the same problem. The re-formation of the Newtonian like flow can also be observed from this figure. Fig. 6 demonstrates the streamline contours of non-newtonian flow computed by using Giesekus constitutive relation with the parameters β = 0.5, α = 0.1, δ = 0.5 and Re = 5. It can be observed from this figure that the reversed flow spot is formed at the right bottom corner of the flow domain while this time it can not cover the whole domain by growing into it. It gets bigger up to We = 2 and from this point on, it gets smaller and smaller again. The Newtonian like flow which is rotating in the counterclockwise direction covers nearly 95 % of the flow region at We = 20. The reversed flow can not cover the whole flow domain in this problem since the Reynolds number is increased to 5. Given the fact that the Reynolds number is the ratio of the inertial forces to the viscous forces, the increase in this number results in enhanced centrifugal forces which can not be overcome by the elastic forces. However the final shapes of the flow domains are similar for the problems in Fig. 4 and Fig. 6. The velocity plots of the flow of 270 Korea-Australia Rheology Journal
7 Rotating-surface-driven non-newtonian flow in a cylindrical enclosure depending on the Weissenberg and Reynolds numbers (up to El = 0.1) can be found in the study of Moroi et al. (2001). 5. Conclusions High Weissenberg number non-newtonian flows of rotating surface driven cylindrical container problem is investigated by generating a robust and efficient matrixfree Inexact Newton-Krylov solver (IN-GMRES). UCM, Oldroyd B and Giesekus constitutive relations are used to model the non-newtonian behavior and the results obtained are compared with the data obtained by POLY- FLOW. Additionally, the selection of the constitutive relation, the effects of the Weissenberg number and effects of the Reynolds number are studied. The simulations point out that the generated algorithm is capable of solving high Weissenberg number problems (up to the Elasticity number limit of 130) when compared to previous studies. The development of the reversed flow spot due to the counterbalancing of the centrifugal forces by elastic forces (with the increasing Weissenberg number) is observed for all of the simulated cases. Additionally, with the further increase in Weissenberg number re-formation of the Newtonian like flow is monitored. Furthermore, it is shown that the increase in the Reynolds number alters the centrifugal forces and changes the flow structure. Acknowledgements This work was funded by the Scientific and Technological Research Council of Turkey (Project: 107M390). Notation Fig. 7. Velocity plots of the Giesekus flow with Re = 5, β = 0.5, α = 0.1 and δ = 0.5. Giesekus constitutive relation with the parameters β = 0.5, α = 0.1, δ = 0.5 and Re = 5 can be seen in Fig. 7, in which the flow formations can also be examined. The comparison of Figs. 4 and 6 reveals that the change in Reynolds number effects the flow field by altering the centrifugal forces. The parametric study of the formation of the flow D Non-Newtonian part of the stress tensor El Elasticity number H Height J Jacobian matrix M Diagonal scaling matrix r Distance in the radial direction R Radius Re Reynolds number S Non-Newtonian part of the stress tensor s Linear step solution vector u u component of the velocity v v component of the velocity We Weissenberg number w w component of the velocity α Mobility factor β Viscosity ratio Γ Circulation δ Aspect ratio ε Matrix-free perturbation parameter ξ Vorticity Korea-Australia Rheology Journal December 2010 Vol. 22, No
8 Yalin Kaptan, Ali Ecder and Kunt Atalik η Inexact Newton parameter λ Relaxation time τ Stress Tensor τ Stream function ψ Angular velocity References Bird, R.B., R.C. Armstrong and O. Hassager, 1987, Dynamics of polymeric liquids, Vol. 1: Fluid mechanics, 2 nd ed., Wiley, New York. Brown, P.N. and Y. Saad, 1990, Hybrid Krylov methods for nonlinear-systems of equations, SIAM Journal on Scientific and Statistical Computing, 11, Crochet, M.J., A.R. Davies and K. Walters, 1991, Numerical simulation of non-newtonian flow, Elsevier, New York. Itoh, M., M. Suzuki and T. Moroi, 2006, Swirling flow of a viscoelastic fluid in a cylindrical casing, Trans. ASME, 128, Kawabata, N., T. Motoyoshi and A. Isao, 1990, A numerical simulation of viscoelastic fluid flow in a two-dimensional channel: application of Lax's Scheme to the constitutive equation [in Japanese], Transactions of the Japan Society of Mechanical Engineers B, 56, Lopez, J.M., 1996, Flow between a stationary and a rotating disk shrouded by a co-rotating cylinder, Phys. Fluids, 8(10), Moroi, T., M. Itoh, K. Fujita and H. Hamasaki, 2001, Viscoelastic flow due to a rotating disk enclosed in a cylindrical casing (Influence of Aspect Ratio), JSME International Series B, 44(3), Pao, H.P., 1970, A numerical computation of a confined rotating flow, Trans. ASME, J. Applied Mechanics, 37, Pao, H.P., 1972, Numerical solution of the navier-stokes equations for flows in the disk-cylinder system, Physics of Fluids, 15, 1, Qin, N., D.K. Ludlow and S.T. Shaw, 2000, A matrix-free preconditioned Newton/GMRES method for unsteady Navier Stokes equations, Int. J. Numer. Meth. Fl., 33, Saad, Y. and M.H. Schultz, 1986, GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM Journal on Scientific and Statistical Computing, 7(3), Saci, R. and P.G. Bellamy-Knights, 1998, Diffusion driven rotating flow in a cylindrical container, Acta Mechanica, 126, Xue, S.C., N. Phan-Thien and R.I. Tanner, 1999, Fully threedimensional, time-dependent numerical simulations of Newtonian and viscoelastic swirling flows in a confined cylinder Part I. Method and steady flows, J. Non-Newtonian Fluid Mech., 87, Korea-Australia Rheology Journal
Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationOutline. Motivation Governing equations and numerical methods Results: Discussion:
Bifurcation phenomena in strong extensional flows (in a cross-slot geometry) F. A. Cruz 1,*, R. J. Poole 2, F. T. Pinho 3, P.J. Oliveira 4, M. A. Alves 1 1 Departamento de Engenharia Química, CEFT, Faculdade
More informationNumerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationA new numerical framework to simulate viscoelastic free-surface flows with the finitevolume
Journal of Physics: Conference Series PAPER OPEN ACCESS A new numerical framework to simulate viscoelastic free-surface flows with the finitevolume method Related content - Gravitational collapse and topology
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationRAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8
HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationFLOW CHARACTERISTICS IN A VOLUTE-TYPE CENTRIFUGAL PUMP USING LARGE EDDY SIMULATION
FLOW CHARACTERISTICS IN A VOLUTE-TYPE CENTRIFUGAL PUMP USING LARGE EDDY SIMULATION Beomjun Kye Keuntae Park Department of Mechanical & Aerospace Engineering Department of Mechanical & Aerospace Engineering
More informationMonolithic FEM multigrid techniques for the simulation of viscoelastic flow
Monolithic FEM multigrid techniques for the simulation of viscoelastic flow A. Ouazzi, H. Damanik, S. Turek, J. Hron Institute of Applied Mathematics, LS III, TU Dortmund http://www.featflow.de European
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationDirect Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids
Direct Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids by P. Y. Huang 1, J. Feng 2, H. H. Hu 3 and D. D. Joseph 1 1 Department of Aerospace Engineering
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationApplication Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus
Application Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus T.S.L Radhika, Aditya Vikram Singh Abstract In this paper, the flow of an incompressible non Newtonian
More informationQuasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section
Article Nihon Reoroji Gakkaishi Vol.34, No.2, 105~113 (Journal of the Society of Rheology, Jaan) 2006 The Society of Rheology, Jaan Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationViscoelastic Flow in Rotating Curved Pipes
Mechanical Engineering Faculty Publications Mechanical Engineering 8-2006 Viscoelastic Flow in Rotating Curved Pipes Yitung Chen University of Nevada, Las Vegas, yitung.chen@unlv.edu Huajun Chen University
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationEFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER
ISTP-1, 5, PRAGUE 1 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA EFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER S. S. Leong School of Mechanical and Manufacturing Engineering
More informationStudy on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time
World Journal of Mechanics, 2018, 8, 301-310 http://www.scirp.org/journal/wjm ISSN Online: 2160-0503 ISSN Print: 2160-049X Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration
More informationNumerical description of start-up viscoelastic plane Poiseuille flow
Korea-Australia Rheology Journal Vol. 1, No. 1, March 009 pp. 47-58 Numerical description of start-up viscoelastic plane Poiseuille flow Kwang Sun Park and Young don Kwon* Department of Chemical Engineering,
More informationHEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS
HEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS Paulo J. Oliveira Departamento de Engenharia Electromecânica, Universidade da Beira Interior Rua Marquês D'Ávila e Bolama, 600
More informationA Numerical Study of Several Viscoelastic Fluid Models
A Numerical Study of Several Viscoelastic Fluid Models Corina Putinar January 3, 6 Abstract Viscoelastic fluids are a type of fluid with a solvent and immersed elastic filaments which create additional
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationLaminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid
Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering
More informationMagnetohydrodynamic Flow of a Liquid Metal in a Curved Circular Duct subject to the Effect of an External Magnetic Field
Paper 85 Civil-Comp Press, 2012 Proceedings of the Eighth International Conference on Engineering Computational Technology, B.H.V. Topping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Magnetohydrodynamic
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationSimulation of Die-swell Flow for Oldroyd-B Model with Feedback Semi-implicit Taylor Galerkin Finite Element Method
KMUTNB Int J Appl Sci Technol, Vol.8, No.1, pp. 55-63, (2015) Simulation of Die-swell Flow for Oldroyd-B Model with Feedback Semi-implicit Taylor Galerkin Finite Element Method Nawalax Thongjub and Vimolrat
More informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
More informationLecture 2: Constitutive Relations
Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt
More information[Supplementary Figures]
[Supplementary Figures] Supplementary Figure 1 Fabrication of epoxy microchannels. (a) PDMS replica is generated from SU-8 master via soft lithography. (b) PDMS master is peeled away from PDMS replica
More informationExperiments at the University of Minnesota (draft 2)
Experiments at the University of Minnesota (draft 2) September 17, 2001 Studies of migration and lift and of the orientation of particles in shear flows Experiments to determine positions of spherical
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationTWO-DIMENSIONAL MAGMA FLOW *
Iranian Journal of Science & Technology, Transaction A, Vol. 34, No. A2 Printed in the Islamic Republic of Iran, 2010 Shiraz University TWO-DIMENSIONAL MAGMA FLOW * A. MEHMOOD 1** AND A. ALI 2 1 Department
More informationTHE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES
THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES Kazuya Yokomizo 1, Makoto Iwamura 2 and Hideki Tomiyama 1 1 The Japan Steel Works, LTD., Hiroshima Research Laboratory,
More informationOn the exact solution of incompressible viscous flows with variable viscosity
Advances in Fluid Mechanics IX 481 On the exact solution of incompressible viscous flows with variable viscosity A. Fatsis 1, J. Statharas, A. Panoutsopoulou 3 & N. Vlachakis 1 1 Technological University
More informationAxisymmetric Hopf bifurcation in a free surface rotating cylinder flow
ANZIAM J. 50 (CTAC2008) pp.c251 C265, 2008 C251 Axisymmetric Hopf bifurcation in a free surface rotating cylinder flow S. J. Cogan 1 G. J. Sheard 2 K. Ryan 3 (Received 13 August 2008; revised 24 October
More informationPeristaltic Flow of Non-Newtonian Fluids through Curved Channels: a Numerical Study
ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 1, 13 Peristaltic Flow of Non-Newtonian Fluids through Curved Channels: a Numerical Study Alireza Kalantari 1, Kayvan Sadeghy 1, and Soheil Sadeqi
More informationNumerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries
Korea-Australia Rheology Journal Vol. 19, No. 4, December 2007 pp. 211-219 Numerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries Youngdon Kwon*
More informationTURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS
HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,
More informationvery elongated in streamwise direction
Linear analyses: Input-output vs. Stability Has much of the phenomenology of boundary la ut analysis of Linearized Navier-Stokes (Cont.) AMPLIFICATION: 1 y. z x 1 STABILITY: Transitional v = (bypass) T
More informationSimulation of Variable Viscosity and Jeffrey Fluid Model for Blood Flow Through a Tapered Artery with a Stenosis
Commun. Theor. Phys. 57 (2012) 133 140 Vol. 57 No. 1 January 15 2012 Simulation of Variable Viscosity and Jeffrey Fluid Model for Blood Flow Through a Tapered Artery with a Stenosis Noreen Sher Akbar 1
More informationON VARIABLE LAMINAR CONVECTIVE FLOW PROPERTIES DUE TO A POROUS ROTATING DISK IN A MAGNETIC FIELD
ON VARIABLE LAMINAR CONVECTIVE FLOW PROPERTIES DUE TO A POROUS ROTATING DISK IN A MAGNETIC FIELD EMMANUEL OSALUSI, PRECIOUS SIBANDA School of Mathematics, University of KwaZulu-Natal Private Bag X0, Scottsville
More informationCONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 6-8 December 1999 CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationStatistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders
Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Toshiki Morita 1, Takashi Watanabe 2 and Yorinobu Toya 3 1. Graduate School of Information
More informationFlow analysis in centrifugal compressor vaneless diffusers
348 Journal of Scientific & Industrial Research J SCI IND RES VOL 67 MAY 2008 Vol. 67, May 2008, pp. 348-354 Flow analysis in centrifugal compressor vaneless diffusers Ozturk Tatar, Adnan Ozturk and Ali
More informationUnsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation
Unsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation Mohamed S. Ebeida Carnegie Mellon University, Pittsburgh, PA 15213 Roger L. Davis and Roland W. Freund
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More information2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationPerformance characteristics of turbo blower in a refuse collecting system according to operation conditions
Journal of Mechanical Science and Technology 22 (2008) 1896~1901 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-008-0729-6 Performance characteristics
More informationReverse flow and vortex breakdown in a shearthinning
Journal of Physics: Conference Series Reverse flow and vortex breakdown in a shearthinning fluid To cite this article: C Cabeza et al 2011 J. Phys.: Conf. Ser. 296 012020 View the article online for updates
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationExact Solutions for the Incompressible Viscous Fluid of a Rotating Disk Flow
Progress in Applied Mathematics Vol. 1, No. 1, 211, pp. 9-97 www.cscanada.org ISSN 1925-251X [Print] ISSN 1925-2528 [Online] www.cscanada.net Exact Solutions for the Incompressible Viscous Fluid of a Rotating
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationLaminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
More informationCOMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
Suranaree J. Sci. Technol. Vol. 20 No. 4; October - December 2013 257 COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationResearch Article Study of the Transient Natural Convection of a Newtonian Fluid inside an Enclosure Delimited by Portions of Cylinders
Research Journal of Applied Sciences, Engineering and Technology 7(5): 3069-3074, 04 DOI:0.906/rjaset.7.644 ISSN: 040-7459; e-issn: 040-7467 04 Maxwell Scientific Publication Corp. Submitted: August 7,
More informationSTUDY OF THE TEMPERATURE DISTRIBUTION IN DISC BRAKES BY THE METHOD OF ORDER-OF-MAGNITUDE ANALYSIS
Regional Conference on Engineering Mathematics, Mechanics, Manufacturing & Architecture (EM 3 ARC) 007 007 Mathematical Sciences in Engineering Editors: A. Zaharim et al. STUDY OF THE TEMPERATURE DISTRIBUTION
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationInertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes
J. Non-Newtonian Fluid Mech., 78 (1998) 27 45 Inertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes Yuriko Renardy *, David O. Olagunju 1 Department of Mathematics and ICAM,
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationNUMERICAL STUDY OF THE EFFECT OF BLADE SZE ON PUMPING EFFECTIVENESS OF A PADDLE IMPELLER IN AN UNBAFFLED MIXING VESSEL
Third International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2003 NUMERICAL STUDY OF THE EFFECT OF BLADE SZE ON PUMPING EFFECTIVENESS OF A PADDLE
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationNumerical Simulation of the Evolution of Reynolds Number on Laminar Flow in a Rotating Pipe
American Journal of Fluid Dynamics 2014, 4(3): 79-90 DOI: 10.5923/j.ajfd.20140403.01 Numerical Simulation of the Evolution of Reynolds Number on Laminar Flow in a Rotating Pipe A. O. Ojo, K. M. Odunfa,
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
More informationNumerical study of vortex shedding in viscoelastic flow past an unconfined square cylinder
Korea-Australia Rheology Journal, 27(3), 213-225 (August 2015) DOI: 10.1007/s13367-015-0022-z www.springer.com/13367 Numerical study of vortex shedding in viscoelastic flow past an unconfined square cylinder
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationLATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY
LATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY M. Y. Gokhale, Ignatius Fernandes Maharashtra Institute of Technology, Pune 4 38, India University of Pune, India Email : mukundyg@yahoo.co.in,
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY. David O. Olagunju. University of Delaware. Newark, DE 19716
ELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY David O. Olagunju Department of Mathematical Sciences University of Delaware Newark, DE 19716 June 15, 1995 Abstract Consider the axisymmetric,
More informationRheometry. II.1 Introduction
II Rheometry II.1 Introduction Structured materials are generally composed of microstructures dispersed in a homogeneous phase [30]. These materials usually have a yield stress, i.e. a threshold stress
More informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More informationNumerical Prediction Of Torque On Guide Vanes In A Reversible Pump-Turbine
Journal of Multidisciplinary Engineering Science and Technology (JMEST) ISSN: 3159 Vol. 2 Issue 6, June - 215 Numerical Prediction Of Torque On Guide Vanes In A Reversible Pump-Turbine Turbine and pump
More informationImplementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows
Short Training Program Report Implementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows Julian Koellermeier RWTH Aachen University Supervisor: Advisor: Prof. Thierry Magin von
More informationIndex. Boundary integral method, 27 Boundary location method, 255 Brinkman number, 14, 158
Index ABFIND,275 Adaptive control, 304, 312 Adaptive process model, 315 Air drag, 221 Alternating Direction Implicit (ADI) technique, 30 Axisymmetric flow, 48 die entry, 58 Boundary integral method, 27
More information5.1 Fluid momentum equation Hydrostatics Archimedes theorem The vorticity equation... 42
Chapter 5 Euler s equation Contents 5.1 Fluid momentum equation........................ 39 5. Hydrostatics................................ 40 5.3 Archimedes theorem........................... 41 5.4 The
More informationInvestigation of reattachment length for a turbulent flow over a backward facing step for different step angle
MultiCraft International Journal of Engineering, Science and Technology Vol. 3, No. 2, 2011, pp. 84-88 INTERNATIONAL JOURNAL OF ENGINEERING, SCIENCE AND TECHNOLOGY www.ijest-ng.com 2011 MultiCraft Limited.
More information