Shear flow of a Newtonian fluid over a quiescent generalized Newtonian fluid

The flow of an upper shear-driven Newtonian fluid above an otherwise still non-Newtonian fluid is considered. The lower fluid is modelled as a generalized Newtonian fluid and set into motion by interfacial shear. By means of similarity transformations, the governing partial differential equations for the two-fluid problem transform exactly into two sets of ordinary differential equations coupled only at the interface. The successful transformation of the two-fluid problem is applied to the particular case when the lower fluid obeys power-law rheology. The resulting three-parameter problem is solved numerically for some different parameter combinations by means of a direct integration approach with the density ratio fixed to unity. We observed that the interfacial velocities decreased with increasing values of the power-law index n in the range from 0.6 to 1.4 whereas the shear-induced motion of the lower fluid penetrates far deeper into a shear-thinning (n < 1) than into a shear-thickening (n > 1) fluid. This phenomenon is ascribed to a corresponding increase of the non-linear viscosity function with lower n-values.     

LCP method for a planar passive dynamic walker based on an event-driven scheme

The main purpose of this paper is to present a linear complementarity problem (LCP) method for a planar passive dynamic walker with round feet based on an event-driven scheme. The passive dynamic walker is treated as a planar multi-rigid-body system. The dynamic equations of the passive dynamic walker are obtained by using Lagrange’s equations of the second kind. The normal forces and frictional forces acting on the feet of the passive walker are described based on a modified Hertz contact model and Coulomb’s law of dry friction. The state transition problem of stick-slip between feet and floor is formulated as an LCP, which is solved with an event-driven scheme. Finally, to validate the methodology, four gaits of the walker are simulated: the stance leg neither slips nor bounces; the stance leg slips without bouncing; the stance leg bounces without slipping; the walker stands after walking several steps.     

Numerical simulation of impact forces by a falling sphere in air using ALE finite element method

ABSTRACT

A numerical method is developed to simulate the process that a falling rigid sphere hits rigid ground and bounces back in air. The problem is treated as fluid-structure interaction problem based on the ALE finite element flow analysis. In order to introduce the numerical process of impact into the present staggered fluid-structure time marching algorithm, the impact force is applied to the equation of motion of the sphere. The magnitude of the impact force is determined by iteration so that the velocity of the sphere after impact converges to zero. Application of the impact force at a single time instant causes unphysical pressure oscillation. This has been suppressed by applying the impact force smoothly over multiple short time steps. In the present method impulse is evaluated instead of impact force. Computations with different density ratio of the sphere to air showed effect of the air on the sphere motion.     

Turbulent Couette–Poiseuille flow with zero wall shear

A particular pressure-driven flow in a plane channel is considered, in which one of the walls moves with a constant speed that makes the mean shear rate and the friction at the moving wall vanish. The Reynolds number considered based on the friction velocity at the stationary wall (u_τ,S) and half the channel height (h) is Re_τ,S = 180. The resulting mean velocity increases monotonically from the stationary to the moving wall and exhibits a substantial logarithmic region. Conventional near-wall streaks are observed only near the stationary wall, whereas the turbulence in the vicinity of the shear-free moving wall is qualitatively different from typical near-wall turbulence. Large-scale-structures (LSS) dominate in the center region and their spanwise spacing increases almost linearly from about 2.3 to 4.2 channel half-heights at this Re_τ,S. The presence of LSS adds to the transport of turbulent kinetic energy from the core region towards the moving wall where the energy production is negligible. Energy is supplied to this particular flow only by the driving pressure gradient and the wall motion enhances this energy input from the mean flow. About half of the supplied mechanical energy is directly lost by viscous dissipation whereas the other half is first converted from mean-flow energy to turbulent kinetic energy and thereafter dissipated.     

A nth-order shear deformation theory for the free vibration analysis on the isotropic plates

In the present paper, a nth-order shear deformation theory is used to perform the free vibration analysis of the isotropic plates. The present nth-order shear deformation theory satisfies the zero transverse shear stress boundary conditions on the top and bottom surface of the plate. Reddy??s third order theory can be considered as a special case of present nth-order theory (n=3). The governing equations and boundary conditions are derived by the principle of virtual work. The governing differential equations of the isotropic plates are solved by the meshless radial point collocation method based on the thin plate spline radial basis function. The effectiveness of the present theory is demonstrated by applying it to free vibration problem of the square and circular isotropic plate.     

Global Solution to the Generalized Riemann Problem in the Presence of a Boundary and Contact Discontinuities

This work is a continuation of our previous work. In the present paper we study the global structure stability of the Riemann solution $u=U(\frac{x}{t})$ containing only contact discontinuities for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary. We prove the existence and uniqueness of a global piecewise C ¹ solution containing only contact discontinuities to a class of the generalized Riemann problems for general n×n quasilinear hyperbolic systems of conservation laws in a half space. Our result indicates that this kind of Riemann solution $u=U(\frac{x}{t})$ mentioned above for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary possesses a global nonlinear structure stability. Some applications to quasilinear hyperbolic systems of conservation laws occurring in physics and other disciplines, particularly to the system describing the motion of the relativistic string in Minkowski space R ^1?+?n , are also given.     

Axisymmetric simulation of the interaction of a rising bubble with a rigid surface in viscous flow

The interaction between a rising deformable gas bubble and a solid wall in viscous liquids is investigated by direct numerical simulation via an arbitrary-Lagrangian–Eulerian (ALE) approach. The flow field is assumed to be axisymmetric. The bubble is driven by gravity only and the motion of the gas inside the bubble is neglected. Deformation of the bubble is tracked by a moving triangular mesh and the liquid motion is obtained by solving the Navier–Stokes equations in a finite element framework. To understand the mechanisms of bubble deformation as it interacts with the wall, the interaction process is studied as a function of two dimensionless parameters, namely, the Morton number (Mo) and Bond number (Bo). We study the range of Bo and Mo from (2, 6.5 × 10⁻⁶) to (16, 0.1). The film drainage process is also considered in this study. It is shown that the deformation of a bubble interacting with a solid wall can be classified into three modes depending on the values of Mo and Bo.     

New integrable problems of motion of a rigid body with a particle oscillating or bouncing in it

Some qualitative aspects of the problem of motion about a fixed point of a rigid body with a particle moving in it in a prescibed (sinusoidal) way was treated in [1–3]. The mechanical system comprised of a rigid body containing an internal mass that moves along a fixed line in the body was considered in several works [4–5]. Recently, an integrable case of this system was found, in which the body is dynamically axisymmetric and moves under no external forces while the particle moves on the axis of dynamical symmetry under the action of Hooke's force to the fixed point [5].In the present note we introduce a more general integrable case in which the particle moves on the axis of dynamical symmetry and is subject to an arbitary conservative force that depends only on the distance from the fixed point. Separation of variables is accomplished and the solution is reduced to quadratures. As a special version of this problem, the case when the particle bounces elastically between two points is briefly discussed.     

Stability of plane oscillations of a satellite in a circular orbit

We study motions of a rigid body (a satellite) about the center of mass in a central Newtonian gravitational field in a circular orbit. There is a known particular motion of the satellite in which one of its principal central axes of inertia is perpendicular to the orbit plane and the satellite itself exhibits plane pendulum-like oscillations about this axis. Under the assumption that the satellite principal central moments of inertia A, B, and C satisfy the relation B = A + C corresponding to the case of a thin plate, we perform rigorous nonlinear analysis of the orbital stability of this motion.In the plane of the problem parameters, namely, the oscillation amplitude ε and the inertial parameter, there exist countably many domains of orbital stability of the satellite oscillations in the linear approximation. Nonlinear orbital stability analysis was carried out in thirteen of these domains. Isoenergetic reduction of the system of equations of the perturbed motion is performed at the energy level corresponding to the unperturbed periodic motion. Further, using the algorithm developed in [1], we construct the symplectic mapping generated by the equations of the reduced system, normalize it, and analyze the stability. We consider resonance and nonresonance cases. For small values of the oscillation amplitude, we perform analytic investigations; for arbitrary values of ε, numerical analysis is used.Earlier, numerical analysis of stability of plane pendulum-like motions of a satellite in a circular orbit was performed in several special cases in [1–4].     

Particle behavior in a turbulent pipe flow with a flat bed

Particle behavior in a turbulent flow in a circular pipe with a bed height h = 0.5R is studied at Re_b = 40,000 and for two sizes of particles (5 μm and 50 μm) using large eddy simulation, one-way coupled with a Lagrangian particle tracking technique. Turbulent secondary flows are found within the pipe, with the curved upper wall affecting the secondary flow formation giving rise to a pair of large upper vortices above two smaller vortices close to the pipe floor. The behavior of the two sizes of particle is found to be quite different. The 50 μm particles deposit forming irregular elongated particle streaks close to the pipe floor, particularly at the center of the flow and the pipe corners due to the impact of the secondary flows. The deposition and resuspension rate of the 5 μm particles is high near the center of the floor and at the pipe corners, while values for the 50 μm particles are greatest near the corners. Near the curved upper wall of the pipe, the deposition rate of the 5 μm particles increases in moving from the wall center to the corners, and is greater than that for the larger particles due to the effects of the secondary flow. The maximum resuspension rate of the smaller particles occurs above the pipe corners, with the 50 μm particles showing their highest resuspension rate above and at the corners of the pipe.     

The hydrodynamic wall effect for a disperse system

The flow of a highly dilute suspension of spheres (radius α) between two parallel ridid planes (distance L) in slow shearing motion is studied. Even for the limiting situation, (α/L) small but finite, there is a layer-one sphere diameter thick—immediately adjacent to the wall in which bulk quantities are so complicated functionals of the parameters of the microstructure that evaluating them seems out of the question. Nevertheless, it is still simple to obtain average bulk quantifies (e.g. apparent viscosity) and even the evaluation of local bulk quantities far away from the wall poses no problem. The reason being that the customary continuum constitutive equation for the bulk stress can be utilized, although a slip velocity has to supplement it. This applies to any disperse system and can be applied to different flows, too. For the spherical suspension at hand an explicit expression for this slip velocity is obtained.     

Forced Convection Flow of Power-Law Fluids Over a Flat Plate Embedded in a Darcy-Brinkman Porous Medium

The steady forced convection flow of a power-law fluid over a horizontal plate embedded in a saturated Darcy-Brinkman porous medium is considered. The flow is driven by a constant pressure gradient. In addition to the convective inertia, also the “porous Forchheimer inertia” effects are taken into account. The pertinent boundary value problem is investigated analytically, as well as numerically by a finite difference method. It is found that far away from the leading edge, the velocity boundary layer always approaches an asymptotic state with identically vanishing transverse component. This holds for pseudoplastic (0 < n < 1), Newtonian (n = 1), and dilatant (n > 1) fluids as well. The asymptotic solution is given for several particular values of the power-law index n in an exact analytical form. The main flow characteristics of physical and engineering interest are discussed in the paper in some detail.     

Transient laminar forced convection to power-law fluids inside ducts resulting from a sudden change in wall temperature

The transient laminar forced convection to power-law fluids in thermally developing, hydrodynamically developed flow inside parallel-plate ducts and circular tubes resulting from a sudden change in wall temperature is studied. The generalized integral transform and the Laplace transform techniques are employed to develop approximate analytic solutions. The local Nusselt number and average fluid temperature are presented over the range of the dimensionless axial coordinate Z varying from 10^?4 to 10^?1 for several dimensionless times. Three different values of the power-law index are considered in the study includedn=1/3,n=1 andn=3 corresponding to, respectively, the pseudoplastic, Newtonian and dilatant fluids.     

Unsteady wave motions of a fluid above an inclined floor

The propagation of unsteady waves above a flat inclined floor within the framework of a linear dispersion model was first studied in [1]. This paper shows how to solve the three-dimensional problem for the case = /4, where is the angle of inclination of the floor plane to the free surface. The two-dimensional problem was studied in [2–4]. In articles [2, 3] asymptotic solutions were found for the Cauchy-Poisson problem for certain values of . In [4], a method is proposed for solving the problem of the wave motion of a fluid due to the displacement of a section of the floor of the basin. However, the complicated structure of the expression obtained by reducing the problem to an inhomogeneous functional equation makes it impossible to study the solution. The aim of the present work is to obtain some exact solutions for the two- and three-dimensional problems of unsteady waves above an inclined floor, which are suitable for calculations and asymptotic estimates.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 65–70, November–December, 1984.     

FInite torsion of slightly compressible rubberlike circular cylinders∗

A regular perturbation procedure and the Rayleigh-Ritz method are used to study the finite torsion of cylinders made of slightly compressible neo-Hookean materials. Two cases are considered: 1. the length of the cylinder is not permitted to change during torsion and 2. the net axial force on the cylinder vanishes. The perturbation procedure fails for certain constitutive relations whereas, in principle, the Rayleigh-Ritz method has general applicability. When it works, the success of the perturbation procedure depends on prior knowledge of the problem for an incompressible material (the zeroth order nonlinear problem). The solution of problem 2. is considerably more complicated than the solution of problem 1. since the complete approximation of order n for problem 2. requires extensive work on the approximation of order (n + 1).     

Solution of a mixed dynamic problem on the displacements given on the boundary of a half-plane lying on the surface of an isotropic homogeneous elastic half-space

We consider the problem on the motion of an isotropic elastic body occupying the half-space z ≥ 0 on whose boundary, along the half-plane x ≥ 0, the horizontal components of displacement are given, while the remaining part of the boundary is stress-free. We seek the solution by the method of integral Laplace transforms with respect to time t and Fourier transforms with respect to the coordinates x, y; the problem is reduced to a system of Wiener-Hopf equations, which can be solved by the methods of singular-integral equations and circulants. We invert the integral transforms and reduce the solution to the Smirnov-Sobolev form. We calculate the tangential stress intensity coefficients near the boundary z = 0, x = 0, |y| < ∞ of the half-plane. The circulant method for solving the Wiener-Hopf system was proposed in [1]. A static problem similar to that considered in the present paper was solved earlier. The Hilbert problem was reduced to a system of Fredholm integral equations in [2]. In the present paper, we solve the above problem by reducing the solution to quadratures and a quasiregular system of Fredholm integral equations. We give a numerical solution of the Fredholm equations and calculate the integrals for the tangential stress intensity coefficients.     

Secular perturbations of quasi-elliptic orbits in the restricted three-body problem with variable masses

In this work we consider the satellite version of the restricted three-body problem when masses of the primary bodies P₀, P₁ vary isotropically with different rates, and their total mass changes according to the joint Meshcherskii law. Equations of motion of the body P₂ of infinitesimal mass are obtained in terms of the osculating elements of the aperiodic quasi-conical motion about the body P₀. Doubly averaging these equations and using the Hill approximation, we have obtained the differential equations determining the secular perturbations of the orbital elements and determined the domains of possible values of the system parameters for which their analytical solutions are expressed in terms of elementary or elliptic functions. The bodies P₀, P₁ mass variation laws for which the corresponding differential equations are integrable, have been found.     

Gravity-driven motion of a deformable drop or bubble near an inclined plane at low Reynolds number

The creeping motion of a three-dimensional deformable drop or bubble in the vicinity of an inclined wall is investigated by dynamical simulations using a boundary-integral method. We examine the transient and steady velocities, shapes, and positions of a freely-suspended, non-wetting drop moving due to gravity as a function of the drop-to-medium viscosity ratio, λ, the wall inclination angle from horizontal, θ, and Bond number, B, the latter which gives the relative magnitude of the buoyancy to capillary forces. For fixed λ and θ, drops and bubbles show increasingly pronounced deformation in steady motion with increasing Bond number, and a continued elongation and the possible onset of breakup are observed for sufficiently large Bond numbers. Unexpectedly, viscous drops maintain smaller separations and deform more than bubbles in steady motion at fixed Bond number over a large range of inclination angles. The steady velocities of drops (made dimensionless by the settling velocity of an isolated spherical drop) increase with increasing Bond number for intermediate-to-large inclination angles (i.e. 45° ? θ ? 75°). However, the steady drop velocity is not always an increasing function of Bond number for viscous drops at smaller inclination angles.     

Axisymmetric spreading of a thin liquid drop with suction or blowing at the horizontal base

The axisymmetric spreading under gravity of a thin liquid drop on a horizontal plane with suction or blowing of fluid at the base is considered. The thickness of the liquid drop satisfies a non-linear diffusion equation with a source term. For a group invariant solution to exist the normal component of the fluid velocity at the base, v_n, must satisfy a first-order quasi-linear partial differential equation. The general form of the group invariant solution for the thickness of the liquid drop and for v_n is derived. Two particular solutions are considered. Each solution depends essentially on only one parameter which can be varied to yield a range of models. In the first solution, v_n is proportional to the thickness of the liquid drop. The base radius always increases even for suction. In the second solution, v_n is proportional to the gradient of the thickness of the liquid drop. The thickness of the liquid drop always decreases even for blowing. A special case is the solution with no spreading or contraction at the base which may have application in ink-jet printing.     

Evolution of deviations from the spherical shape of a vapor bubble in supercompression

This paper considers the evolution of small deviations of a cavitation bubble from a spherical shape during its single compression under conditions of experiments on acoustic cavitation of deuterated acetone. Vapor motion in the bubble and the surrounding liquid is defined as a superposition of the spherical component and its non-spherical perturbation. The spherical component is described taking into account the nonstationary heat conductivity of the liquid and vapor and the nonequilibrium nature of the vaporization and condensation on the interface. At the beginning of the compression process, the vapor in the bubble is considered an ideal gas with a nearly uniform pressure. In the simulation of the high-rate compression stage, realistic equations of state are used. The non-spherical component of motion is described taking into account the effect of liquid viscosity, surface tension, vapor density in the bubble, and nonuniformity of its pressure. Estimates are obtained for the amplitude of small perturbations (in the form of harmonics of degree n = 2, 3, ... with the wavelength λ = 2πR/n, where R is the bubble radius) of the spherical shape of the bubble during its compression until reaching extreme values of pressure, density, and temperature. These results are of interest in the study of bubble fusion since the non-sphericity of the bubble prevents its strong compression.