Research Papers

A Lumped-Parameter Model of Multiscale Dynamics in Steam Supply Systems

[+] Author and Article Information
Hikaru Hoshino

Department of Electrical Engineering,
Kyoto University,
Katsura, Nishikyo,
Kyoto 615-8510, Japan
e-mail: hoshino@dove.kuee.kyoto-u.ac.jp

Yoshihiko Susuki

Department of Electrical Engineering,
Kyoto University,
Katsura, Nishikyo,
Kyoto 615-8510, Japan
e-mail: susuki@ieee.org

Takashi Hikihara

Department of Electrical Engineering,
Kyoto University,
Katsura, Nishikyo,
Kyoto 615-8510, Japan
e-mail: hikihara.takashi.2n@kyoto-u.ac.jp

1Present address: Department of Electrical and Information Systems, Osaka Prefecture University, 1-1 Gakuencho, Nakaku, Sakai, Osaka 599-8531, Japan and JST CREST, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received January 18, 2016; final manuscript received August 7, 2016; published online September 16, 2016. Assoc. Editor: Stefano Lenci.

J. Comput. Nonlinear Dynam 11(6), 061018 (Sep 16, 2016) (9 pages) Paper No: CND-16-1023; doi: 10.1115/1.4034491 History: Received January 18, 2016; Revised August 07, 2016

This paper focuses on multiscale dynamics occurring in steam supply systems. The dynamics of interest are originally described by a distributed-parameter model for fast steam flows over a pipe network coupled with a lumped-parameter model for slow internal dynamics of boilers. We derive a lumped-parameter model for the dynamics through physically relevant approximations. The derived model is then analyzed theoretically and numerically in terms of existence of normally hyperbolic invariant manifold in the phase space of the model. The existence of the manifold is a dynamical evidence that the derived model preserves the slow–fast dynamics, and suggests a separation principle of short-term and long-term operations of steam supply systems, which is analog to electric power systems. We also quantitatively verify the correctness of the derived model by comparison with brute-force simulation of the original model.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Geidl, M. , Koeppel, G. , Favre-Perrod, P. , Klöckl, B. , Andersson, G. , and Fröhlich, K. , 2007, “ Energy Hubs for the Future,” IEEE Power Energy Mag., 5(1), pp. 24–30. [CrossRef]
O'Malley, M. , and Kroposki, B. , 2013, “ Energy Comes Together: The Integration of All Systems,” IEEE Power Energy Mag., 11(5), pp. 18–23. [CrossRef]
Hoshino, H. , Susuki, Y. , and Hikihara, T. , 2014, “ A Nonlinear Dynamical Model of Two-Sites Electricity and Heat Supply System,” International Symposium on Nonlinear Theory and Its Applications, pp. 482–485.
Hoshino, H. , and Susuki, Y. , 2015, “ Graph-Based Modeling and Analysis of Dynamic Flows in Steam Supply Networks,” 54th IEEE Conference on Decision and Control, pp. 1358–1363.
IEA, 2011, “ Co-Generation and Renewables: Solutions for a Low-Carbon Energy Future,” International Energy Agency, Paris, France.
Rezaie, B. , and Rosen, M. A. , 2012, “ District Heating and Cooling: Review of Technology and Potential Enhancements,” Appl. Energy, 93, pp. 2–10. [CrossRef]
Shinji, T. , Sekine, T. , Akisawa, A. , Kashiwagi, T. , Fujita, G. , and Matsubara, M. , 2008, “ Reduction of Power Fluctuation by Distributed Generation in Micro Grid,” Electr. Eng. Jpn., 163(2), pp. 22–29. [CrossRef]
Mueller, S. , Tuth, R. , Fischer, D. , Wille-Haussmann, B. , and Wittwer, C. , 2014, “ Balancing Fluctuating Renewable Energy Generation Using Cogeneration and Heat Pump Systems,” Energy Technol., 2(1), pp. 83–89. [CrossRef]
Åström, K. J. , and Bell, R. D. , 2000, “ Drum-Boiler Dynamics,” Automatica, 36(3), pp. 363–378. [CrossRef]
Kim, T. S. , Lee, D. K. , and Ro, S. T. , 2000, “ Dynamic Behavior Analysis of a Heat Recovery Steam Generator During Start-Up,” Int. J. Energy Res., 24(2), pp. 137–149. [CrossRef]
Wen, C. , and Ydstie, B. , 2009, “ Passivity Based Control of Drum Boiler,” American Control Conference, pp. 1586–1591.
Bujak, J. , 2009, “ Optimal Control of Energy Losses in Multi-Boiler Steam Systems,” Energy, 34(9), pp. 1260–1270. [CrossRef]
Osiadacz, A. J. , 1987, Simulation and Analysis of Gas Networks, E. F. N. Spon, London.
Alobaid, F. , Postler, R. , Ströhle, J. , Epple, B. , and Kim, H.-G. , 2008, “ Modeling and Investigation Start-Up Procedures of a Combined Cycle Power Plant,” Appl. Energy, 85(12), pp. 1173–1189. [CrossRef]
Liu, Q. , Zhang, Z. , Pan, J. , and Guo, J. , 2009, “ A Coupled Thermo-Hydraulic Model for Steam Flow in Pipe Networks,” J. Hydrodyn., Ser. B, 21(6), pp. 861–866. [CrossRef]
Verma, M. P. , 2013, “ Steam Transport Simulation in a Geothermal Pipeline Network Constrained by Internally Consistent Thermodynamic Properties of Water,” Rev. Mex. Cienc. Geol., 30(1), pp. 210–221.
Kevorkian, J. , and Cole, J. , 1996, Multiple Scale and Singular Perturbation Methods, Springer-Verlag, Berlin.
Wiggins, S. , 1994, Normally Hyperbolic Invariant Manifolds in Dynamical Systems, Springer-Verlag, Berlin.
Nishiura, Y. , 2002, Far-From-Equilibrium Dynamics, American Mathematical Society, Providence, RI.
Iri, M. , 1969, Network Flow, Transportation and Scheduling: Theory and Algorithms, Academic Press, New York.
Pozrikidis, C. , 2014, An Introduction to Grids, Graphs, and Networks, Oxford University Press, Oxford, UK.
Machowski, J. , Bialek, J. W. , and Bumby, J. R. , 2008, Power System Dynamics: Stability and Control, 2nd ed., Wiley, New York.
Landau, L. D. , and Lifshitz, E. M. , 1959, Fluid Mechanics, Pergamon, Oxford, UK.
Traviss, D. P. , Baron, A. G. , and Rohsenow, W. M. , 1971, “ Forced-Convection Condensation Inside Tubes,” Heat Transfer Laboratory, Massachusetts Institute of Technology, Cambridge, MA, Technical Report No. 74.
JSME, 1979, Hydraulic Losses in Pipes and Ducts, Maruzen, Tokyo, Japan (in Japanese).
Bergman, T. , Lavine, A. , Incropera, F. , and Dewitt, D. , 2011, Fundamentals of Heat and Mass Transfer, 7th ed., Wiley, New York.
Müller, B. , 1998, “ Low-Mach-Number Asymptotics of the Navier–Stokes Equations,” J. Eng. Math., 34(1), pp. 97–109. [CrossRef]
Principe, J. , and Codina, R. , 2009, “ Mathematical Models for Thermally Coupled Low Speed Flows,” Adv. Theor. Appl. Mech., 2, pp. 93–112.
Felaco, E. , and Gasser, I. , 2013, “ Modelling, Asymptotic Analysis and Simulation of the Gas Dynamics in a Chimney,” J. Math. Ind., 3, pp. 1–20. [CrossRef]
Hirsch, M. W. , Pugh, C. C. , and Shub, M. , 1970, “ Invariant Manifolds,” Bull. Am. Math. Soc., 76(5), pp. 1015–1019. [CrossRef]
Brayton, R. K. , and Moser, J. , 1964, “ A Theory of Nonlinear Networks—I,” Q. Appl. Math., 22(1), pp. 1–33.
Hoshino, H. , Susuki, Y. , and Hikihara, T. , 2014, “ Stability of Two-Sites Regional Power and Heat Supply System—Numerical Simulations Under Step Change of Heat Transfer Reference,” Trans. ISCIE, 27(11), pp. 452–460 (in Japanese). [CrossRef]
Holmgren, M. , 2007, “ X Steam, Thermodynamic properties of water and steam,” Matlab Central, Natick, MA.


Grahic Jump Location
Fig. 1

Overview of steam supply systems with multiple CHP plants that we consider in this paper. It shows an example of three-site system.

Grahic Jump Location
Fig. 2

Schematic diagram of the two-site steam supply. Each block in the figure shows (a) components of the two-site system and (b) components of a boiler.

Grahic Jump Location
Fig. 3

Schematic diagram of the steam transporting pipe l

Grahic Jump Location
Fig. 4

Numerical values of (a) ei(pi) (solid line) and the values of terms on the right-hand side of Eq. (2) and (b) hc(pis(pi)

Grahic Jump Location
Fig. 5

Responses of state variables (p1, p2, u) and heat output rates Qoi′:=msi′hc(pi) from the boilers. The solid lines show the responses of Eq. (14), and the points (×) show the sample plot of responses of Eq. (19). They are initiated by a step changeofthe parameters from (Q1′,Q2′)=(5MJ/s, 5MJ/s ) to. (Q1′,Q2′)=(6MJ/s, 4MJ/s ).

Grahic Jump Location
Fig. 6

Trajectories of the derived model (14) with (Q1′,Q2′)=(6 MJ/s,4 MJ/s). The trajectory starting from the point (p1,p2,u) = (800kPa,800kPa,0m/s) corresponds to the time response presented in Fig. 5 with subsequent long-term response. The trajectory which is parallel to (p1+p2)/2-axis shows an invariant manifold located with direct numerical integration of Eq. (14).

Grahic Jump Location
Fig. 7

Long-term dynamics for the periodic change of the parameter Q′1. Numerical simulation of (a) time responses of the state variables and (b) trajectory in the phase space are shown.

Grahic Jump Location
Fig. 8

Responses of p1, p2, u0:=u(t,0), and uL:=u(t,L) by the original coupled equations. The points (×) show the sequence of time responses of the derived lumped-parameter model (14).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In