An approach to constructing models of multiphase elastic porous media

A novel approach to describing the behaviour of multiphase elastic porous media is proposed. The average values of the physical quantities needed to describe the motions of porous media are formulated using an integral relation. The validity of this relation is taken as the fundamental hypothesis. The integral definition of the average values enables integral relations to be devised for the average values from the integral laws of conservation of mass, momentum and energy and the increase in entropy. Along with the average values, the integral relations contain new variables that can be identified with generalized thermodynamic forces, which can be used to take into account the phase interaction in a porous medium. The integral relations are used to derive differential equations for the rate of entropy change and Gibbs relations for a porous medium as a basis for obtaining the constitutive relations. Relationships between the thermomechanical parameters of the model are established from the Gibbs relations under additional assumptions. The equation for the rate of entropy change can be used to establish relations between the generalized thermodynamic forces and fluxes. A complete system of differential equations in the defining parameters, which describes the motion of multiphase elastic porous media, is finally obtained.     

非饱和土本构关系的混合物理论(Ⅰ)——非线性本构方程和场方程

以混合物理论为基础建立了非饱和土非线性本构方程和场方程.把非饱和土作为3种组分构成的饱和混合物来研究.首先根据土力学成果提出了非饱和土混合物的基本假设,推导出适用于非饱和土混合物的熵不等式;然后采用混合物理论处理本构问题的常规方法得出了非饱和土非线性本构方程;最后把非线性本构方程代入混合物组分动量守恒定律,获得了非饱和土各组分运动的非线性场方程;并且给出了非饱和土混合物的能量守恒方程,从而形成了解决非饱和土混合物热力学过程的完备方程组.     

Boundary value problems in the theory of thermoporoelasticity for materials with double porosity

In this paper the linear quasi-static theory of thermoelasticity for solids with double porosity is considered. The system of equations of this theory is based on the equilibrium equations for solids with double porosity, conservation of fluid mass, constitutive equations, Darcy's law for materials with double porosity and Fourier's law for heat conduction. The basic internal and external boundary value problems (BVPs) of steady vibrations are formulated. The uniqueness and existence theorems for classical solutions of the above mentioned BVPs are proved by means of the potential method (boundary integral equation method) and the theory of singular integral equations. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of variable viscosity and viscous dissipation on the Hagen‐Poiseuille flow and entropy generation

In this article, the steady‐state flow of a Hagen‐Poiseuille modelin a circular pipe is considered and entropy generation due tofluid friction and heat transfer is examined. Because of variationin fluid viscosity, the entropy generation in the flow varies. Inhis model, Arrhenius law is applied for temperature equation‐dependent viscosity, and the influence of viscosity parameters on the entropy generation number and distribution of temperature and velocity is investigated. The governing momentum and energy equations, which are coupled due to the dissipative term in the energy equation, were solved by analytical techniques. The solutions of equations via perturbation method and homotopy perturbation method are obtained and then compared with those of numerical solutions. It is found that the fluid viscosity influences considerably the temperature distribution in the fluid close to the pipe wall, and increasing pipe wall temperature enhances the rate of entropy generation. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 529–540, 2011     

Entropy equation and absolute temperature

In this paper we consider the equivalence between the heat and the entropy balance laws. These two equations are related by an integrating factor, which defines the absolute temperature. This result is obtained applying the thermodynamic laws to a perfect fluid. So that, by means of the entropy equation we introduce the Second Law of Thermodynamics. Two particular cases of constitutive equations, both for the internal energy and for the heat flux, are considered and their corresponding differential equations, useful to study the behaviour of these materials, are also given.     

Thermodynamics of hypoelasticity

A thermodynamic framework for hypoelasticity is constructed based on a modified Gibbs Function which depends on the stress deviator and specific volume. In the hypoelastic equations considered here, the stress deviator is obtained from a rate equation involving the corotational derivative and coefficients which depend on the invariants of the stress deviator as well as specific volume. Non-negative entropy production is enforced, leading to a non-dissipative condition in the sense that entropy is produced only by heat transfer. The resulting conditions lead to relations among the coefficients in the rate equation. These relations underdetermine the coefficients, so that there is leeway in specifying these coefficients. An example of a set of rate equations which satisfy the thermodynamics is presented.     

The thermodynamically consistent equations of a thermoelastic saturated porous medium

A phenomenological model of a porous medium saturated with fluid is considered with in the framework of the hypothesis of interpenetrating continua. Assuming that there are no phase transitions, that the contribution of pulsations to the stress tensor and kinetic energy is small, and the components of the medium are in thermal equilibrium, mass, momentum and energy equations and a law of conservation of compatibility of the deformations and velocities are formulated. Using a representation of the force of interaction of the components in the form of the sum of equilibrium and dissipative components, a new form of inequality is obtained for the rate of entropy production. A definition of a thermoelastic saturated porous medium is given. The symmetry group of such a medium is considered as a set of two groups, corresponding to the symmetry of the skeleton and the fluid. It is shown that, in the class of thermoelastic porous media with an arbitrary type of symmetry of the skeleton, the saturating fluid can only be an ideal fluid, while the thermodynamic potentials and the porosity, stresses and entropies determined by them do not depend on the temperature gradient and the relative fluid velocity. It is found that the condition of incompressibility of only one of the components of the medium leads to the elimination of the porosity from the governing relations, rather than to kinematic limitations. The limitations imposed on the governing relations by the principle of thermodynamic consistency and the requirement of independence of the choice of the frame of reference are investigated. A form of the governing relations, necessary and sufficient to satisfy these principles, is obtained. It is shown that the Biot equations are one of the forms of thermodynamically consistent governing relations. A thermodynamic validation of the effective-stress tensor is given.     

Micromechanical Analysis of Inelastic Deformation of Unidirectional Fibrous Composites under Multiaxial and Shear Loadings

An algorithm is proposed for numerically solving nonlinear 3D problems of micromechanics of a unidirectionally reinforced composite with a regular structure. For the matrix, equations of the deformation theory of plasticity and relations of reduced rigidity in its failure zones are used, whereas the fibers are elastic and indestructible. According to the method of local approximation, fields of microstresses and microstrains are determined in a structural fragment containing nine periodic cells. Boundary conditions of the fragment correspond to an arbitrary combination of longitudinal, transverse, and shear microstresses occurring in the central part of the fragment. The solution to the nonlinear 3D problem is sought by the method of superposition with an iterational refinement based on the successive solution of an antiplane problem and a problem on a generalized plain strain state of the structural segment. Special features of the iteration procedure are considered. The calculated deformation diagrams and ultimate strengths of a unidirectional glass-epoxy composite are presented for several loading trajectories.     

Lp-Lq time decay estimate to the solution of the Cauchy problem of the system of equations for nonlocal model of hyperbolic thermoelasticity theory

The aim of this paper is to present a new system of equations describing nonlocal model of hyperbolic thermoelasticity theory. We used the Papkin and Gurtin approach based on the constitutive relations for internal energy e(x), and heat flux q(x), with integral terms. Such system of equations describes the propagation of thermal perturbation with finite velocity. Using the modified Cagniard–de Hoop's method we constructed the matrix of fundamental solutions for this system of equations in three–dimensional space. Basing on the constructed matrix of fundamental solutions in the explicit formula we represent the solution of the Cauchy problem to this system of equations in the form of some kind of convolutions. Next, applying the method of Sobolev spaces, we obtain the L^p−L^q time decay estimate to the solution of the Cauchy problem. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework

We provide a thermodynamic basis for the development of models that are usually referred to as ??phase-field models?? for compressible, incompressible, and quasi-incompressible fluids. Using the theory of mixtures as a starting point, we develop a framework within which we can derive ??phase-field models?? both for mixtures of two constituents and for mixtures of arbitrarily many fluids. In order to obtain the constitutive equations, we appeal to the requirement that among all admissible constitutive relations that which is appropriate maximizes the rate of entropy production (see Rajagopal and Srinivasa in Proc R Soc Lond A 460:631?C651, 2004). The procedure has the advantage that the theory is based on prescribing the constitutive equations for only two scalars: the entropy and the entropy production. Unlike the assumption made in the case of the Navier?CStokes?CFourier fluids, we suppose that the entropy is not only a function of the internal energy and the density but also of gradients of the partial densities or the concentration gradients. The form for the rate of entropy production is the same as that for the Navier?CStokes?CFourier fluid. As observed earlier in Heida and Málek (Int J Eng Sci 48(11):1313?C1324, 2010), it turns out that the dependence of the rate of entropy production on the thermodynamical fluxes is crucial. The resulting equations are of the Cahn?CHilliard?CNavier?CStokes type and can be expressed both in terms of density gradients or concentration gradients. As particular cases, we will obtain the Cahn?CHilliard?CNavier?CStokes system as well as the Korteweg equation. Compared to earlier approaches, our methodology has the advantage that it directly takes into account the rate of entropy production and can take into consideration any constitutive assumption for the internal energy (or entropy).     

Thermoelasticity with thermal relaxation: An alternative formulation

The theory of thermoelasticity with thermal relaxation for homogeneous materials is formulated upon the basis of the law of balance of energy and the law of balance of entropy, proposed by Green and Naghdi [5]. The non-linear theory is formulated first; then the linearized theory is deduced. The uniqueness of solution of a typical initial, mixed boundary value problem is established.     

重建极性连续统理论的基本定律和原理(Ⅴ)——极性热力连续统

在对传统的微极热弹性理论和热压电弹性理论已进行过再研究的基础上重建极性热力连续统的较为完整的基本均衡方程和边界条件.从较为完整的虚功率原理推导出微极热弹性理论的运动方程和局部能率均衡方程.从较为完整的Hamilton原理通过全变分自然地推导出运动方程,熵均衡方程以及所有边界条件.给出的新的动量均衡方程和局部能率均衡方程与现有理论的结果存在本质的差异.通过过渡和归结可从微极热弹性理论分别得到微态热弹性理论的和偶应力热弹性动力学的结果.最后,按照上述思路直接给出微极热压电弹性理论的结果.     

The heat and moisture transfer balance theory of garment simulation

To establish the human body model to analyze the heat and moisture transfer on body surface, a new explicit definition of rational L-recursion surface is given and the L-recursion surfaces, in Grassmann spaces, are constructed by using blossom method of the homogeneous normal pyramid form. Based on our human body model, the balance theory of garment simulation, the heat and moisture transfer balance equations, called ICAD-balance equations are obtained. The balance theory of garment simulation integrally studies the complex system of human body-fabric-environment. At the same time, the method of obtaining the heat and moisture transfer balance equations is also based on the mass conservation law, the energy conservation law and the Fish law of capillarity. A finite volume method is employed to solve the ICAD-balance equations.     

Radiative processes and non-equilibrium thermodynamics

With the assumption of an elementary physical concept meteorologically effective radiative processes (absorption-emission, scattering) can be included consistently in nonequilibrium thermodynamics of irreversible phenomena. Analogously to the usual Gibbs relations a fundamental equation was formulated for monochromatic light rays as the nucleus of the theory.Using the methods of classical irreversible theory, a complete entropy balance equation is derived in which the entropy variations of the mass as well as of the radiation field are explicitly represented. The resulting entropy source strength function through its analytical structure reveals the dynamical character of the irreversible variation terms. The-expression being positive according to the second law of thermodynamics is found to have a bilinear form as a function of the irreversible fluxes representing the entropy generating radiative processes and their conjugated thermodynamic forces. The mathematical structure and the positive sign of, following the usual line of reasoning, motivate the assumption of constitutive relations for the irreversible radiative processes. These equations developed from purely thermodynamical reasoning turn out to be equivalent to the usual radiative transfer equation which is founded on a very different theoretical concept. A very fundamental relationship can be deduced in this context from the entropy production function. It provides a direct thermodynamical proof that in nonscattering media the definition of a local temperature is necessarily accompanied by the validity of the Kirchhoff law.     

The heat and moisture transfer balance theory of garment simulation

To establish the human body model to analyze the heat and moisture transfer on body surface, a new explicit definition of rational L-recursion surface is given and the L-recursion surfaces, in Grassmann spaces, are constructed by using blossom method of the homogeneous normal pyramid form. Based on our human body model, the balance theory of garment simulation, the heat and moisture transfer balance equations, called ICAD-balance equations are obtained. The balance theory of garment simulation integrally studies the complex system of human body–fabric–environment. At the same time, the method of obtaining the heat and moisture transfer balance equations is also based on the mass conservation law, the energy conservation law and the Fish law of capillarity. A finite volume method is employed to solve the ICAD-balance equations.     

Wave Dispersion Decoupling in Micropolar Thermoelasticity

For a system of field equations of micropolar thermoelasticity we derive a propagation condition for thermoelastic disturbance in a form of monochromatic plane wave in deformation, micro rotation and temperature. The corresponding dispersion relation is given in an explicit form, together with dependence of characteristic coefficients on principal material constants forming the constitutive tensor of isothermal macro and micro elasticity, phenomenological heat conductivity and coupled macroscopic thermoelasticity. It is shown that due to the centro–symmetric nature of microelasticity and particular form of temperature coupling in a free energy function, the separation between the optical and acoustical branch of dispersion relation is inherent. For such systems dispersion relations due to the micropolar fields on one side and macroscopic thermoelastic fields on the other side are completely independent without any cross–coupling. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A higher-order Godunov method for modeling finite deformation in elastic-plastic solids

In this paper we develop a first-order system of conservation laws for finite deformation in solids, describe its characteristic structure, and use this analysis to develop a second-order numerical method for problems involving finite deformation and plasticity. The equations of mass, momentum, and energy conservation in Lagrangian and Eulerian frames of reference are combined with kinetic equations of state for the stress and with caloric equations of state for the internal energy, as well as with auxiliary equations representing equality of mixed partial derivatives of the deformation gradient. Particular attention is paid to the influence of a curl constraint on the deformation gradient, so that the characteristic speeds transform properly between the two frames of reference. Next, we consider models in rate-form for isotropic elastic-plastic materials with work-hardening, and examine the circumstances under which these models lead to hyperbolic systems for the equations of motion. In spite of the fact that these models violate thermodynamic principles in such a way that the acoustic tensor becomes nonsymmetric, we still find that the characteristic speeds are always real for elastic behavior, and essentially always real for plastic response. These results allow us to construct a second-order Godunov method for the computation of three-dimensional displacement in a one-dimensional material viewed in the Lagrangian frame of reference. We also describe a technique for the approximate solution of Riemann problems in order to determine numerical fluxes in this algorithm. Finally, we present numerical examples of the results of the algorithm.     

水动力-热动力学的极值定律

本文对水动力学和更普通性的连续介体动力学中以连续方程与运动方程所表达的现有诸经典守恒定律以外,提出另一最大能量消散率定律.这一定律的推论就是应用水力学中培纶格-波丝最小储存能学说. 凡在运动中消散了的机械能皆转化成为热能,储存在物体里.能量之消散当一定时刻一定温度都使产熵增加.所以,从最大能量消散率可引出热力学第二定律的一个新概念,即机械运动的产熵率也总是一个可能的最大值. 文中建议的这个连续介体极值定律,可从变分原理推导出来,重订热力学第二定律则可藉微观分析加以证明.两者合成水动力-热动力学极值定律     

Finite element methods for micropolar models of granular materials

We present a thermodynamically based finite element scheme for rate-independent materials and demonstrate its application in modelling the rheological behaviour of granular materials. Starting from the laws of thermodynamics, we have recently developed a new class of micropolar-type constitutive relations for two-dimensional densely packed granular media. This class of constitutive laws is expressed in terms of particle-scale properties, thus providing a direct link between observed macroscopic behaviour and the underlying particle–particle interactions. Here, we demonstrate how the connection to the underlying physics can be maintained and carried through to the finite element implementation phase of the modelling process via the same thermodynamical principles used to construct the constitutive laws. Notably, the study indicates that while the traditional Galerkin-FEM method admits a range of weighting functions, the proposed formulation provides an additional constraint that narrows the choice of admissible weighting functions via the second law of thermodynamics. Additionally, this paper presents insights into the finite element implementation of micropolar models deemed to be appropriate for modelling several classes of heterogeneous media (e.g. granular materials, cellular composites and biological materials). As the kinematics and kinetics of micropolar continua are enriched by the addition of rotational degrees of freedom to each material point, the equations governing boundary value problems for such materials differ from those of other continuum models both from the viewpoint of the constitutive law and the governing conservation laws. Analysis of elastoplastic deformation of micropolar continua is presented.