A phenomenological model of electrorheological fluids

The fundamental assumption of the paper is that the extra stress tensor of an electrorheological fluid is an isotropic tensor valued function of the rate of strain tensor D and the vector n (which characterizes the orientation and length N of the fibers formed by application of an electric field). The resulting constitutive equation for is supplemented by the solution of the previously studied time evolution equation for n. Plastic behavior for the shear and normal stresses is predicted. Anticipating that the action of increasing shear rate is i) to orient the fibers more and more in the direction of flow and ii) simultaneously to break up the fibers leads to the conclusion that for the same behavior is encountered as without an electric field. Using realistically possible approximation formulas for the dependence of and N on leads to the Bingham behavior for and power law behavior for large shear rates.

Boundary layer development on a continuous moving surface with a parallel free stream due to impulsive motion

The development of the momentum and thermal boundary layers over a semi-infinite flat plate has been studied when the external stream as well as the plate are impulsively moved with constant velocities. At the same time the temperature of the wall is suddenly raised from T_∞, the temperature of the surrounding fluid, toT _w and maintained at this temperature. The problem has been formulated in a new system of scaled coordinates such that fort?=0 it reduces to Rayleigh type of equation and fort? → ∞ it reduces to Blasius or Sakiadis type of equation. A new scale of time has been used which reduces the region of integration from an infinite region to a finite region which reduces the computational time considerably. The governing singular parabolic partial differential equations have been solved numerically using an implicit finite difference scheme. For some particular cases, analytical solutions have been obtained. The results show that there is a smooth transition from Rayleigh solution to Blasius or Sakiadis solution as the dimensionless timeξ increases from zero to one. The shear stress at the wall is negative for the friction parameterλ<0.5, positive forλ>0.5 and zero forλ=0.5. The zero shear stress at the wall does not imply separation but corresponds to the parallel flow. The surface heat transfer is strongly dependent on the Prandtl numberPr and increases with it. Also forPr<Pr ₀, the surface heat transfer is enhanced as the friction parameterλ increases, but forPr>Pr ₀ it get reduced.     

On immersions of constant mean curvature: Compactness results and finiteness theorems for Plateau's problem

Assuming stability and integral conditions we show that a sequence of immersed surfaces of constant mean curvatureH converges to an immersedH-surface. The latter theorem depends on an oscillation estimate forH-surfaces based on an isoperimetric inequality. These compactness results are utilized to prove that certain Jordan curves only bound finitely many stable and unstable, immersed, smallH-surfaces.     

On the Influence of the Kinetic Energy¶on the Stability of Equilibria¶of Natural Lagrangian Systems

Natural Lagrangian systems (T,Π) on R ² described by the equation are considered, where is a positive definite quadratic form in and Π(q) has a critical point at 0. It is constructively proved that there exist a C ^∞ potential energy Π and two C ^∞ kinetic energies T and such that the equilibrium q(t)≡ 0 is stable for the system (T,Π) and unstable for the system . Equivalently, it is established that for C ^∞ natural systems the kinetic energy can influence the stability. In the analytic category this is not true. Accepted: October 20, 1999     

Derivatives of the Stretch, Rotation and Exponential Tensors in n-Dimensional Vector Spaces

We present a solution for the tensor equation TX + XT ^T = H, where T is a diagonalizable (in particular, symmetric) tensor, which is valid for any arbitrary underlying vector space dimension n. This solution is then used to derive compact expressions for the derivatives of the stretch and rotation tensors, which in turn are used to derive expressions for the material time derivatives of these tensors. Some existing expressions for n = 2 and n = 3 are shown to follow from the presented solution as special cases. An alternative methodology for finding the derivatives of diagonalizable tensor-valued functions that is based on differentiating the spectral decomposition is also discussed. Lastly, we also present a method for finding the derivatives of the exponential of an arbitrary tensor for arbitrary n.     

Viscosity Solutions of a Degenerate Parabolic-Elliptic System Arising in the Mean-Field Theory of Superconductivity

In a Type‐II superconductor the magnetic field penetrates the superconducting body through the formation of vortices. In an extreme Type‐II superconductor these vortices reduce to line singularities. Because the number of vortices is large it seems feasible to model their evolution by an averaged problem, known as the mean-field model of superconductivity. We assume that the evolution law of an individual vortex, which underlies the averaging process, involves the current of the generated magnetic field as well as the curvature vector. In the present paper we study a two‐dimensional reduction, assuming all vortices to be perpendicular to a given direction. Since both the magnetic field H and the averaged vorticity ω are curl‐free, we may represent them via a scalar magnetic potential q and a scalar stream function ψ, respectively. We study existence, uniqueness and asymptotic behaviour of solutions (ψ, q) of the resulting degenerate elliptic‐parabolic system (with curvature taken into account or not) by means of viscosity and weak solutions. In addition we relate (ψ, q) to solutions (ω, H) of the mean‐field equations without curvature. Finally we construct special solutions of the corresponding stationary equations with two or more superconducting phases. (Accepted August 8, 1997)     

Stability Balls in the Rotating Bénard System

Nonlinear stability in the standard rotating Bénard system with free boundaries is investigated. The perturbations are assumed to be three-dimensional and to be periodic in the horizontal directions. Below the critical value of the Rayleigh numberRa there is conditional stability, i.e. there are nonvanishing stability balls such that perturbations with initial values (measured in a suitable norm) in these balls decay exponentially in time. We give here explicit bounds to these stability balls from below in terms of the parameters of the system, i.e.Ra, the Taylor numberT and the Prandtl numberPr as well as the size of the periodicity cell of the perturbation. The bounds are valid in the entire parameter space; in particular, forPr<1 and for arbitrarily large values ofT. They provide a qualitative explanation for the experimental observation of subcritical instabilities in the rangePr<1. The method is based on a mode expansion of the perturbation equations and explicit estimates of the semigroup operator as well as of the nonlinearity.     

Controllable infinitesimal deformations in homogeneously deformed compressible elastic materials

A controllable static deformation is a deformation that may be maintained in all materials of a given class under the action of surface forces alone. For compressible, homogeneous, isotropic elastic materials the only controllable deformations are homogeneous. However, it is known that there are solutions of the static equations of finite elasticity, linearized about a finite homogeneous deformation, which do not correspond to homogeneous deformations. These approximate solutions are investigated here. Three cases arise, depending on whether none, two, or three of the basic principal stretches are equal.Nomenclature A arbitrary vector potential - a ₁, a ₂, a ₃ bounding coordinates of body - B, B _ij left Cauchy-Green tensor - C, C _ijpq elasticity tensor - c, c ₁, c ₂, c ₃ arbitrary constants - N ₀, N ₁, N ₂ elastic response functions - n vector normal to surface of body - T ₁, T ₂, T ₃ surface tractions - t ₁, t ₂, t ₃ surface tractions - t, t _ij Cauchy stress tensor - t ⁰, t _ij ⁰ Cauchy stress corresponding to basic homogeneous deformation - u, u _i infinitesimal displacement from basic homogeneous deformation - X, X _i position vector in reference state - x, x _i position vector - arbitrary function - _ij Kronecker delta - , ₁, ₂, ₃ principal stretches - arbitrary function - arbitrary function - arbitrary function - I, II, III principal invariants of B     

A second-order approximation to natural convection for large Rayleigh numbers and small Prandtl numbers

The problem under investigation is that of fluid flow within an enclosed rectangular cavity. It is assumed that one wall is maintained at a constant temperature T₁ (hot wall) and the other wall is maintained at a constant temperature T₀ (cold wall). At the remaining walls, two separate cases are studied. In the first, an adiabatic boundary condition is assumed. That is, the normal derivative of the temperature function is assumed to be 0. In the second, it is assumed the temperature varies linearly from T₀ to T₁. The purpose of this paper is the application of a second order numerical technique to the problem of fluid flow within a heated closed cavity. The method is a modification of a method developed by Shay¹ and applied to the driven cavity problem. In order to test the viability of this technique, it was decided to extend the technique to the problem of natural convection in a square. Jones² proposed that this problem is suitable for testing techniques that may be applied to a wide range of practical problems such as reactor insulation, cooling of radioactive waste containers, solar energy collection and others.³ The technique makes use of second-order finite difference approximations to all derivatives in the governing equations. Furthermore, second-order approximations are also used to determine boundary vorticities and, when the adiabatic boundary condition is used, for the boundary temperatures as well. In some works, where second-order approximations are used at interior points, second-order boundary approximations have been sacrificed in favour of a more stable, but first-order boundary approximation. The current approximations are generated by writing the unknown value of a function at a given interior node as a linear combination of unknown function values at all of the neighbouring nodes. Then the function values at these neighbouring nodes are expanded in a Taylor series about the given node. Through appropriate regrouping of terms and the use of the equations to the solved, constraints are imposed on the coefficients of the linear combination to yield a second-order approximation. As it turns out, there are more unknowns than constraints and, as a result, we are left with some freedom in choosing coefficients. In this work this freedom was used to choose coefficients in such a way as to maximize stability of the resulting system of equations. In other words, the approximations to the governing partial differential equation are individually determined at each point dependent on the direction of flow in order to generate the best possible stability. This idea is analogous to that used in the derivation of the upwind method. However, the current method is second-order accurate where the upwind method is only first-order accurate. Thus, what is generated is an easily implemented second-order method that yields a system of equations that has proved easy to solve. The system of equations is solved via the method of successive overrelaxation. The stability of the method is shown in the convergence for a wide range of Rayleigh numbers, Prandtl numbers and mesh sizes. Level curves of the stream, vorticity and temperature functions are provided for Rayleigh numbers (Ra) as large as 100,000, Prandtl numbers (Pr) as small as 0.0001, and mesh sizes as small as 0.0125. Values of the Nusselt number have also been calculated through the use of Simpson's rule, and a second order approximation to the normal derivative of the temperature along the cold wall. Comparisons are made with other current works to aid in the verification of this methods' accuracy and also with the first-order upwind method to demonstrate superiority over the first-order method.     

Further Study on Pure Shear

It is known that the Cauchy stress tensor T is a pure shear when trT = 0. An elementary derivation is given for a coordinate system such that, when referred to this coordinate system, the diagonal elements of T vanish while the off-diagonal elements τ ₁, τ ₂, τ ₃, are the pure shears. The structure of τ _i (i = 1, 2, 3) depends on one non-dimensional parameter q = 54(detT)² / [tr(T ²)]³, 0 ≤ q ≤ 1. When q = 0, one of the three τ _i vanishes. A coordinate system can be chosen such that the remaining two have the same magnitude or one of the remaining two also vanishes. When q = 1, all three τ _i have the same magnitude. However, there is a one-parameter family of coordinate systems that gives the same three τ _i . For q ≠ 0 or 1, none of the three τ _i vanishes and the three τ _i in general have different magnitudes. Nevertheless, a coordinate system can be chosen such that two of the three τ _i have the same magnitude. ^★Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

Extended Polar Decompositions for Plane Strain

In a finite deformation at a particle of a continuous body, a triad of infinitesimal material line elements is said to be “unsheared” when the angles between the three pairs of line elements of the triad suffer no change. In a previous paper, it has been shown that there is an infinity of unsheared oblique triads. With each oblique unsheared triad may be associated an “extended polar decomposition” F = QG = HQ of the deformation gradient F, in which Q is a rotation tensor, and G, H are not symmetric. Both G and H have the same real eigenvalues which are the stretches of the elements of the triad. In this paper, a detailed analysis of extended polar decompositions is presented in the case when the finite deformation is that of plane strain. Then, we may deal with a 2 × 2 deformation gradient F′ = Q′G′ = H′Q′ instead of the full 3 × 3 tensor F. In this case, the extended polar decompositions are associated with “unsheared pairs,” i.e., pairs of infinitesimal material line elements in the plane of strain which suffer no change in angle in the deformation. If one arm of an unsheared pair is chosen in the plane of strain, then, in general, its companion in the plane is determined. It follows that all possible extended polar decompositions may then be described in terms of a single parameter, the angle that the chosen arm makes with a coordinate axis in the plane. Explicit expressions for G′ and H′ are obtained, and various special cases are discussed. In particular, we note that the expressions for G′ and H′ remain valid even when the chosen arm is along a “limiting direction,” that is the direction of a line element which has no companion element in the plane forming an unsheared pair with it. The results are illustrated by considering the cases of simple shear and of pure shear.Dedicated to Professor Piero Villaggio as a symbol of our friendship and esteem.     

THE PROBLEMS OF THE NONLINEAR UNSYMMETRICAL BENDING FOR CYLINDRICALLY ORTHOTROPIC CIRCULAR PLATE（II）

ＴＨＥＰＲＯＢＬＥＭＳＯＦＴＨＥＮＯＮＬＩＮＥＡＲＵＮＳＹＭＭＥＴＲＩＣＡＬ．ＢＥＮＤＩＮＧＦＯＲＣＹＬＩＮＤＲＩＣＡＬＬＹＯＲＴＨＯＴＲＯＰＩＣＣＩＲＣＵＬＡＲＰＬＡＴＥ（ＩＩ）ＨｕａｎｇＪｉａｙｉｎ（黄家寅）;ＱｉｎＳｈｅｎｇｌｉ（秦圣立）;Ｘｉ...     

On Shearing, Stretching and Spin

An analysis is presented of stretching, shearing and spin of material line elements in a continuous medium. It is shown how to determine all pairs of material line elements at a point x, at time t, which instantaneously are not subject to shearing. For a given pair not subject to shearing, a formula is presented for the determination of a third material line element such that all three form a triad not subject to shearing, instantaneously. It is seen that there is an infinity of such triads not subject to shearing. A new decomposition of the velocity gradient L is introduced. In place of the classical decomposition of Cauchy and Stokes, L=d+w, where d is the stretching tensor and w is the spin tensor, the new decomposition is L=?+ℳ, where ?, called the ldquo;modified” stretching tensor, is not symmetric, and ℳ, called the “modified” spin tensor, is skew-symmetric – the tensor ? being chosen so that it has three linearly independent real right (and left) eigenvectors. The physical interpretation of this decomposition is that the material line elements along the three linearly independent right eigenvectors of ? instantaneously form a triad not subject to shearing. They spin as a rigid body with angular velocity μ (say) associated with ℳ. Also, for each decomposition L=?+ℳ, there is a decomposition L=? ^T +, where is also skew-symmetric. The triad of material line elements along the right eigenvectors of ? ^T (the set reciprocal to the right eigenvectors of ?) is also instantaneously not subject to shearing and rotates with angular velocity (say) associated with . It is seen that the vorticity vector ω is the mean of the two angular velocities μ and , ω =(μ+)/2. For irrotational motion, ω =0, so that μ=-; any triad of material line elements suffering no shearing rotates with angular velocity equal and opposite to that of the reciprocal triad of material line elements. It is proved that provided d is not spherical, there is an infinity of choices for ? and ℳ in the decomposition L=?+ℳ. Two special types of decompositions are introduced. The first type is called “CCS-decomposition” (where CCS is an abbreviation for Central Circular Section). It is associated with the infinite family of triads (not subject to shearing) with a common edge along the normal to one plane of central circular section of an ellipsoid ? associated with the stretching tensor, and the two other edges arbitrary in the other plane of central circular section of ?. There are two such CCS-decompositions. The second type is called “triangular decomposition”, because, in a rectangular cartesian coordinate system, ? has three off-diagonal zero elements. There are six such decompositions. Received 14 November 2000 and accepted 2 August 2001     

Plane strain problems in second-order elasticity theory

The equations of second-order elasticity are developed in polar coordinates R, θ for plane strain deformations of incompressible isotropic elastic materials. By considering a ‘displacement function’ the second-order problem is reduced to the solution of an equation of the form ⁴ψ = g(R, Θ) where ² is Laplace's differential operator and g(R, Θ) depends only on the first-order solution. The displacement function technique is then applied to obtain a second-order solution to the problem of an elastic body contained between two concentric rigid circular boundaries, when the outer boundary is held fixed and the inner is subjected to a rigid body translation.     

On the existence of common fixed points for a pair of lipschitzian mappings in banach spaces

ＩｎｔｒｏｄｕｃｔｉｏｎＬｅｔＫｂｅａｎｏｎｅｍｐｔｙｓｕｂｓｅｔｏｆａＢａｎａｃｈｓｐａｃｅＸ .ＴｈｅｎａｍａｐｐｉｎｇＴ :Ｋ→ＫｉｓｓａｉｄｔｏｂｅａＬｉｐｓｃｈｉｔｚｉａｎｍａｐｐｉｎｇｉｆ,ｆｏｒｅａｃｈｉｎｔｅｇｅｒｎ≥ 1 ,ｔｈｅｒｅｅｘｉｓｔｓａｃｏｎｓｔａｎｔｋｎ >0ｓｕｃｈｔｈａｔ‖Ｔｎｘ-Ｔｎｙ‖ ≤ｋｎ‖ｘ-ｙ‖　ｆｏｒａｌｌｘ ,ｙ∈Ｋ .ＡＬｉｐｓｃｈｉｔｚｉａｎｍａｐｐｉｎｇＴｉｓｓａｉｄｔｏｂｅｕｎｉｆｏｒｍｌｙｋ＿Ｌｉｐｓｃｈｉｔｚｉａｎｉｆｋｎ =ｋｆｏｒａｌｌｎ ≥ 1 ;ｎｏ…     

Flow and heat transfer due to impulsive motion of a cone in a viscous fluid

Numerical solution has been obtained for the development of the flow over a cone which is impulsively set into motion. Initially the flow is described by the solution of Rayleigh and then it tends to the ultimate steady state solution of Falkner-Skan equation. But due to the leading edge effect the semi-similar equation describing the transient flow changes its character after certain time and the solution depends also on the ultimate steady state solution of the Falkner-Skan equation. A second-order upwind difference scheme has been used for discretisation. The temperature distribution and heat transfer has also been obtained for constant wall temperature as well as for constant heat flux at the wall. With the increase ofm, Falkner-Skan parameter, the magnitude of skin friction and wall heat transfer increases. It has been found that form≥?0.275 flow separation does not occur.     

Anisotropic polyconvex energies on the basis of crystallographic motivated structural tensors

In large strain elasticity the existence of minimizers is guaranteed if the variational functional to be minimized is sequentially weakly lower semicontinuous (s.w.l.s.) and coercive. Therefore, polyconvex functions which are always s.w.l.s. are usually considered. For isotropic as well as for transversely isotropic and orthotropic materials constitutive functions that are polyconvex already exist. The main goal of this contribution is to provide a new method for the construction of polyconvex hyperelastic models for more general anisotropy classes. The fundamental idea is the introduction of positive definite second-order structural tensors G=HH^T encoding the anisotropies of the underlying crystal. These tensors can be viewed as a push-forward of a cartesian metric of a fictitious reference configuration to the real reference configuration. Here the driving transformations H in the push-forward operation are mappings of the cartesian base vectors of the fictitious configuration onto crystallographic motivated base vectors. Restrictions of this approach are based on the polyconvexity condition as well as on the usage of second-order structural tensors and pointed out in detail.     

Group classification of one-dimensional nonisentropic equations of fluids with internal inertia

A systematic application of the group analysis method for modeling fluids with internal inertia is presented. The equations studied include models such as the nonlinear one-velocity model of a bubbly fluid (with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski Zhurnal Prikladnoj Mekhaniki i Tekhnitheskoj Fiziki 3, 102–111, 1960; Kogarko Dokl. AS USSR 137, 1331–1333, 1961; Wijngaarden J. Fluid Mech. 33, 465–474, 1968), and the dispersive shallow water model (Green and Naghdi J. Fluid Mech. 78, 237–246, 1976; Salmon 1988). These models are obtained for special types of the potential function W(r,[(r)\dot],S){W(\rho,\dot \rho,S)} (Gavrilyuk and Teshukov Continuum Mech. Thermodyn. 13, 365–382, 2001). The main feature of the present paper is the study of the potential functions with W _S  ≠ 0. The group classification separates these models into 73 different classes.