Exact solutions of the equations of motion for an incompressible viscoelastic Maxwell medium

The unsteady plane-parallel motion of a incompressible viscoelastic Maxwell medium with constant relaxation time is considered. The equations of motion of the medium and the rheological relation admit an extended Galilean group. The class of solutions of this system which are partially invariant with respect to the subgroup of the indicated group generated by translation and Galilean translation along one of the coordinate axes is studied. The system does not have invariant solutions, and the set of partially invariant solutions is very narrow. A method for extending the set of exact solutions is proposed which allows finding solutions with a nontrivial dependence of the stress tensor elements on spatial coordinates. Among the solutions obtained by this method, the solutions describing the deformation of a viscoelastic strip with free boundaries is of special interest from a point of view of physics. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 16–23, March–April, 2009.     

Mathematical model of an incompressible viscoelastic Maxwell medium

Nonstationary motions of incompressible viscoelastic Maxwell continuum with a constant relaxation time are considered. Because in an incompressible continuous medium, pressure is not a thermodynamic variable but coincides with the stress-tensor trace to within a factor, it follows that, separating the spherical part from this tensor, one can assume that the remaining part of the stress tensor has zero trace. In the case of an incompressible medium, the equations for the velocity, pressure, and stress tensor form a closed system of first-order equations which has both real and complex characteristics, which complicates the formulation of the initial-boundary-value problem. Nevertheless, the resolvability of the Cauchy problem can be proved in the class of analytic functions. Unique resolvability of the linearized problem was established in the classes of functions of finite smoothness. The class of effectively one-dimensional motions for which the subsystem of three equations is a hyperbolic one was studied. The results of an asymptotic analysis of the latter imply the possible formation of discontinuities during the evolution of the solution. The general system of equations of motion admits an infinite-dimensional Lie pseudo-group which contains an extended Galilean group. The theorem of the invariance of the conditions on the a priori unknown free boundary was proved to obtain exact solutions of free-boundary problems. The problem of deformation of a viscoelastic strip subjected to tangential stresses applied to the free boundary is considered as an example of application of this theorem. In this problem, a scale effect of short-wave instability caused by the absence of diagonal dominance of the stress tensor deviator was found.     

Nonlinear waves in incompressible viscoelastic Maxwell medium

In this paper we study two-dimensional flows of incompressible viscoelastic Maxwell media with Jaumann corotational derivative in the rheological constitutive law. In the general case, due to the incompressibility condition, the equations of motion have both real and complex characteristics. Group properties of this system are studied. On this basis, two submodels of the Maxwell model are selected, which can be reduced to hyperbolic ones. More precisely, we consider plane shear flow between two parallel planes and Couette type flow caused by the inertial cylinder rotation. As a result, we obtain the closed systems of three equations of mixed type, which describe nonlinear transverse waves in an incompressible Maxwell fluid. It is demonstrated that discontinuities can develop in elastic media even from smooth initial data. Stability of shocks in the Maxwell fluid with and without retardation time is discussed.     

A series solution for the effective properties of incompressible viscoelastic media

This paper presents a series solution for the homogenization problem of a linear viscoelastic periodic incompressible composite. The method uses the Laplace transform and the correspondence principle which are combined with the classical expansion along Neumann series of the solution of the periodic elasticity problem in Fourier space. The terms of the Neumann series appear as decoupled, containing geometry dependent terms and viscoelastic properties dependent terms which are polynomial fractions whose inverse Laplace transforms are provided explicitly.     

Strain localization in an oscillating Maxwell viscoelastic cylinder

The transient rotation responses of simple, axisymmetric, viscoelastic structures are of interest for interpretation of experiments designed to characterize materials and closed structures such as the brain using magnetic resonance techniques. Here, we studied the response of a Maxwell viscoelastic cylinder to small, sinusoidal displacement of its outer boundary. The transient strain field can be calculated in closed form using any of several conventional approaches. The solution is surprising: the strain field develops a singularity that appears when the wavefront leaves the center of the cylinder, and persists as the wavefront reflects to the outer boundary and back to the center of the cylinder. The singularity is alternately annihilated and re-initiated upon subsequent departures of the wavefront from the center of the cylinder until it disappears in the limit of steady state oscillations. We present the solution for this strain field, characterize the nature of this singularity, and discuss its potential role in the mechanical response and evolved morphology of the brain.     

Nonlinear motion equations for a non-Newtonian incompressible fluid in an orthogonal coordinate system

This paper describes the work carried out to investigate the pore pressures occurring in secondary consolidation. A theoretical approach and an experimental technique was developed in order to conduct the study. By considering compression to occur only due to water leaving the soil it was possible to derive an expression for the dissipation of pore pressure in the secondary phase. By further simplified assumptions which are based on experimental observations, the above general solution was reduced to a simple formula which predicted the observed behaviour of pore water pressures during secondary consolidation.     

Nonlinear motion equations for a Non-Newtonian incompressible fluid in an orthogonal coordinate system

Summary Starting with an assumed relationship between the stress tensor, the non-Newtonian viscosity, and the strain rate tensor, the nonlinear equations of motion are developed for use in any orthogonal coordinate system. The resulting equations are written in terms of the scalar velocities, the non-Newtonian viscosity, the metric coefficients, and their derivatives.The non-Newtonian viscosity is assumed to be a scalar function of the strain rate tensor, and so depends upon the invariants of the strain rate tensor. For convenience, the necessary invariants are written out in complete form for use in any orthogonal coordinate system, in terms of the scalar velocities, the metric coefficients, and their derivatives.Using the resulting motion equations and a model of this type of viscosity, theOstwald-de Waele model, an example of time dependent flow is solved using a continuous time-discrete space method programmed on an analog computer. e _ij strain rate tensor - body force density, dynes/cm³ - F ₁,F ₂,F ₃ components of body force density, dynes/cm³ - g acceleration of gravity - H function of time - h ₁,h ₂,h ₃ metric coefficients - I ₁,I ₂,I ₃ invariants - m constant - P pressure, dynes/cm² - r radius, cm - t time, sec - velocity vector, cm/sec - v ₁,v ₂,v ₃ velocities in thex ₁,x ₂ andx ₃ directions, respectively, cm/sec - v _n (t) velocity of thenth node, cm/sec - x ₁,x ₂,x ₃ coordinate directions - z coordinate, cm - unit tensor - _ij Kronecker delta - _ij 2e _ij - nabla - _ijk alternating unit tensor - non-Newtonian viscosity, dynes/cm² - ₀, ₁ constant viscosities, dynes sec/cm², dynes sec ^m /cm² - angle, radians - v ₀,v ₁ constant kinematic viscosities, cm²/sec, cm² sec ^m-2 - density, g/cm³ - _ij stress tensor - fluid dilation With 3 figures     

Minimum principles for incompressible viscoelastic solids

Summary Four minimum principles are established, governing the quasi-static deformations of linear, isotropic, incompressible viscoelastic solids.

---

Sommario Si stabiliscono quattro principi di minimo che governano le deformazioni quasi statiche di solidi viscoelastici, lineari, isotropici e incompressibili.

     

Controllability of viscoelastic stresses for nonlinear Maxwell models

We investigate a class of models for viscoelastic fluids, in which the elastic stress is determined by a conformation tensor, and the conformation tensor is linked to the velocity field by a system of ordinary differential equations. We study the question which values of the conformation tensor can be reached in a homogeneous flow, subject to a given initial condition and arbitrary velocity fields. This problem is a special “easy” case for the question of controllability of viscoelastic flows. For a class of models, we show that constraints on the values of the conformation tensor are given by lower and/or upper bounds on its determinant. The behavior of seemingly similar models, e.g. the PTT, Giesekus and Peterlin dumbbell models, turns out to be surprisingly different.     

饱和黏弹性多孔介质中的平面波及能量耗散

研究了流体饱和不可压黏弹性多孔介质中的非均匀平面波及其能量流和能量耗散规律. 在流 相和固相物质微观不可压、固相骨架宏观服从积分型本构关系和小变形的假定下，利用 Helmholtz分解，得到了饱和黏弹性多孔介质中非均匀平面波的一般解以及纵波、横波相速 度和衰减率等的解析表达式，分析了平面波传播矢量和衰减矢量之间的关系. 数值结果表明 孔隙流体与固相骨架间的相互作用以及固相骨架的黏性对波的相速度、衰减率等有着显著的 影响. 同时，得到了饱和黏弹性多孔介质的能量方程，给出了能量流矢量和能量耗散率. 对 非均匀平面纵波和横波，推导了平均能量流矢量和平均能量耗散率的解析表达式.     

Unsteady flow of ground-water in an aquifer system with viscoelastic properties

The transient flow behaviour of groundwater in an aquifer-aquitard system with viscoelastic properties is studied. On the basis of previous work (Hantush, Neuman, Brutsaert, Corapcioglu), new partial-differential integral equations are assembled. The well-known equations (Hantush, Brutsaert) are special cases of these equations, which describe the flow of slightly compressible groundwater in aquifer-aquitard layers with viscoelastic properties.Analytical solutions of the coupled partial differential-integral equations are obtained by Laplace transform. The viscoelastic properties enhance the heterogeneities of an elastic aquifer system with delay and feed. Numerical inversion of the Laplace transformed drawdown is in good agreement with the analytical solutions obtained.     

Modulation equations for a mixture of gas bubbles in an incompressible liquid

Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 86–92, March–April, 1989.     

Numerical study around the corotational Maxwell model for the viscoelastic fluid flows

We present numerical results for the FEM (finite element method) presented in [Comput. Methods Appl. Mech. Engrg. 191 (2002) 5045–5065]. This method is devoted to the approximation of fluid flows obeying the Oldroyd model. A particularity of this method, is to take into account the purely viscoelastic case, the so-called Maxwell model, important in practice. Numerical results are given for a fluid flowing in an abrupt plane 4 to 1 contraction. We use the corotational Maxwell model as benchmark in the choice of our computations. Results are also given for the upper convected Maxwell model. Interesting effects appear on the velocity profile: a phenomenon of quasi slip at the downstream wall.     

Characteristic properties of the system of equations of a shear flow with nonmonotonic velocity profile

Institute of Hydrodynamics, Russian Academy of Sciences; Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizka, Vol. 36, No. 3, pp. 53–59, May–June, 1995.     

On the smallest scale for the incompressible Navier-Stokes equations

We prove that, for solutions to the two- and three-dimensional incompressible Navier-Stokes equations, the minimum scale is inversely proportional to the square root of the Reynolds number based on the kinematic viscosity and the maximum of the velocity gradients. The bounds on the velocity gradients can be obtained for two-dimensional flows, but have to be assumed in three dimensions. Numerical results in two dimensions are given which illustrate and substantiate the features of the proof. Implications of the minimum scale result, to the decay rate of the energy spectrum are discussed.Research was supported in part by the National Acronautics and Space Administration under NASA Contract No. NAS1-18107, while the second author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. Additional support was provided by the National Science Foundation under Grant DMS-8312264 and the Office of Naval Research under Contract N-00014-83-K-0422.     

Implementation of the Lagrange-Galerkin method for the incompressible Navier-Stokes equations

This paper is concerned with the implementation of Lagrange-Galerkin finite element methods for the Navier-Stokes equations. A scheme is developed to efficiently handle unstructed meshes with local refinement, using a quad-tree-based algorithm for the geometric search. Several difficulties that arise in the construction of the right-hand side are discussed in detail and some useful tricks are proposed. The resulting method is tested on the lid-driven square cavity and the vortex shedding behind a rectangular cylinder and is found to give satisfactory agreement with previous works. A detailed analysis of the effect of time discretization is included.