Hu, Howard H.
Email Address
ORCID
Disciplines
relationships.isProjectOf
relationships.isOrgUnitOf
Position
Introduction
Research Interests
Collection
18 results
Search Results
Now showing 1 - 10 of 18
Publication Rheology of a Suspension of Elastic Particles in a Viscous Shear Flow(2011-11-01) Hu, Howard H; Gao, Tong; Ponte-Castañeda, PedroIn this paper we consider a suspension of elastic solid particles in a viscous liquid. The particles are assumed to be neo-Hookean and can undergo finite elastic deformations. A polarization technique, originally developed for analogous problems in linear elasticity, is used to establish a theory for describing the finite-strain, time-dependent response of an ellipsoidal elastic particle in a viscous fluid flow under Stokes flow conditions. A set of coupled, nonlinear, first-order ODEs is obtained for the evolution of the uniform stress fields in the particle, as well as for the shape and orientation of the particle, which can in turn be used to characterize the rheology of a dilute suspension of elastic particles in a shear flow. When applied to a suspension of cylindrical particles with initially circular cross-section, the theory confirms the existence of steady-state solutions, which can be given simple analytical expressions. The two-dimensional, steady-state solutions for the particle shape and orientation, as well as for the effective viscosity and normal stress differences in the suspension, are in excellent agreement with direct numerical simulations of multiple-particle dispersions in a shear flow obtained by using an arbitrary Lagrangian–Eulerian (ALE) finite element method (FEM) solver. The corresponding solutions for the evolution of the microstructure and the rheological properties of suspensions of initially spherical (three-dimensional) particles in a simple shear flow are also obtained, and compared with the results of Roscoe (J. Fluid Mech., vol. 28, 1967, pp. 273–293) in the steady-state regime. Interestingly, the results show that sufficiently soft elastic particles can be used to reduce the effective viscosity of the suspension (relative to that of the pure fluid).Publication Direct Simulation of Initial Value Problems for the Motion of Solid Bodies in a Newtonian Fluid. Part 2. Couette adn Poiseuille Flows.(1994-05-11) Feng, James; Hu, Howard H.; Joseph, Daniel D.This paper reports the results of a two-dimensional finite element simulation of the motion of a circular particle in a Couette and a Poiseuille flow. The size of the particle and the Reynolds number are large enough to include fully nonlinear inertial effects and wall effects. Both neutrally buoyant and non-neutrally buoyant particles are studied, and the results are compared with pertinent experimental data and perturbation theories. A neutrally buoyant particle is shown to migrate to the centreline in a Couette flow, and exhibits the Segré-Silberberg effect in a Poiseuille flow. Non-neutrally buoyant particles have more complicated patterns of migration, depending upon the density difference between the fluid and the particle. The driving forces of the migration have been identified as a wall repulsion due to lubrication, an inertial lift related to shear slip, a lift due to particle rotation and, in the case of Poiseuille flow, a lift caused by the velocity profile curvature. These forces are analysed by examining the distributions of pressure and shear stress on the particle. The stagnation pressure on the particle surface are particularly important in determining the direction of migration.Publication Direct Numerical Simulation of the Sedimentation of Solid Particles with Thermal Convection(2003-01-01) Gan, Hui; Chang, Jianzhong; Feng, James J.; Hu, Howard H.Dispersed two-phase flows often involve interfacial activities such as chemical reaction and phase change, which couple the fluid dynamics with heat and mass transfer. As a step toward understanding such problems, we numerically simulate the sedimentation of solid bodies in a Newtonian fluid with convection heat transfer. The two-dimensional Navier–Stokes and energy equations are solved at moderate Reynolds numbers by a finite-element method, and the motion of solid particles is tracked using an arbitrary Lagrangian–Eulerian scheme. Results show that thermal convection may fundamentally change the way that particles move and interact. For a single particle settling in a channel, various Grashof-number regimes are identified, where the particle may settle straight down or migrate toward a wall or oscillate laterally. A pair of particles tend to separate if they are colder than the fluid and aggregate if they are hotter. These effects are analysed in terms of the competition between the thermal convection and the external flow relative to the particle. The mechanisms thus revealed have interesting implications for the formation of microstructures in interfacially active two-phase flows.Publication The Dynamics of Two Spherical Particles in a Confined Rotating Flow: Pedalling Motion(2008-03-25) Mukundakrishnan, Karthik; Hu, Howard H.; Ayyaswamy, Portonovo S.We have numerically investigated the interaction dynamics between two rigid spherical particles moving in a fluid-filled cylinder that is rotating at a constant speed. The cylinder rotation is about a horizontal axis. The particle densities are less than that of the fluid. The numerical procedure employed to solve the mathematical formulation is based on a three-dimensional arbitrary Larangian–Eulerian (ALE), moving mesh finite-element technique, described in a frame of reference rotating with the cylinder. Results are obtained in the ranges of particle Reynolds number, 1Publication Lubricated Pipelining: Stability of Core-Annular Flow. Part 2(1989-01-05) Hu, Howard H.; Joseph, Daniel D.In this paper, we study the linearized stability of three symmetric arrangements of two liquids in core-annular Poiseuille flow in round pipes. Deferring to one important application, we say oil and water when we mean more viscous and less viscous liquids. The three arrangements are (i) oil is in the core and water on the wall, (ii) water is in the core and oil is outside and (iii) three layers, oil inside and outside with water in between. The arrangement in (iii) is our model for lubricated pipelining when the pipe walls are hydrophobic and i t has not been studied before. The arrangement in (ii) was studied by Hickox (1971) who treated the problem as a perturbation of long waves, effectively suppressing surface tension and other essential effects which are necessary to explain the flows observed, say, in recent experiments of W. L. Olbricht and R. W. Aul. The arrangement in (i) was studied in Part 1 of this paper (Preziosi, Chen & Joseph 1987). We have confirmed and extended their pseudo-spectral calculation by introducing a more efficient finite-element code. We have calculated neutral curves, growth rates, maximum growth rate, wavenumbers for maximum growth and the various terms which enter into the analysis of the equation for the evolution of the energy of a small disturbance. The energy analysis allows us to identify the three competing mechanisms underway : interfacial tension, interfacial friction and Reynolds stress. Many results are presented.Publication Particle Motion in a Liquid Film Rimming the Inside of a Partially Filled Rotating Cylinder(2003-07-24) Joseph, Daniel D.; Wang, J.; Bai, R.; Yang, B. H.; Hu, Howard H.Both lighter- and hydrophobic heavier-than-liquid particles will float on liquid–air surfaces. Capillary forces cause the particles to cluster in typical situations identified here. This kind of clustering causes particles to segregate into islands and bands of high concentrations in thin liquid films rimming the inside of a slowly rotating cylinder partially filled with liquid. A second regime of particle segregation, driven by secondary motions induced by off-centre gas bubbles in a more rapidly rotating cylinder at higher filling levels, is identified. A third regime of segregation of bidisperse suspensions is found in which two layers of heavier-than-liquid particles that stratify when there is no rotation, segregate into alternate bands of particles when there is rotation.Publication Migration of a sphere in tube flow(2005-10-10) Yang, B. H; Wang, J.; Joseph, D. D; Hu, Howard H; Pan, T. -W.; Glowinski, R.The cross-stream migration of a single neutrally buoyant rigid sphere in tube flow is simulated by two packages, one (ALE) based on a moving and adaptive grid and another (DLM) using distributed Lagrange multipliers on a fixed grid. The two packages give results in good agreement with each other and with experiments. A lift law L=CUs (Ωs— Ωse) analogous to L=ρUΓ which was proposed and validated in two dimensions is validated in three dimensions here; C is a constant depending on material and geometric parameters, Us is the slip velocity and it is positive, Ωs is the slip angular velocity and Ωse is the slip angular velocity when the sphere is in equilibrium at the Segré–Silberberg radius. The slip angular velocity discrepancy Ωs— Ωse is the circulation for the free particle and it changes sign with the lift. A method of constrained simulation is used to generate data which is processed for correlation formulas for the lift force, slip velocity, and equilibrium position. Our formulae predict the change of sign of the lift force which is necessary in the Segré–Silberberg effect. Our correlation formula is compared with analytical lift formulae in the literature and with the results of two-dimensional simulations. Our work establishes a general procedure for obtaining correlation formulae from numerical experiments. This procedure forms a link between numerical simulation and engineering practice.Publication Stability of Core-Annular Flow with Very Small Viscosity Ratio(1990-07-24) Hu, Howard H.; Lundgren, Thomas S.; Joseph, Daniel D.It is known that the stability problem for core-annular flow of very viscous crude oil and water is singular, the water annulus appears to be inviscid with boundary layers at the pipe wall and at the interface. In the present paper, this singular problem is treated by the method of matched asymptotic expansions using ∈ = m/ℝα as a small parameter. There are two cases of instability corresponding to different positions of the critical point in the annulus. One case is when the critical point is far away from the interface, the other is when the critical point is close to the interface within a distance of order ∈1/3. In both cases, the equations for the eigenvalues are derived, and the explicit forms for the neutral curves are given. The stability problem is also treated by the modified finite element code used by Hu and Joseph [J. Fluid Mech. 205, 359 ( 1989); Phys. Fluids A 1, 1659 ( 1989)], taking into account the boundary layers at the pipe wall and at the interface. The results of the two methods agree where they overlap, but the finite element technique goes further.Publication Direct Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids(1997-02-04) Huang, P. Y.; Feng, James J; Hu, Howard H.; Joseph, Daniel D.This paper reports the results of direct numerical simulation of the motion of a two-dimensional circular cylinder in Couette flow and in Poiseuille flow of an Oldroyd-B fluid. Both neutrally buoyant and non-neutrally buoyant cylinders are considered. The cylinder's motion and the mechanisms which cause the cylinders to migrate are studied. The stable equilibrium position of neutrally buoyant particles varies with inertia, elasticity, shear thinning and the blockage ratio of the channel in both shear flows. Shear thinning promotes the migration of the cylinder to the wall while inertia causes the cylinder to migrate away from the wall. The cylinder moves closer to the wall in a narrower channel. In a Poiseuille flow, the effect of elastic normal stresses is manifested by an attraction toward the nearby wall if the blockage is strong. If the blockage is weak, the normal stresses act through the curvature of the inflow velocity profile and generate a lateral force that points to the centreline. In both cases, the migration of particles is controlled by elastic normal stresses which in the limit of slow flow in two dimensions are compressive and proportional to the square of the shear rate on the body. A slightly buoyant cylinder in Couette flow migrates to an equilibrium position nearer the centreline of the channel in a viscoelastic fluid than in a Newtonian fluid. On the other hand, the same slightly buoyant cylinder in Poiseuille flow moves to a stable position farther away from the centreline of the channel in a viscoelastic fluid than in a Newtonian fluid. Marked effects of shear thinning are documented and discussed.Publication Direct Simulation of Initial Value Problems for the Motion of Solid Bodies in a Newtonian Fluid Part 1. Sedimentation(1994) Feng, James; Hu, Howard H.; Joseph, Daniel D.This paper reports the result of direct simulations of fluid-particle motions in two dimensions. We solve the initial value problem for the sedimentation of circular and elliptical particles in a vertical channel. The fluid motion is computed from the Navier-Stokes equations for moderate Reynolds numbers in the hundreds. The particles are moved according to the equations of motion of a rigid body under the action of gravity and hydrodynamic forces arising from the motion of the fluid. The solutions are as exact as our finite-element calculations will allow. As the Reynolds number is increased to 600, a circular particle can be said to experience five different regimes of motion: steady motion with and without overshoot and weak, strong and irregular oscillations. An elliptic particle always turn its long axis perpendicular to the fall, and drifts to the centreline of the channel during sedimentation. Steady drift, damped oscillation and periodic oscillation of the particle are observed for different ranges of the Reynolds number. For two particles which interact while settling, a steady staggered structure, a periodic wake-action regime and an active drafting- kissing-tumbling scenario are realized at increasing Reynolds numbers. The nonlinear effects of particle-fluid, particle-wall and interparticle interactions are analysed, and the mechanisms controlling the simulated flows are shown to be lubrication, turning couples on long bodies, steady and unsteady wakes and wake interactions. The results are compared to experimental and theoretical results previously published.