Non-isothermal spherical Couette flow of Oldroyd-B fluid

The present paper is concerned with non-isothermal spherical Couette flow of Oldroyd-B fluid in the annular region between two concentric spheres. The inner sphere rotates with a constant angular velocity while the outer sphere is kept at rest. The viscoelasticity of the fluid is assumed to dominate the inertia such that the latter can be neglected in the momentum and energy equations. An approximate analytical solution is obtained through the expansion of the dynamical variable fields in power series of Nahme number. Non-homogeneous, harmonic for axial- velocity and temperature equations and biharmonic for stream function equations, have been solved up to second order approximation. In comparison of the present work with isothermal case; [1,2], two additional terms; a first order velocity and a second order stream function are stem as a result of the interaction between the fluid viscoelasticity and temperature profile. These contributions prove to be the most important results for rheology in this work.     

Non-isothermal spherical Couette flow of Oldroyd-B fluid

The present paper is concerned with non-isothermal spherical Couette flow of Oldroyd-B fluid in the annular region between two concentric spheres. The inner sphere rotates with a constant angular velocity while the outer sphere is kept at rest. The viscoelasticity of the fluid is assumed to dominate the inertia such that the latter can be neglected in the momentum and energy equations. An approximate analytical solution is obtained through the expansion of the dynamical variable fields in power series of Nahme number. Non-homogeneous, harmonic for axial- velocity and temperature equations and biharmonic for stream function equations, have been solved up to second order approximation. In comparison of the present work with isothermal case; [1,2], two additional terms; a first order velocity and a second order stream function are stem as a result of the interaction between the fluid viscoelasticity and temperature profile. These contributions prove to be the most important results for rheology in this work.     

On the boundary value problem of the biharmonic operator on domains with angular corners

The paper is concerned with boundary singularities of weak solutions of boundary value problems governed by the biharmonic operator. The presence of angular corner points or points at which the type of boundary condition changes in general causes local singularities in the solution. For that case the general theory of V. A. Kondrat'ev provides a priori estimates in weighted Sobolev norms and asymptotic singular representations for the solution which essentially depend on the zeros of certain transcendental functions. The distribution of these zeros will be analysed in detail for the biharmonic operator under several boundary conditions. This leads to sharp a priori estimates in weighted Sobolev norms where the weight function is characterized by the inner angle of the boundary corner. Such estimates for “negative” Sobolev norms are used to analyse also weakly nonlinear perturbations of the biharmonic operator as, for instance, the von Kármán model in plate bending theory and the stream function formulation of the steady state Navier-Stokes problem. It turns out that here the structure of the corner singularities is essentially the same as in the corresponding linear problem.     

Numerical study of spherical Couette flows for certain zenith-angle-dependent rotations of boundary spheres at low Reynolds numbers

The numerical method with splitting of boundary conditions developed previously by the first and third authors for solving the stationary Dirichlet boundary value problem for the Navier-Stokes equations in spherical layers in the axisymmetric case at low Reynolds numbers and a corresponding software package were used to study viscous incompressible steady flows between two con-centric spheres. Flow regimes depending on the zenith angle ?? of coaxially rotating boundary spheres (admitting discontinuities in their angular velocities) were investigated. The orders of accuracy with respect to the mesh size of the numerical solutions (for velocity, pressure, and stream function in a meridional plane) in the max and L ₂ norms were studied in the case when the velocity boundary data have jump discontinuities and when some procedures are used to smooth the latter. The capabilities of the Richardson extrapolation procedure used to improve the order of accuracy of the method were investigated. Error estimates were obtained. Due to the high accuracy of the numerical solutions, flow features were carefully analyzed that were not studied previously. A number of interesting phenomena in viscous incompressible flows were discovered in the cases under study.     

A parametric method of constructing Poincaré mappings in hydrodynamic systems

The plane-parallel motion of the particles of an incompressible medium reduces to an investigation of a Hamilton system. The stream function is a Hamilton function. A Hamilton function, which depends periodically on time and corresponds to the agitation of an incompressible medium in a domain which varies periodically with time, is considered. This agitation of the medium is due to dynamic chaos. The transition to dynamic chaos is described by investigating the location of the Lagrangian particles over time intervals which are multiples of the period (Poincaré points (PP)). The set of PP can be obtained using a Poincaré mapping in the phase flow. The method which has been developed is used to investigate the plane-parallel motion of the particles in an incompressible fluid in a thin layer, bounded by a flat bottom, rectilinear side walls and an upper boundary which is deformed according to a specified periodic law. The motion of the particles is determined from Hamilton's system of equations. The Hamiltonian (the stream function) is found in the thin-layer approximation and depends on two dimensionless parameters: the amplitude of deformation and the tangential velocity in the deforming boundary. The characteristic boundary, which separates the domain of the chaotic motion of the PP from the domain of ordered motion, is determined numerically in the domain of the two parameters. The topological structure of the phase trajectories up to the transition to chaotic conditions is analysed using the Poincaré mapping, found with an accuracy up to the third order with respect to the amplitude. The phase trajectories of the PP, found analytically, turn out to be close to the trajectories of the initial Hamilton system, determined numerically. The mapping found in the domain of the two dimensionless parameters enables one to determine, qualitatively, the boundary of the transition to chaos.     

About singularities at angular points of a trailing edge under the Joukowskii–Kutta Condition

The linear problem for the velocity potential around a slightly curved thin finite wing is considered under the Joukowskii–Kutta hypothesis. The exponents of possible singularities of solutions at angular points on wing's trailing edge are expressed in terms of eigenvalues of mixed boundary value problems for the Beltrami–Laplace operator on the hemisphere and the semicircle. These singularities have a structure such that the circulation function turns out to be continuous in interior angular points of the trailing edge. In the case of trapezoidal shape of the wing ends there occur square-root singularities of the velocity field at the trailing edge endpoints and the same singularities, of course, are extended along the lateral sides of the wake behind the wing. It is proved that for any angular point on the trailing edge the exponents of all above-mentioned singularities form a countable set in the upper complex half-plane with the only accumulation point at infinity. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

The optimal convergence rates of non-isentropic subsonic Euler flows through the infinitely long three-dimensional axisymmetric nozzles

This paper is concerned about the optimal convergence rates of non-isentropic subsonic flows at far fields in three-dimensional infinitely long axisymmetric nozzles. By using the stream function formulation for the compressible Euler equations, the subsonic Euler flows are equivalent to a quasilinear elliptic equation of the stream function. The key points to prove the convergence rates of subsonic flows at far fields are the choice of compared functions and the maximum principles.     

Creeping flow of a second grade fluid in a corner

A variant of Taylor’s (1962) [23] scraper problem, in which, the lower plate rotates is considered. The non-linear partial differential equations governing the flow of a second grade fluid are modeled and solved by using the domain perturbation technique considering the angular velocity of the rotating plate as a small parameter. Also the rheology of the second grade fluid is examined by depicting the profiles of the velocity, stream function, pressure and stress fields.     

Flow past a porous sphere at small Reynolds number

low of an incompressible viscous fluid past a porous sphere has been discussed. The flow has been divided in three regions. The Region-I is the region inside the porous sphere in which the flow is governed by Brinkman equation with the effective viscosity different from that of the clear fluid. In Regions II and III clear fluid flows and Stokes and Oseen solutions are respectively valid. In all the three regions Stokes stream function is expressed in powers of Reynolds number. Stream function of Region II is matched with that of Region I at the surface of the sphere by the conditions suggested by Ochao-Tapia and Whitaker and it is matched with that of Oseen’s solutions far away from the sphere. It is found that the drag on the sphere reduces significantly when it is porous and it decreases with the increase of permeability of the medium.Received: February 7, 2002; revised: April 8, 2003 / June 9, 2004     

Stabilization and optimal decay rate for a non-homogeneous rotating body-beam with dynamic boundary controls

In this paper, we consider the boundary stabilization of a flexible beam attached to the center of a rigid disk. The disk rotates with a non-uniform angular velocity while the beam has non-homogeneous spatial coefficients. To stabilize the system, we propose a feedback law which consists of a control torque applied on the disk and either a dynamic boundary control moment or a dynamic boundary control force or both of them applied at the free end of the beam. By the frequency multiplier method, we show that no matter how non-homogeneous the beam is, and no matter how the angular velocity is varying but not exceeding a certain bound, the nonlinear closed loop system is always exponential stable. Furthermore, by the spectral analysis method, it is shown that the closed loop system with uniform angular velocity has a sequence of generalized eigenfunctions, which form a Riesz basis for the state space, and hence the spectrum-determined growth condition as well as the optimal decay rate are obtained.     

Stokes流和理想流体轴对称问题的解析函数法

本文利用流函数解的完备性和共轭势函数的概念,导出了轴对称Stokes流和理想流体完备的速度和压力的解析函数表达式解.作为它的应用,我们求出关于球的缓慢绕流问题的解.     

Two-link billiard trajectories: extremal properties and stability

Two-link periodic trajectories of a plane convex billiard, when a point mass moves along a segment which is orthogonal to the boundary of the billiard at its end points, are considered. It is established that, if the caustic of the boundary lies within the billiard, then, in a typical situation, there is an even number of two-link trajectories and half of them are hyperbolic (and, consequently, unstable) and the other half are of elliptic type. An example is given of a billiard for which the caustic intersects the boundary and all of the two-link trajectories are hyperbolic. The analysis of the stability is based on an analysis of the extremum of a function of the length of a segment of a convex billiard which is orthogonal to the boundary at one of its ends.     

Optimal rigid body motions,part 2: Minimum time solutions

The time-optimal control of rigid-body angular rates is investigated in the absence of direct control over one of the angular velocity components. The existence of singular subarcs in the time-optimal trajectories is explored. A numerical survey of the optimality conditions reveals that, over a large range of boundary conditions, there are in general several distinct extremal solutions. A classification of extremal solutions is presented, and domains of existence of the extremal subfamilies are established in a reduced parameter space. A locus of Darboux points is obtained, and global optimality of the extremal solutions is observed in relation to the Darboux points. The continuous dependence of the optimal trajectories with respect to variations in control constraints is noted, and a procedure to obtain the time-optimal bang-bang solutions is presented.This work was supported in part by DARPA under Contract No. ACMP-F49620-87-C-0016, by SDIO/IST under Contract No. F49620-87-C0088, and by Air Force Grant AFOSR-89-0001.     

Surface waves on steady perfect‐fluid flows with vorticity

This is a theory of two‐dimensional steady periodic surface waves on flows under gravity in which the given data are three quantities that are independent of time in the corresponding evolution problem: the volume of fluid per period, the circulation per period on the free stream line, and the rearrangement class (equivalently, the distribution function) of the vorticity field. A minimizer of the total energy per period among flows satisfying these three constraints is shown to be a weak solution of the surface wave problem for which the vorticity is a decreasing function of the stream function. This decreasing function can be thought of as an infinite‐dimensional Lagrange multiplier corresponding to the vorticity rearrangement class being specified in the minimization problem. (Note that functional dependence of vorticity on the stream function was not specified a priori but is part of the solution to the problem and ensures the flow is steady.) To illustrate the idea with a minimum of technical difficulties, the existence of nontrivial waves on the surface of a fluid flowing with a prescribed distribution of vorticity and confined beneath an elastic sheet is proved. The theory applies equally to irrotational flows and to flows with locally square‐integrable vorticity. © 2011 Wiley Periodicals, Inc.     

A priori bounds and well-posedness of a system associated with unsteady boundary layer flows

An integral equation with singularities is introduced to characterize unsteady laminar boundary layer flows and some properties of solutions of this integral equation are investigated. Utilizing these properties, a priori bounds are obtained for the skin friction function and the similarity stream function and the well-posedness of solutions is proved.     

Exterior Stokes flows with stick-slip boundary conditions

Steady two-dimensional creeping flows induced by line singularities in the presence of an infinitely long circular cylinder with stick-slip boundary conditions are examined. The singularities considered here include a rotlet, a potential source and a stokeslet located outside a cylinder and lying in a plane containing the cylinder axis. The general exterior boundary value problem is formulated and solved in terms of a stream function by making use of the Fourier expansion method. The solutions for various singularity driven flows in the presence of a cylinder are derived from the general results. The stream function representation of the solutions involves a definite integral whose evaluation depends on a non-dimensional slip parameter l₁\lambda_1. For extremal values, l₁ = 0\lambda_1 = 0 and l₁ = 1\lambda_1 = 1, of the slip parameter our results reduce to solutions of boundary value problems with stick (no-slip) and perfect slip conditions, respectively.¶The slip parameter influences the flow patterns significantly. The plots of streamlines in each case show interesting flow patterns. In particular, in the case of a single rotlet/stokeslet (with axis along y-direction) flows, eddies are observed for various values of l₁\lambda_1. The flow fields for a pair of singularities located on either side of the cylinder are also presented. In these flows, eddies of different sizes and shapes exist for various values of l₁\lambda_1 and the singularity locations. Plots of the fluid velocity on the surface show locations of the stagnation points on the surface of the cylinder and their dependencies on l₁\lambda_1 and singularity locations.     

Conservation laws for laminar axisymmetric jet flows with weak swirl

The conservation laws for laminar axisymmetric jet flows with weak swirl are studied here. The multiplier approach is used to derive the conservation laws for the system of three boundary layer equations for the velocity components governing flow in laminar axisymmetric jet flows with weak swirl. Conservation laws for the system of two partial differential equations for the stream function are also derived.     

Accumulation of tangent points with the boundary and Lagrangian manifolds in problems with phase constraints

The phenomenon of the existence of optimal trajectories having an infinite number of tangent points with the boundary of the phase constraint is studied. The optimal synthesis and the structure of the corresponding Lagrangian manifolds are obtained. The main tool is the blowing up procedure for the resolution of the singularities of the Poincaré mapping of the constraint boundary on itself. A singularity appears at the point of accumulation of the tangent points with the boundary.     

Numerical study of the basic stationary spherical couette flows at low Reynolds numbers

Previously developed iterative numerical methods with splitting of boundary conditions intended for solving an axisymmetric Dirichlet boundary value problem for the stationary Navier-Stokes system in spherical layers are used to study the basic spherical Couette flows (SCFs) of a viscous incompressible fluid in a wide range of outer-to-inner radius ratios R/r (1.1 ≤ R/r ≤ 100) and to classify such SCFs. An important balance regime is found in the case of counter-rotating boundary spheres. The methods converge at low Reynolds numbers (Re), but a comparison with experimental data for SCFs in thin spherical layers show that they converge for Re sufficiently close to Re_cr. The methods are second-order accurate in the max norm for both velocity and pressure and exhibit high convergence rates when applied to boundary value problems for Stokes systems arising at simple iterations with respect to nonlinearity. Numerical experiments show that the Richardson extrapolation procedure, combined with the methods as applied to solve the nonlinear problem, improves the accuracy up to the fourth and third orders for the stream function and velocity, respectively, while, for the pressure, the accuracy remains of the second order but the error is nevertheless noticeably reduced. This property is used to construct reliable patterns of stream-function level curves for large values of R/r. The possible configurations of fluid-particle trajectories are also discussed.     

Thin Sets and Boundary Behavior of Solutions of the Helmholtz Equation

The Martin boundary for positive solutions of the Helmholtz equation in n-dimensional Euclidean space may be identified with the unit sphere. Let v denote the solution that is represented by Lebesgue surface measure on the sphere. We define a notion of thin set at the boundary and prove that for each positive solution of the Helmholtz equation, u, there is a thin set such that u/v has a limit at Lebesgue almost every point of the sphere if boundary points are approached with respect to the Martin topology outside this thin set. We deduce a limit result for u/v in the spirit of Nagel–Stein (1984).