- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 524718, 7 pages
Nonholonomic Geometry of Viscoanelastic Media and Experimental Confirmation
1Department of Mathematics and Computer Science, University of Messina, Viale F. Stagno d'Alcontres 31, 98166 Messina, Italy
2Department of Mathematics, University of Salerno, Via Ponte Don Melillo, 84084 Fisciano, Italy
Received 16 July 2013; Accepted 9 August 2013
Academic Editor: Cristian Toma
Copyright © 2013 Armando Ciancio and Carlo Cattani. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
A thermodynamical model for viscoanelastic media is analyzed using the nonholonomic geometry. A 27-dimensional manifold is introduced, and the differential equations for the geodetics are determined and analytically solved. It is shown that, in this manifold, the best specific entropy is a harmonic function. In the linear case the propagation of transverse acoustic waves is studied, and the theoretical results are compared with some experimental data from a polymeric material (polyisobutylene).
From macroscopic point of view the most popular mathematical approaches to nonequilibrium thermodynamics are based both on the Caratheodory theory [1–3] which involves Pfaff equations and on the contact structure of thermodynamic state space [4–6]. If these equations are completely integrable, then the thermodynamics is called holonomic otherwise nonholonomic.
The theory of nonequilibrium holonomic thermodynamics for mechanical phenomena in continuous media was developed in the 90s (see e.g.,  and references therein). The phenomenological equations were derived [8, 9] by introducing the tensorial internal variables () which occur in the entropy production.
In particular, if one linearizes this theory by neglecting the cross effects among the irreversible phenomena (heat flow, mechanical viscosity, and anelastic deformations), the following rheological equations for distortional phenomena are obtained : where and are the deviators of the stress () and the strain () tensors, respectively.
The quantities are given by in which is the relaxation time and are the state coefficients, while and are the phenomenological coefficients related to the following physical phenomena: In this paper we reconsider the theory from the point of view of the nonholonomic geometry.
Finally, in Section 5 we study the transverse waves and we will show the connection between complex numbers and the shear complex modulus. By applying these results to a polymeric material, as polyisobutylene we will also show that the expected results of the theoretical model are in agreement with the experimental data.
2. Gibbs-Pfaff Equation of Viscoanelastic Media
Let us define the space of states where is the specific entropy, is the absolute temperature, is the specific internal energy, is the symmetric equilibrium stress tensor, is the total symmetric strain tensor, and is the symmetric affinity stress tensor conjugate to the anelastic tensor (the symmetric tensorial thermodynamic internal variable).
Let be the specific volume related to the mass density (homogeneous media) by . The Gibbs-Pfaff equation of viscoanelastic media is This equation defines in the nonholonomic contact distribution of dimension so that the highest dimension of integral manifold is .
The representation of this integral manifold (of maximum dimension) is usually given as follows: This parametrization is based on the arbitrary -function . So that we have a family of integral manifolds of dimension indexed on the arbitrary function .
The -representation is a generic element in the -jet space converted into In order to fix a representative for the specific entropy , we need a supplementary condition as follows.
Theorem 1. If is an homogeneous function of order one, then
Proof. The condition of homogeneity of order one gives the PDE Then the entropy (solution of this PDE) appears as the potential (9).
Corollary 2. If the specific entropy is an homogeneous function of order one, then(1)the variables are conjugated to the intensive variables ;(2)the variables are not essential parameters (because they are not independent).
3. Specific Entropy via Least Squares Lagrangian
The most convenient way to fix a representative of the specific entropy is to look at (5) as a partial derivative evolution equation and to build the least squares method Lagrangian and the functional where , , and .
The extremals are solutions of the Euler-Lagrange PDE where is the total derivative with respect to the variable .
In our case, we get So that, by replacing the partial derivatives of in (17), there follows the Laplace equation for the entropy Consequently, we have the following.
Theorem 3. The best entropy for the nonholonomic nonequilibrium thermodynamics is an harmonic function.
Any curve in the distribution (5) is described by In order to be a geodesic, this curve must minimize the energy functional In short, we must solve the problem To solve this problem, we use the method of Lagrange multipliers. For this we defined the constrained Lagrangian where is the Lagrangian multiplier.
Theorem 4. The geodesics of the nonholonomic viscoanelastic distribution are solutions of the Euler-Lagrange ODEs for the Lagrangian (23) where are the generalized coordinates: So that explicitly from (23) and (24) we have
Let be the given (constant) initial values.
By some explicit computation we can easily show that the following.
Proof. Let us first solve the simplest equation: from (26)3 we have From (26)1 it is and deriving (26)2 By replacing with the previous expression we get that is This is a linear nonhomogeneous second-order (harmonic) equation whose solution is With this function, from the expression we can compute also : that is, The solution is Concerning (26)4,5,6,7 we can notice that it is enough to solve (26)4,7 since their solutions are formally equal to the solutions of (26)5,6 since these equations coincide with (26)4,7 apart from the substitutions: Thus from (26)7 it is that is, If we put this expression in (26)4 we get and by some manipulation we get the (vectorial) linear second order harmonic equation for : The solution is so that by integrating (41) and using the previous equation, we can easily get the expression for With similar computations we get also the last two equations of (28).
It should be noticed that the the projection of the geodesics (28) into different planes gives rise to well known curves. For instance, in the plane (28) are the parametric equations of a cycloid. Moreover, by assuming that all the initial values are positive, the asymptotic limits give which are in agreement with physical consideration especially for the entropy which is upper bounded. Analogously we have similar bounded asymptotic limits for the vectorial functions.
5. Experimental Approach to the Linear Response Theory
In this section we study another aspect of the transversal waves propagation in viscoanelastic media, and we apply the theoretical results to a polymeric material as the polyisobutylene.
We consider, with respect to the Cartesian orthogonal axes (), the following displacement being and the complex wave number, so that is the phase velocity and is connected with the attenuation of the waves.
As , from (1) one obtains where The complex shear velocity is  where , (storage modulus), and (loss modulus) are, respectively, linked with the nondissipative and dissipative phenomena, and their experimental curves are plotted in Figure 1.
By putting where and , (52) becomes For the Polyisobutilene we have  the characteristic values which are and the graphics confirm (see Figure 2) with experimental data the validity of the model proposed for viscoanelastic phenomena in continuous media.
From the viewpoint of nonholonomic irreversible thermodynamics, it is shown that the best specific entropy is an harmonic function in a 27-dimensional manifold. The differential equations of geodetics are obtained and the corresponding curves are explicitly computed. In the linearized theory it is shown that the theoretical results are in agreement with the experimental data in the case of polymeric material (Polyisobutilene) (Figure 2).
Conflict of Interests
The authors declare that there is no conflicts of interests regarding the publication of this paper.
- C. Caratheodory, Undersuchungen über die Grunlagen der Thermodynamik, vol. 2, Gesammelte Mathematische Werke, Munuch, Germany, 1995.
- R. Hermann, Geometry, Physics, and Systems, vol. 18, Marcel Dekker, New York, NY, USA, 1973.
- R. Mrugała, “Submanifolds in the thermodynamic phase space,” Reports on Mathematical Physics, vol. 21, no. 2, pp. 197–203, 1985.
- C. Udrişte, O. Dogaru, and I. Ţevy, Extrema with Nonholonomic Constraints: Selected Papers, vol. 4 of Monographs and Textbooks, Geometry Balkan Press, Bucharest, Romania, 2002, Selected papers.
- C. Stamin and C. Udrişte, “Nonholonomic geometry of gibbs contact structure,” UPB Scientific Bulletin A, vol. 72, no. 1, pp. 153–170, 2010.
- J. Merker and M. Krüger, “On a variational principle in thermodynamics,” Continuum Mechanics and Thermodynamics, 2012.
- G. D. C. Kuiken, Thermodynamics of Irreversible to Diffusion and Rheology, Wiley Tutorial Series in Theoretical Chemistry, John Wiley & Sons, Chichester, UK, 1994.
- G. A. Kluitenberg and V. Ciancio, “On linear dynamical equations of state for isotropic media—I: general formalism,” Physica A, vol. 93, no. 1-2, pp. 273–286, 1978.
- V. Ciancio and G. A. Kluitenberg, “On linear dynamical equations of state for isotropic media—II. Some cases of special interest,” Physica A, vol. 99, no. 3, pp. 592–600, 1979.
- V. Ciancio, An Introduction to Thermomechanical of Continuous Media. Rheology, vol. 10 of Monographs and Textbooks, Geometry Balkan Press, Bucharest, Romania, 2009.
- A. Ciancio, “An approximate evaluation of the phenomenological and state coefficients for visco-anelastic media with memory,” UPB Scientific Bulletin A, vol. 73, no. 4, pp. 3–14, 2011.
- N. G. McCrum, B. F. Read, and G. Williams, Anelastic and Dielectric Effects in Plymeric Solids, John Wiley & Sons, London, UK, 1967.
- P. Hedvig, Dielectric Spectroscopy of Polymers, Adam Hilger, Bristol, UK, 1977.
- W. Graham, Comprehensive Polymer Science: the Synthesis, Characterization, Reactions & Applications of Polymers, vol. 2, chapter 16, edited by G. Allen and J.C. Bevington, Pergamon Press, 1989.
- I. M. Ward and D. W. Hadley, Mechanical Properties of Solid Polymers, John Wiley & Sons, 1993.