On Infinitesimal Shear

The setting for this note is the theory of infinitesimal strain in the context of classical linearized elasticity. As a body is subjected to a deformation the angle between a pair of material line elements through a typical point P is changed. The decrease in angle is called the shear of this pair of elements. Here, we determine all pairs of material line elements at P which are unsheared in a deformation. It is seen, in general, that corresponding to any given material line element in a given plane through P, there is one corresponding “companion” material line element such that the given element and its conjugate are unsheared in the deformation. There are two exceptions. If the plane through P is a plane of central circular section of the strain ellipsoid, then every material line element through P in this plane has an infinity of companion elements in this plane – all pairs of material line elements in the plane(s) of central circular section of the strain ellipsoid are unsheared. If the plane through P is not a plane of central circular section of the strain ellipsoid, then there are two exceptional material line elements through P such that neither of them has a companion material line element forming an unsheared pair with it. The directions of these exceptional elements in the plane are called “limiting directions”. It is seen that it is the pair of elements along the limiting directions in a plane which suffer the maximum shear in that plane. A geometrical construction is presented for the determination of the extensional strains along the pairs of elements which are unsheared. Also, it is shown that knowing one unsheared pair in a plane and their extensions is sufficient to determine the principal extensions and the principal axes in this plane. Expressions for all unsheared pairs in a given plane are given in terms of the normals to the planes of central circular sections of the strain ellipsoid. Finally, for a given material line element, a formula is derived for the determination of all other material line elements which form an unsheared pair with the given element.     

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.     

In General There Are No Unsheared Tetrads

Suppose the principal stretches are all different at a point P in a deformed body. In this case, it has been shown [1] that generally there is an infinity of non coplanar infinitesimal material line elements at P which remain unsheared following the deformation – that is, the angle between the arms of each pair of material line elements forming the triad remains unchanged. Here it is shown that in this case when all three principal stretches at P are different, there is no set of four infinitesimal material line elements, no three of which are coplanar, and such that the angle between each pair of the six pairs of material line elements is unchanged following the deformation. It is only when all three principal stretches at P are equal to each other, that there are unsheared tetrads at P, and in that case all tetrads are unsheared. This revised version was published online in June 2006 with corrections to the Cover Date.     

Similarities Between Rotation and Stratification Effects on Homogeneous Shear Flow

Linear theory is applied to examine rotation and buoyancy effects on homogeneous turbulent shear flows with given vertical velocity shear, S=d/dx ₃. In the rotating shear case (where the rotation vector is perpendicular to the plane of the mean flow, Ω _i =Ωδ _{i 2}), general solutions for the Fourier components of the fluctuating velocity are proposed. These solutions are compared with those proposed in the literature for the Fourier components of the fluctuating velocity and density in the case of a homogeneous stratified shear flow with vertical density gradient, S _ρ=d/dx ₃. It is shown that from the normal mode stability stand point the Bradshaw parameter B=2Ω/S(1+2Ω/S) (in the rotating shear case) and the Richardson number R _i (in the statified shear case) play similar roles in identifying the stability for all the wave components except in the case where Ω·k=0, for which rotation has no effects on the flow. Analysis of the long-time behavior of the non-dimensional spectral density of energy, e ^g , is carried out. In the stable case, e ^g has decaying oscillations or undergoes a power law decay in time. Analytical solutions for the streamwise two-dimensional energy ℰ _ii ¹/2 (i.e. the limit at k ₁=0 of the one-dimensional energy spectra) are proposed. At large time, ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (3R _i +1)/(4R _i ) except at R _i =1 it stays constant in time. Similar behavior for ℰ _ii ¹(t)/ℰ _ii ¹(0) is also observed in the rotating shear case (ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (1+4B)/(4B)). Due to the behavior of the dimensionless spectral density of energy in both flow cases, the turbulent kinetic energy, /2, the production rate, ?, and the rate due to the buoyancy forces, ℬ, are split into two parts, , ?=?₁+?₂, ℬ=ℬ₁+ℬ₂ (in the stratified shear case, both ?₁ and ℬ₁ vanish when R _i >?, while in the rotating shear case one has ℬ=0). It is shown that when rotation is “cyclonic” (i.e. Ω/S>0), part reaches maximum magnitudes at St ≈2, independent of the B value, and the first time to a zero crossing of ?₂ occurs at this particular value. When rotation is “anticyclonic” (i.e. Ω/S<0) one finds St ≈1.6 instead of St ≈2. In the stratified shear case, both ?₂ and ℬ₂ cross zero at Nt=St ≈2, and part reaches maximum magnitudes at this particular value. These results and in particular those for the turbulent kinetic energy are compared with previous direct numerical simulation (DNS) results in homogeneous stratified shear flows. Received 30 July 2001 and accepted 19 February 2002     

Identities in Finite Strain

First of all the deformation is considered of two infinitesimal material line elements lying along vectors M,N emanating from a particle at X in a body. For all M,N lying in a given plane, an identity is derived relating the stretches along M,N and the angles of the pair of infinitesimal material line elements before and after deformation. Then, the deformation is considered of three non-coplanar infinitesimal material line elements lying along vectors M,N,P emanating from a particle at X in a body. An identity is derived relating the stretches along M,N,P and the angles between the three pairs of infinitesimal material line elements before and after deformation. The identity is factored leading to easy interpretation. The special case of infinitesimal strain is considered.      

On Anisotropic Elastic Materials for which Young''s Modulus E ( n ) is Independent of n or the Shear Modulus G ( n , m ) is Independent of n and m

It is shown that, among anisotropic elastic materials, only certain orthotropic and hexagonal materials can have Young modulus E(n) independent of the direction n or the shear modulus G(n,m) independent of n and m. Thus the direction surface for E(n) can be a sphere for certain orthotropic and hexagonal materials. The structure of the elastic compliance for these materials is presented, and condition for identifying if the material is orthotropic or hexagonal is given. We also study the case in which n of E(n) and n, m of G(n,m) are restricted to a plane. When E(n) is a constant on a plane so are G(n,m) and Poisson's ratio ν(n,m). The converse, however, does not necessarily hold. A plane on which E(n) is a constant can exist for all anisotropic elastic materials. In particular, existence of such a plane is assured for trigonal, hexagonal and cubic materials. In fact there are four such planes for a cubic material. For these materials, not only E(n) is a constant, two other Young's moduli, the three shear moduli and the six Poisson's ratio on the plane are also constant.     

Asymptotic Stability of Riemann Solutions for a Class of Multidimensional Systems of Conservation Laws with Viscosity

We prove the asymptotic stability of two-state nonplanar Riemann solutions for a class of multidimensional hyperbolic systems of conservation laws when the initial data are perturbed and viscosity is added. The class considered here is those systems whose flux functions in different directions share a common complete system of Riemann invariants, the level surfaces of which are hyperplanes. In particular, we obtain the uniqueness of the self-similar L^∞ entropy solution of the two-state nonplanar Riemann problem. The asymptotic stability to which the main result refers is in the sense of the convergence as t→∞ in L_loc¹ of the space of directions ξ = x/t. That is, the solution u(t, x) of the perturbed problem satisfies u(t, tξ)→R(ξ) as t→∞, in L_loc¹(ℝⁿ), where R(ξ) is the self-similar entropy solution of the corresponding two-state nonplanar Riemann problem.     

Effects of Large Scales Motion in Unstable Stratified Shear Flows

Homogeneous turbulence under unstable uniform stratification (N ² < 0) and vertical shear is investigated by using the linear theory (or the so-called rapid distortion theory, RDT) for an initial isotropic turbulence over a range −∞ ≤ R _i =N ²/S ² ≤ 0. The initial potential energy is zero and P _r =1 (i.e. the molecular Prandtl number).One-dimensional (streamwise) k ₁−spectra, especially Θ₃₃(k ₁) (i.e., that of the vertical kinetic energy, are investigated. In agreement with previous experiments, it is found that the unstable stratification affects the turbulence quantities at all scales. A significant increase of the vertical kinetic energy is observed at low wavenumbers k ₁ (i.e. large scales) due to an increase of the stratification . The effect of the shear (S) is appreciable only at high wavenumbers k ₁ (i.e. small scales).Based on the importance of the spectral components with k ₁ = 0, the asymptotic forms of Θ _ij (k ₁ = 0) or equivalently the so-called “two-dimensional” energy components (2DEC) are analyzed in detail. The asymptotic form for the ratio of 2DEC is compared to the long-time limit of the ratio of real energies. In the unstable stratified shearless case (S=0,N ² ≠ 0) the variances and the covariances of the velocity and the density are derived analytically in terms of the Weber functions, while when S ≠ 0 and N ² ≠ 0 they are obtained numerically (−100 ≤ R _i <0 and . The results are discussed in connection to previous experimental results in unstable stratified open channel flows cooled from above by Komori et al. Phy Fluids 25, 1539–1546 (1982).It is shown that the Richardson number dependence of the long-time limit of the ratios of real energies is well described by this “simple” model (i.e. the dependence of the long-time limit of 2DEC on R _i ). For example, the ratio of the potential energy to the kinetic energy (q ²/2), approaches −R _i /(1−R _i ), the ratio of turbulent energy production by buoyancy forces to production by shearing forces (i.e. the flux Richardson number, R _if ), approaches R _i . Also, the Richardson number dependence of the principal angle (β) of the Reynolds stress tensor and the angle (β_ρ) of the scalar flux vector is fairly predicted by this model .On the other hand, it is found that the above ratios are insensitive to viscosity, while the ratios ɛ /q ² and , depend on the viscosity and they evolve asymptotically like t ⁻¹. The turbulent Froude number, F _rt =(L _Oz /L _E )^2/3, where L _Oz and L _E are the Ozmidov length scale and the Ellison length scale, respectively, evolves asymptotically like t ^−1/3.     

On Determinants having Complex Conjugate Columns or Rows

It is proved that the determinant, det D, of an N × N matrix D having n ≤ N/2 pairs of complex conjugate columns (or rows), while all other elements are real-valued, is given by det D =(−2i)ⁿdet S,n = 0,1,2,... in which S is a certain residual matrix having real-valued elements. Thus, det D is either real-valued or pure imaginary according as n is even (including n = 0 ) or odd, respectively. The general theorem is illustrated in an example. This revised version was published online in August 2006 with corrections to the Cover Date.     

On Aerodynamic Forces for Viscous Compressible Flow

The paper is aimed at reviewing and adding some new results to our recent work on a force theory for viscous compressible flows around a finite body. It has been proposed to analyze aerodynamic forces directly in terms of fluid elements of nonzero vorticity and density gradient. Let ρ denote the density, u the velocity, and ω the vorticity. It is demonstrated that for largely separated flows about bluff bodies, there are two major source elements: R _e(x) =−?u ²∇ρ·∇ϕ and V _e(x) =−u×ω·∇ϕ, where ϕ is an acyclic potential, generated by the solid body moving with unit velocity in the negative direction of the force considered. In particular, under mild conditions, the (unique) choice of ϕ enforces that the elements R _e(x) and V _e(x) decay rapidly away from the body. Four kinds of finite body are considered: a circular cylinder, a sphere, a hemi sphere-cylinder, and a delta wing of elliptic section—all in transonic-to-supersonic regimes. From an extensive numerical study carried out for these bodies, it is found that these two elements contribute to 95% or more of the total drag or lift for all the cases under consideration. Moreover, R _e(x) due to density gradient becomes progressively important relative to V _e(x) due to vorticity as the Mach number increases. The present method of force analysis enables effective analysis and assessment of relative importance of aerodynamics forces, contributed from individual flow structures. The analysis could therefore be very much useful in view of the rapid growth in numerical fluid dynamics; detailed (either local or global) flow information is often available. The paper is dedicated to Sir James Lighthill in honor of his great contributions to aeronautics on the occasion of the publication of his collected works. Received 3 January 1997 and accepted 11 April 1997     

Scherk-Type Capillary Graphs

This paper concerns the regularity of a capillary graph (the meniscus profile of liquid in a cylindrical tube) over a corner domain of angle α. By giving an explicit construction of minimal surface solutions previously shown to exist (Indiana Univ. Math. J. 50 (2001), no. 1, 411–441) we clarify two outstanding questions. Solutions are constructed in the case α = π/2 for contact angle data (γ₁, γ₂) = (γ, π − γ) with 0 < γ < π. The solutions given with |γ − π/2| < π/4 are the first known solutions that are not C² up to the corner. This shows that the best known regularity (C^{1, ∈}) is the best possible in some cases. Specific dependence of the H?lder exponent on the contact angle for our examples is given. Solutions with γ = π/4 have continuous, but horizontal, normal vector at the corners in accordance with results of Tam (Pacific J. Math. 124 (1986), 469–482). It is shown that our examples are C^{0, β} up to and including the corner for any β < 1. Solutions with |γ − π/2| > π/4 have a jump discontinuity at the corner. This kind of behavior was suggested by numerical work of Concus and Finn (Microgravity sci. technol. VII/2 (1994), 152–155) and Mittelmann and Zhu (Microgravity sci. technol. IX/1 (1996), 22–27). Our explicit construction, however, allows us to investigate the solutions quantitatively. For example, the trace of these solutions, excluding the jump discontinuity, is C^2/3.     

Symmetry of Mountain Pass Solutions of Some Vector Field Equations

We prove radial symmetry (or axial symmetry) of the mountain pass solution of variational elliptic systems − AΔu(x) + ∇ F(u(x)) = 0 (or − ∇.(A(r) ∇ u(x)) + ∇ F(r,u(x)) = 0,) u(x) = (u ₁(x),...,u _N (x)), where A (or A(r)) is a symmetric positive definite matrix. The solutions are defined in a domain Ω which can be , a ball, an annulus or the exterior of a ball. The boundary conditions are either Dirichlet or Neumann (or any one which is invariant under rotation). The mountain pass solutions studied here are given by constrained minimization on the Nehari manifold. We prove symmetry using the reflection method introduced in Lopes [(1996), J. Diff. Eq. 124, 378–388; (1996), Eletron. J. Diff. Eq. 3, 1–14].     

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     

Flow in micellar cubic crystals

 We investigate and compare the mechanical properties of two micellar cubic crystals. Both are obtained from aqueous triblock-copolymer solutions of similar block structure, (PEO)₇₆(PPO)₂₉(PEO)₇₆ and (PEO)₁₂₇(PPO)₄₈(PEO)₁₂₇, designated Pluronic F68 and F108 respectively. Extensive small angle X-ray scattering experiments under shear and common rheometry provide the following results: (i) the solid phases of the two polymeric solutions are constituted of micelles that are ordered in a bcc (F68) and fcc (F108) structure. The overall appearance of both unsheared samples is polycrystalline; (ii) SAXS showed that both systems undergo a structural reorientation under shear, with the dense planes arranged parallel to the velocity - vorticity plane, thus facilitating flow; (iii) F68 showed a clear transition in the flow curve, associated with a textural change, but the F108 system exhibited a continuous evolution; (iv) the bcc crystal appears to have no measurable yield stress and flows even for low applied stresses. This was in contrast to the fcc sample which showed a clear yield stress, separating creep and flow regimes. Received: 30 November 1999/Accepted: 30 November 1999     

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.     

Viscoelastic properties of ultra-high viscosity alginates

Ultra-high viscosity alginates were extracted from the brown seaweeds Lessonia nigrescens (UHV_N, containing 61% mannuronate (M) and 2% guluronate (G)) and Lessonia trabeculata (UHV_T, containing 22% M and 78% G). The viscoelastic behavior of the aqueous solutions of these alginates was determined in shear flow in terms of the shear stress σ ₂₁, the first normal stress difference N ₁, and the shear viscosity η in isotonic NaCl solutions (0.154 mol/L) at T = 298 K in dependence of the shear rate [(g)\dot]\dot{\gamma} for solutions of varying concentrations and molar masses (3–10 × 10⁵ g/mol, homologous series was prepared by ultrasonic degradation). Data obtained in small-amplitude oscillatory shear (SAOS) experiments obey the Cox–Merz rule. For comparison, a commercial alginate with intermediate chemical composition was additionally characterized. Particulate substances which are omnipresent in most alginates influenced the determination of the material functions at low shear rates. We have calculated structure–property relationships for the prediction of the viscosity yield, e.g., η–M _w–c–[(g)\dot]\dot{\gamma} for the Newtonian and non-Newtonian region. For the highest molar masses and concentrations, the elasticity yield in terms of N ₁ could be determined. In addition, the extensional flow behavior of the alginates was measured using capillary breakup extensional rheometry. The results demonstrate that even samples with the same average molar mass but different molar mass distributions can be differentiated in contrast to shear flow or SAOS experiments.     

Elastic properties of a dilute surfactant solution in the shear induced state

The first normal stress differences N ₁ of a highly dilute cationic surfactant solution are investigated in a cone-and-plate rheometer. In continuation of a previous paper (Nowak 1998), where the buildup of a shear induced structure in such a solution was attained after a reduced deformation, the N ₁ turned out to be in proportion to the square of the shear rate γ˙ reduced by a critical value γ˙ _c in a first range above γ˙ _c . At higher shear rates the N ₁ tend to lower values than predicted by this relation. Relaxation experiments were performed in the same geometry to determine the characteristic time scales of the shear induced state's decay. In the lower range above &γdot; _c the stress decay is a monoexponential process, while a second time constant has to be introduced to describe the relaxation in that range, where the N ₁ deviate from the parabolic dependence of the reduced shear rate. Received: 10 May 1999 Accepted: 15 November 2000     

Comparison of the vortex shedding behavior of a single plate and a plate array

This study compares the shedding behavior around and downstream of a single plate positioned in a flow field alone with the shedding behavior around and downstream of the same plate positioned in an array of identical plates. The shedding frequencies and corresponding Strouhal numbers based on chord [S _r (c)] and based on thickness [S _r (t)] are obtained using a hot-wire anemometer. In comparison with the plate positioned as a single plate, the same plate placed in a plate array shows increases in S _r (c) of up to 55.5% and produces a dominant peak in the power spectra that is wider by a factor of 3.5. In contrast to the single-plate results, which exhibit step changes in S _r (c) of about 0.6 at c/t ≈ 6, 8 and 11, the plate positioned in an array shows only one abrupt transition at c/t ≈ 4. Received: 26 January 1999/Accepted: 7 February 2000     

Lagrangian Time-Scales in Homogeneous Non-Gaussian Turbulence

Lagrangian time-scales in homogeneous non-Gaussian turbulence were studied using a one-dimensional Lagrangian Stochastic Model. The existence of two time-scales τ _L and T _L , one typical of the inertial subrange and the other which is an integral property, is outlined. Variations of the ratio T _L /τ _L in the plane skewness-flatness (S, F) are shown and a connection with the statistical constraint F ≥S ² + 1 is evidenced. The Lagrangian autocorrelation function ρ(t) of particle velocity was computed for some values of (S, F). It is shown that for small times, say t < T _L , the influence of non-Gaussianity is negligible and ρ(t) presents the same behaviour as in the Gaussian case regardless of variations in (S, F).As the time increases, departures from Gaussianity are observed and autocorrelation turns out to be always larger than in the Gaussiancase. This is supported by some considerations in terms of information entropy, which is shown to decrease with increasing departures from Gaussianity. Spectral analysis of Lagrangian velocity shows that non-Gaussianity is relevant only to large scales of the stochastic process and that the expected inertial subrange decay ω⁻² is attained by spectra of all simulations, except for one case in which the model probability density function is bimodal, due to the vicinity to the statistical limit. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.