Journal of Physics: Conference Series PAPER OPEN ACCESS Control of Stirling engine. Simplified, compressible model To cite this article: P I Plotnikov et al 216 J. Phys.: Conf. Ser. 722 1232 View the article online for updates and enhancements. This content was downloaded from IP address 37.44.192.194 on 29/12/217 at 17:
Journal of Physics: Conference Series 722 (216) 1232 doi:1.188/1742-6596/722/1/1232 Control of Stirling engine. Simplified, compressible model P I Plotnikov 1,2, J Soko lowski 3 and A Żochowski 4 1 Lavrentyev Institute of Hydrodynamics, Novosibirsk, Russia 2 Novosibirsk State University, Novosibirsk, Russia 3 Université de Lorraine, Institut Élie Cartan, Nancy, France 4 Instytut Badan Systemowych Polskiej Akademii Nauk, Warszawa, Polska E-mail: plotnikov@hydro.nsc.ru Abstract. A one-dimensional free boundary problem on a motion of a heavy piston in a tube filled with viscous gas is considered. The system of governing equations and boundary conditions is derived. The obtained system of differential equations can be regarded as a mathematical model of an exterior combustion engine. The existence of a weak solution to this model is proved. The problem of maximization of the total work of the engine is considered. 1. Problem formulation The paper deals with a free boundary problem which describes the motion of a heavy piston in a tube filled with viscous gas. It is assumed that the viscous gas occupies the interval (, 1). The piston is modeled by a material point a (, 1). Its motion is described by a function x = a(t), t (, T ). We also assume that the dynamics of the piston is governed by the equation mä(t) + G(ȧ(t), a(t)) = F, (1) where F is the hydrodynamical force acting on the piston, G is a given function. This means that the piston is connected with some mechanical device which dynamics is described by the function G. The state of the gas is completely described by the density ϱ(x, t), the velocity u(x, t), and the temperature ϑ(x, t). Among the other thermodynamic quantities there are the internal energy e, the pressure p, and the entropy s. The three quantities are assumed to be functions of state variable ϱ and θ. Typically, the pressure formula admits the representation The internal energy is defined by the relation which yields the representation p = p e (ϱ) + p θ (ϱ)θ. e ϱ = p e(ϱ) ϱ 2, e = P e (ϱ) + Q(θ), P e (ϱ) = ϱ 1 z 2 p e (z)dz, Content from this work may be used under the terms of the Creative Commons Attribution 3. licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Published under licence by Ltd 1
Journal of Physics: Conference Series 722 (216) 1232 doi:1.188/1742-6596/722/1/1232 where Q is an arbitrary function of θ. The entropy s(ϱ, θ) is defined by the relations s = θ 1 Q (z) ϱ dz P θ (ϱ), P θ (ϱ) = z 1 p θ (z) z 2 dz. In the simple case of barotropic gas the pressure is a function of the density, p = p(ϱ), and the entropy s equals zero. This eliminates the temperature from the equations. In this framework, the problem can be formulated as follows. Denote by Ω the rectangular (, 1) (, T ) and set Ω = {(x, t) : < x < a(t), < t < T }, Ω + = {(x, t) : a(t) < x < 1, < t < T } The problem is to find the velocity, the density, the temperature, and the trajectory of the piston satisfying the following equations and boundary conditions (ϱu) + t ( ( u 2 ϱ t ϱ t + x (ϱu) = in Ω±, (2) x (ϱu2 ) + x (p νu x) = in Ω ±, (3) )) 2 + e + x ( u (ϱ u2 2 + ϱe + p θ (κ(θ) x x ) ν x (uu x) = in Ω ±, u(a(t), t) = ȧ(t) in (, T ), (5) mä(t) + G(ȧ(t), a(t)) = (T + (a(t), t) T (a(t), t)) in (, T ), (6) ( κ θ ) + ( (a(t), t) κ θ ) (a(t), t) = (T (t) T + (t))u(a(t), t) in (, T ), (7) x x )) θ + (a(t), t) = θ (a(t), t)). (8) Here κ(θ) is a given positive function, mostly a constant. The stress is defined by the equality We also use the notation T ± (t) := T = νu x p. lim T(x, t), x a(t)± ϑ± (t) := lim ϑ(x, t). (9) x a(t)± The system of equations (2)-(8) should be supplemented with the boundary conditions. Choose functions U ± : (, T ) R, Θ ± : (, T ) R +, and R ± : (, T ) R +. Denote by Σ ± in [, T ] the inlets and by Σ ± out [, T ] the outlets defined by U, We set the boundary conditions Σ in Σ + in = {t (, T ) : U(, t) > }, = {t (, T ) : U(1, t) < }, (1) Σ ± out = (, T ) \ Σ± in u(, t) = U (t), u(1, t) = U + (t) for t (, T ), (11) ϱ(, t) = R (t) on Σ in, (12) ϱ(1, t) = R + (t) on Σ + in, (13) θ(, t) = Θ (t), θ(1, t) = Θ + (t) for t (, T ). (14) (4) 2
Journal of Physics: Conference Series 722 (216) 1232 doi:1.188/1742-6596/722/1/1232 At the initial moment t = the velocity and the density distributions are prescribed u(x, ) = u (x), ϱ(x, ) = R (x), θ(, x) = Θ (x) for x (, 1). (15) The initial position of the piston is given a() = a. (16) The equations, boundary, and initial conditions (2) (16) form a closed system of equations for the velocity, the density, the temperature, and the trajectory of the piston. This system can be regarded as a mathematical model of an exterior combustion engine often called the Stirling engine. An intensive treatment of one-dimensional compressible Navier-Stokes equations starts with pioneering book by Antontsev, Kazhikhov, and Monakhov [1] on the global theory for onedimensional homogeneous problems. Although the theory of strong solutions is satisfactory in what concerns the specific boundary conditions, many issues of global behavior of solutions for large data in general boundary conditions are far from being understood. A global theory of weak solutions was developed by P.L. Lions. We refer the reader to the books by Lions [3], Feireisl [2], and Novotný & Straškraba [4] for the state of the art in the domain. The theory of weak solutions to inhomogeneous boundary value problems for compressible Navier-Stokes equations was developed in the book by Plotnikov & Sokolowski [5]. The problem on the motion of the piston was studied in the paper by Shelukhin [6] under the assumptions that U ± = and G =. In this paper we are focused on the case when the thermodynamical process is adiabatic and the dynamics of the piston is conservative. This means that p = cϱ γ, γ > 1, G(ȧ, a) = V (a), (17) where V : (, 1) R is a mechanical potential. In this case equations (2) (16) read (ϱu) t ϱ t + x (ϱu) = in Ω±, (18) + x (ϱu2 ) + x (p νu x) = in Ω ±, (19) u(a(t), t) = ȧ(t) in (, T ), (2) mä(t) + V (a(t)) = (T + (a(t), t) T (a(t), t)) in (, T ). (21) u(, t) = U (t), u(1, t) = U + (t) for t (, T ), (22) ϱ(, t) = R (t) on Σ in, (23) ϱ(1, t) = R + (t) on Σ + in, (24) u(x, ) = u (x), ϱ(x, ) = R (x) >, for x (, 1), a() = a. (25) Furthermore, we assume that the potential V conditions H.1 V : (, 1) R + is a smooth positive function such that and the given data satisfy the following lim V (a) = lim V (a) =. (26) a + a 1 3
Journal of Physics: Conference Series 722 (216) 1232 doi:1.188/1742-6596/722/1/1232 H.2 The functions U ± belong to the class C 2 [, T ]. Every set Σ ± in consists of finite number of open intervals H.3 The functions R ± L (, T ) are nonnegative. The function R L (, 1) is strictly positive. Moreover, u L (, 1), < a < 1. Condition (26) prevents the destruction of the mechanical system. 2. Weak solution. Existence theorem. The triplet (ϱ, u) : Ω R + R, a : (, T ) (, 1) is a weak solution to problem (18) (25) if the following conditions are satisfied: The state variables satisfy The integral identity ϱ, p, ϱ u 2 L (, T ; L 1 (Ω)), u L 2 (, T ; W 1,2 (, 1)), a W 1,2 (, 1) (27) Ω u(, t) = U (t), u(1, t) = U + (t) on (, T ). (28) 1 ϱuζ t + ϱu 2 ζ x + (p νu x )ζ x dxdt + +mu (a )ζ(a ) + holds for all functions ζ C 1 (Ω) such that The integral identity (ϱ t ψ + ϱu ψ) dxdt = Ω ϱ (x)u (x)ζ(x, )dx (29) ( mȧ(t) 2 x ζ(a(t), t) V (a(t))ζ(a(t), t) ) dt = ζ(x, T ) =, ζ(, t) = ζ(1, t) =. holds for all functions ψ C 1 (Ω) satisfying ψr + U + dt ψr U dt 1 ψ = in a neighborhood of Σ ± out {(, 1) {T }}, R ψ(, ) dx (3) which means that the support of the test function ψ meets the boundary of the rectangular Ω at the inlet Σ ± in and the bottom (, 1) {}. The first main result of the study is the following Theorem 1 Assume that Conditions (H.1) (H.3) are fulfilled. Then problem (18) (25) has a weak solution which satisfies the estimates u L 2 (,T ;W 1,2 (,1)) + ϱ u 2 L (,T ;L 1 (,1)) + p L (,T ;L 1 (,1)) c, sup ȧ(t) 2 + V (a(t)) c, t (,T ) where the constant c depends only on T as well as the boundary and initial data. The proof is based on the methods developed in [5]. Following [5] we use the multistep regularization. The weak solutions is obtained as a limit of a sequence of approximate solutions to regularized equations. The main difficulty is the derivation of estimate (31). (31) 4
Journal of Physics: Conference Series 722 (216) 1232 doi:1.188/1742-6596/722/1/1232 3. The optimization problem. Recall that the power of the hydrodynamical forces acting on the piston is equal to (T + (a(t), t) T (a(t), t))ȧ(t). Assume that the functions U ± and u, R are fixed and satisfy conditions (H.1) (H.3). In this case the total work W of hydrodynamical forces acting on the piston becomes the function of R ±, W (R ± ) = (T + (a(t), t) T (a(t), t))ȧ(t)dt (32) In view of Theorem 1 the functional W is well defined on the set of bounded nonnegative functions. Fix an arbitrary N > and consider the set of admissible data M = {R ± L (, T ) : R ± L (,T ) N}. The natural question is to find the boundary data which maximize the total work of the engine. This problem can be formulated as follows: to find an admissible couple R ± M such that W (R ± ) = sup W (R ± ), (33) R ± M where the cost function is given by (32). The following theorem on solvability of variational problem (33) is the second main result of the study. Theorem 1 Assume that conditions (H.1) (H.3) are fulfilled. Then problem (33) has at least one solution R ± M. The proof is based on Theorem 1. The key observation is that the functional W is star weakly continuous on M. References [1] Antontsev S N, Kazhikhov A V and Monakhov V N 199 Boundary Value Problems in Mechanics of Nonhomogeneous Fluids (Amsterdam: North-Holland) [2] Feireisl E 24 Dynamics of Viscous Compressible Fluids, (Oxford: Oxford Univ. Press) [3] Lions P L 1998 Mathematical Topics in Fluid Dynamics; Compressible Models (Oxford: Clarendon Press). [4] Novotný A and Straškraba I 24 Introduction to the Mathematical Theory of Compressible Flow (Oxford: Oxford Univ. Press) [5] Plotnikov P I and Sokolowski J 212 Compressible Navier-Stokes equations. Theory and shape optimization (Basel: Springer) [6] Shelukhin V V 1982 Dynamica sploshnoi sredy [in Russian] 57 131 52 5