arxiv: v1 [math.na] 8 Feb 2017

Transcription

1 Variational iscretization of te nonequilibrium termoynamics of simple systems arxiv:17.594v1 [mat.na] 8 Feb 17 François Gay-Balmaz Hiroaki Yosimura CNRS - LMD - IPSL Scool of Science an Engineering Ecole Normale Supérieure e Paris Wasea University 4 Rue Lomon 755 Paris, France Okubo, Sinjuku, Tokyo , Japan francois.gay-balmaz@lm.ens.fr yosimura@wasea.jp 7 November 18 Abstract In tis paper, we evelop variational integrators for te nonequilibrium termoynamics of simple close systems. Tese integrators are obtaine by a iscretization of te Lagrangian variational formulation of nonequilibrium termoynamics evelope in Gay-Balmaz an Yosimura [17a], an tus exten te variational integrators of Lagrangian mecanics, to inclue irreversible processes. In te continuous setting, we erive te structure preserving property of te flow of suc systems. Tis property is an extension of te symplectic property of te flow of te Euler-Lagrange equations. In te iscrete setting, we sow tat te iscrete flow solution of our numerical sceme verifies a iscrete version of tis property. We also present te regularity conitions wic ensure te existence of te iscrete flow. We finally illustrate our iscrete variational scemes wit te implementation of an example of a simple an close system. Contents 1 Introuction Nonequilibrium termoynamics of simple systems 3.1 Funamental laws of nonequilibrium termoynamics Variational formulation for nonequilibrium termoynamics Discretization of te variational formulation Variational integrators in Lagrangian mecanics Variational integrators for te termoynamics of simple systems Structure preserving properties Termoynamics of simple systems - continuous case Termoynamics of simple systems - iscrete case Examples Variational iscretization scemes Example: a mass-spring-friction system moving in an ieal gas Numerical tests

2 1 Introuction 1 Introuction Nonequilibrium termoynamics is a penomenological teory wic aims to ientify an escribe te relations among te observe macroscopic properties of a pysical system an to etermine te macroscopic ynamics wit te elp of funamental laws of termoynamics e.g. Stueckelberg an Sceurer [1974]. Te fiel of nonequilibrium termoynamics naturally inclues macroscopic isciplines suc as classical mecanics, flui ynamics, elasticity, an electromagnetism. It is well known tat te equations of motion of classical mecanics, i.e., te Euler- Lagrange equations can be erive from Hamilton s variational principle applie to te action functional associate to te Lagrangian of te mecanical system. One of te many features of te variational formulation is tat it amits a iscrete version wic allows te erivation of structure preserving numerical scemes for te system. Suc scemes, calle variational integrators, see, Wenlant an Marsen [1997], Marsen an West [1], Lew, Marsen, Ortiz an West [4] are obtaine via a iscrete version of Hamilton s principle an are originally base on Moser-Veselov iscretizations see, Veselov [1988], Veselov [1991], Moser an Veselov [1991]. Several extensions of tis meto ave been evelope, for example to treat te case of force mecanical systems Kane, Marsen, Ortiz, an West [] or nonolonomic mecanical systems Cortés an Martínez [1], McLaclan an Perlmutter [6]. In Gay-Balmaz an Yosimura [17a,b], we ave evelope a Lagrangian variational formulation for nonequilibrium termoynamics wic extens Hamilton s principle of classical mecanics by allowing te inclusion of irreversible penomena in bot iscrete an continuum systems, i.e., systems wit finite an infinite egrees of freeom. Te irreversibility is encoe into a nonlinear nonolonomic constraint given by te expression of te entropy prouction associate to all te irreversible processes involve. From a matematical point of view, te variational formulation of Gay-Balmaz an Yosimura [17a,b] may be regare as a nonlinear generalization of te Lagrange- Alembert principle use in nonolonomic mecanics, see e.g., Bloc [3]. In orer to formulate te nonolonomic constraint, to eac irreversible process is associate a variable calle te termoynamic isplacement tat generalizes te termal isplacement introuce in Green an Nagi [1991], following von Helmoltz [1884]. Te introuction of suc variables allows te efinition of a corresponing variational constraint. In te present paper, we evelop variational integrators for nonequilibrium termoynamics by iscretizing te Lagrangian variational formulation evelope in Gay-Balmaz an Yosimura [17a]. Te resulting numerical scemes are tus extensions of te variational integrators of Lagrangian mecanics tat enable to inclue irreversible penomena. In te present paper, we restrict our iscussions to te case of simple close systems, i.e., close systems in wic one termal scalar variable an a finite set of mecanical variables are sufficient to escribe entirely te state of te system, toug we will be able to evelop our iscrete teory to anle more general cases incluing te nonequilibrium termoynamics of continuum systems. A key property of variational integrators in Lagrangian mecanics is teir symplecticity, meaning tat te iscrete flow, similarly to te flow of te continuous system, preserves a symplectic form. Tis ensures an excellent long-time energy beavior, see Hairer, Lubic, an Wanner [6]. Wen irreversible effects are consiere in te ynamics, te symplecticity of te flow may be lost at te continuous level, so tere is no ope to iscretize te system wit a symplectic integrator, in general. In te paper, we sall present a property of te flow F t of simple close systems in termoynamics, wic reuces to te symplecticity

3 Nonequilibrium termoynamics of simple systems 3 of te flow in absence of termal effects. Tis property as te form F t Ω Ω = t F s ωs, for all t, were Ω is an entropy-epenent symplectic form, ω is a one-form encoing te effects of friction an temperature, an is te exterior erivative. We ten sow tat our numerical sceme verifies a iscrete version of tis formula an terefore it reuces to a symplectic integrator in absence of termal effects. Te paper is organize as follows. In Section we review te funamental laws governing te nonequilibrium termoynamics of macroscopic systems by following te axiomatic formulation of Stueckelberg an Sceurer [1974]. Ten, we also review te Lagrangian variational formulation of nonequilibrium termoynamics evelope in Gay-Balmaz an Yosimura [17a], wic is an extension of Hamilton s principle of classical mecanics tat allows te inclusion of irreversible penomena. In Section 3, after recalling some basic facts about variational integrators in Lagrangian mecanics, we propose a iscrete version of te variational formulation for nonequilibrium termoynamics of simple close systems an euce a variational integrator for tese systems. In Section 4, we present a property of te flow F t of simple close systems in termoynamics, wic reuces to te symplecticity of te flow in absence of termal effects. Ten we sow tat te iscrete flow of our variational integrator verifies a iscrete version of tis property. We also stuy te regularity conitions wic ensure te existence of te iscrete flow. Finally, in Section 5, we illustrate te implementation of our integrator wit an example of a simple system. Nonequilibrium termoynamics of simple systems In tis section we first review te funamental laws governing te nonequilibrium termoynamics of macroscopic systems. We follow te axiomatic formulation of termoynamics evelope by Stueckelberg aroun 196 see, for instance, Stueckelberg an Sceurer [1974], wic is well suite for te stuy of nonequilibrium termoynamics as a general macroscopic ynamic teory tat extens classical mecanics to account for irreversible processes. Neeless to say, it is important to point out tat tis axiomatic formulation inclues te escription of systems out of equilibrium an is not restricte to te treatment of equilibrium states an transition from one equilibrium state to anoter. Ten, we review te Lagrangian variational formulation of nonequilibrium termoynamics from Gay-Balmaz an Yosimura [17a] wic is an extension of Hamilton s principle of classical mecanics to allow te inclusion of irreversible penomena. For brevity, in tis paper, we will restrict to te case of simple an close systems..1 Funamental laws of nonequilibrium termoynamics For te macroscopic escription of nonequilibrium termoynamics, we ave te following laws, see Stueckelberg an Sceurer [1974]: I First law: For every system, tere exists an extensive scalar state function E, calle energy, wic satisfies t Et = P ext W t + PH ext t + PM ext t, were t enotes time, PW ext t is te power ue to external forces acting on te mecanical variables of te system, PH ext ext t is te power ue to eat transfer, an PM t is te power ue to matter transfer between te system an te exterior.

4 . Variational formulation for nonequilibrium termoynamics 4 II Secon law: For every system, tere exists an extensive scalar state function S, calle entropy, wic obeys te following two conitions. a Evolution part: If te system is aiabatically close, te entropy S is a non-ecreasing function wit respect to time, i.e., St = It, t were It is te entropy prouction rate of te system accounting for te irreversibility of internal processes. b Equilibrium part: If te system is isolate, as time tens to infinity te entropy tens towars a finite local maximum of te function S over all te termoynamic states ρ compatible wit te system, i.e., lim St = max S[ρ]. t + ρ compatible In tis context, a system is sai to be close if tere is no excange of matter between te system an te exterior, i.e., PM ext t = ; a system is sai to be aiabatically close if it is close an tere is no eat excanges between te system an te exterior, i.e., PM ext ext t = PH t = ; an a system is sai to be isolate if it is aiabatically close an tere is no mecanical power excange between te system an te exterior, i.e., PM extt = ext t = P t =. P ext H W By efinition, te evolution of an isolate system is sai to be reversible if It =, namely, te entropy is constant. In general, te evolution of a system is sai to be reversible, if te evolution of te total isolate system wit wic it interacts is reversible. In tis paper, we only consier simple an close systems. By efinition, a simple system 1 is a system were one scalar termal variable S an a finite set q i, q i of mecanical variables are sufficient to escribe entirely te state of te system, an we assume tat tere is no power ue to matter transfer PM ext between te system an te exterior since te system is close.. Variational formulation for nonequilibrium termoynamics Consier a simple close system escribe by a mecanical variable q Q an one entropy variable S R. Let L = Lq, q, S : T Q R R be te Lagrangian of te system, F ext : T Q R T Q te external force, F fr : T Q R T Q te friction force, an PH ext te power ue to eat transfer between te system an te exterior. Te forces are fiber preserving maps, i.e., F fr q, q, S, F ext q, q, S Tq Q, were Tq Q enotes te cotangent space to Q at q. A common form for te Lagrangian is Lq, q, S := K mec q, q Uq, S, were K mec : T Q R enotes te kinetic energy of te mecanical part of te system assume to be inepenent of S an U : Q R R enotes te potential energy, wic is a function of bot te mecanical variable q an te entropy S. Te variational formulation for te termoynamics of simple close systems is efine as follows; see Def..1 in Gay-Balmaz an Yosimura [17a]. 1 In Stueckelberg an Sceurer [1974] tey are calle élément e système Frenc. We coose to use te Englis terminology simple system instea of system element. See also Gay-Balmaz an Yosimura [17a].

5 . Variational formulation for nonequilibrium termoynamics 5 A curve qt, St Q R, t [, T ] R is a solution of te variational formulation if it satisfies te variational conition δ T Lq, q, St + T F ext q, q, S, δq t =, Variational Conition.1 for amissible variations δqt an δst subject to te constraint L S q, q, SδS = F fr q, q, S, δq, Variational Constraint. an if it satisfies te nonlinear nonolonomic constraint L S q, q, SṠ = F fr q, q, S, q P ext H. Penomenological Constraint.3 Taking variations of te integral in.1, integrating by part an using δq = δqt =, it follows T L q L t q + F ext, δq + L S δs t =, were te variations δq an δs ave to satisfy te variational constraint.. Now, replacing L S δs by te virtual work expression F fr q, q, S, δq accoring to. an using te penomenological constraint, te curve qt, St satisfies te following evolution equations for te termoynamics of te simple close system L t q L q = F ext q, q, S + F fr q, q, S, L S Ṡ = F fr q, q, S, q P ext H..4 Notice tat te explicit expression of te constraint.3 involves penomenological laws for te friction force F fr ; tis is te reason wy we refer to it as a penomenological constraint. Te constraint. is calle a variational constraint since it is a conition on te variations to be use in.1, wic follows from.3 by formally replacing te velocity by te corresponing virtual isplacement, an by removing te contribution from te exterior of te system. Suc a simple corresponence between te penomenological an variational constraints still ols for te general class of termoynamical systems consiere in Gay- Balmaz an Yosimura [17a,b]; Gay-Balmaz [17]. Energy balance law. Te energy associate wit L : T Q R R is te function L E : T Q R R efine by Eq, q, S := q, q Lq, q, S. Using te system.4 an efining PW ext := F ext, q, we obtain te energy balance law : t E = L t q L q, q L S Ṡ = P W ext + PH ext, wic is consistent wit te first law of termoynamics. Notice tat energy is preserve wen te system is isolate, i.e., wen PW ext = P H ext = P M ext =, consistently wit te first law of termoynamics.

6 3 Discretization of te variational formulation 6 Entropy prouction. Te temperature is given by minus te partial erivative of te Lagrangian wit respect to te entropy, T = L S, wic is assume to be positive. So te secon equation in.4 reas T Ṡ = P ext H F fr q, q, S, q. Accoring to te secon law of termoynamics, for aiabatically close systems, i.e., wen P H ext = P M ext =, entropy is increasing. So te friction force F fr must be issipative, tat is F fr q, q, S, q, for all q, q, S T Q R. For te case in wic te force is linear in velocity, i.e., F fr q, q, S = λq, S q,, were λq, S : T q Q T q Q R is a two covariant tensor fiel, tis implies tat te symmetric part λ sym of λ as to be positive. For a simple system, te internal entropy prouction as te form It = 1 T F fr q, q, S, q. Recovering Hamilton s principle. In absence of te entropy variable an te external force, te constraints isappear an ence te variational formulation given in equations.1.3 reuces to Hamilton s principle of Lagrangian mecanics T δ Lq, qt =,.5 for variations δqt vanising at te enpoints, i.e., δq = δqt =. Reversibility. As we recalle earlier, te evolution of an isolate system is sai to be reversible if te entropy is constant. In te case of an isolate simple system, in view of te secon equation in.4 tis means tat te evolution is suc tat F fr qt, qt, St, qt =. 3 Discretization of te variational formulation In tis section we first make a brief review of some basic facts about variational integrators in Lagrangian mecanics. Ten we propose a iscrete version of te variational formulation for nonequilibrium termoynamics of simple close systems an euce a variational integrator for tese systems. We also present a conition wic ensures te existence of te flow of te integrator an we make several comments on te construction of te constraint. 3.1 Variational integrators in Lagrangian mecanics Variational integrators are numerical scemes tat arise from a iscrete version of Hamilton s variational principle.5; see, for instance, Wenlant an Marsen [1997] an Marsen an West [1]. Let Q be te configuration manifol of a mecanical system an let L : T Q R be a Lagrangian. Suppose tat a time step as been fixe, enote by {t k = k k =,..., N} te sequence of times iscretizing [, T ], an by q : {t k } N k= Q, q k := q t k te corresponing iscrete curve. A iscrete Lagrangian L : Q Q R is an approximation of te time integral of te continuous Lagrangian between two consecutive configurations q k an q k+1 L q k, q k+1 tk+1 t k Lqt, qtt, 3.1

7 3. Variational integrators for te termoynamics of simple systems 7 were q k = qt k an q k+1 = qt k+1. Equippe wit suc a iscrete Lagrangian, one can now formulate a iscrete version of Hamilton s principle.5 accoring to L q k, q k+1 =, N 1 δ k= for variations δq k vanising at te enpoints. Tus, if we enote D i te partial erivative wit respect to te i t variable, tree consecutive configuration variables q k 1, q k, q k+1 must verify te iscrete analogue of te Euler-Lagrange equations: D L q k 1, q k + D 1 L q k, q k+1 =. 3. Tese iscrete Euler-Lagrange equations efine, uner appropriate conitions, an integration sceme wic solves for q k+1, knowing te two previous configuration variables q k 1 an q k. A iscrete Lagrangian L is calle regular if te following maps, calle iscrete Legenre transforms, are local iffeomorpisms: F + L : Q Q T Q, F + L q, q 1 = q 1, D L q, q 1 T q 1 Q F L : Q Q T Q, F + L q, q 1 = q, D 1 L q, q 1 T q Q. 3.3 In fact it is enoug to prove tat one of tese maps is a local iffeomorpism. Tis turns out to be equivalent to te invertibility of te matrix D 1 D L q, q 1 for all q, q 1. Uner te regularity ypotesis, te sceme 3. yiels a well-efine iscrete flow : Q Q Q Q; q k 1, q k q k, q k+1 tat is symplectic: F L F L Ω L = Ω L, 3.4 were te symplectic form Ω L := F ± L Ω can is efine wit respect to eiter F + L or FL. External forces can be ae using a iscrete version of te Lagrange- Alembert principle in a similar manner, see Marsen an West [1]. 3. Variational integrators for te termoynamics of simple systems Let us first exten te concept of iscrete Lagrangian 3.1 from mecanics to te nonequilibrium termoynamics of simple close systems escribe by a mecanical variable q Q an one entropy variable S R. Definition 3.1. Consier a simple close system wit Lagrangian L = Lq, q, S : T Q R R, suppose tat a time step as been fixe, an enote by {t k = k k =,..., N} te sequence of times iscretizing [, T ]. A iscrete Lagrangian is a function L : Q Q R R R, wic is an approximation of te time integral of L between two consecutive states q k, S k an q k+1, S k+1 : L q k, q k+1, S k, S k+1 tk+1 were qt i = q i an St i = S i, for i = k, k + 1. t k Lqt, qt, Stt,

8 3. Variational integrators for te termoynamics of simple systems 8 One example of suc a iscrete Lagrangian, wen Q is a vector space, may be given by L q k, q k+1, S k, S k+1 := [ L q k, q k+1 q k, S k + L q k+1, q ] k+1 q k, S k+1. Similarly, we efine te iscrete analogue of external an friction forces as follows. Definition 3.. Consier an external force F ext : T Q R T Q an a friction force F fr : T Q R T Q, wic are fiber preserving maps, i.e., F fr q, q, S, F ext q, q, S T q Q. We efine iscrete friction forces an iscrete exterior forces to be maps F fr, F fr+, F ext, F ext+ : Q Q R R T Q, suc tat te following approximation ols F fr q k, q k+1, S k, S k+1, δq k + F fr+ q k, q k+1, S k, S k+1, δq k+1 tk+1 t k F fr qt, qt, St, δqt, similarly for F ext±, F ext, were qt i = q i, St i = S i, δqt i = δq i, an δst i = δs i, for i = k, k + 1. Tese iscrete forces are require to be fiber preserving in te sense tat π Q F fr± = π ± Q, π Q F ext± = π ± Q, were π Q : T Q Q is te cotangent bunle projection an π Q, π+ Q : Q Q R R Q are efine by π Q q, q 1, S, S 1 = q an π + Q q, q 1, S, S 1 = q 1. Construction of te constraint. For te case of nonolonomic mecanics wit linear constraint, te iscrete constraint can be constructe from a finite ifference map, see Cortés an Martínez [1] an McLaclan an Perlmutter [6]. We sall exten tis construction to our nonlinear situation an wit te entropy variable. Following McLaclan an Perlmutter [6], a finite ifference map ϕ Q on a manifol Q is a iffeomorpism ϕ Q : N Q T Q, were N Q is a neigboroo of te iagonal Q in Q Q an T Q is a neigboroo of te zero section of T Q, wic satisfies te following conitions: 1. ϕ Q Q is te zero section of T Q;. τϕ Q N Q = Q; 3. τϕ Q q, q = q. All tree conitions can be equivalently escribe as: ϕ Q q, q = q. Definition 3.3. Taking two finite ifference maps ϕ Q : N Q T Q an ϕ R : N R T R, we efine te finite ifference map ϕ = ϕ Q ϕ R : N Q R T Q R by ϕq k, q k+1, S k, S k+1 = ϕ Q q k, q k+1, ϕ R S k, S k+1, 3.5 were te neigboroos are N Q R = N Q N R, Q R T Q R = T Q T R. = Q R, an

9 3. Variational integrators for te termoynamics of simple systems 9 Recall tat, in te continuous setting, te penomenological constraint is te subset C K T Q R efine by q, q, S, Ṡ C K L S q, q, SṠ F fr q, q, S, q =, 3.6 were we assume PH ext = for simplicity. Notice tat for any pysically relevant Lagrangian see also Assumption II in 4.3 below, te function P : T Q R R efine by L P q, q, S, Ṡ := S q, q, SṠ F fr q, q, S, q, 3.7 is a submersion, since P = L Ṡ S q, q, S, being minus te temperature. Tus C K is a coimension one submanifol of T Q R. Notice also tat te zero section is inclue in C K. In orer to formulate te iscrete version of te penomenological constraint, we nee to efine a iscrete version CK Q Q R R of te submanifol C K T Q R. Suc a iscrete version is written wit te elp of a function P : Q Q R R R as C K = {q, q 1, S, S 1 Q Q R R P q, q 1, S, S 1 = }. 3.8 In te efinition below, we present a way to construct CK from a given finite ifference map. We will ten sow ow to construct bot CK an L in a consistent way. Definition 3.4. Given te constraint in 3.6 an a finite ifference map ϕ in 3.5, te associate iscrete constraint is efine by C K := ϕ 1 C K T Q R Q Q R R. 3.9 In tis case te function P in 3.8 is obtaine by composing te function P : T Q R R in 3.7 wit te finite ifference map ϕ. It is possible to construct bot te iscrete penomenological constraint an te iscrete Lagrangian in a consistent way. Inee, suppose tat a finite ifference map ϕ : Q Q R R T Q R is given, ten one can construct C K as in 3.9 an L as L := L π ϕ, 3.1 were we recall is te time step an π : T Q R = T Q T R T Q R is te canonical projection. Tis formula can be interprete in two ways. On one an, as L = L ϕ, were L := π L is te lifte Lagrangian on T Q R, wile it can be written as L = L Ψ, were we efine te iscretizing map Ψ by Ψ := π ϕ : Q Q R R T Q R. Remark 3.5. We will sow tat te construction of bot L an CK from a unique finite ifference map ϕ is not neee to obtain te structure preserving properties in 5. One can coose a finite ifference map ϕ an a iscretizing map Ψ wic are not necessarily relate troug Ψ = π ϕ. For example, in nonolonomic mecanics linear case, tere are examples of integrators in wic CK an L are not constructe from te same finite ifference mapping, but wic perform extremely well, see 4.18 in McLaclan an Perlmutter [6]. Definition 3.6. By analogy wit te continuous variational constraint., we efine te iscrete variational constraint by imposing te following constraint on δq k an δs k as D 3 L q k, q k+1, S k, S k+1 δs k + D 4 L q k, q k+1, S k, S k+1 δs k+1 = F fr q k, q k+1, S k, S k+1, δq k + F fr+ q k, q k+1, S k, S k+1, δq k

10 3. Variational integrators for te termoynamics of simple systems 1 Definition 3.7 Discrete variational formulation for te nonequilibrium termoynamics of simple systems. Given a iscrete Lagrangian L, iscrete friction forces F fr±, external forces F ext±, an a iscrete penomenological constraint CK, a iscrete curve q, S = {q k, S k } N k= is a solution of te variational formulation if it satisfies te iscrete variational conition L q k, q k+1, S k, S k+1 N 1 δ k= + N 1 k= F ext q k, q k+1, S k, S k+1, δq k + F ext+ q k, q k+1, S k, S k+1, δq k+1 =, for variations satisfying te iscrete variational constraint 3.11 an were te iscrete curve q, S = {q k, S k } N k= is subject to te iscrete penomenological constraint q k, q k+1, S k, S k+1 C K. A irect application of tis variational formulation yiels te following result. Teorem 3.8. A iscrete curve q, S = {q k, S k } N k= is a solution of te variational formulation if an only if it satisfies te following iscrete evolution equations: D 1 L q k, q k+1, S k, S k+1 + D L q k 1, q k, S k 1, S k +F fr + F ext q k, q k+1, S k, S k+1 +F fr+ + F ext+ 3.1 q k 1, q k, S k 1, S k =, q k, q k+1, S k, S k+1 CK. Discrete flow map. By applying te implicit function teorem, we see tat if te following matrix [ D D 1 L r + D F r D 4D 1 L r + D 4 F r ] 3.13 D P r D 4 P r is invertible for all r := q, q 1, S, S 1, were we set F 3.1 yiels a well-efine iscrete flow := F fr + F ext, ten te sceme F L : q k, q k+1, S k, S k+1 C K q k+1, q k+, S k+1, S k+ C K It is easy to ceck tat a matrix of te form 3.13 is invertible if an only if D 4 P r an te matrix D D 1 L r + D F r 1 D4 D 1 L r + D 4 F D 4 P r r D P r 3.15 is invertible. Tis criteria generalizes to te case of termoynamics, te regularity criteria of te iscrete Lagrangian of iscrete mecanics, namely te conition tat D D 1 L q, q 1 is invertible for all q, q 1, see 3.3, wic may be recovere from 3.15 wen te entropy variable an te forces are absent.

11 4 Structure preserving properties 11 4 Structure preserving properties In tis Section, we present a property of te flow F t of a simple an close system wic reuces to te symplecticity of te flow in absence of termal effects. Ten we sow tat te iscrete flow of our numerical integrator verifies a iscrete version of tis property. 4.1 Termoynamics of simple systems - continuous case It is well known tat wen te Lagrangian L : T Q R of a mecanical system is regular, ten te flow F t : T Q T Q of te Euler-Lagrange equations preserves te symplectic form Ω L = FL Ω can calle te Lagrangian two-form on T Q: F t Ω L = Ω L. 4.1 In orer to formulate te extension of tis property to te case of te termoynamics of simple systems, we first make below some efinitions an assumptions concerning te Lagrangian function in termoynamics. Regularity an assumptions on te Lagrangian. R, te Legenre transform is efine by FL : T Q R T Q, FLq, v, S := Given a Lagrangian L : T Q R q, L v q, v, S. Te only ifference wit te stanar case in mecanics is te epenence on S. By efinition, we say tat te Lagrangian Lq, v, S is regular if an only if for eac S fixe, te map q, v T Q FLq, v, S T Q is a local iffeomorpism. One easily cecks tat tis is equivalent to te invertibility of te matrix L v i v q, v, S for all q, v, S. j We efine te following two Lagrangian forms on T Q R, namely, te Lagrangian oneform Θ L q, v, S := FL Θ can = L q, v, Sq v an te Lagrangian two-form wic reas locally Ω L q, v, S = L v i q j q, v, Sqi q j + Ω L := Θ L Ω T Q R, L v i v j q, v, Sqi v j + L v i S q, v, Sqi S. In absence of te entropy variable, tese forms recover te usual Lagrangian forms on T Q efine in Lagrangian mecanics. We now write two pysical assumptions mae on te Lagrangian L : T Q R R. Assumption I: A first pysical restriction on te Lagrangian is te following assumption L v i S =, 4.

12 4.1 Termoynamics of simple systems - continuous case 1 wic means tat te temperature T = L S oes not epen on v or, equivalently, te momentum p = L v oes not epen on S. In oter wors, L L q, v, S = S S q, S an L v L q, v, S = q, v. v It follows from Assumption I 4. tat te Lagrangian is necessarily of te form Lq, v, S = Kq, v Uq, S, for two functions K : T Q R an U : Q R R. Lagrangian two-form reas Uner Assumption I, te Ω L q, v, S = L v i q j q, v, Sqi q j + L v i v j q, v, Sqi v j. In tis case, Ω L can be seen as a S-epenent two-form on T Q. Moreover, Ω L is symplectic on T Q, for eac S fixe, if an only if te Lagrangian is regular. Assumption II: Any pysical Lagrangian must verify te conition L q, v, S <, for all q, v, S, 4.3 S since L S = T is ientifie wit minus te temperature. If Assumption I 4. is verifie, ten Assumption II 4.3 reas simply U q, S >, S for all q, S. Structure preserving property. Recall tat given a Lagrangian L : T Q R R, an te forces F fr, F ext : T Q R T Q, te evolution equations are given by te system.4, rewritten ere for te curve qt, vt, St as L L qt, vt, St qt, vt, St t q q = F ext qt, vt, St + F fr qt, vt, St, 4.4 L S qt, vt, StṠt = F fr qt, vt, St vt PH ext t, qt = vt, were qt = q t. We assume tat te Lagrangian is regular an tat te pysical assumptions 4. an 4.3 are verifie. In tis case, one observes tat 4.4 gives a well-efine first orer orinary ifferential equation for te curve qt, vt, St an, terefore, a wellefine flow F t. Let us ientify T Q R wit te space of solution of 4.4 by using te corresponence q, v, S T Q R F t q, v, S = qt, vt, St Solutions of 4.4, were F t : T Q R T Q R is te flow of te system 4.4.

13 4.1 Termoynamics of simple systems - continuous case 13 We efine te orizontal one-forms ω fr, ω ext Ω 1 T Q R associate to te friction an external forces by ω fr q, v, S δq, δv, δs := F fr q, v, S, δq, ω ext q, v, S δq, δv, δs := F ext q, v, S, δq, were δq, δv, δs T q,v,s T Q R. We also efine te one-form ω τ := T S on T Q R by ω τ q, v, S δq, δv, δs := T q, SδS = L q, v, SδS. S In orer to erive te structure preserving property, we sall exten te argument use in Marsen an West [1, 1..3]. Let us efine te restricte action map as Ŝ : T Q R R, Ŝq, v, S := Te erivative of tis map reas T L F t q, v, S t. Ŝq, v, S δq, δv, δs T L = q qt, qt, St L L t=t qt, qt, St, δqt t + qt, qt, St, δqt t q q t= + T L qt, qt, StδStt S T = F fr+ext qt, qt, St, δqt T T qt, StδStt T = Ft ω fr+ext q, v, S δq, δv, δs t t + Θ Lqt, qt, St δqt, δvt, δst T + FT Θ L Θ L q, v, S δq, δv, δs Ft ω τ q, v, S δq, δv, δs t, were we use te notations F fr+ext := F fr +F ext an ω fr+ext := ω fr +ω ext. Tus we obtain te relation T Ŝ = F T Θ L Θ L Ft ω fr+ext+τ t. as one-forms on T Q R, were ω fr+ext+τ := ω fr +ω ext +ω τ. By taking te exterior erivative of tis equality, we obtain te following result. Teorem 4.1. Consier a simple termoynamic system an assume tat te Lagrangian Lq, q, S is regular an te pysical assumptions 4. an 4.3 are verifie. Ten 4.4 efines a well-efine flow F t on T Q R. Tis flow verifies te following generalization of te symplectic property 4.1 of te flow in classical mecanics: t=t t= Note tat we can write tis property as T FT Ω L = Ω L Ft ω fr+ext+τ t. 4.5 F T Ω L = Ω L T F t ω fr+ext+τ t.

14 4. Termoynamics of simple systems - iscrete case Termoynamics of simple systems - iscrete case In tis section, we will sow tat te iscrete flow of our variational integrator satisfies a iscrete analogue of te property 4.5. We assume tat te iscrete termoynamical system satisfies te regularity criteria Tis ensures te existence of te iscrete flow F L : CK C K : qk, qk+1, Sk, Sk+1 CK q k+1, q k+, S k+1, S k+ CK, obtaine by solving te numerical sceme 3.1, namely, D 1 L q k, q k+1, S k, S k+1 + D L q k 1, q k, S k 1, S k +F fr + F ext q k, q k+1, S k, S k+1 +F fr+ + F ext+ q k 1, q k, S k 1, S k =, q k, q k+1, S k, S k+1 CK. 4.6 In orer to formulate te property of te iscrete flow, we nee to efine te following iscrete forms on Q Q R R. Definition 4.. Given a iscrete termoynamical system wit iscrete Lagrangian L an iscrete friction an external forces F fr± an F ext±, we efine te iscrete one-forms Θ L,F q, q 1, S, S 1 := D 1 L q, q 1, S, S 1 q F q, q 1, S, S 1 q, Θ + L,F q, q 1, S, S 1 := D L q, q 1, S, S 1 q 1 + F + q, q 1, S, S 1 q 1, were F ± := F fr± + F ext±, an te iscrete one-forms ω fr q, q 1, S, S 1 : = F fr q, q 1, S, S 1 q + F fr+ q, q 1, S, S 1 q 1, ω ext q, q 1, S, S 1 : = F ext q, q 1, S, S 1 q + F ext+ q, q 1, S, S 1 q 1, ω τ q, q 1, S, S 1 : = D 3 L q, q 1, S, S 1 S D 4 L q, q 1, S, S 1 S 1, wic are te iscrete analogue of te one-forms ω fr, ω ext, ω τ efine in te continuous case earlier. Te one-forms Θ ± L,F Ω 1 Q Q R R are relate to te canonical one-form Θ can Ω 1 T Q as Θ ± L,F = FL ± F Θ can, were FL ± F : Q Q R R T Q are te iscrete Legenre transforms wit force efine by FL F q, q 1, S, S 1 := q, D 1 L q, q 1, S, S 1 F q, q 1, S, S 1 = q, p, FL + F q, q 1, S, S 1 := q 1, D L q, q 1, S, S 1 + F + q, q 1, S, S 1 = q 1, p 1. Tese one-forms are te natural extensions of te one-forms for te iscrete Euler-Lagrange equations wit external forces consiere in Marsen an West [1]. We sow below tat te iscrete flow 3.14 satisfies a iscrete analogue of te property 4.5 of te continuous flow obtaine in Teorem 4.1. To obtain tis result, we exten te argument use in Marsen an West [1, 1.3.]. Similarly wit te continuous case earlier, we ientify te space of solutions of 3.1 wit te space of initial conitions q, q 1, S, S 1 an we efine te restricte iscrete action map Ŝ q, q 1, S, S 1 := N 1 k= L q k, q k+1, S k, S k+1,

15 4. Termoynamics of simple systems - iscrete case 15 were on te rigt an sie, te iscrete action functional is evaluate on te solution of 3.1 wit initial conitions q, q 1, S, S 1. Teorem 4.3. Consier te numerical sceme 4.6 arising from te iscrete variational formulation of Definition 3.7 for te nonequilibrium termoynamic of a simple system. Assume tat te regularity criteria 3.13 is verifie. Ten te sceme 4.6 inuces a well-efine iscrete flow F L : CK C K : q k, q k+1, S k, S k+1 C K q k+1, q k+, S k+1, S k+ C K. Moreover, tis flow verifies te following property N 1 Ω F + L L,F Ω L,F = N 1 k= F k L ω fr+ext+τ, 4.7 wic is a iscrete version of te property 4.5 of te flow of a simple an close system. Tis property is also an extension to nonequilibrium termoynamics of te symplectic property 3.4 of te flow of a variational integrator in classical mecanics. Proof. Te first part of te teorem as been proven above. We now prove formula 4.7. Using te notations L k := L q k, q k+1, S k, S k+1 an F k± := F fr± +F ext± q k, q k+1, S k, S k+1, we compute te erivative of Ŝ as Ŝq, q 1, S, S 1 δq, δq 1, δs, δs 1 = N 1 k= D 1L k δq k + D L k δq k+1 + D 3L k δs k + D 4L k δs k+1 N 1 = D 1L δq + D 1L k + D L k 1 δq k + D L N 1 δq N + k=1 N 1 = D 1L δq k=1 F k = D 1L δq + F δq + F k 1+ δq k + D L N 1 δq N + N 1 k= F k δq k + F k+ δq k+1 + F N 1+ δq N + D L N 1 δq N + = Θ L,F q, q 1, S, S 1δq, δq 1, δs, δs 1 N 1 k=1 N 1 k= N 1 k= N 1 k= D 3L k δs k + D 4L k δs k+1 D 3L k δs k + D 4L k δs k+1 D 3L k δs k + D 4L k δs k+1 ω fr+ext q k, q k+1, S k, S k+1 δq k, δq k+1, δs k, δs k+1 + Θ + L,F q N 1, q N, S N 1, S N δq N 1, δq N, δs N 1, δs N ω τ q k, q k+1, S k, S k+1 δq k, δq k+1, δs k, δs k+1. N 1 k= By using te notation F k L := F L... F L, }{{} k

16 5 Examples 16 we can write tis ifferential as Ŝ = F N 1 L Θ + L,F Θ L,F N 1 k= F k L ω fr+ext+τ. Taking te exterior erivative of tis relation, we ave te result. 5 Examples In tis section, we evelop several numerical scemes base on te variational integrator for te nonequilibrium termoynamics erive in Section 3 by consiering several stanar iscretizations of a given Lagrangian. Ten, we illustrate our scemes wit te example of te mass-spring-friction system moving in an ieal gas. 5.1 Variational iscretization scemes We consier tree stanar types of approximation of te time integral of a given Lagrangian. Tis leas to numerical scemes wic are extensions of te Verlet sceme, of te variational mipoint rule sceme as well as of te symmetrize Lagrangian variational integrator. Let us assume Q = R n. Variational sceme 1. Let us first coose te finite ifference map ϕ : Q Q R R T Q R as ϕq k, q k+1, S k, S k+1 = q k, S k, q k+1 q k, S k+1 S k. For a given Lagrangian Lq, q, S, te iscrete Lagrangian in 3.1 tus reas L q k, q k+1, S k, S k+1 = L q k, q k+1 q k, S k an te iscrete penomenological constraint 3.9 is given ere by L S q k, S k S k+1 S k = F fr q k, q k+1 q k qk+1 q k, S k. Te natural iscretization of te forces F fr an F ext associate to tis iscretization of te Lagrangian may be given as follows see Marsen an West [1, 3..5]: F fr q k, q k+1, S k, S k+1 = F fr F ext q k, q k+1, S k, S k+1 = F ext q k, q k+1 q k q k, q k+1 q k, S k, F fr+ q k, q k+1, S k, S k+1 =, Te first equation in 3.1 tus becomes [ 1 L q k, q k+1 q k, S k L q k 1, q ] k q k 1, S k 1 L v v q = F fr q k, q k+1 q k, S k + F ext q k, q k+1 q k, S k., S k, F ext+ q k, q k+1, S k, S k+1 =. q k, q k+1 q k, S k

17 5.1 Variational iscretization scemes 17 For te stanar Lagrangian were v = q, we obtain te following numerical sceme: Lq, v, S = 1 m v Uq, S, 5.1 Sceme 1: m q k+1 q k + q k 1 + U q q k, S k = F fr q k, q k+1 q k, S k + F ext q k, q k+1 q k, S k, U S q k, S k S k+1 S k = F fr q k, q k+1 q k qk+1 q k, S k. Tis is an extension of te Verlet sceme to nonequilibrium termoynamics. Te matrix 3.13 for Sceme 1 as te entries: A 11 = m + F v, F := F fr + F ext, A 1 =, 5. A 1 = 1 F fr q 1 q + 1 v F fr, A = 1 U S, were F ext = F ext q, q1 q, S, F fr = F fr q, q1 q, S an U = Uq, S. Te regularity criteria 3.13 is tus satisfie if an only if U S q, S 1 an m fr F + q, q 1 q, S. v Te first conition is always satisfie uner te pysical assumption 4.3. Te secon conition is satisfie for all friction forces tat are linear in velocity. Variational sceme. More generally, we can coose a finite ifference map of te form ϕq k, q k+1, S k, S k+1 = 1 αq k + αq k+1, q k+1 q k, 1 αs k + αs k+1, S k+1 S k for some parameter α [, 1]. For α = 1, we ave qk + q k+1 L q k, q k+1, S k, S k+1 := L, q k+1 q k, S k + S k+1. Te natural iscretization of te force F fr associate to tis iscretization of te Lagrangian is see Marsen an West [1, 3..5] F fr q k, q k+1, S k, S k+1 = 1 F fr qk + q k+1, q k+1 q k, S k + S k+1 = F fr+ q k, q k+1, S k, S k+1, similarly for F ext. Te iscrete penomenological constraint 3.9 is given ere by L qk + q k+1, S k + S k+1 Sk+1 S k S = F fr qk + q k+1, q k+1 q k, S k + S k+1 qk+1 q k.

18 5.1 Variational iscretization scemes 18 Te first equation in 3.1 is [ 1 L qk + q k+1 v 1 [ L qk 1 + q k q = 1 F fr qk + q k F ext qk + q k+1, q k+1 q k, S k + S k+1 L qk 1 + q k v, q k q k 1, S k 1 + S k + L qk + q k+1 q, q k+1 q k, S k + S k+1, q k+1 q k, S k + S k F fr qk 1 + q k + 1 F ext qk 1 + q k, q k q k 1, S ] k 1 + S k, q k+1 q k, S ] k + S k+1, q k q k 1, S k 1 + S k, q k q k 1, S k 1 + S k Te corresponing expressions for arbitrary α [, 1] are erive similarly. For te stanar Lagrangian 5.1 we obtain te following numerical sceme: Sceme : m q k+1 q k + q k = 1 F fr qk + q k F ext qk + q k+1 [ U qk 1 + q k, S k 1 + S k q, q k+1 q k, S k + S k+1, q k+1 q k, S k + S k+1 [ 1 U qk 1 + q k, S k 1 + S k + U qk + q k+1 S S + U qk + q k+1, S ] k + S k+1 q + 1 F fr qk 1 + q k + 1 F ext qk 1 + q k, S ] k + S k+1 Sk+1 S k, q k q k 1, S k 1 + S k, q k q k 1, S k 1 + S k., = F fr qk + q k+1, q k+1 q k, S k + S k+1 qk+1 q k. Tis is an extension to nonequilibrium termoynamics of te variational mipoint rule sceme. Te matrix 3.13 as te entries: A 11 = m U q + F 4 q + 1 F v, F := F fr + F ext, A 1 = 4 A 1 = 1 U S q U S q + F 4 S, S 1 S + 1 F fr q + 1 F fr q1 q + 1 v F fr, 5.3 were A = 1 U S 1 S + 1 U S S + 1 F fr S q + q 1 U = U, S + S 1 q 1 q, an F fr = F fr q + q 1, q 1 q, S + S 1.

19 5.1 Variational iscretization scemes 19 Variational sceme 3. We can coose to approximate te Lagrangian by te symmetrize iscrete Lagrangian as follows L q k, q k+1, S k, S k+1 = 1 L q k, q k+1 q k, S k + 1 L q k+1, q k+1 q k, S k+1. Te associate natural coice of iscrete forces is given by F fr q k, q k+1, S k, S k+1 = 1 F fr q k, q k+1 q k, S k, F fr+ q k, q k+1, S k, S k+1 = 1 F fr q k+1, q k+1 q k, S k+1, similarly for F ext± ; see Kane, Marsen, Ortiz, an West [, p. 9] α =. A natural iscrete penomenological constraint is given ere by [ L S q k, S k + L ] S q Sk+1 S k k+1, S k+1 = [F fr q k, q k+1 q k, S k + F fr q k+1, q ] k+1 q k qk+1 q k, S k+1. Te first equation in 3.1 is 1 1 [ L v q k, q k+1 q k 1 [ L q + L v = 1 F fr q k, q k+1 q k, S k q k, q k q k 1 q k+1, q k+1 q k + L q L v q k, q k+1 q k + F fr q k, q k q k 1, S k + 1 F ext q k, q k+1 q k, S k q k 1, q k q k 1 ] + F ext q k, q k q k 1, S k. L v q k, q k q k 1 ] For te stanar Lagrangian 5.1 we obtain te following numerical sceme: Sceme 3: m q k+1 q k + q k [ + U q q k, S k = 1 [ F ext q k, q k+1 q k, S k F fr q k, q k+1 q k, S k + F ext q k, q k q ] k 1, S k, + F fr q k, q k q ] k 1, S k [ U S q k, S k + U ] S q Sk+1 S k k+1, S k+1 = [F fr q k, q k+1 q k, S k + F fr q k+1, q k+1 q k ] qk+1 q k, S k+1.

20 5. Example: a mass-spring-friction system moving in an ieal gas Tis is a symmetrize Lagrangian variational integrator applie to nonequilibrium termoynamics. Te matrix 3.13 for Sceme 3 as te entries: A 11 = m + 1 F v, fr F := F + F ext, A 1 = 1 F S, [ A 1 = U 1 S 1 S 1 + S q F fr v + 1 F fr 1 v A = U 1 S 1 S U + S S + U1 1 S ] fr F + 1 q1 q + F fr + F1 fr, q + F 1 fr S q 1 q, 5.4 were F = F q, q1 q, S, F fr = F fr q, q1 q, S, F fr 1 = F fr q 1, q1 q, S 1, U = Uq, S an U 1 = Uq 1, S Example: a mass-spring-friction system moving in an ieal gas We consier te example of a mass-spring-friction system moving in a close room fille wit an ieal gas. Tis system, enote Σ, is illustrate in Fig.5.1. We refer to Ferrari an Gruber [1] for te erivation of te equations of motion for tis system from Stueckelberg s point of view. We consier tis simple example since we can take avantage of te fact tat te equations of evolution for tis system can be explicitly solve. Tis allows us to easily estimate te numerical valiity of te sceme to simulate te entropy, temperature, an internal energy beaviors. Continuous setting. Te Lagrangian L for te system Σ is given by Lx, ẋ, S = 1 mẋ Ux, S, were x, ẋ, S enotes te state of te system, m is te mass of te soli, an k is te spring constant, an were Ux, S = 1 kx + US is te potential energy. Te internal energy of te ieal gas is given by U = cnrt, were c is te gas constant, N is te number of moles, R is te universal gas constant, an T is te temperature. Note tat U may be rewritten as a function 1 US, N, V = U e 1 S cr N S N N c +1 1 V c, N V were U inicates te initial value of te internal energy, N is te initial mole number of te ieal gas, an V is te volume of te room wit te ieal gas, wic is assume to be constant, i.e., V = V. We assume tat te friction force is given by a viscous force as F fr x, ẋ, S = λẋ, were λ is te penomenological coefficient etermine experimentally an also tat te system Σ is subject to an external force F ext exerte from te exterior Σ ext. We also assume tat te system is aiabatically close, so te power ue to eat transfer between te system an te exterior is zero, i.e., PH ext = an tere is no cange in te number of moles of te gas, i.e., N = N. Wit te above coice of Lagrangian an force, te general equations.4 arising from te variational formulation.1.3 yiel te time evolution equations of te couple mecanical an termal system as mẍ = kx λẋ + F ext x, ẋ, S, Ṡ = 1 T λẋ, 3 For simplicity, we neglect te internal energy of te soli. It is given by U s = 3N srt, were N s is number of mole of te soli, an it can be easily inclue in our iscussion.

21 5. Example: a mass-spring-friction system moving in an ieal gas 1 gas k F fr F k m x ext F ext Figure 5.1: A mass-spring-friction system in a room wit gas were T = U S S = U cn R e 1 crn S S = T e 1 crn S S. Te total energy, given by Ex, ẋ, S = 1 mẋ +Ux, S, verifies te energy balance equation t E = F ext x, ẋ, S, ẋ. Exact solutions. Consier te special case in wic tere is no external force, i.e., F ext =. In tis case te time evolution equations are given by mẍ = kx λẋ an T Ṡ = λẋ, 5.5 an te total energy is preserve. Tese equations can be easily solve explicitly, see Ferrari an Gruber [1]. were Setting x = x an v = ẋ, te solution of te first equation in 5.5 is xt = e κt x cosωt + v + κx sinωt, 5.6 ω κ = λ k m, ω = m, ω = ω κ, κ < ω. In orer to solve te secon equation in 5.5, we note tat t U = T Ṡ = λẋ an from wic we obtain te evolution of te temperature as In te above, an ence ft = λ = t ẋ s 1 mv + 1 kx T t = T + 1 cnr ft, t U = cnr T, 5.7 were ft := λ t ẋt = e v κt cosωt κ v + κx + x ω sinωt ω 1 4km λ e λ m t 4kmmv + λv x + kx ẋ ss. 5.8 λvλm + 4v mkx + λkx cosωt λmv kx 4km λ 1/ sinωt.

22 5. Example: a mass-spring-friction system moving in an ieal gas Using 5.8 an 5.7, we get te explicit evolution of te entropy as T t St = S + cnr ln. 5.9 T Variational iscretizations. We now apply our variational iscretization scemes 1,, 3 to tis example. Sceme 1: Te first sceme yiels m x k+1 x k + x k 1 U S S k S k+1 S k + kx k = λ x k+1 x k, xk+1 x k = λ. Te secon equation can be restate as were S k+1 = λ T k xk+1 x k + S k, T k = U S S k = U cn R e 1 crn S k S = T e 1 crn S k S. For tis example, te matrix 5. of te variational iscretization sceme 1 reas m λ 1 T x, S Tus, since m >, λ, an T >, te iscrete flow of te extene Verlet sceme is well-efine.. Sceme : Te secon sceme yiels m x k+1 x k + x k k xk 1 + x k + x k + x k+1 = 1 λ xk+1 x k + x k x k 1, [ 1 U S Sk 1 + S k + U S ] Sk + S k+1 Sk+1 S k = λ xk+1 x k. For tis example, te matrix 5.3 of te variational iscretization sceme reas m λ λ x 1 x T S1 S cn R + 1 were T = U S+S 1 S. Tus, since m >, λ, T >, an cn R >, te iscrete flow of te variational mipoint rule sceme is well-efine.,

23 5.3 Numerical tests 3 Sceme 3: Te tir example yiels m x k+1 x k + x k 1 + kx k = 1 λ xk+1 x k + x k x k 1, [ 1 U S S k + U ] S S Sk+1 S k xk+1 x k k+1 = λ. For tis example, te matrix 5.4 of te variational iscretization sceme 3 reas m λ 4λ x 1 x T 1 cn R S1 S + T + T 1 1 were T = U S S an T 1 = U S S 1. Tus, since m >, λ, T 1 >, T > an cn R >, te iscrete flow of te symmetrize Lagrangian variational integrator is well-efine. 5.3 Numerical tests We illustrate te beavior of te tree variational scemes for te mass-spring-friction system, by consiering two cases of pysical parameters an various values for te friction coefficient, namely λ =,., 5, an 1 [N s/m]. For eac of te five values of λ we isplay te evolutions of te position, entropy, total energy Eq k, q k+1, S k, relative energy errors Eq k,q k+1,s k E E, internal energy, an temperature. Eac figure sows te results for te tree scemes as well as te exact solution, troug 1 5 time steps. Case 1. For te first series of numerical tests, we coose te time step = 1 3 [s] an we set te parameters of te system Σ as follows: m = 5 [kg], N = 1 [mol], k = 5 [N/m], V = [m 3 ]. Te initial conitions are x =.3 [m], x 1 =.3 [m], T = 3 [K], S = [J/K]. Te results for λ =, see Figures 5. an 5.3, consistently recover te beavior obtaine troug a usual variational iscretization of te Euler-Lagrange equations for te conservative mass-spring system in classical mecanics. In particular, for eac sceme te internal energy US k = U is preserve an te temperature, given by T k = U S S k, remains a constant, see Figure 5.4. Exactly as in te continuous case, in absence of friction in an isolate simple system, te entropy an temperature stay constant, te system is reversible, an te ynamics is completely escribe by te Euler-Lagrange equations. For all te cases wit friction, λ =., 5, 1, te numerical solutions of te position, entropy, internal energy, an temperature reprouce te correct beaviors for all te tree scemes, as we see from a irect comparison wit te exact solutions, see Figures 5.5, 5.7, 5.8, 5.1, 5.11, For Sceme 1, te relative energy error is boune by 1 8 for all values of λ, an ecreases in time, wereas for Sceme an 3, te relative energy error is boune by 1 11 for all values of λ, see Figures 5.6, 5.9, 5.1, an stays constant in time.,

24 5.3 Numerical tests Sceme Entropy [J/K] Position [m] Sceme 1 Sceme Sceme 3 Sceme 1 Sceme 3 Figure 5.: Time evolutions of position an entropy Case 1: λ = Sceme Sceme 3 Energy Errors Total Energy [kj] Sceme 1 Sceme 1 Sceme 3 Sceme Sceme 1 Sceme 3 Sceme Temperature [K] Internal Energy [J] Figure 5.3: Total energy an relative energy error Case 1: λ = Sceme Sceme 1 Sceme 3 Figure 5.4: Internal energy an temperature Case 1: λ = 4

25 5.3 Numerical tests Sceme Position [m] Entropy [J/K] Sceme 1 Sceme Sceme 3 Sceme 1 Sceme 3 Figure 5.5: Time evolutions of position an entropy Case 1: λ =. Sceme 3 Sceme Energy Errors Total Energy [kj] Sceme 1 Sceme 1 Sceme Sceme 3 Figure 5.6: Total energy an relative energy error Case 1: λ =. Sceme 1 Sceme Sceme 3 Temperature [K] Internal Energy [J] Sceme 3 Sceme Sceme 1 Figure 5.7: Internal energy an temperature Case 1: λ =. 5

26 5.3 Numerical tests 6 Position [m] Sceme 1 Sceme Sceme 3 Entropy [J/K] Sceme 3 Sceme Sceme 1 Figure 5.8: Time evolutions of position an entropy Case 1: λ = 5 Total Energy [kj] Sceme 1 Sceme 3 Sceme Energy Errors Sceme 1 Sceme Sceme 3 Figure 5.9: Total energy an relative energy error Case 1: λ = 5 Sceme 1 Internal Energy [J] Sceme 3 Sceme Temperature [K] Sceme 3 Sceme Sceme 1 Figure 5.1: Internal energy an temperature Case 1: λ = 5

27 5.3 Numerical tests 7 Position [m] Sceme 1 Sceme Sceme 3 Entropy [J/K] Sceme Sceme 3 Sceme 1 Figure 5.11: Time evolutions of position an entropy Case 1: λ = 1 Total Energy [kj] Sceme 1 Sceme 3 Sceme Energy Errors Sceme Sceme 3 Sceme 1 Figure 5.1: Total energy an relative energy error Case 1: λ = 1 Internal Energy [J] Sceme Sceme 1 Sceme 3 Temperature [K] Sceme 3 Sceme Sceme 1 Figure 5.13: Internal energy an temperature Case 1: λ = 1

28 5.3 Numerical tests 8 Case. For te secon series of numerical tests, we coose te time step = 1 3 [s] an we set te parameters of te system Σ as follows: m = 1 [kg], N = [mol], k = [N/m], V = [m 3 ]. Te initial conitions are x =.1 [m], x 1 =.1 [m], T = 3 [K], S = [J/K]. Te results for λ = are sown in Figures 5.14 an 5.15, an, similarly to Case 1, consistently recover te beavior obtaine troug a usual variational iscretization of te Euler-Lagrange equations. In particular, te internal energy US k = U is preserve an te temperature, given by T k = U S S k, remains a constant, see Figure For all te cases wit friction, λ =., 5, 1, te numerical solutions of te position, entropy, internal energy, an temperature reprouce te correct beaviors for all te tree scemes, as we see from a irect comparison wit te exact solutions, see Figures 5.17, 5.19, 5., 5., 5.3, 5.5. For Sceme 1, te relative energy error is boune by 1 6 for all values of λ, an ecreases in time, wereas for Sceme an 3, te relative energy error is boune by 1 9 for all values of λ, see Figures 5.18, 5.1, 5.4, an stays constant in time. Sceme 1 Sceme Position [m] Entropy [J/K] Sceme Sceme 3 Sceme 1 Sceme 3 Figure 5.14: Time evolutions of position an entropy Case : λ = Sceme 1 Total Energy [kj] Sceme 1 Sceme 3 Sceme Energy Errors Sceme 3 Sceme Figure 5.15: Total energy an relative energy error Case : λ =

29 5.3 Numerical tests Sceme 3 Sceme Temperature [K] Internal Energy [J] Sceme 1 Sceme Sceme 1 Sceme 3 Figure 5.16: Internal energy an temperature Case : λ = Sceme 1 Sceme Position [m] Entropy [J/K] Sceme Sceme3 Sceme 3 Sceme1 Figure 5.17: Time evolutions of position an entropy Case : λ = Sceme Sceme Sceme Sceme Energy Errors Total Energy [kj] Sceme 3 Sceme 1 Figure 5.18: Total energy an relative energy error Case : λ =. 9