Ordered Rate Constitutive Theories for Non-Classical Thermofluids Based on Convected Time Derivatives of the Strain and Higher Order Rotation Rate Tensors Using Entropy Inequality

This paper considers non-classical continuum theory for thermoviscous fluids without memory incorporating internal rotation rates resulting from the antisymmetric part of the velocity gradient tensor to derive ordered rate constitutive theories for the Cauchy stress and the Cauchy moment tensor based on entropy inequality and representation theorem. Using the generalization of the conjugate pairs in the entropy inequality, the ordered rate constitutive theory for Cauchy stress tensor considers convected time derivatives of the Green’s strain tensor (or Almansi strain tensor) of up to orders nε as its argument tensors and the ordered rate constitutive theory for the Cauchy moment tensor considers convected time derivatives of the symmetric part of the rotation gradient tensor up to orders nΘ. While the convected time derivatives of the strain tensors are well known the convected time derivatives of higher orders of the symmetric part of the rotation gradient tensor need to be derived and are presented in this paper. Complete and general constitutive theories based on integrity using conjugate pairs in the entropy inequality and the generalization of the argument tensors of the constitutive variables and the representation theorem are derived and the material coefficients are established. It is shown that for the type of non-classical thermofluids considered in this paper the dissipation mechanism is an ordered rate mechanism due to convected time derivatives of the strain tensor as well as the convected time derivatives of the symmetric part of the rotation gradient tensor. The derivations of the constitutive theories presented in the paper is basis independent but can be made basis specific depending upon the choice of the specific basis for the constitutive variables and the argument tensors. Simplified linear theories are also presented as subset of the general constitutive theories and are compared with published works.


Introduction
The conservation and balance laws for non-classical thermofluids incorporating internal rotation rates due to antisymmetric part of the velocity gradient tensor have been presented by . These consist of usual conservation and balance laws of classical continuum mechanics: conservation of mass (CM), balance of linear momenta (BLM), balance of angular momenta (BAM), first and second laws of thermodynamics (FLT, SLT), but require some modification and/or re-derivations and additionally a new balance law, balance of moment of moment (BMM), is needed due to presence of new physics associated with the rotation rates. Yang et al. [5] derived BMM balance law based on static considerations. Surana et al. [6] showed that a balance law must be based on rate considerations and presented a consistent derivation in Lagrangian and Eulerian descriptions.
Surana et al. [3] presented an ordered rate constitutive theory for the Cauchy stress tensor using convected time derivatives of up to orders n of the Green's or Almansi strain tensors. However, the constitutive theory for the Cauchy moment tensor was limited to the convected time derivative of order one of the symmetric part of the rotation gradient tensor. In this paper convected time derivatives of the symmetric part of the rotation gradient tensor of up to orders n Θ are derived. These are then utilized as argument tensors of the Cauchy moment tensor in conjunction with conjugate pairs in entropy inequality. The ordered rate constitutive theories for the Cauchy stress tensor (up to orders n ε ) and the Cauchy moment tensor (up to order n Θ ) are derived based on theory of isotropic tensors using basis independent approach. The basis independent constitutive theories so derived can be easily made basis specific by basis dependent choices of constitutive variables and their argument tensors. The most general derivations of the constitutive theories for the Cauchy stress and Cauchy moment tensors are based on representation theorem and integrity established using their argument tensors. Constitutive theory for the heat vector is also derived using representation theorem and integrity based on the argument tensors of heat vector. Material coefficients are derived for all ordered rate constitutive theories presented in this paper. Simplified linear ordered rate constitutive theories are presented for the Cauchy stress tensor, Cauchy moment tensor and the heat vector.
The notations used in this paper have been used by the authors in the current literature, nonetheless some description and their use in deriving conservation and balance laws are presented in the following. Over bar is used on quantities to express quantities in the current configuration in Eulerian description, that is, all quantities with over bars are functions of current coordinatesx x x and time t. The mass density, Helmholtz free-energy density, temperature, and entropy density of the fluid in the current configuration are denotedρ,Φ,θ, andη, respectively, and all are functions ofx x x. The Cauchy stress tensor in Eulerian description in contravariant basis is denoted byσ σ σ (0) . The superscript "0" is used to signify that it is rate of order zero and the lowercase parenthesis distinguish it from the second Piola-Kirchhoff stress tensor σ σ σ [0] used in Lagrangian description. Dot on any quantity refers to the material derivative. As explained above undeformed and deformed configurations can be used in the derivatives as long as the final equations from the conservation and balance laws containx i and t and do not have displacements and strains in them as these are not considered for fluent continua.

Motivation and New Contributions in the Present Work
In order to illustrate the motivation for this work and the new contribution in this paper for non-classical fluent continua, we must begin with thermoviscous fluids based on classical continuum mechanics (CCM). It is well known that for such fluids the commonly used constitutive theories for deviatoric Cauchy stress tensor dσ σ σ is Newton's law of viscosity. From the CCM point of view this constitutive theory in Equation (1) is a much simplified form of the constitutive theory based on integrity derived using representation theorem in which symmetric part of the velocity gradient tensor and temperature are argument tensors of the deviatoric Cauchy stress tensor in Equation (2): We shall show that dσ σ σ in this case is basis independent. It is easy to show [7] That is, [D] is basis independent, hence dσ σ σ in Equation (1) which is entirely dependent on [D] is also basis independent. Constitutive theory in Equation (1) is what is used in fluid mechanics writings.
It can be shown that for finite strain rates the constitutive theory resulting from Equation (4) is valid, where as the constitutive theory due to Equation (5) is only valid for small strain rates [7]. Inclusion of γ γ γ (i) ; i = 2, 3, . . . , n as argument tensors of deviatoric Cauchy stress dσ σ σ (0) makes it basis dependent.

Remark 1.
(1) Inclusion of γ (i) ; i = 2, 3, . . . , n as additional argument tensors of deviatoric Cauchy stress dσ σ σ (0) may appear non-justifiable based on the conjugate pair tr ([σ] [L]) = tr ([σ] [D]) in the entropy inequality (used in fluid mechanics writings), but if we realize that the rate of work used in the derivation of energy equation only considers velocity and force, then it is not surprising to see the absence of other conjugate terms that would support inclusion of γ (i) ; i = 2, 3, . . . , n as additional argument tensors of dσ σ σ (0) .
(2) From continuum mechanics point of view, we do not seek for applications of the constitutive theory based on γ (i) ; i = 2, 3, . . . , n but rather seek if their inclusion as additional argument tensors of dσ σ σ (0) would lead to a consistent constitutive theory which it does.
(3) This work described above has been published by Surana et al. in various papers, some are cited in this paper. (4) A clear explanation of the use of higher order convected time derivatives of the strain tensors as argument tensors of deviatoric Cauchy stress for thermoviscous fluids based on CCM and consistency of the resulting constitutive theory and the rationale supporting such derivations is essential in explaining the use of these concepts in non-classical fluent continua incorporating internal rotation rates. (5) We remark that the presence of [D] or γ (1) in the constitutive theory for dσ σ σ (0) is responsible for dissipation mechanism in the fluent continua. Thus, each γ (i) ; i = 1, 2, . . . , n contributes to dissipation mechanism, hence the constitutive theory could be referred to as ordered rate constitutive theory for dσ σ σ (0) . (6) In classical continuum mechanics the constitutive theories for polymeric fluids makes use of convected time derivatives of the strain tensor [8]. The Maxwell model only uses first convected time derivative of Green's strain tensor, however the Oldroyd-B model uses first as well as second convected time derivatives of the Green's strain tensor for deriving constitutive theory for the deviatoric Cauchy stress tensor. This suggest that a constitutive theory for non-classical polymeric fluid for Cauchy moment tensor could possibly benefit from the higher order convected rotation rates derived in this paper. This work is currently in progress.
The non-classical continuum mechanics considered here for thermofluids incorporates antisymmetric part of the velocity gradient tensor (rotation rates) in the conservation and balance laws. This physics is neglected in classical continuum mechanics. The resistance offered to these by the deforming continua results in moments. One could show that upon using Cauchy principle, we can derive the concept of Cauchy moment tensor (contravariant or covariant). The derivation of balance based on this new physics shows that Cauchy moment tensor is symmetric [3] and that Cauchy moment tensor and the symmetric part of the gradients of rotation rate tensor are rate of work conjugate (entropy inequality). The similarity of these conjugate quantities with the conjugate pair dσ (0) and [D] is obvious. Replacing [D] = γ (1) with γ (i) ; i = 1, 2, . . . , n as the argument tensors of dσ σ σ (0) suggests, that perhaps we must also consider a parallel and consistent derivation of the constitutive theory for Cauchy moment tensor in which higher order convected time derivatives of rotation rate tensors are considered as its argument tensors. The first concern is that will this lead to a consistent constitutive theory? We show in the paper that it does. The second concern may be, where do we use this theory. We make some comments in the following to address these.

Comments
(1) Just like dσ σ σ (0) and each γ (i) ; i = 1, 2, . . . , n provides a mechanism of dissipation in CCM, Cauchy moment tensor and each convected time derivative of the rotation rate tensor will also result in is additional mechanism of dissipation that is absent in classical mechanics. (2) We shall observe that the conjugate pairs in the entropy inequality provide limited information regarding the argument tensors of the constitutive variables. Initial selection of the argument tensors of the constitutive variables based on conjugate pairs in the entropy inequality is augmented in the present work to include additional desired physics based on the higher order convected time derivatives of the Green's strain tensor for the deviatoric Cauchy stress tensor and higher order convected rotation rates for Cauchy moment tensor. Inclusion of higher order convected time derivatives in both deviatoric Cauchy stress tensor and the Cauchy moment tensor provide hierarchical and ordered rate dissipation mechanism. (3) Do we have applications that require constitutive theory for Cauchy moment tensor with higher order convected time derivatives of the rotation rate tensor? Maybe not at the moment, but the motivation for the present work is to derive more complete and consistent theories for the future that allow more complex physics than what we may be used to at the present. This is the motivation for the work presented here. (4) This paper presents new constitutive theory for non-classical fluent continua, applications are not the thrust of the paper.

Preliminaries, Rotation Rates, Gradients of Rotation Rates, Convected Time Derivatives of Strain Tensor and Rotation Gradient Tensor
In the following, a brief explanation of notations is necessary as some of the notations are new. x i andx i denote the position coordinates of a material point in the reference and current configurations, respectively, in a fixed frame (x-frame).x Now we can define [W]: The non-classical continuum theory in this paper considers [L] in its entirety, that is, both [D] and [W] and their convected time derivatives are considered in the derivation of the conservation and balance laws and the constitutive theories.

Conservation and Balance Laws
The conservation and balance laws for non-classical continuum mechanics in Eulerian description have been presented in References [1-3]. As shown by Surana et al. [7], the contravariant Cauchy stress tensorσ σ σ (0) is physical as it is derived using contravariant components of stresses acting on the faces of the deformed tetrahedron in the current configuration. This choice necessitates that we use convected time derivatives of the Green's strain tensor as rate of work conjugate quantities toσ σ σ (0) . On the other hand if we decide to choose covariant Cauchy stress tensorσ σ σ (0) as a stress measure, then the convected time derivatives of Almansi strain tensor must be used as rate of work conjugate quantities toσ σ σ (0) . It is well known (see Reference [7]) that the choice ofσ σ σ (0) as a stress measure is only meaningful when x i ∼ = x i . This precludes its use for finite strain rate processes. Jaumann stress measure (0)σ σ σ J that is average ofσ σ σ (0) andσ σ σ (0) obviously suffers from the same limitations asσ σ σ (0) . We shall show in a later section that we have the following conjugate pairs: in which γ (i) , γ (i) are convected time derivatives of the Green's and Almansi strain tensors and (i) γ J are Jaumann rates. These conjugate pairs are clearly basis dependent. That is, if we choose a stress tensor in a preferred basis then the basis for the strain measure and its convected time derivatives is automatically established based on conjugate pair(s). Instead of choosing conjugate pairs from Equation (14) we consider the basis independent measures in the derivation of the constitutive theories. Let the following be conjugate pair(s): Similarly, for the Cauchy moment tensor we consider similar basis independent conjugate pairs in the derivation of the constitutive theory for the Cauchy moment tensor: in which (i) γ Θ ; i = 1, 2, . . . , n θ are the convected time derivatives of the rotation gradient tensor. The derivations of these is presented in a following section. Choosing desired pairs from Equation (14) and using them in Equation (15) will yield the constitutive theory for Cauchy stress tensor in the desired basis. The same will hold for Equation (16). We summarize the conservation and balance laws for non-classical continuum mechanics (NCCM) based on internal rotation rates in whichē,η,θ,q q q,ḡ g g, (1) γ γ γ Θ , and (1)Θ Θ Θ are specific internal energy, entropy density, absolute temperature, heat flux, temperature gradients, gradients of the first convected internal rotation rates, and internal rotation rates, respectively.

Remark 2.
(1) Equations (17)-(21) are a system of eight partial differential equations (CM (1), BLM (3), BAM (3), (6)), thus we need fifteen equations for the mathematical model defined by Equations (17)- (21) to have closure. These are obtained from the constitutive theory. We note thatφ,ē, and,η are not dependent variables as these can be expressed in terms of the variables already considered. (2) While in fluid mechanics and gas dynamics the entropy inequality seems to play a very limited role based on the writing on the subject, primarily in inferring whether a process is irreversible or not or whether the calculated solutions violate the SLT, which could happen if the calculated solutions are in error. In view of the fact that all fluids must possess viscosity, irreversibility in fluid mechanics and gas dynamics is intrinsic hence is a given fact. Thus, it appears that the role of entropy inequality is rather peripheral.
In continuum mechanics (classical or non-classical) the role of entropy inequality is rather central in the development of the constitutive theories: (i) determination of constitutive variables and the initial determination of their argument tensors based on the conjugate pairs in one of the most significant and crucial conclusions we are able to draw from the entropy inequality. This gives us a starting point in the derivation of the constitutive theories. (ii) Helmholtz free energy density in the entropy inequality facilitates constitutive theories for equilibrium stress or total stress for elastic solid continua depending upon the physics considered. (iii) The role of Helmholtz free energy density and the rate of work conjugate pairs in entropy inequality are vital in the derivation of the constitutive theories for solid continua as well as fluent continua.
(3) We remark that once the constitutive theories are derived using the conjugate pairs in the entropy inequality, the entropy inequality condition is identically satisfied by the solution obtained using the balance laws incorporating these constitutive theories. When using numerical methods such as finite element method L 2 -norm of the residuals over the discretization [9,10] can be used to ensure that the calculated solution is indeed the solution of the mathematical model consisting of the conservation and balance laws incorporating constitutive theories derived using entropy inequality. Hence, every calculated solution with sufficient accuracy always satisfies SLT provided the conditions resulting from the conjugate pairs are not violated in the derivation of the constitutive theories. (4) Thus, in continuum mechanics the entropy inequality plays a central role as without it thermodynamically constitutive theories are not possible. FLT and SLT need to be recast in slightly different forms so that the conjugate pairs helpful in deriving constitutive theories can be established.

Constitutive Theories
First we consider some preliminary material to present alternate forms of some terms in the FLT and the SLT. Let (1) γ Θ represent the gradient of the rotation rate tensor (basis independent), then We note that (1) γ Θ is also the first convected time derivative of the rotation gradient tensor. We consider decomposition of (0)σ , (0)m , [L], and (1) γ Θ into symmetric and antisymmetric tensors (indicated by back subscripts s and a): We note the following: We can also write Equation (25) as Similarly, We can also write Equation (28) as With simple calculations we can show that Using Equations (26), (29), and (30) in FLT and SLT (Equations (21) and (22)) we can obtain Equations (17)- (20) and (31), (32) constitute the complete mathematical model in Eulerian description.

Constitutive Variables and Their Argument Tensors
From Equations (17)- (20) and (31), (32), we determine thatφ,η, sm m m, andq q q are constitutive variables. At this stage the argument tensors of , and the conjugate pairq iḡī θ . From conservation of mass in Lagrangian description ρ 0 = |J| ρ, that is, compressibility is due to |J| = ρ 0 ρ , hence it is fitting to consider 1 ρ as an argument tensor in Eulerian description as opposed toρ of the constitutive variables. At a later stage dependence on 1 ρ can be replaced by dependence onρ (using calculus). Furthermore we considerθ as argument tensor of the constitutive variables (thermoviscous fluid). Thus, we can write q q q =q q q 1 ρ ,ḡ g g,θ .
However, we have not used principle of equipresence for (0)σ σ σ, (0)m m m, andq q q as the conjugate pairs specifically dictate their argument tensors as used in Equations (33)-(35).

Convected Time Derivatives of Strain Tensors and Convected Rotation Rates
The argument tensors of the constitutive variables in Equations (33)-(37) can be enhanced to permit more comprehensive constitutive theories. Convected time derivatives of the Green's strain tensor and the Almansi strain tensor have been presented in Reference [7]. We summarize these in the following.
(1) We note that Thus the first convected time derivatives of ε [0] and ε [0] are basis independent.

Gradients of the Higher Order Convected Rotation Rates
To our knowledge, the gradients of the higher order convected rotation rates presented in the following is its first presentation. Consider and∇ ∇ ∇ ×v v v = e e e 1 (Θ x 1 ) (1) + e e e 2 (Θ x 2 ) (1) + e e e 2 (Θ x 3 ) (1) .
Thus, we have in which γ Θ (1) and γ Θ (1) are gradients of the rotation rates (of order one) in co-and contra-variant bases. The symmetric and the antisymmetric components of γ Θ (1) and γ Θ (1) are given by and We can define the second Piola-Kirchhoff tensors corresponding to s γ Θ (1) and s γ Θ (1) by Following the details in Reference [7] and Sections 4.2.1 and 4.2.2 we can write the following for s γ Θ (k) ; k = 2, 3, . . . , n Θ , the symmetric parts of the gradient of the higher order convected rotation rates in covariant basis Likewise, for s γ Θ (k) ; k = 2, 3, . . . , n Θ , the symmetric parts of the gradient of the higher order convected rotation rates in contravariant basis can be defined by Symmetric parts of the gradient of the convected rotation rate tensors called Jaumann tensors can be defined as In order to make the derivation of the constitutive theories basis independent we use the notation (k) s γ Θ ; k = 1, 2, . . . , n Θ for the symmetric parts of the gradient of convected rotation rates. Depending upon the choice of the basis for the Cauchy moment tensor, that is, contra-or co-variant or Jaumann basis we can choose the corresponding conjugates as s γ Θ (k) , s γ Θ (k) , and (k) s γ J Θ , respectively.

Entropy Inequality and Constitutive Theories
Consider entropy inequality Equation (32) and the argument tensors ofφ in Equation (60). Using the chain rule . φ = Dφ Dt can now be obtained: From the continuity equation Using Equation (62) in Equation (61) We note that ∂φ Using Equation (64) in Equation (63) .

Incompressible Matter
For incompressible matter density is constant, henceρ = ρ 0 , ∂φ ∂ρ = 0, and e sσ σ σ for incompressible matter cannot be derived using ∂φ ∂ρ , that is, using Equation (84). The problem arises due to the fact that the incompressibility condition is not present in the entropy inequality, but can be incorporated in it by simply adding the following to the entropy inequality Equation (83)p θ (1) γ lk δ lk = 0, in whichp θ is an arbitrary Lagrange multiplier. This yields: Using ∂φ ∂ρ = 0 and regrouping terms in Equation (89), Since e  p θ is called mechanical pressure. Sincep θ is an arbitrary Lagrange multiplier,p θ cannot be determined from the deformation field. The entropy inequality reduces to the same as in the compressible case, that is, Equation (86).

Constitutive Theories for d
q q q =q q q ρ,ḡ g g,θ .
The final form of the entropy inequality is given bȳ

Remark 6.
(1) The constitutive theory Equation (104) is based on integrity, hence utilizes complete basis, but unfortunately requires determination of too many material coefficients. This constitutive theory utilizes convected time derivatives of the strain tensor up to orders n ε . (2) We can derive simplified constitutive theory from Equation (104) by only retaining the desired terms on the right hand side of Equation (104). In the following we consider a constitutive theory that is linear in convected time derivatives of the strain tensor in which higher degree and all product terms are neglected. Introducing more familiar notation for material coefficients: This constitutive theory contains each strain rate and its trace and is obviously nonlinear due to the fact that (i) γ ; i = 2, 3, . . . , n ε are nonlinear functions of the components of the velocity gradient tensor. Each convected time derivative of the strain tensor requires two material coefficients (µ i , λ i ).
This constitutive theory is basis independent as γ (1) = γ (1) = [D] = (1) γ J , hence d We substitute m α i ; i = 0, 1, . . . , N m in Equation (113) and collect coefficients (those defined in Ω configuration) of the terms defined in the current configuration and introduce new coefficients for the terms defined in the known configuration Ω. Then, we obtain the following: In Equation (116), m a j , m b i , m c ij , m d i , and m α tm ; i = 1, 2, . . . , N m ; j = 1, 2, . . . , M m are material coefficients defined in the known configuration Ω. The material coefficients are functions ofρ| Ω , θ Ω , and m I j Ω ; j = 1, 2, . . . , M m . This constitutive theory requires (M m + 2N m + N m M m + 1) material coefficients.

Remark 7.
(1) The constitutive theory Equation (116) is based on integrity, hence utilizes complete basis, but unfortunately requires determination of too many material coefficients. This constitutive theory utilizes symmetric part of the gradients of the convected rotation rates up to orders n Θ . (2) It is possible to derive simplified constitutive theory for s γ Θ ; i = 1, 2, . . . , n Θ in which the higher degree terms of (i) s γ Θ as well as all product terms are neglected: This constitutive theory contains sm m m is given by: in which tr (i) s γ Θ = 0, hence in this constitutive theory only one material coefficient is associated with the rate of order one. (4) It is instructive to consider a constitutive theory for Hence,

If we choose
(1) then the resulting constitutive theory is given by Equation (116) with N m = 2, M m = 3. The generators and invariants are defined in Equation (122).

Constitutive Theory forq q q
Consider the argument tensors ofq q q in Equation (96), that is, q q q =q q q ρ,ḡ g g,θ .
Using representation theorem (following Reference [7]) we can derive the following constitutive theory, keeping in mind that here the generator isḡ g g and the invariant isḡ g g ·ḡ g g: q q q = −kḡ g g − k 1 ḡ g g ·ḡ g g − (ḡ g g ·ḡ g g)| Ω ḡ g g − k 2 θ −θ Ω ḡ g g, in which k, k 1 , and k 2 can be functions ofρ,θ Ω , and (ḡ g g ·ḡ g g)| Ω . A linear constitutive theory forq q q is given byq q q = −kḡ g g.
This is the Fourier heat conduction law.

Complete Mathematical Model
The equations resulting from the conservation and balance laws and constitutive theories, based on integrity and simplified are presented in this section. It is shown that this mathematical model has closure. ∂ρ ∂t Constitutive theories for e sm , andq q q: (a) Based on integrity: q q q = −kḡ g g − k 1 ḡ g g ·ḡ g g − (ḡ g g ·ḡ g g)| Ω ḡ g g − k 2 θ −θ Ω ḡ g g.

Summary and Conclusions
In this section, we present a summary of the work presented in the paper, highlight new and meritorious features and present some conclusions.
(1) Role of entropy inequality in the derivation of the constitutive theories for Cauchy stress tensor, Cauchy moment tensor, and heat vector are discussed in Section 3, remark (2). (2) The paper considers the conservation and balance laws of non-classical mechanics [1][2][3] in Eulerian description in which internal rotation rates naturally arise due to anti-symmetric part of the velocity gradient tensor. (3) Surana et al. [6] have shown that in this NCCM, the conservation and balance laws of CCM need modifications and a new balance law 'balance of moment of moments' is essential to ensure dynamic equilibrium of the deforming fluent continua. (4) It has been shown [3] that, in this physics, a more complete constitutive theory for the symmetric part of the Cauchy stress tensor can be derived using convected time derivatives of the Green's or Almansi strain tensor. However, the constitutive theories for the symmetric Cauchy moment tensor have been limited to the symmetric part of the gradients of rotation rates. (5) In this paper, we show that parallel to the convected time derivatives of the strain tensors (up to orders n ε ), it is also possible to establish symmetric gradient tensors of the higher order rotation rates. The derivation presented in this paper derives a recursive form through which rotation rates of any order (n Θ ) can be obtained. (6) The constitutive theories for Cauchy stress and Cauchy moment tensors are derived by: (i) first establishing their argument tensors using conjugate pairs in the entropy inequality; (ii) then augmenting these by the higher order convected time derivatives; (iii) consideration of the argument tensors ofφ using principle of equipresence and using . φ in the entropy inequality; (iv) decompositions of Cauchy stress and Cauchy moment tensors into equilibrium and deviatoric parts.
(7) Constitutive theory for equilibrium stress tensor is established for compressible and incompressible fluent continua. (8) Constitutive theories for deviatoric part of the symmetric Cauchy stress and Cauchy moment tensor and heat vector are derived using representation theorem. These theories are based on integrity (complete basis), hence are complete and are naturally nonlinear in terms of their argument tensors. Simplified constitutive theories are also presented.