Z. angew. Math. Phys. 53 (2002) 923 948
0044-2275/02/060923-26 Zeitschrift f�r angewandte
� 2002 Birkh�user Verlag, Basel Mathematik und Physik ZAMP
Modeling of the microwave drying process of
aqueous dielectrics
K. R. Rajagopal and L. Tao
Dedicated to the memory of Eugen Soós
Abstract. We propose a thermodynamic framework for describing the microwave drying pro-
cess of aqueous dielectrics based on Maxwell-Lorentz field equations and mixture theory. Several
issues are discussed such as the form of entropy equation; the constitutive relations for the macro-
scopic electric polarization vectors, Cauchy stresses, heat fluxes, internal momentum supplies,
etc., for each component of the mixture: porous solid, water and gas in different regions; and the
interfacial jump conditions between different regions in the mixture. A brief examination of the
status of material frame indifference within the context of our framework is presented.
Mathematics Subject Classification (2000). 76A99, 76S05, 76T30, 78A25, 78A40, 80A17,
80A20.
Keywords. Aqueous dielectric drying, frame indifference, Maxwell s equations, microwave heat-
ing, mixture theory.
1. Introduction
We derive most of our energy from electromagnetic radiation from our environ-
ment. The complete electromagnetic spectrum spans the gamut from exceedingly
short wave lengths associated with ł -rays of the order of 10-5 �m , to X-rays
between 10-4 �m to 10-2 �m , the ultraviolet to the visible spectrum onto in-
frared rays of wave lengths between 1 �m and 102 �m and microwave radiation
with wave lengths greater than 102 �m . This electromagnetic radiation that is
received by us is converted into energy in various forms: chemical, mechanical
and thermal, and the nature of the conversion depends on the body that is being
irradiated as well as the wave length of the electromagnetic radiation. Energy
associated with certain wave lengths is not absorbed by certain bodies (molecules)
while that associated with other wave lengths is, in fact this is the basis for time
resolved spectroscopy. In this paper we are interested in the conversion of electro-
magnetic radiation of wave lengths in the range of microwaves into energy in the
924 K. R. Rajagopal and L. Tao ZAMP
thermal form, commonly referred to as heat.
Microwave provides an effective means to dry aqueous dielectric materials due
to its ability to transfer energy directly to the interior of the materials which
generates internal heating through polarization. To date there is no comprehensive
model in place to describe this phenomenon and thus it offers an interesting and
challenging problem concerning the modeling of this drying process. The present
work attempts to propose a thermodynamic formulation regarding the modeling
within the context of Maxwell-Lorentz field equations and mixture theory under
the premise that such a continuum treatment is appropriate.
Any attempt at developing a model of the process of drying an aqueous di-
electric through microwave radiation has to come to grips with how to model the
supply of the electromagnetic radiation to the body. The mechanism by means
of which energy in the electromagnetic form is converted to heat is not yet well
understood for even reasonably simple molecules, though some understanding has
been achieved for simple monatomic gases. While dealing with continua we at
the very least need to provide a phenomenological continuum description of the
radiation received by the body, i.e., we need to provide a constitutive equation
for the radiant energy received or absorbed by a body. We realize this need for a
constitutive description for we recognize that for the class of  black bodies and
for electromagnetic radiation in a certain range of wave lengths we do have such a
specification. However, for bodies other than this special class and for electromag-
netic radiation of different wave lengths, we know little concerning the constitutive
form.
In most thermodynamic studies in classical continuum thermomechanics, ra-
diation is essentially ignored, and even if it is not ignored it is eliminated from
consideration in the entropy equation in virtue of the following practice. Let us
consider the somewhat common scheme in classical continuum thermomechanics
where the energy equation is given for a single continuum in the local form
dU
� =tr (� L) -"�q+�r,
dt
where � denotes the density, U the specific internal energy, � the Cauchy
stress, L the velocity gradient, q the heat flux vector, r the radiant heating and
" denotes the eulerian spatial gradient. Let us suppose that the second law is
given through the Clausius-Duhem inequality ([13])
ds q r
� + "� - = � e" 0, (")
dt T T
where s , T and � are, respectively, the entropy density, absolute temperature
and the internal source of the entropy besides r . We can now eliminate r from
the above equation by substituting for it from the energy equation to obtain what
is popularly referred to as the dissipation inequality, namely
ds 1 dU 1 q�"T
� = � - + tr (� L) - e" 0. ("")
2
dt T dt T �T
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 925
Next, the standard practice in continuum thermodynamics is to require that all
processes that the body can undergo meet the above inequality. This usually
places restrictions on the forms that the constitutive equations for the internal
energy, entropy, stress and heat flux can take. There are some problems with
regard to such an approach. First, (**) places no restrictions on the constitutive
representation for the radiant heating, though (*) indicates that it cannot be
arbitrary. Of course, not all forms of the radiant heating, when substituted into
the energy equation will lead to well-posed problems, but the way, in which the
second law is brought into play, namely (**), provides no restrictions with regard
to the entropy production due to radiant heating, directly. Second, it is common
to exploit (**) by assuming that (**) holds for arbitrary processes. However, it is
unreasonable to expect (**) to hold for any real material for all arbitrary processes
for the constitutive assumptions that we start with are only expected to hold for
a certain class of processes. For example, if we are trying to model rubber as an
elastic material, we cannot really subject it to arbitrary deformations, it will either
start behaving inelastically or fracture and fall apart. While mathematically we
can define an ideal elastic body, restrictions obtained from such a definition will
not have its intended meaning for a real body. To be more explicit, if (**) is
required to hold only for a class of processes, then the restrictions imposed are
far less severe and thus a larger class of models would be allowable. This point
cannot be overemphasized. Constitutive relations for bodies have meaning only
for a specified class of processes. To overcome the first difficulty we offer to include
the the radiant heating as a part of the internal source of entropy by redefining �
which will allow the second law to place rather strigent restrictions on the form
of the heating. A similar scheme has also been adopted by Eringen and Maugin
[3], though they still employ a term r/T in the entropy equation to account for
some other forms of heating, which is not necessary based on the argument offered
above. The second problem can be avoided by seeking and exploring constitutive
relations that satisfy the entropy inequality automatically. Furthermore, there are
other important issues concerning the second law that are far from clear. For
instance, the exact form of the second law for open systems is a matter of much
debate. Also, whether the second law should be enforced as a global inequality for
a body (or for a system) or whether it could be enforced locally is far from obvious
especially as one deals, in the classical context, with sources and sinks. We shall
not get into a discussion of these issues here as they are beyond the intended scope
of this work. Here, we shall assume that the second law can be enforced pointwise.
We are yet faced with the onerous task of assuming a constitutive form for
the microwave heating in an aqueous dielectric. When a body is supplied with
the radiant energy, the vibrational energy, the rotational energy or the translation
energy of the molecules or the atoms can be changed. To try to capture the
manner in which the energy in its electromagnetic form converts into the energy
in its thermal form needs fundamentally an evaluation at the quantum mechanical
level. This is not the spirit of this presentation as we are primarily interested in
926 K. R. Rajagopal and L. Tao ZAMP
trying to capture the form for the radiation from a phenomenological perspective.
Unfortunately, there is not sufficient amount of experimental data to do this. We
fill this lacuna with a judicious guess work for the radiation, following a statistical
scheme for the basis of the Maxwell-Lorentz field equations.
We are interested in describing microwave drying and we thus have to deal with
electromagnetic fields. Any attempt at incorporating electromagnetic effects nec-
essarily requires us to include Maxwell s equations into the set of motion equations.
Thus, we need to decide whether we will enforce the requirement of material frame
indifference for this problem. We are aware of the controversy surrounding this
topic (Edelen and McLennan [1]; Wang [15]) and we will not pursue a thorough
discussion of it. We require Galilean invariance, instead of the more restrictive ma-
terial frame indifference, within the framework of classical non-relativistic theory,
for the sake of consistency to put the constitutive relations on the electromag-
netic and the thermodynamic quantities on the same footing. Also, this choice is
motivated by a short analysis to be presented in the Appendix.
We should mention that we do not pursue a scheme of volume-averaging here to
help close the mathematical model of microwave drying (for example, see Whitaker
[16] for details). This type of averaging scheme requires one to know the specific
constitutive relations for the motion of each constituent as a physical continuum
in the sub-domain it occupies in the mixture. With this information, an initial-
boundary value problem can be formally set up, the interactions between the
constituents are supposedly accounted for through some interfacial jump condi-
tions, and then the volume averaging can be carried out. There is some merit
to this scheme for suggesting forms of balance type equations, but it does not
resolve the issue of furnishing the constitutive relations necessary to make the
balance equations closed due to the following reasons. First, there are molecu-
lar interactions among the constituents that can significantly change the behavior
of each constituent and we may not have comprehensive information concerning
these molecular interactions. Second, each constituent will be treated as a math-
ematical continuum occupying the whole domain of the mixture, instead of the
sub-domain the constituent itself occupies physically. This treatment makes the
constitutive relations for the constituent as a mathematically homogenized con-
tinuum differ from the constitutive relations associated with the constituent in
the sub-domain it occupies, since they involve different types of field distributions.
Additional assumptions have to be introduced to overcome this problem in the
volume-averaging scheme, which is tantamount to constitutive assumptions. Fur-
thermore, a conventional scheme of volume averaging to develop the equations of
motion has some (conceptual) difficulties of its own. One difficulty occurs when
the spatial point is near a bounded boundary where a constant representative el-
ementary volume (REV) is not applicable. (A non-constant REV could be used
with complexities.) Another problem is the prescription of boundary conditions
consistent with the volume averaging scheme. An alternative is to seek the scheme
of ensemble averaging. And this issue has been discussed rather thoroughly in the
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 927
studies of turbulence modeling.
In the following section, we shall formulate, within the context of Maxwell-
Lorentz field equations and mixture theory, the electromagnetic field equations
and the equations of motion on the porous solid, water and gas in an aqueous
dielectric for the process of microwave drying. The microwave heating makes it
possible for the presence of several mixture regions inside an aqueous dielectric
body: mixture of a porous solid and water; mixture of a porous solid, water and
gas; and mixture of a porous solid and gas. And therefore, proper constitutive
relations of the macroscopic electric polarization vectors, Cauchy stresses, heat
fluxes, internal momentum supplies and so on have to be provided to characterize
the thermomechanical behavior of these mixtures, which is carried out on the basis
of continuum thermodynamics. The interfacial jump conditions between these
different mixture regions are treated in section 3. In the appendix we discuss
briefly an issue related to the requirement of material frame indifference. The
framework presented here is tentative and much more needs to be done to fix the
specific forms of the constitutive relations based on experimental data and the
solving of specific drying problems.
2. Formulation
Let us consider the microwave drying process of an aqueous dielectric. We have
a porous dielectric solid infused with water that is heated by an electromagnetic
wave of frequency ( � ) in the range of about 300MHz to 3GHz. The mechanisms of
heating are supposedly orientation polarization and Maxwell-Wagner polarization
due to an alternating electric field (Hasted [6]; Metaxas and Meredith [9]) (these
polarizations are akin to a phenomenological description). In the initial stage
of the heating process, we have a mixture of a porous solid and water. As the
c
temperature of the mixture increases to some critical value T , water starts to
evaporate (phase change) and a mixture of three constituents, a porous solid, water
and gas (dry steam), is present. Finally water will be converted into gas (dry
steam) completely, and we have a mixture of a porous solid and gas. These three
mixtures can co-exist simultaneously in different regions of the aqueous dielectric
body due to the non-uniform motion and heating. In the following, we will denote
the quantities associated with the three constituents by the superscripts s , w ,
and g , respectively.
Introducing the electric field E and the magnetic induction field B and fol-
lowing a statistical scheme on the basis of the Maxwell-Lorentz field equations
(Eringen and Maugin [3]; de Groot and Suttorp [5]), we can obtain the Maxwell
field equations for the mixture,
"�E = -"� Pą, (1)
928 K. R. Rajagopal and L. Tao ZAMP
"E "Pą
-c-1 + "�B = c-1 + "�(Pą � vą) , (2)
"t "t
"�B =0, (3)
"B
c-1 + "�E = 0, (4)
"t
where Pą and vą are, respectively, the macroscopic electric polarization vector
and the velocity associated with constituent ą , and c is the speed of light in
vacuum. In deriving these equations, we have assumed that the effect of magnetic
polarization is negligible and all the constituents are electrically neutral. Also,
these equations are formulated within the framework of classical, non-relativistic
theory and they are Galilean invariant up to the order of c-1 . That is, the
equations retain the same form under the transformation
x = x + Vt, E = E - c-1V�B, B = B + c-1V�E,
Pą = Pą, ą = s, w, g (5)
where V is constant. The structure of Pą will be given by constitutive as-
sumptions. We should remark that the above equations hold in all three kinds
of mixtures and the summation in the equations shall take proper terms in each
mixture.
We now formulate the equations of motion for the three mixtures, respectively,
within the context of mixture theory and on the basis of the scheme mentioned
above while ignoring the multi-poles other than dipoles.
2.1. Case 1: Mixture of a porous solid infused with water
When the temperature of the aqueous dielectric body is lower than the critical
c
temperature T , we have a mixture of a porous solid infused with water. The
balance of mass equations for constituents s and w are
"xs
�sdetFs = �sR, Fs := , (6)
"Xs
"�w
+ "�(�wvw) =0, (7)
"t
where �ą is the effective mass density in the present configuration, Xs is a
reference configuration for the solid and �sR is the effective mass density of the
solid constituent in the reference configuration.
The balance of linear momentum equations, for ą = s , w , are
"
ą ą
(�ąvą) +"�(vą�ąvą) ="�� + fem + fą� + �ąg, (8)
"t
� =ą
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 929
ą
fem := c-1 c ("E)�Pą +("B)�(Pą�vą)
"
+ (Pą�B) +"�(vąPą�B) , (9)
"t
fsw + fws = 0. (10)
ą ą
Here � is the Cauchy stress, fem is the electromagnetic force acting on the
constituent ą , fą� is the internal momentum supply to the constituent ą from
ą
the constituent � , and g is the gravitational acceleration. The stress � is
asymmetric and meets
ą ą
� - (� )T = (Eą"Pą - Pą"Eą) , Eą := E + c-1vą�B. (11)
Let Uą be the internal energy density associated with the constituent ą , then
we have for ą = s , w ,
" 1 1
�ą |vą|2 + Uą + "� vą�ą |vą|2 + Uą
"t 2 2
ą
= "�(� vą - qą) +hą + fą� �vą + hą� + �ąg�vą, (12)
em
� =ą
"Pą "B
hą := �E - c-1 (Pą�vą)� + "�(vąPą�E) , (13)
em
"t "t
fsw �vs + hsw + fws�vw + hws =0, (14)
where qą is the diffusive flux, hą is the energy source due to the electromag-
em
netic field, and hą� is the internal energy supply to the constituent ą from the
constituent � . Assuming that the solid constituent and the water have the same

temperature field T , we can deduce from the above equations, with the help of
the balance of linear momentum equations, that
dą ą dą ą
- �ą Uą + E � P -"�[qą +PąEą�(vą -vs)]
dt dt
w w
ą
+ tr [(E� + Pą�Eą1) Lą] + fsw +("Pw) Ew - "P E
+"�(Pw "Ew) �(vw - vs) =0, (15)

Different temperature fields can be associated with the solid and the water, and the con-
vective heat transfer between the two constituents can be accounted for through a constitutive
relation. We do not get into this issue in the formulation here, both for the sake of simplicity
and due to the many uncertainties associated with the drying of the aque- ous dielectric.
930 K. R. Rajagopal and L. Tao ZAMP
where
dą " ą ą
ą
:= + vą�", Lą := vi ,j, E := FsT Eą, P := (Fs)-1 Pą,
ij
dt "t
ą ą w
� := � + Pą"Eą, [("Pw) Ew]i := P Ewj. (16)
E j,i
ą ą
Here, � is symmetric and ensures that � satisfies (11).
E
A balance equation for the entropy for the mixture as a whole is adopted,
together with the non-negativity of the internal entropy production,
"
(�ąsą) +"�(�ąsąvą) -"�Ją = �ą =� e"0 (17)
"t
where są , Ją and �ą are, respectively, the entropy density, the entropy flux
and the entropy production associated with the constituent ą . This relation will
ą
yield certain restrictions on the constitutive forms of Pą , � , fą� , Uą , qą ,
są and Ją . To this end, we will assume that the solid constituent behaves like
an elastic body and
ą
{Uą, są} = {Uą, są} P , Fs, �w, T , ą = s, w. (18)
We do not include the porosity � explicitly in the above relation on the basis that
the true mass density of water �wR is supposedly a function of T only. Note that
� = �w/�wR .
Following Liu s scheme [8] which treats the above balance equations as con-
straints that can be appended to (17), we have
"są "Uą dąT
0 d" � = �ą - Z - "�{Ją+ Z [qą+ PąEą�(vą - vs)]}
"T "T dt
ą "są "Uą dą ą
+ ZE + �ą ą - Z � P
ą
dt
"P "P
1 "są "Uą T
s
+Ztr � + Ps�Es1 + Fs �ą - Z Ls
E
Z "Fs "Fs
1 "są "Uą
w
+Ztr � + Pw �Ew1 - �w �ą - Z 1 Lw
E
Z "�w "�w
w w
+Z fsw +("Pw) Ew - "P E +"�(Pw "Ew)
1 1 "sw "Uw T
+ EwPw �"Z + �w("Fs): - Z
Z Z "Fs "Fs
1 "ss "Us
- �s - Z "�w �(vw - vs) + qą�"Z. (19)
Z "�w "�w
Here Z is the Lagrange multiplier related to the balance of energy equation (15).
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 931
Next, pick
"są "Uą
- Z =0
"T "T
or
1 "Uą "są -1
= .
Z "T "T
Further, if the structures of Uą and są are such that
"Uą "są "Aą
= T , or są = - , (20)
"T "T "T
with Aą := Uą - Tsą being the Helmholtz potential of constituent ą , then,
1
Z = , (21)
T
and relation (19) can be recast as
qą + �ąNąuą + Eą�(vą - vs) Pą
0 d" � = - "� Ją +
T
1 ą "Aą dą ą
+ E - �ą ą � P
T dt
"P
1 �s�w
s
+ tr � + Ps�Es1 + (Ns - Nw)
E
T �
"Aą T
-Fs �ą Ls
"Fs
1 �w�s
w
+ tr � + Pw �Ew1 + (Nw - Ns)
E
T �
"Aą
+�w �ą 1 Lw
"�w
1 w w
+ fsw +("Pw) Ew - "P E +"�(Pw "Ew)
T
"Aw T "As
-�w ("Fs): + �s "�w
"Fs "�w
�w�s
+"� (Nw - Ns) �(vw - vs)
�
�w�s 1
+ qą+Ew �(vw -vs)Pw + (Nw -Ns)(vw -vs) �" . (22)
� T
Here
� := �ą, �v := �ąvą, uą := vą - v,
932 K. R. Rajagopal and L. Tao ZAMP
and Ną is a second order tensor (with Ną - N� being symmetric) to be deter-
mined that has the same unit as that of Aą . We notice that there are other forms
ą
to represent the entropy equation: For instance, we may adopt P = �sR/�s
(Fs)-1 Pą as in electromagnetic elastic solids, provided that corresponding mod-
ifications are made in the terms involving Lą , vą - v� and so on in the above
equation. Certainly, these various forms are equivalent, but they may yield differ-
ent consequences in practice, as we may choose to meet these equivalent represen-
tations by different constitutive structures that are sufficient to meet the equivalent
forms. It is rather difficult to judge at this moment which form is better suited
to the problem under consideration, and we will proceed by adopting the form of
(22) and examine its consequences.
In what follows, we construct a model on the basis of the above relation. For
the entropy flux, we choose
qą + �ąNąuą + Eą�(vą - vs) Pą
Ją = - , ą = s, w, (23)
T
to eliminate the divergence term from (22).
In the case of dielectric polarizations, we take (de Groot and Mazur [4])
dą ą ą "Aą
P = ą E - �ą ą , ą = s, w, (24)
dt
"P
ą
where  is positive definite and symmetric. It can be checked that this equa-
tion is Galilean invariant up to the order of c-1 . Furthermore, we partition Aą
according to
1 ą ą
Aą = Aą(Cs, �w, T) + Hą(Cs, �w, T)P �P , ą = w, s, (25)
0
2
where Hą is positive definite and symmetric and Cs := FsT Fs due to the re-
quirement of Galilean invariance on Aą (Tao, Humphrey and Rajagopal [12]),
and we get
dą ą ą ą
P = ą E - �ąHąP , ą = s, w. (26)
dt
ą
We are able to determine Hą and  under the condition of small solid defor-
mations by resorting to the theory of dielectric polarizations in which P = � E
is employed under E = Emaxexp(i�t) (Hasted [6]; Hippel [7]). Here, � is the
complex electric susceptibility of a dielectric, ( � is denoted as [re( � )]-i[im( � )]
with re( � ) and im( � ) being real), � is the frequency of the electromagnetic
wave field, and the material is at rest. Suppose that the susceptibility of the aque-
M M
ous dielectric is obtained experimentally as � = � (�s, �w, T, �) (some
data concerning susceptibility can be found in the book of Metaxas and Mered-
s s s
ith [9]), we may assume that Ps = � E with � estimated by � =
M w w M M
� (�s, 0, T, �) and Pw = � E with � = � - � (�s, 0, T, �) .
Substituting these expressions for Ps and Pw into (26), along with Fs H" 1 ,
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 933
vą = 0 and E = Emaxexp(i�t) results in, for ą = s, w ,
ą = � [im (�ą)] + � [re (�ą)] [im (�ą)]-1 [re (�ą)] , (27)
�
Hą = [ą]-1 [re (�ą)] [im (�ą)]-1 . (28)
�ą
The Cauchy stresses, the internal momentum supply and the diffusive flux are
taken as
�s�w "Aą T
s
� = -Ps�Es1 - (Ns - Nw) +Fs �ą ,
E
� "Fs
�w�s "Aą
w
� = -Pw �Ew1 - (Nw - Ns) - �w �ą 1, (29)
E
� "�w
w w "Aw T
fsw = - ("Pw) Ew+ "P E -"�(Pw "Ew)+�w("Fs):
"Fs
"As �w�s
-�s "�w -"� (Nw - Ns) + bsw(vw - vs)+Bsw"T, (30)
"�w �
and
�w�s
qą = -Ew �(vw - vs) Pw - (Nw - Ns)(vw - vs)
�
-�sw"T + Ksw (vw - vs) . (31)
Here bsw is positive definite and symmetric, accounting for the effect of drag and
sw
� is the thermal conductivity of the mixture which is also positive definite and
symmetric. Both Bsw and Ksw are due to the cross effects associated with the
relative velocity and the temperature gradient.
With these constitutive relations, the entropy equation (22) reduces to
1 ą "Aą ą "Aą
0 d" � = E - �ą ą �ą E - �ą ą
T
"P "P
1
+ [bsw (vw - vs) +Bsw"T ]�(vw - vs)
T
1
+[-�sw"T + Ksw (vw - vs)]�" . (32)
T
sw
This inequality will impose restrictions on bsw , Bsw , � and Ksw . Certainly,
more comprehensive models can be suggested that include more cross effects.
2.2. Case 2: Mixture of a porous solid with gas
When the temperature of the aqueous dielectric is higher than the critical temper-
c
ature T , there is a mixture of a porous solid with gas (dry steam). The balance
equations for constituents s and g and the mixture can be obtained as in the
934 K. R. Rajagopal and L. Tao ZAMP
above case of the mixture of a porous solid with water and they are recorded be-
low for the sake of completeness. One exception is that the porosity � will be
included in formulating Uą and są since the true mass density of gas �gR is not
determined completely by T . Then,
"xs
�sdetFs = �sR, Fs := , (33)
"Xs
"�g
+ "�(�gvg) =0, (34)
"t
"
ą ą
(�ąvą) +"�(vą�ąvą) ="�� + fem + fą� + �ąg, (35)
"t
� =ą
ą
fem := c-1 c ("E)�Pą +("B)�(Pą�vą)
"
+ (Pą�B) +"�(vąPą�B) , (36)
"t
fsg + fgs = 0. (37)
ą ą
� - (� )T = (Eą"Pą - Pą"Eą) , (38)
" 1 1
�ą |vą|2 + Uą + "� vą�ą |vą|2 + Uą
"t 2 2
ą
= "�(� vą - qą) +hą + fą� �vą + hą� + �ąg�vą, (39)
em
� =ą
"Pą "B
hą := �E - c-1 (Pą�vą)� + "�(vąPą�E) , (40)
em
"t "t
fsg �vs + hsg + fgs�vg + hgs =0, (41)
dą ą dą ą
- �ą Uą+ E � P -"�[qą +PąEą�(vą -vs)]
dt dt
g g
ą
+ tr [(E� + Pą�Eą1) Lą] + fsg +("Pg) Eg - "P E
+"�(Pg " Eg) �(vg - vs) =0, (42)
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 935
qą + �ąNąuą + Eą�(vą - vs) Pą
0 d" � = - "� Ją +
T
1 ą "Aą dą ą 1 "Aą ds
+ E - �ą ą � P - �ą �
T dt T "� dt
"P
1 �s�g
s
+ tr � + Ps�Es1 + (Ns - Ng)
E
T �
"Aą T
-Fs �ą Ls
"Fs
1 �g�s
g
+ tr � + Pg �Eg1 + (Ng - Ns)
E
T �
"Aą
+�g �ą 1 Lg
"�g
1 g g
+ fsg +("Pg) Eg - "P E +"�(Pg "Eg)
T
"Ag T "As "Ag
-�g ("Fs): + �s "�g - �g "�
"Fs "�g "�
�g�s
+"� (Ng - Ns) �(vg - vs)
�
�g�s 1
+ qą+ Eg �(vg - vs) Pg + (Ng - Ns)(vg - vs) �" , (43)
� T
ą
{Uą, są, Aą} = {Uą, są, Aą} P , Fs, �g, �, T . (44)
Based on relation (43), we choose the following model for the mixture of the
porous solid and gas. We take
qą + �ąNąuą + Eą�(vą - vs) Pą
Ją = - , ą = s, g, (45)
T
ds s s "As g
P = s E - �s s , P = 0. (46)
dt
"P
As indicated above, we partition As in the form of
1 s s
As = As(Cs, �g, �, T) + Hs(Cs, �g, �, T)P �P , (47)
0
2
and we get
ds s s s
P = s E - �sHsP . (48)
dt
936 K. R. Rajagopal and L. Tao ZAMP
s
We can estimate Hs and  in the same way as we did in Case 1,
s = � [im (�s)] + � [re (�s)] [im (�s)]-1 [re (�s)] ,
�
Hs = [s]-1 [re (�s)] [im (�s)]-1 , (49)
�s
s M
with � = � (�s, 0, T, �) .
The evolution of porosity, the Cauchy stresses, the internal momentum supply
and the heat flux are taken to be
ds� "Aą
= -ł �ą , (50)
dt "�
�s�g "Aą T
s
� = -Ps�Es1 - (Ns - Ng) +Fs �ą , (51)
E
� "Fs
�g�s "Aą
g
� = -Pg �Eg1 - (Ng - Ns) - �g �ą 1, (52)
E
� "�g
g g "Ag T
fsg = - ("Pg) Eg + "P E -"�(Pg "Eg) +�g("Fs):
"Fs
"As "Ag �g�s
-�s "�g + �g "� -"� (Ng -Ns)
"�g "� �
+bsg (vg - vs) +Bsg"T, (53)
and
�g�s
qą = -Eg �(vg - vs) Pg - (Ng - Ns)(vg - vs)
�
-�sg"T + Ksg (vg - vs) , (54)
sg
where ł e" 0, bsg and � are positive definite and symmetric. The entropy
equation (43) reduces to
1 s "As s "As ł "Aą 2
0 d" � = E - �s s �s E - �s s + �ą
T T "�
"P "P
1
+ [bsg (vg - vs) +Bsg"T ]�(vg - vs)
T
1
+[-�sg"T + Ksg (vg - vs)]�" . (55)
T
sg
This inequality will impose restrictions on bsg , Bsg , � and Ksg .
2.3. Case 3: Mixture of a porous solid with water and gas
Suppose that the temperature of the aqueous dielectric in some region reaches
c
the critical temperature T due to the microwave heating, then water starts to
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 937
evaporate and gas (dry steam) comes into being. (It may also happen that the
gas condenses into water when the temperature drops.) There is now a mixture
of three constituents, the porous solid, the water and the gas, and the above
formulation needs to be extended to account for the mass conversion between the
water and the gas constituents due to the phase change. We will assume that the
c
constituents have the same temperature field T which depends on the motion of
the mixture.
The balance equations for mass, linear momentum and angular momentum for
constituents s , w and g and the mixture as a whole are
"xs
�sdetFs = �sR, Fs := , (56)
"Xs
"�ą
+ "�(�ąvą) =Mą, (57)
"t
Mą =0, Ms =0, (58)
where Mw and Mg characterize the mass conversion between the water and the
gas constituents, and
"
ą ą
(�ąvą) +"�(vą�ąvą) ="�� + fem + fą� + Mąvą + �ąg, (59)
"t
� =ą
ą
fem := c-1 c ("E)�Pą +("B)�(Pą�vą)
"
+ (Pą�B) +"�(vąPą�B) , (60)
"t
fsw + fws = 0, fsg + fgs = 0,
fwg + Mwvw + fgw + Mgvg = 0, (61)
ą ą
� - (� )T = (Eą"Pą - Pą"Eą) . (62)
We have adopted three separate constraints in (61) to account for the interactions
between the three constituents and this might be relaxed by adopting only one
constraint for the mixture as a whole.
c
Since the temperature field of the mixture T during the process of phase
change is supposedly known, the balance of energy for the mixture will be used to
determine the mass conversion rate Mg (or Mw ). Let Ucą denote the internal
938 K. R. Rajagopal and L. Tao ZAMP
energy density of constituent ą , we have, for ą = s , w , g ,
" 1 1
�ą |vą|2 + Ucą + "� vą�ą |vą|2 + Ucą
"t 2 2
1
ą
= "�(� vą - qą) +hą + Mą Ucą + |vą|2
em
2
+ fą� �vą + hą� + �ąg�vą, (63)
� =ą
"Pą "B
hą := �E - c-1 (Pą�vą)� + "�(vąPą�E) , (64)
em
"t "t
fsw �vs + hsw + fws�vw + hws =0, fsg �vs + hsg + fgs�vg + hgs =0,
1
fwg �vw + hwg + Mw Ucw + |vw|2
2
1
+fgw �vg + hgw + Mg Ucg + |vg|2 =0. (65)
2
Assuming that
ą
c
Ucą = Ucą P , Fs, �w, �g, �, T , ą = s, w, g, (66)
summarizing equation (63) and eliminating the kinetic energy terms with the help
of (59), we obtain the expression for the mass conversion rate
1 "Ucą "Ucą -1
Mg = |vw - vg|2 + Ucg - Ucw + �ą -
2 "�g "�w
"Ucą dąTc
� - �ą c -"�[qą +PąEą�(vą -vs)]
"T dt
ą "Ucą dą ą "Ucą dą�
+ E - �ą ą � P - �ą
dt "� dt
"P
"Ucą T
s
+tr � + Ps�Es1 - Fs �ą Ls
E
"Fs
"Ucą
w
+tr � + Pw �Ew1 + �w �ą 1 Lw
E
"�w
"Ucą
g
+tr � + Pg �Eg1 + �g �ą 1 Lg
E
"�g
w w
+ fsw +("Pw) Ew - "P E +"�(Pw "Ew)
"Ucw T "Ucs
-�w ("Fs): + �s "�w �(vw - vs)
"Fs "�w
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 939
g g
+ fsg +("Pg) Eg - "P E +"�(Pg "Eg)
"Ucg T "Ucs
-�g ("Fs): + �s "�g �(vg - vs)
"Fs "�g
"Ucw "Ucg
+ fgw - �w "�g + �g "�w �(vw - vg) . (67)
"�g "�w
We adopt the following entropy equation for the mixture
"
0 d" � = (�ąscą) +"�(�ąscąvą) -"�Ją
"t
dą
= �ą scą -"�Ją +Mg (scg - scw) , (68)
dt
where scą , Ją are, respectively, the entropy density and the entropy flux associ-
ated with constituent ą . Substituting (67) into the above equation and assuming
that
ą
c
scą = scą P , Fs, �w, �g, �, T , ą = s, w, g,
(69)
results in
"scą 1 "Ucą dąTc qą+PąEą�(vą -vs)
0 d" � = �ą c - -"� Ją+
c
"T � "T dt �
1 ą "scą "Ucą dą ą
+ E + �ą � - � P
ą ą
� dt
"P "P
1 "scą "Ucą dą�
+ �ą � -
� "� "� dt
1 "scą "Ucą T
s
+ tr � + Ps�Es1 + Fs �ą � - Ls
E
� "Fs "Fs
1 "scą "Ucą
w
+ tr � + Pw �Ew1 - �w �ą � - 1 Lw
E
� "�w "�w
1 "scą "Ucą
g
+ tr � + Pg �Eg1 - �g �ą � - 1 Lg
E
� "�g "�g
1 w w
+ fsw +("Pw) Ew - "P E +"�(Pw "Ew)
�
"scs "Ucs "scw "Ucw T
-�s � - "�w +�w("Fs): � -
"�w "�w "Fs "Fs
940 K. R. Rajagopal and L. Tao ZAMP
1
+� EwPw �" �(vw - vs)
�
1 g g
+ fsg +("Pg) Eg - "P E +"�(Pg "Eg)
�
"scs "Ucs "scg "Ucg T
-�s � - "�g + �g("Fs): � -
"�g "�g "Fs "Fs
1
+� EgPg �" �(vg - vs)
�
1 "scw "Ucw
+ fgw + �w � - "�g
� "�g "�g
"scg "Ucg 1
-�g � - "�w �(vw - vg) + qą�" , (70)
"�w "�w �
where
"scą "scą -1 1
� := scg - scw + �ą - |vw - vg|2
"�g "�w 2
"Ucą "Ucą
+Ucg - Ucw + �ą - . (71)
"�g "�w
To construct a model that satisfies the above relation, we take
"scą 1 "Ucą
- =0,
c
"T � "Tc
and assume that the structures of scą and Ucą are such that
"Ucą c "scą "Acą
= T , or scą = - , (72)
"Tc "Tc "Tc
c
with Acą := Ucą - T scą being the Helmholtz potential associated with con-
c
stituent ą . Then we have � = T ,
qą+�ąNąuą+ Eą�(vą - vs) Pą
0 d" � = - "� Ją+
Tc
1 ą "Acą dą ą 1 "Acą ds�
+ E - �ą ą � P - �ą
c
Tc dt T "� dt
"P
1 �s
s
+ tr � + Ps�Es1 + [�w (Ns - Nw) +�g (Ns -Ng)]
E
c
T �
"Acą T
-Fs �ą Ls
"Fs
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 941
1 �w
w
+ tr � + Pw �Ew1 + [�s (Nw - Ns) +�g (Nw -Ng)]
E
c
T �
"Acą
+�w �ą 1 Lw
"�w
1 �g
g
+ tr � + Pg �Eg1 + [�s (Ng - Ns) +�w (Ng -Nw)]
Tc E �
"Acą
+�g �ą 1 Lg
"�g
1 w w
+ fsw +("Pw) Ew - "P E +"�(Pw "Ew)
Tc
"Acs "Acw T "Acw
+�s "�w - �w ("Fs): - �w "�
"�w "Fs "�
�w�s
+"� (Nw - Ns) �(vw - vs)
�
1 g g
+ fsg +("Pg) Eg - "P E +"�(Pg "Eg)
Tc
"Acs "Acg T "Acg
+�s "�g - �g ("Fs): - �g "�
"�g "Fs "�
�g�s
+"� (Ng - Ns) �(vg - vs)
�
1 "Acw "Acg
+ fgw - �w "�g + �g "�w
Tc "�g "�w
�w�g
+"� (Nw - Ng) �(vw - vg)
�
+ qą + Ew �(vw - vs) Pw + Eg �(vg - vs) Pg
�w�s
+ (Nw - Ns)(vw - vs)
�
�g�s
+ (Ng - Ns)(vg - vs)
�
�w�g 1
+ (Nw - Ng)(vw - vg) �" , (73)
� Tc
942 K. R. Rajagopal and L. Tao ZAMP
and
"scą "scą -1 1
c
T = scg - scw + �ą - |vw - vg|2
"�g "�w 2
"Ucą "Ucą
+Ucg - Ucw + �ą - . (74)
"�g "�w
On the basis of (73), we will simply choose, for ą = s, w, g ,
qą + �ąNąuą + Eą�(vą - vs) Pą
Ją = - , (75)
c
T
g
P = 0, (76)
ds s s "Acs
P = s E - �s s , (77)
dt
"P
dw w w "Acw
P = w E - �w w , (78)
dt
"P
ds� "Acą
= -ł �ą , (79)
dt "�
�s
s
� = -Ps�Es1 - [�w (Ns - Nw) +�g (Ns -Ng)]
E
�
"Acą T
+Fs �ą , (80)
"Fs
�w
w
� = -Pw �Ew1 - [�s (Nw - Ns) +�g (Nw -Ng)]
E
�
"Acą
-�w �ą 1, (81)
"�w
�g
g
� = -Pg �Eg1 - [�s (Ng - Ns) +�w (Ng -Nw)]
E
�
"Acą
-�g �ą 1, (82)
"�g
w w "Acs
fsw = - ("Pw) Ew + "P E -"�(Pw"Ew) -�s "�w
"�w
"Acw T "Acw �w�s
+�w("Fs): + �w "� -"� (Nw -Ns)
"Fs "� �
c
+bsw (vw - vs) +Bsw"T , (83)
g g "Acs
fsg = - ("Pg) Eg + "P E -"�(Pg "Eg) -�s "�g
"�g
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 943
"Acg T "Acg �g�s
+�g("Fs): + �g "� -"� (Ng -Ns)
"Fs "� �
c
+bsg (vg - vs) +Bsg"T , (84)
"Acw "Acg �w�g
fgw = �w "�g - �g "�w -"� (Nw -Ng)
"�g "�w �
c
+bgw (vw - vg) +Bgw"T , (85)
qą = -Ew �(vw - vs) Pw - Eg �(vg - vs) Pg
�w�s �g�s
- (Nw - Ns)(vw - vs) - (Ng - Ns)(vg - vs)
� �
�w�g
c
- (Nw - Ng)(vw - vg) - �swg"T
�
+Ksw (vw - vs) +Ksg (vg - vs) +Kgw (vw - vg) . (86)
Then, equation (73) reduces to
1 ą "Acą ą "Acą ł "Acą 2
0 d" � = ą E -�ą ą � E -�ą ą + �ą
c c
T T "�
"P "P
ą=s,w
1
c
+ [bsw (vw - vs) +Bsw"T ]�(vw - vs)
Tc
1
c
+ [bsg (vg - vs) +Bsg"T ]�(vg - vs)
c
T
1
c
+ [bgw (vw - vg) +Bgw"T ]�(vw - vg)
c
T
c
+[-�swg"T + Ksw (vw - vs) +Ksg (vg - vs)
1
+Kgw (vw - vg)]�" , (87)
Tc
which will impose restrictions on the coefficients.
As in the case of the mixture of s and w , we will adopt
1 ą ą
c c
Acą = Acą(Cs, �w, �g, �, T ) + Hcą(Cs, �w, �g, �, T )P �P , ą = s, w, (88)
0
2
then, equations (77) and (78) become
dą ą ą ą
P = ą E - �ąHcąP , ą = s, w. (89)
dt
Ignoring the polarization of the gas constituent, we may also estimate Hcą and
ą M M c s M c
 as above, by taking � = � (�s, �w, T , �) , � = � (�s, 0, T , �)
w M M c
and � = � - � (�s, 0, T , �) ,
ą =�[im (�ą)] + � [re (�ą)] [im (�ą)]-1 [re (�ą)] ,
�
Hcą = [ą]-1 [re (�ą)] [im (�ą)]-1 . (90)
�ą
944 K. R. Rajagopal and L. Tao ZAMP
3. Interfacial jump conditions
To have a determinate model of the microwave drying process of an aqueous di-
electric body, we need to furnish the interfacial jump conditions. These interfaces
include the boundary of the body (relative to its environment) and the bound-
ary between any two mixture regions of {s, w} , {s, w, g} and {s, g} inside the
body. Certainly, we notice the fuzzy nature related to the latter boundaries, and
we will idealize them as singular surfaces without any surface property attached
to them for the sake of simplicity. The jump conditions are obtained by following
the scheme illustrated in [2] and [14]. We may attach the velocity of the solid
constituent vs to the fields of E and B in deriving these conditions, though it
is not essential.
Let S(t) be an interface with normal n and move with the velocity of � and
let us define
�ł
�łĆ�ł := Ć+ - Ć-.
�ł
�ł �ł
(91)
Here Ć+ and Ć- are the values of Ć from positive and negative sides of n of
S(t) . Then, the jump conditions associated with the Maxwell field equations (1)
through (4) are
�ł �ł
�ł �ł
�łE + Pą�ł�n =0,
�ł �ł
(92)
�ł �ł
�ł �ł
�ł �ł
�łc-1 E + Pą ׽ + B + c-1 vą�Pą�ł�n = 0,
�ł �ł
�ł �ł (93)
�ł �ł
�ł
�łB�ł�n =0,
�ł
�ł �ł
(94)
�ł �ł
�łE - c-1B׽�ł�n = 0.
�ł �ł
�ł �ł (95)
The balance of mass for the mixtures yields
�ł �ł
�ł �ł
�ł
�ł �ł
�ą (vą - �)�ł�n =0. (96)
�ł �ł
To derive the jump conditions associated with the balances of linear momenta and
energy for the mixtures, we rewrite the balance equations, with the help of (1) to
(4), in the forms
"
ą
�ąvą + c-1E�B = "� (-vą�ąvą + � + Pą"Eą)
"t
1
- |E|2 + |B|2 1 + E"E + B"B + �ąg, (97)
2
" 1 1
�ą |vą|2 + Uą + |E|2 + |B|2
"t 2 2
1
ą
= "� - vą�ą |vą|2 + Uą + � vą - qą + vąPą�E
2
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 945
- cE� B - c-1 Pą�vą + �ąvą�g. (98)
Correspondingly, we will get the following jump conditions
�ł
�ł
�ł
�ł
�ł
ą
�ł
�ąvą"(vą - �) - (� )T - Eą"Pą - c-1 (E�B)"�
�ł
�ł
�ł
�ł
�ł
�ł
�ł
�ł
1
�ł
+ |E|2 + |B|2 1 - E"E - B"B�łn = 0, (99)
�ł
�ł
2 �ł
�ł
�ł
�ł
�ł
�ł
1
�ł
ą
�ł
�ł �ą |vą|2 + Uą (vą - �) - � vą + qą - Pąvą�E
�ł
�ł 2
�ł
�ł
�ł
�ł
�ł
1
�ł
- |E|2 + |B|2 � + cE�B�ł�n =0. (100)
�ł
�ł
2 �ł
�ł
Several more sets of interfacial jump conditions can be generated. One results
from the continuity of displacement ws (:= xs - Xs) and the balance of mass
for the solid constituent
�ł �ł
�łws�ł
�ł �ł
= 0, (101)
�ł �ł
�ł�s (vs - �)�ł�n =0.
�ł �ł
�ł �ł (102)
However, we may not have the freedom to obtain similar jump conditions asso-
ciated with the balances of linear momenta and energy for the solid constituent
due to the problem associated with identifying the momentum and energy of the
electromagnetic field related to the solid constituent and the ambiguity of whether
to decompose some quantities into surface flux terms or source terms. The relation
(102) is met automatically when the interface is a part of the boundary "D of the
aqueous dielectric body since � = vs|"D on one side and �s = 0 on the other
side. Also, the requirement can be imposed that the temperature is continuous
across the interfaces between the region of {s, w} and the region of {s, w, g} and
c
between {s, w, g} and {s, g} and this temperature equals T . The other set of
conditions follows from the entropy equation for the mixture as a whole [12],
�ł �ł
�ł �ł
�ł
�ł �ł
�ąsą (vą - �) - Ją�ł�n =0, (103)
�ł �ł
which is consistent with the assumption that the interface does not possess intrinsic
surface properties of its own, like surface mass density, surface entropy and so on.
4. Summary
We have presented, within the context of Maxwell-Lorentz field equations and mix-
ture theory, the electromagnetic field equations and the equations of motion for the
946 K. R. Rajagopal and L. Tao ZAMP
porous solid, water and gas in an aqueous dielectric corresponding to the process
of microwave drying. Three mixtures are considered: mixture of a porous solid
and water; mixture of a porous solid, water and gas; and mixture of a porous solid
and gas. On the basis of continuum thermodynamics, constitutive relations are
proposed for the macroscopic electric polarization vectors, Cauchy stresses, heat
fluxes, internal momentum supplies and so on to characterize the thermomechan-
ical behaviors of these mixtures. The interfacial jump conditions between these
different mixture regions were derived. This framework is tentative at best, some
revisions might be needed upon further studies, and much more needs to be done
with regard to fixing the specific forms of the constitutive relations based on ex-
perimental data and by solving specific drying problems. The appendix discusses
briefly an issue related to the requirement of material frame indifference.
Appendix
Though a subject of controversy (Edelen and McLennan [1]; Wang [15]), the re-
quirement of material frame indifference (MFI) plays a major role in constitutive
theory (Noll [10]; Oldroyd [11]; Truesdell and Toupin [14]). Part of the disagree-
ment lies in whether MFI should be adopted or a less restrictive requirement, the
principle of Galilean invariance, should be used to formulate constitutive relations.
Confusion has also been caused in averaged turbulence modeling and mixture the-
ory, wherein the situation is worse in that there seems no sharp criterion to classify
some quantities as material properties, due to the averaged fields involving large
spatial and temporal scales. In this appendix, we examine a feature of MFI briefly
from a physical perspective, a perspective related to observers measuring of length
and time interval necessary in obtaining constitutive relations, not merely talking
about these material properties in a purely abstract and conceptual sense. And
the intended purpose is to provide some justification for our adoption of Galilean
invariance in this work.
Let {x, t} be an inertial frame, and let a body in motion be represented in
this frame. Consider another frame x ,
x = Q(t)x + a(t), QQT = QT Q = 1
and assume that an observer is located at x0 in the frame of x . We have to
restrict x0 , Q and a such that
Ł
Ł
QQT (x0 - a) + a <c
(104)
where c is the speed of light in vacuum and  <<1 . This is essential in order
that we may neglect relativistic effects. Otherwise, the observer s measuring scales
of length and time would differ nontrivially from that of the observer s located in
the frame of x , and we have to go beyond the framework of Newtonian mechanics
which is needed to establish the corresponding relations between the quantities
Vol. 53 (2002) Modeling of the microwave drying process of aqueous dielectrics 947
Ł
Ł
associated with these two frames. Hence, there is the constraint that Q , a , a
and x0 have to meet.
We now take a = C (t + t0) with C and t0 being constant in order to
simplify the analysis and we have
Ł
QQT (x0 - C (t + t0)) + C < c.
(105)
The independence between Q , x0 and C which is necessary for applying MFI
to formulate constitutive relations results in
|C| <c (106)
Ł
by picking Q = 0 ;
Ł
QQT x0 <c
(107)
by picking C = 0 ; and
Ł
QQT C (t + t0) < 3c
(108)
from (105) to (107) and the triangle inequality. Furthermore, as a proper estima-
tion, we can infer from (106) and (108) that
"
Ł
QQT |t + t0| d"3 2.
(109)
"
Ł
That is, the product of the rotating speed QQT / 2 and time t + t0 is bounded
Ł
above. Specifically, if we choose t0 = 12 (hours) (or both t0 = 0 and QQT
constant) and t " [0, 12] (hours), the rotating speed will be smaller than the spin
of Earth. This implies essentially that Q should be treated as constant or the
frame x should be approximated as not rotating, because we neglect the effect of
Earth s spin in measuring material properties in most cases.
Acknowledgement
We thank the National Science Foundation and the National Institutes for Health
for support of this work.
References
[1] D. G. B. Edelen and J. A. McLennan, Material indifference: a principle or a convenience,
Int. J. Engng. Sci. 11 (1973), 813-817.
[2] A. C. Eringen, Mechanics of Continua. Wiley, New York, 1980.
[3] A. C. Eringen and G. A. Maugin, Electrodynamics of Continua I. Springer-Verlag, New
York, 1989.
[4] S. R. de Groot and P. Mazur, Non-Equilibrium Thermodynamics. Dover Publications, Inc.,
New York, 1984.
948 K. R. Rajagopal and L. Tao ZAMP
[5] S. R. de Groot and L. G. Suttorp, Foundations of Electrodynamics. North-Holland Publish-
ing Co., Amsterdam, 1972.
[6] J. B. Hasted, Aqueous Dielectric. Chapman and Hall, London, 1973.
[7] A. R. von Hippel, Dielectrics and Waves. The MIT Press, Cambridge, Massachusetts, 1954.
[8] I-Shin Liu, Method of Lagrange multipliers for exploitation of the entropy principle, Arch.
Rat. Mech. Anal., 46 (1972), 131-148.
[9] A. C. Metaxas and R. J. Meredith, Industrial Microwave Heating. Peter Peregrinus Ltd,
London, 1983.
[10] W. Noll, On the continuity of the solid and fluid state, J. Ratl. Mech. Anal. 4 (1955), 13-81.
[11] J. G. Oldroyd, On the formulation of rheological equations of state, Proc. Roy. Soc. A200
(1950), 523-541.
[12] L. Tao, L., J. D. Humphrey and K. R. Rajagopal, A mixture theory for heat-induced alter-
ations in hydration and mechanical properties in soft tissues, Int. J. Engng. Sci. 39 (2001),
1535-1556.
[13] C. Truesdell, Rational Thermodynamics. Springer-Verlag, New York, 1984.
[14] C. Truesdell and R. Toupin, The Classical Field Theories. In Handbuch der Physik Vol.
III/1, W. Flugge, ed, Springer-Verlag, New York, 1960.
[15] C. C. Wang, On the concept of frame-indifference in continuum mechanics and in the kinetic
theory of gases, Arch. Ratl. Mech. Anal. 58 (1975), 381-393.
[16] S. Whitaker, Coupled transport in multiphase systems: A theory of drying. In Advances in
Heat Transfer, Vol. 31, 1-104. Academic Press, New York, 1998.
K. R. Rajagopal and L. Tao
Department of Mechanical Engineering
Texas A&M University
College Station, Texas 77843-3123
USA
e-mail: [e-mail]
To access this journal online:
http://www.birkhauser.ch