Modified DLM method for finite-volume simulation of particle flow
|
|
- Joel Parker
- 6 years ago
- Views:
Transcription
1 Modified DLM method for finite-volume simulation of particle flow A. M. Ardekani, S. Dabiri, and R. H. Rangel Department of Mechanical and Aerospace Engineering, University of California, Irvine, CA , USA A Distributed-Lagrange-Multiplier(DLM)-based method is implemented for the simulation of particulate flow. Initially, we show that our fluid-particle solver produces results which are in good agreement with numerical studies of benchmark viscous flow problems. Subsequently, the bouncing motion of a solid sphere onto a solid plate in an ambient fluid is considered. Comparing results for the coefficient of restitution for different Stokes numbers shows very good agreement with the experimental results. Finally, the effect of the coefficient of restitution for a dry collision on the vortex dynamics associated with a colliding sphere onto a wall is discussed. I. Introduction The interaction between fluid and particles in a dispersed flow occurs in several natural and industrial applications such as sedimentation, combustion systems, crystal growth, and microfluidic devices. For example, the design and utilization of microfluidic devices for biomedical purposes requires a better understanding of particulate flow in low and intermediate Reynolds number flows. Particulate flow simulations under the Stokes-flow assumption have been conducted using Stokesiandynamics techniques. 1 For dilute suspensions, the unsteady motion of two solid spheres has also been analyzed. 2 However, for intermediate Reynolds numbers the use of numerical simulation is generally unavoidable. The flow around a single object can be calculated using a fixed body-fitted grid in a non-inertial reference frame. For multiple objects moving relative to each other, remeshing is required along the moving interfaces. The flow field is usually solved using finite-element 3 or boundary-element 4 methods. An alternative is to use a fixed cartesian grid which has the advantage of permitting the use of fast solvers. Calhoun 5 and separately Russell et al. 6 utilized a cartesian grid using a stream function-vorticity formulation in order to simulate irregular shapes and multiple objects. The non boundary-fitted methods are usually simpler and more efficient and have been used by several researchers Another promising method developed in the last decade is a Distributed-Lagrange-Multiplier (DLM) technique using a fictitious domain as presented by Glowinski et al., 9 who used a finite-element method with a fixed structured grid thus eliminating the need of remeshing required for unstructured boundary-fitted grids. In this method, the entire domain is treated as fluid but the fluid inside the particle domain satisfies a rigidity constraint by using Lagrange multipliers. Patankar et al. 11 formulated their DLM solution by forcing the deformation tensor in the particle domain to be zero thus eliminating U and ω as the variables from the coupled system of equations, where U and ω are the translational and angular velocity of the particle, respectively. This formulation introduces a stress field in the particle domain similar to the pressure in an incompressible fluid. Sharma et al. 10 presented the formulation of DLM for steady Stokes flow using a control-volume method. Rebound of colliding particles has also been studied in recent years. Davis et al. 12 employed an elastohydrodynamic approach and showed that the pertinent parameter for collision in the fluid is not the Reynolds number Re but the Stokes number St = 1 ρ p 9 ρ f Re where ρ p and ρ f are the particle and fluid densities, respectively. No rebound occurs for St lower than a critical value due to the fact that elastic energy stored by the particle deformation is dissipated in the fluid. The experiments by Joseph et al. 13 show that the Graduate Student Researcher, AIAA student member Professor, AIAA senior member 1 of 13
2 rebound velocity depends strongly on the impact Stokes number and weakly on the elastic properties of the material, where the Stokes number is defined using the approaching velocity of the particles. It has been shown that below a St of 10, no rebound occurs For impact Stokes numbers larger than 500 the coefficient of restitution asymptotes to that for a dry collision. Whereas several experimental studies have been conducted on the influence of the fluid on the collision process, few numerical studies address this issue. In this work, we numerically study the effects of Stokes number on the rebound velocity and obtain results which are in agreement with the experimental results by Gondret et al. 15 We first implement DLM for fixed cylindrical particles in low to intermediate Reynolds number flow. The constraint of zero velocity satisfies the condition of rigid body motion inside the particle. Thus, initially it is enough to set the velocity field inside the particle equal to zero and this can be done by means of a rigidity force field in the particle domain. Even though the advantage of using DLM resides in the use of a fixed grid for moving boundaries, utilizing DLM for the flow around stationary objects provides the opportunity of experimenting with irregular geometries simply by using a rectilinear grid. Furthermore, modification of the flow solver due to the presence of the particle is simple. Subsequently, DLM is implemented for the case of moving particles. The bouncing motion of a solid sphere onto a wall is numerically simulated and the vortex dynamics associated with this problem is discussed. Following the derivation procedure of Patankar et al., 10, 16 we derive the governing equations in Section II. The numerical implication will be described in Section III and the method is verified by applying it to different fluid-particle interaction problems in Sections IV and V. II. Theoretical Development In this section, the DLM method is described for the motion of one particle but the method can be easily extended for additional particles. Let Γ represent the fluid boundary which is not shared with the particle. The solid domain and its boundary are denoted by P and P, respectively. The computational domain is Ω, including both the fluid and the particle. The governing equations in the fluid domain are ρ f Du Dt = σ + ρ f g in Ω \ P (1) u = 0 in Ω \ P (2) u = U i on P (3) σ n = t on P (4) u t=0 = u 0 (x) in Ω \ P (5) in addition to the outer boundary conditions on Γ. In these equations, u is the fluid velocity, g is the acceleration of gravity, ρ f is the fluid density, n is the normal unit vector, U i is the particle velocity, and t is the traction vector on the particle surface. The initial velocity u 0 satisfies the continuity equation, σ = pi + τ is the stress tensor, p is the pressure field, I is the identity tensor, and τ is the viscous stress: The governing equations in the particle domain are τ = 2µD[u] = µ [ u + ( u) T ]. (6) ρ P Du Dt = σ + ρ Pg in P (7) u = 0 in P (8) D[u] = 0 in P (9) u = U i on P (10) σ n = t on P (11) u t=0 = u 0 (x) in P. (12) where ρ P is the particle density. Equation (9) satisfies the continuity equation but in order to facilitate numerical implementation, equation (8) is retained. As pointed out by Sharma et al., 10 the rigidity constraint gives rise to a stress field inside the particle which is a function of three scalar Lagrange multipliers for three dimensional problems. Thus for a two dimensional problem, the stress field can be represented in terms of 2 of 13
3 two scalars: σ = pi + τ + D[λ] (13) where λ represents the Lagrange multipliers and τ is zero inside the particle due to the rigidity constraint. The governing equations in the entire domain can be combined as: ρ Du = σ + ρg + f in Ω (14) Dt u = 0 in Ω (15) u t=0 = u 0 (x) in Ω (16) where ρ = ρ f ρ P in Ω \ P in P (17) in addition to the outer boundary conditions on Γ. In the above equations, f = D[λ] is zero everywhere except in the particle domain and leads to the rigid body motion inside the particle. III. Numerical Implementation A finite-volume method using a staggered grid for incompressible flow is implemented. The SIMPLE algorithm 17 is used to solve the fluid equations with modifications to account for the presence of particles. The Crank-Nicolson scheme is used for time discretization. The discretised momentum equations in the x and y directions are a i,j u i,j = a nb u nb p I,J p I 1,J δx u V u + b i,j + F x (18) a I,j v I,j = a nb v nb p I,J p I,J 1 δx v V v + b I,j + F y (19) Where V u and V v are the volumes of the u-cell and v-cell, respectively; u and v represent the horizontal and vertical components of the velocity field; I and J are the nodes at the center of the main control volumes; while i and j represent the nodes at the center of the control volumes for u and v, respectively as shown in Fig. 1. Subscript nb refers to corresponding neighboring staggered control volumes; b is the momentum source term which includes the gravity term; and F represents the rigidity force which makes the velocity field inside the particle domain satisfy the rigidity constraint. v I-1,j+2 v I,j+2 u i-1,j+1 u i,j+1 u i+1,j+1 v I-1,j+1 v I,j+1 u i-1,j u i,j u i+1,j P I,J v I-1,j v I,j Figure 1: Staggered grid, u control volume and corresponding volume fraction 3 of 13
4 R=0.1 6 Figure 2: Relative position of cylinder and boundaries. Following the SIMPLE algorithm and taking the particles into account, the force F is added as an unknown variable and F = F + F where F is the force predicted at each iteration (the force calculated in the previous iteration) and F is the correction force. At each iteration, the rigidity force is defined as and similarly in the y direction F xi,j = F xi,j + Cρ P A i,j (u i,j u Ri,J ) (20) F yi,j = F yi,j + Cρ P A I,j (v I,j v RI,j ) (21) where C is a dimensional constant which includes an under-relaxation factor; A i,j and A I,j are the area of u-cell and v-cell, respectively; while u R is the velocity vector rigidified inside the particle, equal to u outside the particle domain, and defined as follows u R = (1 φ)u + φu P (22) where φ is the volume fraction occupied by the particle in each control volume, defined separately for u and v-cells. For example, as shown in Fig. 1, φ u = V hatched V u. In the present study, φ is exactly calculated while u P, defined only in the particle domain, can be calculated as follows: u P = U P + ω r (23) where U P and ω are the particle translational and angular velocities. By using conservation of linear and angular momentum for the solid particle, one can calculate the particle translational and angular velocities 10, 16 as follows: M P U P = ρudx and I P ω = r ρudx (24) P Defining a rigidity force as in equations (20) and (21) guarantees that, upon convergence at each time step, u = u R everywhere in the domain. Finally, the modifications due to the presence of the particle can be included by adding F as a source term in the momentum equations. P IV. Results and Discussion A. Flow past a stationary circular cylinder As in the first case, the flow past a stationary circular cylinder is modeled. This case is designed to test the solution process without the added complication of the particle motion. Figure 2 shows the geometry used for this case. The boundary conditions are 4 of 13
5 v u y = 0, y = 0 on the top and bottom boundaries Du Dt u = 0, x = 0 on the right (exit flow) boundary u = U in, v = 0 on the left (inlet flow) boundary For fixed objects we have u p = 0 and u R = u(1 φ) (25) The drag force can be determined by calculating the summation of F in the particle domain. Table 1: Summary of results for 1 < Re < 40 C D (L) Re Takami& Keller (0.25) 2.00 (0.93) 1.72 (1.61) 1.54 (2.32) Dennis& Chang (0.26) 2.05 (0.94) 1.52 (2.35) Fornberg (0.91) 1.50 (2.24) Calhoun (0.91) 1.62 (2.18) present work (0.24) 2.08 (0.90) 1.75 (1.57) 1.55 (2.23) Table 1 shows a summary of results for different Reynolds numbers where C D is the drag coefficient and L represents the ratio of circulation length to the particle diameter. The results compare very favorably with both experimental and previous computational results. For the data presented in Table 1, the mesh size is 1 11 of the particle radius. The geometry used for Re = 1 and Re = 4 is different from the one shown in figure 2. Since the upstream flow is affected by the presence of the cylinder as much as downstream flow, we double the width of the flow domain and the center of the cylinder is moved towards the center of the computational domain C D C L C T C D,C L,C T Time Figure 3: Time dependent lift, drag, and torque coefficient at Re = Unsteady vortex shedding: Re = 100 For higher Reynolds number, instability occurs and the classic oscillatory wake behind the cylinder can be easily captured. Because our algorithm has a slight built-in asymmetry in its details, there is no need to artificially perturb the flow field to initiate the unsteady behavior. In figure 3, we plot lift, drag, and torque coefficients, C L, C D, C T versus nondimensional time tu in D, respectively, at Re = 100 where D is the diameter of particle. The forces which impose zero velocity field in the particle domain are equal to the reaction forces applied on the cylinder to fix it. Therefore, the lift and drag coefficients can be calculated easily just by 5 of 13
6 calculating the summation of F y and F x in the particle domain, respectively. The summary of results is shown in table 2. Generally, our results are well within the range of results reported by other researchers. Table 2: Summary of results for Re = 100 C D C L St Braza et al ± ±0.25 Liu et al ± ± Calhoun ± ± Russell et al ± ± present work 1.36 ± ± (a) Streamlines (b) Vorticity contours Figure 4: Streamlines and vorticity contours for flow past a biperiodic array. Re= 5 and volume fraction is 0.2. Streamline contour levels are 0:0.1:1 and 0.45:0.01:0.55. Vorticity contour levels are -31:2:31. B. Flow past a biperiodic array Flow past an infinite periodic array of cylinders is considered in this section. The purpose of this example is to show that our method works well for low-reynolds-number flows. Flow past a cylinder array in a unit square with periodic boundary conditions in both the horizontal and the vertical directions is considered. The following boundary conditions are utilized: u y = 0, v = 0 on the top and bottom boundaries Periodic boundary conditions on the downstream and upstream velocities A uniform grid is used for Re = 5 and a volume fraction of 0.2. Reynolds number is based on particle diameter and fluid average upstream velocity (U 0 ). Streamlines and vorticity contours are shown in figure 4. The solution was computed on one cylinder and plotted periodically. Comparing our drag coefficient F x µu 0 with the results by Koch and Ladd 23 shows only a 1.2% difference. 6 of 13
7 V. Numerical Implementation For Moving Objects A. Sedimentation of a spherical particle towards a wall In this section, we examine the motion of a spherical particle moving towards a wall with a non-zero dry coefficient of restitution. All dimensional quantities are in the CGS units unless otherwise stated. Specifically, the collision of a steel sphere with density of 7.8 colliding with a glass wall is investigated. The dry coefficient of restitution for this system is equal to 0.97 ± 0.02 according to the data reported by Gondret et al. 15 The sphere is moving in oil with density of and viscosity of 0.1. The sphere diameter is 0.3 and the gravitational acceleration is g = 981. A fixed nonuniform structured grid is employed. The smallest mesh size is near the contact region where higher resolution is needed. The numerical simulation is performed utilizing axisymmetry. The distance h from the bottom of the sphere to the wall and the particle velocity are shown in figure 5. The Reynolds number based on the approach velocity of the sphere towards (a) Distance between particle surface and the wall (b) Vertical velocity Figure 5: Sedimentation of a spherical particle and its collision with the wall. Re = 63, St = 58, h min = 1.2µm, and e dc = the wall is 62.7 and the impact Stokes number is 58. The effective roughness height is h min = µm. As it can be seen, rebound trajectories are non-parabolic since the velocity decreases nonlinearly with time. A jump in velocity occurs when the surface of the particle and the wall come into contact. Subsequently, the particle moves away from the wall and a marked decrease in particle velocity is observed. As the sphere moves away from the wall and the gap widens, the velocity varies more slowly. Figure 6 compares the results using our numerical method with the experimental results by Gondret et al. For larger Stokes number, higher rebound is observed, as expected. It can be seen that the coefficient of restitution we obtain numerically is in very good agreement with the experimental measurements. More detailed results can be found in Ardekani and Rangel. 24 B. Vortex dynamics of a colliding sphere onto a wall When a sphere moves towards a wall and stops upon making contact with it, a secondary vortex ring is generated from the wake vorticity. The evolution of vorticity for this problem is considered numerically by Thompson et al. 25 and experimentally by Eames and Dalziel. 26 In both these studies, the sphere sticks to the wall after the collision. In the present work, we discuss the vortex dynamics associated with the collision of a sphere and a wall for non zero coefficient of restitution. First, we compare our numerical results with experimental results by Eames and Dalziel. 26 In this case, a sphere is set into motion from a distance of L = 7.5D above the wall and moves with constant velocity and a Reynolds number of 850. A nonuniform structured grid and a domain size of 10D 20D with nodes are employed. Figure 7 shows the vorticity contours when the sphere sticks to the wall after collision. 7 of 13
8 e/e dc St Figure 6: Coefficient of restitution normalized by that for dry collision as a function of St. Present results ( ) where h min = 0.7µm. Experimental measurement for different materials by Gondret et al. 15 tungsten carbide ( ), steel ( ), glass ( ), Teflon ( ), Derlin ( ), polyurethane ( ), and Nylon ( ). Roughness in these experimental cases is less than 1µm. Lubrication theory of Davis et al. 12 ( ). In this case, the sphere is moving with constant velocity before the collision. The left and right sides of each frame in the figure show the present numerical results and the experimental results by Eames and Dalziel, 26 respectively. This comparison shows very good agreement. In frames (a,b) of figure 7 the wake vorticity of the sphere approaching the wall can be observed. As the sphere stops, the wake vorticity moves towards the wall due to its inertia. The wake vortex threads over the sphere and generates a secondary vortex ring as seen in frame (c). This coherent structure, composed of the wake and secondary vortices, leaves the sphere and strikes the wall as seen in frame (d). The secondary vortex is stretched by the wake vortex and becomes a sheet-like structure (e-f). Finally, the secondary vortex extends radially and breaks into two parts. One part is advected around the wake vortex while the second component remains trapped close to the point of contact between the sphere and the wall in frame (g). The motion of a sphere falling under gravity is considered next. The sphere is released from a height of L = 5D with a Reynolds number of 510 based on the initial velocity. The Reynolds number increases to 865 when the particle collides with the wall. The vorticity contours are shown in figure 8. The coefficient of restitution for a dry collision is zero in the left half and 0.5 in the right half of each frame of figure 8. In frame (a), both sides correspond to a time before collision and are the same. As time elapses (frame (b)), a secondary vortex ring is generated in both cases while the wake vortex is moving towards the wall due to its inertia. In the right half of frame (b), a secondary vortex is growing and becoming the wake vortex for the up-going sphere. A separated region on the sphere surface with vorticity of the same sign as the primary vortex is observed in both cases. In addition, a region with the opposite sign of the primary vorticiy is generated in both cases at the wall. In left half of the frame (c), a coherent structure leaves the sphere and extends radially whereas in the right side of this frame, since the particle is moving downward, counterclockwise vorticity spreads out below the sphere. The secondary vortex is divided into two portions: one above the sphere on the surface while the larger portion of this secondary vortex remains between sphere and the wall. In frame (d), since the particle is moving downwards, the clockwise vorticity above the sphere is diffused while the counterclockwise vorticity below the sphere grows. The coherent structure which has extended radially and the primary vortex far above the sphere are similar to the case without rebound. In frames (e-f) of figure 8, the vortices have diffused and clockwise or counterclockwise vorticity is observed around the sphere for upward and downward motion of the sphere, respectively. The separation distance versus nondimensional time is shown in figure 9. Figure 10 shows a comparison between zero coefficient of restitution (left half of each frame) and dry coefficient of restitution (right half of each frame). Frames (a-c) are similar to the previous case, except that in this case, a secondary vortex, which is the wake vortex during the upward motion of the sphere, 8 of 13
9 (a) -τ 0 (b) 0 (c) τ 0 (d) 2τ 0 (e) 3τ 0 (f) 4τ 0 (g) 5τ 0 Figure 7: Collision of a sphere onto a wall. The left and right frames correspond to the present numerical results and the experimental results by Eames and Dalziel,26 respectively. Re = 850, hmin = 9.2 µm, and the time difference between adjacent figures is τ 0 = D U. The vorticity contours are shown and the black line is the streakline. 9 of 13
10 (a) before collision (b) 0.76τ (c) 1.56τ (d) 2.36τ (e) 3.16τ (f) 3.39τ Figure 8: Collision of a sphere onto a wall. The coefficient of restitution for a dry collision is zero in the left half and 0.5 in the right half of each frame. Vorticity contours are shown. Re = 865, h min = 9.8 µm, and τ D = g. 10 of 13
11 4 3 h/d Time/τ Figure 9: Nondimensional separation distance. Re = 865, h min = 9.8 µm, and e dc = 0.5 is larger than the one observed at the same time in figure 8. In frame (d-e), the clockwise vorticity below the sphere moves above it due to the change in direction of the motion of the sphere. In frame (f) this clockwise vorticity which remained from the secondary vortex is advected from the sphere and deviates the wake vortex of the downward moving sphere. The coherent structure near the wall is similar to that for the case without rebound. In general, vorticity is stronger for higher dry coefficient of restitution due to the fact that less energy is dissipated during the collision process. The separation distance versus nondimensional time is shown in figure 9. VI. Conclusions We describe a finite-volume algorithm using a Distributed-lagrange-multiplier for solving particulate flow. We show that our fluid-particle solver produces results which are in good agreement with experimental and numerical studies of benchmark viscous flow problems. One advantage of this approach is that the modifications needed to account for the presence of particles can be easily implemented. Comparison of our numerical results for the bouncing motion of a solid sphere onto a wall shows very good agreement with experimental results. The evolution of vorticity due to collision and the effect of the coefficient of restitution are discussed. References 1 Brady, J. F. and Bossis, G., Stokesian dynamics, Annu. Rev. Fluid Mech., Vol. 20, 1988, pp Ardekani, A. M. and Rangel, R. H., Unsteady motion of two solid spheres in Stokes flow, Phys. of Fluids, Vol. 18, Hu, H. H., Patankar, N., and Zhu, M. Y., Direct numerical simulations of fluid-solid systems using the arbitrary Lagrangian-Eulerian technique, J. Comput. Phys., Vol. 169, 2001, pp Pozrikidis, C., Orbiting motion of a freely suspended spheroid near a plane wall, JFM, Vol. 541, 2005, pp Calhoun, D., A cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions, J. Comput. Phys., Vol. 176, 2002, pp Russell, D. and Wang, Z. J., A cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow, J. Comput. Phys., Vol. 191, 2003, pp Sethian, J. A. and Smereka, P., Level set methods for fluid interfaces, Ann. Rev. Fluid Mech., Vol. 35, 2003, pp Esmaeeli, A. and Tryggvason, G., A front tracking method for computations of boiling in complex geometries, Int. J. Multiphase Flow, Vol. 30, 2004, pp Glowinski, R., Pan, T. W., Hesla, T. I., and Joseph, D. D., A distributed Lagrange multiplier fictitious domain method for particulate flows, Int. J. Multiphase Flow, Vol. 25, 1999, pp Sharma, N., Chen, Y., and Patankar, N. A., A distributed Lagrange multiplier based computational method for the simulation of particulate-stokes flow, Comput. Methods Appl. Mech. Engrg, Vol. 194, 2005, pp patankar, N. A., Singh, P., Joseph, D. D., Glowinski, R., and Pan, T. W., A new formulation of the distributed Lagrange multiplier/ fictitious domain method for particulate flows, Int. J. Multiphase Flow, Vol. 26, 2000, pp Davis, R. H., Serayssol, J. M., and Hinch, E. J., The elastohydrodynamic collision of two spheres, J. Fluid Mech., Vol. 163, 1986, pp Joseph, G., Zenit, R., Hunt, M., and Rosenwinkel, A., Particle-wall collisions in a viscous fluid, J. Fluid Mech., Vol. 433, 2001, pp of 13
12 (a) before collision (b) 0.76τ (c) 1.56τ (d) 2.36τ (e) 3.16τ (f) 4.76τ (g) 5.56τ Figure 10: Collision of a sphere onto a wall. The coefficient of restitution for a dry collision is zero in the left half and 1.0 in the right half of each frame. Vorticity contours are shown. Re = 865, h min = 9.8 µm, and τ D = g. 12 of 13
13 Figure 11: Nondimensional separation distance. Re = 865, h min = 9.8 µm, and e dc = Gondret, P., Hallouin, E., Lance, M., and Petit, L., Experiments on the motion of a solid sphere toward a wall: From viscous dissipation to elastohydrodynamic bouncing, Phys. Fluids, Vol. 11, 1999, pp Gondret, P., Lance, M., and Petit, L., Bouncing motion of spherical particles in fluids, Phys. Fluids, Vol. 14, 2002, pp Sharma, N. and Patankar, N. A., A fast computation technique for direct numerical simulation of rigid particulate flows, J. Comput. Phys., Vol. 205, 2005, pp Patankar, S. V., Numerical heat transfer and fluid flow, McGraw-Hill, Takami, H. and Keller, H. B., Steady two-dimensional viscous flow of an incompressible fluid past a circular cylinder, Phys. Fluids Suppl., Vol. 12, 1969, pp. II51 II Dennis, S. C. R. and Chang, G. Z., Numerical solutions for steady flow past a circular cylinder at Reynolds number up to 100, J. Fluid mech., Vol. 42, 1970, pp Fornberg, B., A numerical study of steady viscous flow past a circular cylinder, J. Fluid mech., Vol. 98, 1980, pp Braza, M., Chassaing, P., and Minh, H. H., Numerical study and physical analysis of pressure and velocity fields in the near wake of a circular cylinder, J. Fluid mech., Vol. 165, 1986, pp Liu, C., Zheng, X., and Sung, C. H., Preconditioned multigrid methods for unsteady incompressible flows, J. Comput. Phys., Vol. 139, 1998, pp Koch, D. L. and Ladd, A. J. C., Moderate Reynolds number flows through periodic and random arrays of aligned cylinders, J. Fluids Mech., Vol. 349, 1997, pp Ardekani, A. M. and Rangel, R. H., A computational method for particulate flow with collisions, in preparation. 25 Thompson, M. C., Hourigan, K., Cheung, A., and Leweke, T., Hydrodynamics of a particle impact on a wall, Applied Mathematical Modelling, Vol. 30, 2006, pp Eames, I. and Dalziel, S. B., Dust resuspension by the flow around an impacting sphere, J. Fluid Mech., Vol. 403, 2000, pp of 13
Journal of Computational Physics
Journal of Computational Physics 227 (2008) 10094 10107 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Collision of multi-particle
More informationParticle-wall collision in a viscoelastic fluid
J. Fluid Mech. (2009), vol. 633, pp. 475 483. c 2009 Cambridge University Press doi:10.1017/s0022112009990632 Printed in the United Kingdom 475 Particle-wall collision in a viscoelastic fluid A. M. ARDEKANI
More informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More informationValidation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
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 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 informationStudy of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields
Study of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields Jie Zhang 1, Cheng Wang 1 1. Department of Mechanical and Aerospace Engineering, Missouri University of Science
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 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 informationPairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray
Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray Georges Akiki, Kai Liu and S. Balachandar * Department of Mechanical &
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationNumerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow
Fluid Structure Interaction VII 147 Numerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow H. G. Sung 1, H. Baek 2, S. Hong 1 & J.-S. Choi 1 1 Maritime and Ocean
More informationVelocity fluctuations resulting from an. immersed sphere-wall collision
Under consideration for publication in J. Fluid Mech. 1 Velocity fluctuations resulting from an immersed sphere-wall collision By A. R U I Z - A N G U L O A N D R. Z E N I T Instituto de Investigaciones
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
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 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 informationVORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
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 informationDIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW
DIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW http://www.aem.umn.edu/solid-liquid_flows Sponsored by NSF-Grand Challenge Grant Fluid Mechanics & CFD Computer Scientists D.D. Joseph Y. Saad R. Glowinski
More informationSingle Curved Fiber Sedimentation Under Gravity. Xiaoying Rong, Dewei Qi Western Michigan University
Single Curved Fiber Sedimentation Under Gravity Xiaoying Rong, Dewei Qi Western Michigan University JunYong Zhu, Tim Scott USDA Forest Products Laboratory ABSTRACT Dynamics of single curved fiber sedimentation
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
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 informationA particle flow specific boundary element formulation for microfluidic applications
A particle flow specific boundary element formulation for microfluidic applications Besim BARANOĞLU 1, Barbaros ÇETİN 2,* * Corresponding author: Tel.: +90 (312) 290-2108; Fax: +90 (312) 266-4126; Email:
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 informationDynamics of particle-particle collisions in a viscous liquid
PHYSICS OF FLUIDS 18, 121506 2006 Dynamics of particle-particle collisions in a viscous liquid F.-L. Yang and M. L. Hunt a Division of Engineering and Applied Sciences, California Institute of Technology,
More informationBeyond the Point Particle: LES-Style Filtering of Finite-Sized Particles
ILASS Americas th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 7 Beyond the Point Particle: LES-Style Filtering of Finite-Sized Particles Brooks Moses and Chris Edwards Department
More informationA numerical study on the effects of cavitation on orifice flow
PHSICS OF FLUIDS, A numerical study on the effects of cavitation on orifice flow S. Dabiri, W. A. Sirignano, and D. D. Joseph, University of California, Irvine, California 9697, USA University of Minnesota,
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More information8.012 Physics I: Classical Mechanics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical Mechanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS INSTITUTE
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationNumerical Investigation of Laminar Flow over a Rotating Circular Cylinder
International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:3 32 Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder Ressan Faris Al-Maliky Department of
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More informationFigure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m
1. For the manometer shown in figure 1, if the absolute pressure at point A is 1.013 10 5 Pa, the absolute pressure at point B is (ρ water =10 3 kg/m 3, ρ Hg =13.56 10 3 kg/m 3, ρ oil = 800kg/m 3 ): (a)
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 informationNumerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder
Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)
More informationLarge-scale flow structures in particle wall collision at low Deborah numbers
European Journal of Mechanics B/Fluids 26 (2007) 511 530 Large-scale flow structures in particle wall collision at low Deborah numbers M. Guala a, A. Stocchino b, a Institute for Environmental Engineering,
More informationEnhancement of Heat Transfer by an Electric Field for a Drop Translating at Intermediate Reynolds Number
Rajkumar Subramanian M. A. Jog 1 e-mail: milind.jog@uc.edu Department of Mechanical, Industrial, and Nuclear Engineering, University of Cincinnati, Cincinnati, OH 45221-0072 Enhancement of Heat Transfer
More informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationNumerical Simulation of Elongated Fibres in Horizontal Channel Flow
Martin-Luther-Universität Halle-Wittenberg Mechanische Verfahrenstechnik 4th Workshop on Two-Phase Flow Predictions Halle, 7-0 September 05 Numerical Simulation of Elongated Fibres in Horizontal Channel
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationShape tilting in viscoelastic fluids
Shape tilting in viscoelastic fluids Shape tilting refers to the fact that falling long bodies tend to line up along the longest line in the body. A cube ill fall in such a ay that the line beteen vertices
More informationConvective Vaporization and Burning of Fuel Droplet Arrays
Convective Vaporization and Burning of Fuel Droplet Arrays Guang Wu and William A. Sirignano Mechanical & Aerospace Engineering University of California, Irvine May 19, 2009 Background and Task In many
More informationReynolds number scaling of inertial particle statistics in turbulent channel flows
Reynolds number scaling of inertial particle statistics in turbulent channel flows Matteo Bernardini Dipartimento di Ingegneria Meccanica e Aerospaziale Università di Roma La Sapienza Paolo Orlandi s 70th
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More informationInfluence of Rotation on the Heat and Fluid Flow around a Circular Cylinder
International Journal of Modern Studies in Mechanical Engineering (IJMSME) Volume 4, Issue 4, 2018, PP 16-29 ISSN 2454-9711 (Online) DOI: http://dx.doi.org/10.20431/2454-9711.0404003 www.arcjournals.org
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationNumerical study of the steady state uniform flow past a rotating cylinder
Numerical study of the steady state uniform flow past a rotating cylinder J. C. Padrino and D. D. Joseph December 17, 24 1 Introduction A rapidly rotating circular cylinder immersed in a free stream generates
More informationNumerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement
Numerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement Sourav Kumar Kar a, 1,, Harshit Mishra a, 2, Rishitosh Ranjan b, 3 Undergraduate Student a, Assitant Proffessor
More information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
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 informationSingle curved fiber sedimentation under gravity
Single curved fiber sedimentation under gravity Xiaoying Rong, Dewei Qi, Guowei He, JunYong Zhu, Tim Scott Abstract Dynamics of single curved fiber sedimentation under gravity are simulated by using the
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More information8.012 Physics I: Classical Mechanics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical Mechanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS INSTITUTE
More informationDEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
DRAFT Proceedings of ASME IMECE: International Mechanical Engineering Conference & Exposition Chicago, Illinois Nov. 5-10, 2006 IMECE2006-14867 DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationINTERACTION BETWEEN A PAIR OF DROPS ASCENDING IN A LINEARLY STRATIFIED FLUID
Proceedings of the ASME 2013 Fluids Engineering Division Summer Meeting FEDSM 2013 July 7-11, 2013, Incline Village, Nevada, USA FEDSM 2013-16046 INTERACTION BETWEEN A PAIR OF DROPS ASCENDING IN A LINEARLY
More informationExternal and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationChapter 4: Fluid Kinematics
Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to
More informationA Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles
A Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles Aleksandar Donev Courant Institute, New York University APS DFD Meeting San Francisco, CA Nov 23rd 2014 A. Donev (CIMS)
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationx j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2
Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological
More informationarxiv: v1 [physics.flu-dyn] 16 Nov 2018
Turbulence collapses at a threshold particle loading in a dilute particle-gas suspension. V. Kumaran, 1 P. Muramalla, 2 A. Tyagi, 1 and P. S. Goswami 2 arxiv:1811.06694v1 [physics.flu-dyn] 16 Nov 2018
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 informationNumerical simulation of the flow behind a rotary oscillating circular cylinder
PHYSICS OF FLUIDS VOLUME 10, NUMBER 4 APRIL 1998 Numerical simulation of the flow behind a rotary oscillating circular cylinder Seung-Jin Baek and Hyung Jin Sung a) Department of Mechanical Engineering;
More informationDYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1551-1555 1551 DYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS by Zhan-Hong WAN a*, Zhen-Jiang YOU b, and Chang-Bin WANG c a Department of Ocean
More informationINTRODUCTION OBJECTIVES
INTRODUCTION The transport of particles in laminar and turbulent flows has numerous applications in engineering, biological and environmental systems. The deposition of aerosol particles in channels and
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationTwo-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet
PHYSICS OF FLUIDS 18, 085103 2006 Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet Suman Muppidi and Krishnan Mahesh Aerospace Engineering & Mechanics,
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
More informationTsorng-Whay Pan. phone: (713) Web page: pan/
Tsorng-Whay Pan Department of Mathematics University of Houston Houston, TX 77204 e-mail: pan@math.uh.edu phone: (713) 743-3448 Web page: www.math.uh.edu/ pan/ Education: 1990 Ph. D., Mathematics 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 informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationContribution of inter-particle collisions on kinetic energy modification in a turbulent channel flow
Contribution of inter-particle collisions on kinetic energy modification in a turbulent channel flow Valentina Lavezzo a, Alfredo Soldati a,b a Dipartimento di Energetica e Macchine and b Centro Interdipartimentale
More informationApplication of the Immersed Boundary Method to particle-laden and bubbly flows
Application of the Immersed Boundary Method to particle-laden and bubbly flows, B. Vowinckel, S. Schwarz, C. Santarelli, J. Fröhlich Institute of Fluid Mechanics TU Dresden, Germany EUROMECH Colloquium
More informationWhere does Bernoulli's Equation come from?
Where does Bernoulli's Equation come from? Introduction By now, you have seen the following equation many times, using it to solve simple fluid problems. P ρ + v + gz = constant (along a streamline) This
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationDynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device *
Chin. J. Chem. Eng., 15(3) 333 338 (2007) Dynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device * YU Zhaosheng( 余钊圣 ) a,b, **, SHAO Xueming( 邵雪明 ) a and Roger
More information8.012 Physics I: Classical Mechanics Fall 2008
IT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical echanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ASSACHUSETTS INSTITUTE
More informationTheoretical Developments in Group Combustion of Droplets and Sprays
Theoretical Developments in Group Combustion of Droplets and Sprays William A. Sirignano University of California, Irvine Collaborations: Guang Wu, current student; Randall Imaoka, former student, US Navy;
More informationSimple shear flow of collisional granular-fluid mixtures
Manuscript Click here to download Manuscript: Manuscript_r1.docx 1 Simple shear flow of collisional granular-fluid mixtures 2 3 4 5 D. Berzi 1 1 Department of Environmental, Hydraulic, Infrastructure,
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 informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationDynamics of magnetic chains in a shear flow under the influence of a uniform magnetic field
Dynamics of magnetic chains in a shear flow under the influence of a uniform magnetic field Tae Gon Kang, 1, a) Martien A. Hulsen, 2 and Jaap M.J. den Toonder 2 1) School of Aerospace and Mechanical Engineering,
More informationNUMERICAL STUDY OF TURBULENT CHANNEL FLOW LADEN WITH FINITE-SIZE NON-SPHERICAL PARTICLES
NUMERICAL STUDY OF TURBULENT CHANNEL FLOW LADEN WITH FINITE-SIZE NON-SPHERICAL PARTICLES Luca Brandt and Mehdi Niazi Ardekani Linné FLOW Centre and SeRC KTH Mechanics SE 44, Stockholm, Sweden luca@mech.kth.se
More information