1 Introduction

Thermodynamic variational principles (later “principles”) are now accepted tools to develop kinetic (evolution) equations in systems with constraints among the state variables. It should be mentioned that already Lord Rayleigh [25] set the first steps in 1871 by introduction of a “dissipation function,” for details see Sect. 2.1 of this paper. For sake of completeness, we would like to add that “dissipation potentials” have been used in mechanics of materials since several decades, for an overview see [4]. However, groundbreaking contributions were published by Onsager [23] and Ziegler [29] and [30]. With respect to relevant contributions to the principles, we refer the reader to Hackl and Fischer [15] published in 2008 and to the overview [12] published in 2014.

From the point of view of applications of the principles, we consider Ziegler’s concept [29, 30] as one of the most successful ones for processes not too far from the equilibrium and obeying the laws of linear non-equilibrium thermodynamics, see, e.g., the textbooks by de Groot and Mazur [6], Callen [3] and Gyarmati [14]. Ziegler’s concept maximizes the dissipation with the constraint that the dissipation D as bilinear form of forces and fluxes is equal to a particular dissipation function Q, a physical quantity, e.g., a quadratic form of the fluxes (for explanation of the terms see the next section). Generally, this extremal principle, supplemented by proper conservation laws as well as boundary and contact conditions, leads to a set of partial differential equations for the evolution of the system, which may be rather difficult to handle even for quite simple systems. However, material science requires models describing the evolution of the material system only by a limited number of the most appropriate time dependent state variables. A straightforward way, how one can obtain the time dependent state variables, has been provided by reformulating Ziegler’s concept in terms of discrete state variables as shown by Svoboda and Turek [27].

From this time on, this concept has been denominated as the “Thermodynamic Extremal Principle” (TEP). Since this formulation avoids having to work directly with partial differential equations, the TEP has become an efficient tool to find kinetic equations for the internal variables in a wide variety of fields. Diverse applications to thermodynamic problems from multicomponent diffusion to grain growth have been reported, for references see, e.g., [12]. However, one aspect has not yet been dealt with in detail in the literature, namely the situation, when the sets of kinetic variables in the dissipation D and in the dissipation function Q are not the same ones. Generally, the sets of kinetic variables may be related by functional dependencies, which can be considered as constraints. The goal of this note is now to show how to implement efficiently such functional relations in the maximization procedure.

2 System description

2.1 Variables and Gibbs energy

  • Remark 1: We use the matrix notation with “T” as superscript meaning a transposed matrix.

  • Remark 2: All quantities are described in the actual configuration with “div” and “grad” and the Laplace operator \(\varDelta \) as operators there.

First of all, we introduce a representative volume element (RVE) of a multicomponent system. The RVE is supposed to be a sufficiently small part of the system so that it can be considered as homogeneous but large enough to allow a statistical approach (see also Sect. 4.2).

To define the thermodynamic material system the following quantities are introduced:

  1. 1.

    the external state variables, collected for sake of easier writing in a vector \( \mathbf {X} \),

  2. 2.

    the internal state variables, collected in a vector \( \mathbf {x} \),

  3. 3.

    the Gibbs thermodynamic potential G as function \(G\left( \mathbf {X} , \mathbf {x} ; \mathbf {p} \right) \) defined on the RVE with given material parameters collected in a vector p.

We would like to mention the definition of internal variables and their appearance in the dissipation D worked out in detail in the seminal paper by Coleman and Gurtin [5] 50 years ago. The term “dissipative structures” was explained by Maugin [20]. Recently this author published also an according very interesting overview (a “saga”) [21].

2.2 Dissipation and dissipation function

According to standard continuum thermodynamics, see e.g. [5], the dissipation D per RVE follows, with \( \mathbf {X} \) assumed as fixed, from the well established relation

$$\begin{aligned} D=-\frac{\partial G}{\partial \mathbf {x} } \cdot \dot{ \mathbf {x} }= \mathbf {f} \cdot \dot{ \mathbf {x} }, \end{aligned}$$
(1)

where the vector \(-\partial G / \partial \mathbf {x} \) is denoted as collection of the thermodynamic forces \( \mathbf {f} \) and \(\dot{ \mathbf {x} }\) as collection of the thermodynamic fluxes. We use the term “collection,” since a thermodynamic force can be a scalar or a vectorial quantity. With respect to the terms thermodynamic “forces” and “fluxes,” we refer to Onsager [23] and the book by Baluffi et al. [2]. The dissipation D may be considered as a pure mathematical construct, of course, stemming from the physical entity G.

Therefore, to reflect the internal mechanisms of dissipation, the dissipation function Q is introduced which, in a rather general form, is a function \(Q\left( \mathbf {X} , \mathbf {x} ,\dot{ \mathbf {x} };\tilde{ \mathbf {p} }\right) \) with given material parameters, collected in vector \(\tilde{ \mathbf {p} }\) (note that \( \mathbf {p} \) and \(\tilde{ \mathbf {p} }\) need not to be identical). In many practical cases, however, Q depends on a set \(\left\{ \mathbf {K} \right\} \) of kinetic variables, whose members do not appear as components of the vector \(\dot{ \mathbf {x} }\). The reader may think of diffusion, where the concentrations (or molar fractions) of the components are internal state variables \( \mathbf {x} \) and, therefore, their rates appear in D. However, Q includes the set of fluxes \( \mathbf {j} _{k} \) of the components (counting index k) and not the vector \(\dot{ \mathbf {x} }\). Therefore, let us separate the internal variables in primary ones, \( \mathbf {x} _\mathrm {p} \) (with \(M_\mathrm {p} \) members) whose rates \(\dot{ \mathbf {x} }_\mathrm {p} \) appear also in Q and secondary ones, \( \mathbf {x} _\mathrm {s} \) (with \(M_\mathrm {s} \) members), whose rates \(\dot{ \mathbf {x} }_{s} \) do appear only in D but not in Q. The set of kinetic variables \(\left\{ \mathbf {K} \right\} \) consists of \(M_\mathrm {w} \) scalar members collected in a vector \( \mathbf {w} \) and a subset \(\left\{ \mathbf {W} \right\} \) with \(M_\mathrm {W} \) vector members \( \mathbf {W} _{k} \) (k is the index; think, e.g., on the a/m fluxes). Then Q follows as \(Q\left( \mathbf {X} , \mathbf {x} ,\dot{ \mathbf {x} }_\mathrm {p} ,\left\{ \mathbf {K} \right\} ;\tilde{ \mathbf {p} }\right) \). Since both, D and Q, describe the dissipation process, they must be equal, i.e.,

$$\begin{aligned} D\left( \mathbf {X} , \mathbf {x} ;\dot{ \mathbf {x} }_\mathrm {p} ,\dot{ \mathbf {x} }_\mathrm {s} ; \mathbf {p} \right) = \mathbf {f} \cdot \dot{ \mathbf {x} }=Q\left( \mathbf {X} , \mathbf {x} ,\dot{ \mathbf {x} }_\mathrm {p} ;\left\{ \mathbf {K} \right\} ;\tilde{ \mathbf {p} }\right) . \end{aligned}$$
(2)

The TEP in its original form, see, e.g., [12], involving the same kinetic variables \(\dot{ \mathbf {x} }\) in both, D and Q, states that the dissipation obtains a maximum subjected to the constraint, Eq. (2). The according Lagrangian L follows as

$$\begin{aligned} L=D-\lambda \left( Q-D\right) \end{aligned}$$
(3)

with \(\lambda \) being a Lagrange multiplier. Maximizing of L with respect to the kinetic variables \(\dot{ \mathbf {x} }\) provides a relation between the thermodynamic forces \( \mathbf {f} \) and the kinetic variables \(\dot{ \mathbf {x} }\) as

$$\begin{aligned} \mathbf {f} =\left( \frac{Q}{\left( \partial Q / \partial \dot{ \mathbf {x} } \right) \cdot \dot{ \mathbf {x} } } \right) \frac{\partial Q}{\partial \dot{ \mathbf {x} }} . \end{aligned}$$
(4)

It should be noted again that Q as dissipation function goes back to Lord Rayleigh [25]. Furthermore, it may be of interest that Onsager [23] introduced a dissipation potential \(\phi \left( \mathbf {X} , \mathbf {x} ,\dot{ \mathbf {x} }; \mathbf {q} \right) \) with given material parameters, collected in a vector \( \mathbf {q} \), for the thermodynamic forces as

$$\begin{aligned} \frac{1}{T} \, \mathbf {f} =\frac{\partial \phi }{\partial \dot{ \mathbf {x} }} \end{aligned}$$
(5)

with T as the temperature.

3 The thermodynamic extremal principle (TEP) and constraints

3.1 Dissipation and constraints

The dissipation D, Eq. (2), can be reformulated as

$$\begin{aligned} D= \mathbf {f} _\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {f} _\mathrm {s} \cdot \dot{ \mathbf {x} }_\mathrm {s} , \end{aligned}$$
(6)

with

$$\begin{aligned} \mathbf {f} _\mathrm {p} =-\frac{\partial G}{\partial \mathbf {x} _\mathrm {p} } ,\quad \mathbf {f} _\mathrm {s} =-\frac{\partial G}{\partial \mathbf {x} _\mathrm {s} } . \end{aligned}$$
(7)

The dissipation \(D_\mathrm {system} \) of the total system follows as

$$\begin{aligned} D_\mathrm {system} =\int _{V_\mathrm {system} } D \, \mathrm {d}V=\int _{V_\mathrm {system} } \left( \mathbf {f} _\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {f} _\mathrm {s} \cdot \dot{ \mathbf {x} }_\mathrm {s} \right) \, \mathrm {d}V. \end{aligned}$$
(8)

In many cases, a functional relation exists between \(\dot{ \mathbf {x} }_\mathrm {s} \) and \(\dot{ \mathbf {x} }_\mathrm {p} \) and the set \(\left\{ \mathbf {K} \right\} \) in the form

$$\begin{aligned} \dot{ \mathbf {x} }_\mathrm {s} = \mathbf {A} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {B} \cdot \mathbf {w} + \mathbf {C} \cdot \mathrm {div}\left( \left\{ \mathbf {W} \right\} \right) \end{aligned}$$
(9)

with the following definitions:

  1. 1.

    \( \mathbf {A} \) is a given matrix \(A_{ij} \) with \(M_\mathrm {s} \) rows and \(M_\mathrm {p} \) colums;

  2. 2.

    \( \mathbf {B} \) is a given matrix \(B_{ij} \) with \(M_\mathrm {s} \) rows and \(M_\mathrm {w} \) colums;

  3. 3.

    \( \mathbf {C} \) is a given matrix \(C_{ij} \) with \(M_\mathrm {s} \) rows and \(M_\mathrm {W} \) colums;

  4. 4.

    \(\left\{ \mathbf {W} \right\} \) represents a set with \(M_\mathrm {W} \) vectors \( \mathbf {W} _k\);

  5. 5.

    \(\mathrm {div}\left( \{ \mathbf {W} \}\right) \) represents a vector with \(M_\mathrm {W} \) elements \(\mathrm {div}\left( \mathbf {W} _{k} \right) \).

Finally, \( \mathbf {C} \cdot \mathrm {div}\left( \left\{ \mathbf {W} \right\} \right) \) represents a vector \( \mathbf {c} \) with the elements \(\displaystyle c_{i} =\sum _{k=1}^{M_\mathrm {W} }C_{ik} \, \mathrm {div}\left( \mathbf {W} _{k} \right) , \quad i=1,\ldots ,M_\mathrm {s}.\)

Taking into account the above definitions and inserting of Eq. (9) into the integrand on the r.s. of Eq. (8) yield

$$\begin{aligned} D= \mathbf {f} _\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {f} _\mathrm {s} \cdot \dot{ \mathbf {x} }_\mathrm {s} =\left( \mathbf {f} _\mathrm {p} + \mathbf {A} ^\mathrm {T} \cdot \mathbf {f} _\mathrm {s} \right) \cdot \dot{ \mathbf {x} }_\mathrm {p} +\left( \mathbf {B} ^\mathrm {T} \cdot \mathbf {f} _\mathrm {s} \right) \cdot \mathbf {w} +\sum _{i=1}^{M_\mathrm {s} }f_{\mathrm {s},i} \sum _{k=1}^{M_\mathrm {W} }C_{ik} \; \mathrm {div}\left( \mathbf {W} _{k} \right) . \end{aligned}$$
(10)

If we introduce the vector \(\tilde{ \mathbf {f} }_\mathrm {s} = \mathbf {C} ^{T} \cdot \mathbf {f} _\mathrm {s} \) and integrate the last term on the r.s. of Eq. (10)\({}_{2}\) over the total volume \(V_\mathrm {system} \) and apply the Gauss theorem (assuming that \( \mathbf {W} _{k} \cdot \mathbf {n} =0\) with \( \mathbf {n} \) being the normal to \(\partial V_\mathrm {system} \)), then we can reformulate (10)\({}_{2}\) with \(\tilde{ \mathbf {f} }_\mathrm {p} = \mathbf {f} _\mathrm {p} + \mathbf {A} ^{T} \cdot \mathbf {f} _\mathrm {s} \) as

$$\begin{aligned} D= \mathbf {f} _\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {f} _\mathrm {s} \cdot \dot{ \mathbf {x} }_\mathrm {s} =\tilde{ \mathbf {f} }_\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} +\left( \mathbf {B} ^{T} \cdot \mathbf {f} _\mathrm {s} \right) \cdot \mathbf {w} -\sum _{k=1}^{M_\mathrm {W} }\mathrm {grad}\left( \tilde{f}_{\mathrm {s},k} \right) \cdot \mathbf {W} _{k} . \end{aligned}$$
(11)

We define now \( \mathbf {f} _{\mathrm {w}} = \mathbf {B} ^{T} \cdot \mathbf {f} _\mathrm {s} \) and \( \mathbf {f} _{\mathrm {W},k}=-\mathrm {grad}\left( \tilde{f}_{\mathrm {s},k} \right) \) and can formulate D as

$$\begin{aligned} D=\tilde{ \mathbf {f} }_\mathrm {p} \cdot \dot{ \mathbf {x} }_\mathrm {p} + \mathbf {f} _{\mathrm {w}} \cdot \mathbf {w} +\sum _{k=1}^{M_\mathrm {W} } \mathbf {f} _{\mathrm {W},k} \cdot \mathbf {W} _{k} . \end{aligned}$$
(12)

For sake of easier writing, we denominate the rates as \(\dot{ \varvec{\xi } } = \{\dot{ \mathbf {x} }_\mathrm {p} , \mathbf {w} , \mathbf {W} _{k}\}\) (\(k=1,\ldots ,M_\mathrm {W}\)), so that we can formulate D, Eq. (2), as \(D= \tilde{D}\left( \mathbf {X} , \mathbf {x} ,\; \dot{ \varvec{\xi } };\; \mathbf {p} \right) \) and \(Q= \)\(\tilde{Q}\left( \mathbf {X} , \mathbf {x} ,\; \dot{ \varvec{\xi } };\; \tilde{ \mathbf {p} }\right) \). Furthermore, the thermodynamic forces are denoted as

$$\begin{aligned} \mathbf {f} _\xi = \frac{\partial \tilde{D}}{\partial \dot{ \varvec{\xi } }} = \{\tilde{ \mathbf {f} }_\mathrm {p}, \mathbf {f} _{\mathrm {w}}, \mathbf {f} _{\mathrm {W},1},\ldots , \mathbf {f} _{\mathrm {W},M_\mathrm {W}}\}. \end{aligned}$$
(13)

Finally, we can rewrite Eq. (2) as

$$\begin{aligned} \tilde{D}\left( \mathbf {X} , \mathbf {x} ,\dot{ \varvec{\xi } }; \mathbf {p} \right) =\tilde{Q}\left( \mathbf { X} , \mathbf {x} ,\dot{ \varvec{\xi } }; \tilde{ \mathbf {p} }\right) . \end{aligned}$$
(14)

3.2 The TEP in modified form

The modification of D to \(\tilde{D}\) as above allows applying the TEP in its original form given by Eqs. (2) and (3), which will be replaced now by Eq. (14) and the modified Lagrangian

$$\begin{aligned} \tilde{L}=\tilde{D}-\lambda \left( \tilde{Q}-\tilde{D}\right) . \end{aligned}$$
(15)

Maximizing of \(\tilde{L}\) with respect to the kinetic variables \(\dot{ \varvec{\xi } }\) provides a relation between the thermodynamic forces \( \mathbf {f} _\xi \) and the rates of the kinetic variables \(\dot{ \varvec{\xi } }\) as

$$\begin{aligned} \mathbf {f} _\xi =\left( \frac{\tilde{Q}}{\left( \partial \tilde{Q} / \partial \dot{ \varvec{\xi } } \right) \cdot \dot{ \varvec{\xi } } } \right) \frac{\partial \tilde{Q}}{\partial \dot{ \varvec{\xi } }} . \end{aligned}$$
(16)

Since one needs finally the evolution equations for \(\dot{ \varvec{\xi } }\), it is necessary to invert Eq. (16). If Q is supposed as a homogeneous function in \(\dot{ \varvec{\xi } }\) of order l, yielding \({\partial \tilde{Q} / \partial \dot{ \varvec{\xi } }} \cdot \dot{ \varvec{\xi } }=l \,\tilde{Q}\), it follows

$$\begin{aligned} \mathbf {f} _\xi =\frac{1}{l} \frac{\partial \tilde{Q}}{\partial \dot{ \varvec{\xi } }} . \end{aligned}$$
(17)

If \(\tilde{Q}\) is a quadratic function in \(\dot{ \varvec{\xi } }\), which means \(l=2\), having the form \(\tilde{Q}=\dot{ \varvec{\xi } } \cdot \mathbf {U} \left( \mathbf {X} , \mathbf {x} \right) \cdot \dot{ \varvec{\xi } }\) with \( \mathbf {U} \) being a positive-definite symmetric matrix involving material and structural parameters, then this yields

$$\begin{aligned} \dot{ \varvec{\xi } }= \mathbf {U} ^{-1} \cdot \mathbf {f} _\xi . \end{aligned}$$
(18)

No problem from the mathematical point of view does occur, if the order of homogeneity l is unequal to 2 but equal for all processes. If this is not the case, we refer to the papers by Hackl et al. [16, 17] and the recent paper by Goddard [13] referring to the necessity of introducing a dissipation potential for some distinct problems.

4 Application of the TEP to local and non-local diffusion

One of the most relevant processes in the evolution of material systems is the diffusion process. Let us consider a multicomponent system consisting of M components with the site fractions \(y_{k} ,k=1,\ldots ,M\). To each component a diffusive flux \( \mathbf {j} _{k} \) is assigned. The diffusion is vacancy-mediated (as it is the case, when K substitutional and \(\left( M-K\right) \) interstitial components are diffusing). The vacancies occupy a site fraction \(y_{0} \) and diffuse with a flux \( \mathbf {j} _{0} \). Let us assume that \(K\le M\) substitutional components diffuse. As outlined in detail in Svoboda et al. [26] Sects. 2.1 and 2.2 two basic identities exist: one for the site fractions and the other for the fluxes. Both identities lead to two constraints, namely

$$\begin{aligned} y_{0} =1-\sum _{k=1}^{K}y_{k} ,\quad \mathbf {j} _{0} =-\sum _{k=1}^{K} \mathbf {j} _{k} . \end{aligned}$$
(19)

Note, that constraint Eq. (19)\(_1\) can directly be inserted into the Gibbs energy G, since the \(y_{k}\) are internal variables. Furthermore, constraint Eq. (19)\(_2\) can also directly be inserted into the dissipation function Q as “classical” function of the squares of the fluxes.

The rate, with which the vacancies are generated or annihilated, is denoted as \(\alpha \), \(0\le \left| \alpha \right| \le \infty \), with \(\alpha =0\) for no generation or annihilation and with \(\left| \alpha \right| \Rightarrow \infty \) for ideal generation or annihilation. The mass balance, also treated in Svoboda et al. [26], is expressed for the site fractions as

$$\begin{aligned} \dot{y}_{k} =-y_{k} \alpha -\bar{\varOmega }\mathrm {div}\, \mathbf {j} _{k} ,\quad k=1,\ldots ,M. \end{aligned}$$
(20)

The quantity \(\bar{\varOmega }\) is the molar volume of one mole of lattice sites and depends, therefore, on the individual site fractions \(y_{k}, k=1,\ldots ,M\). The mass balance activates now a third set of constraints, which will be dealt with below.

Going back to Fick’s work from 1855 one needs a constitutive law for the fluxes \( \mathbf {j} _{k} \), \(k=1,\ldots ,M\). Therefore, the TEP is employed to derive such a constitutive (or evolution) law, see again Svoboda et al. [26], Sect. 4 there, and the rather recent paper by Fischer et al. [11].

4.1 Local diffusion

This is the standard case for a deterministic situation, where no stochastic processes are considered. We concentrate only on constraints Eq. (20). The dissipation D, see Eq. (1), has as kinetic variables the \(\dot{y}_{k} \) and \(\alpha \). The dissipation function Q, see Eq. (2), has as kinetic variables the fluxes \( \mathbf {j} _{k} \) and \(\alpha \). Note that we avoid the rate of a diffusion-induced creep strain for sake of an easier formulation. The dissipation D (denoted now as \(\tilde{D}\)) is calculated according to Eqs. (9) to (12) by assigning the actual quantities to the symbols in general Eqs. (11) and (12). The flux \(\dot{ \mathbf {x} }_\mathrm {p}\) consists of only one member \(\dot{x}_\mathrm {p} =\alpha \), \( \mathbf {x} _\mathrm {s} \) corresponds to the vector \( \mathbf {y} \) with \( \mathbf {y} ^{T} =\left( y_{1} ,y_{2} ,\ldots ,y_{M} \right) \) and \(\dot{ \mathbf {x} }_\mathrm {s} \) to \(\dot{ \mathbf {y} }\). With the notation of Eqs. (9, 10) we have A degenerating to a (1x1 matrix) with the only element \(-1\), \( \mathbf {w} = \mathbf {0} , \mathbf {C} =-\bar{\varOmega }\; \mathbf {I} \) (\( \mathbf {I} \) is the unity matrix) and \( \mathbf {W} _{k} = \mathbf {j} _{k} \) with \(k=1,\ldots ,M\) and \(M_{W} =M_\mathrm {s} =M\). Then we find with Eq. (12)

$$\begin{aligned} \tilde{D}=\left( f_{\alpha } - \mathbf {f} _\mathrm {s} \cdot \mathbf {y} \right) \; \alpha +\sum _{k=1}^{M}\left( \mathrm {grad}\left( \bar{\varOmega }f_{\mathrm {s},k} \right) \right) \cdot \mathbf {j} _{k} . \end{aligned}$$
(21)

4.2 Non-local diffusion

Non-local concepts in continuum mechanics are already dealt with since the pioneering work of Mindlin on non-local plate theory, more than half a century ago, see, e.g., Mindlin’s Collected Works in [18]. A general view on non-local continuum fields was published by Eringen [9]. Here, it shall be mentioned that discrepancies have been detected between the Eringen model in differential and integral form, see, e.g., recent paper [10]. The rather popular gradient approach was related to the non-local field theory by Aifantis [1]. Locally varying material properties as the elastic moduli have motivated to introduce non-local theoretical concepts, see, e.g., recent contributions [7, 22]. With respect to an extremal principle recent paper [24] concerning the Principle of Minimum of Potential Energy of an elastic system shall be mentioned. Particularly of increasing relevance are computational algorithms for non-local diffusion, together with constraints as, e.g., volume constraints, taking into account stochastic processes (as jumping of atoms or Brownian motion). Here we refer to papers by Du et al. [8] and Tian et al. [28] and Lehoucq and Rowe [19] concerning the analysis and approximation of non-local diffusion problems.

We offer a problem solution concept for non-local diffusion by an efficient application of the TEP with constraints. For sake of demonstrating the concept in an easy as possible way we consider the case of diffusion of only one interstitial element with the site fraction y, to which a chemical potential \(\mu =\mu \left( y\right) \) is assigned, and the activity of vacancies represented by their generation/annihilation rate \(\alpha \) [since \(\alpha \) is involved in the site fraction balance of y, see Eq. (19)\(_1\)]. Furthermore, the following (simplifying) assumptions are applied:

  1. 1.

    the molar volume \(\bar{\varOmega }\) is assumed as constant quantity in space and time;

  2. 2.

    the role of the stress state is considered as negligible; consequently, diffusion-induced creep is not activated;

  3. 3.

    the volume of RVE, denoted as \(\mathrm {d}V\), is as a constant quantity; therefore, the material derivative \(\mathrm {d}(\mathrm {d}V)/\mathrm {d}t =0\) or, in other words, \(\mathrm {div}\; \mathbf {v} = \mathbf {0} \) with \( \mathbf {v} \) being the material velocity vector;

  4. 4.

    the quantity y and the chemical potential \(\mu \) (see below) are replaced by

    $$\begin{aligned} \dot{y}\Rightarrow \left( 1-\theta \right) \dot{y}+\theta \dot{\bar{y}}=\dot{\tilde{y}},\quad \mu \Rightarrow \left( 1-\theta \right) \mu +\theta \bar{\mu }=\tilde{\mu }. \end{aligned}$$
    (22)

The quantities \(\bar{y}\) and \(\bar{\mu }\) are the non-local averages of the local quantities y and \(\mu \) resp., e.g., over the RVE. Furthermore, a time-independent switcher \(\theta \) is introduced \(0\le \theta \le 1\), which allows investigating the interaction of a pure deterministic behavior (\(\dot{y}, \mu \), local) and pure stochastic behavior (\(\dot{\bar{y}}, \bar{\mu }\), non-local). A short introduction to the definition of a non-local quantity is presented in the “Appendix”.

The dissipation D follows with Eq. (22) as

$$\begin{aligned} D=-\mu \dot{\tilde{y}}=-\mu \left( \left( 1-\theta \right) \dot{y}+\theta \dot{\bar{y}}\right) . \end{aligned}$$
(23)

The dissipation function Q writes as

$$\begin{aligned} Q={ \mathbf {j} ^{2} / A} +K\alpha ^{2} \end{aligned}$$
(24)

with A as a given mobility coefficient and K as a given bulk viscosity; for details see [26], Sects. 4 and 5 there.

As constraint the balance for the site fraction y is expressed in \(\tilde{y}\) and writes with Eq. (20)

$$\begin{aligned} \dot{\tilde{y}}+\tilde{y}\alpha +\bar{\varOmega } \, \mathrm {div}\; \mathbf {j} =0, \end{aligned}$$
(25)

and with Eq. (22) as

$$\begin{aligned} \left( 1-\theta \right) \dot{y}+\theta \dot{\bar{y}}+\left[ \left( 1-\theta \right) y+\theta \bar{y}\right] \alpha +\bar{\varOmega } \,\mathrm {div}\; \mathbf {j} =0. \end{aligned}$$
(26)

The TEP is formulated for the system with Eqs. (2326) as

$$\begin{aligned} L=\int _{V_\mathrm {system} } \left[ D-\lambda \left( D-Q\right) +2\,\eta \left( \dot{\tilde{y}}+\tilde{y}\alpha +\bar{\varOmega } \, \mathrm {div}\; \mathbf {j} \right) \right] \; \mathrm {d}V. \end{aligned}$$
(27)

Now two Lagrange multipliers, namely the constants \(\lambda \) and \(\eta \), appear. Before calculating the Lagrange parameters, it is useful to apply the “duality condition” [Eq. (47) in the “Appendix”] to the expressions \(\mu \dot{\bar{y}}\) in D and \(\eta \dot{\bar{y}}\) in the constraint in Eq. (27), which can be replaced by \(\bar{\mu }\dot{y}\) and \(\bar{\eta }\dot{y}\), resp. Furthermore, the Gauss theorem allows replacing \(\eta \; \mathrm {div}\; \mathbf {j} \) by \(-\mathrm {grad} \, \eta \cdot \mathbf {j} \) in Eq. (27), since we assume a closed system with a zero normal projection of \( \mathbf {j} \) at the surface. These both steps allow to rewrite L, Eq. (27), as

$$\begin{aligned} L=-\left( 1-\lambda \right) \int _{V_\mathrm {system} } \tilde{\mu }\dot{y} \, \mathrm {d}V+\lambda \int _{V_\mathrm {system} } \left( \mathbf {j} ^{2} / A +K\alpha ^{2} \right) \, \mathrm {d}V + 2\int _{V_\mathrm {system} } \left( \tilde{\eta }\dot{y}+\eta \tilde{y}\alpha -\bar{\varOmega } \, \mathrm {grad}\,\eta \cdot \mathbf {j} \right) \, \mathrm {d}V . \end{aligned}$$
(28)

We perform variations with respect to \(\dot{y}, \mathbf {j} \) and \(\alpha \) yielding

$$\begin{aligned} \text {ad} \quad \dot{y}{:}&\qquad -\left( 1-\lambda \right) \tilde{\mu }+2\tilde{\eta }=0, \end{aligned}$$
(29)
$$\begin{aligned} \text {ad} \;\;\; \alpha {:}&\qquad \lambda K\alpha +\eta \tilde{y}=\lambda K\alpha +\tilde{\eta }y=0, \end{aligned}$$
(30)
$$\begin{aligned} \text {ad} \quad \mathbf {j} {:}&\qquad \lambda \, \mathbf {j} / A -\bar{\varOmega } \, \mathrm {grad}\,\eta =0. \end{aligned}$$
(31)

After multiplication of Eq. (29) by \(\dot{y}\), of Eq. (30) by \(\alpha \) and (31) by \( \mathbf {j} \), integration over the system and keeping in mind \(D=Q\), see Eqs. (23, 24), one finds \(\lambda =-1\) and from

$$\begin{aligned} \text {ad} \; (29){:}&\quad \tilde{\eta }=\tilde{\mu }, \end{aligned}$$
(32)
$$\begin{aligned} \text {ad} \; (30){:}&\quad \alpha =\eta \tilde{y}/ K, \end{aligned}$$
(33)
$$\begin{aligned} \text {ad} \; (31){:}&\quad \mathbf {j} =-A\bar{\varOmega }\,\mathrm {grad}\,\eta . \end{aligned}$$
(34)

Obviously, \( \mathbf {j} \) depends on the Lagrange parameter \(\eta \), see Eq. (34), which is related to the chemical potential \(\mu \) with the non-local relation Eq. (32), \(\tilde{\eta }=\tilde{\mu }\).

If \(\theta =0\), which refers to a pure local problem, then immediately the established equation

$$\begin{aligned} \mathbf {j} =-A\bar{\varOmega }\,\mathrm {grad}\,\mu \end{aligned}$$
(35)

appears.

However, if \(\theta \ne 0\), one has to extract \(\eta \) from \(\tilde{\eta }\) together with \(\tilde{\eta }=\tilde{\mu }\). To separate \(\eta \) from \(\tilde{\eta }\) we can use the relation

$$\begin{aligned} \eta =\beta _{0} \bar{\eta } + \beta _{2} \varDelta \bar{\eta }+\text {h.o.t.}, \end{aligned}$$
(36)

derived in the “Appendix”, see Eq. (45). The coefficients \(\beta _{0}\) and \(\beta _{2}\) are given also in the “Appendix”, see Eqs. (44) and (46), respectively. Then Eq. (32), with the assumption of a known function \(\tilde{\mu }=\left( 1-\theta \right) \mu +\theta \bar{\mu }\), allows formulating a partial differential equation based on Eq. (36) as

$$\begin{aligned} \beta _{2} \left( 1-\theta \right) \varDelta \bar{\eta }+\left( \beta _{0} +\left( 1-\beta _{0} \right) \theta \right) \bar{\eta }=\tilde{\mu }. \end{aligned}$$
(37)

The solution of Eq. (37), together with the appropriate boundary condition, see Eqs. (49) or (50) in the “Appendix”, gives \(\bar{\eta }\) from which \(\eta \) is obtained by applying Eq. (36). This closes the formulation.

5 Conclusion

The “Thermodynamic Extremal Principle” (TEP), based on Ziegler’s work, can be formulated in the frame of linear thermodynamics as a maximization problem with only one constraint equating the dissipation in the system (bilinear form in fluxes and forces) and the dissipation function (quadratic form in fluxes accounting for material properties like diffusivities and interface mobilities). This concept of the TEP can also be used in a straightforward way, if some further constraints, e.g., as balance or conservation equations meeting specific mechanisms etc., exist. The respective procedure is outlined in this paper. Then the TEP provides in the current form original constitutive laws implicitly respecting the constraints. Moreover, the TEP allows treating systems for which the dissipation and the dissipation function are expressed in different kinetic variables coupled by linear relations representing further constraints. A proper redefinition of the kinetic variables may significantly reduce the number of kinetic variables and constraints and provide the most effective description of the system. One should keep in mind that the redefinition does not influence the dissipation and thus also the system kinetics. The developed concept is demonstrated on an example for local and non-local multicomponent diffusion.