Performant and Simple Numerical Modeling of District Heating Pipes with Heat Accumulation

This paper compares approaches for accurate numerical modeling of transients in the pipe element of district heating systems. The distribution grid itself affects the heat flow dynamics of a district heating network, which subsequently governs the heat delays and entire efficiency of the distribution. For an efficient control of the network, a control system must be able to predict how “temperature waves” move through the network. This prediction must be sufficiently accurate for real-time computations of operational parameters. Future control systems may also benefit from the accumulation capabilities of pipes. In this article, the key physical phenomena affecting the transients in pipes were identified, and an efficient numerical model of aboveground district heating pipe with heat accumulation was developed. The model used analytical methods for the evaluation of source terms. Physics of heat transfer in the pipe shells was captured by one-dimensional finite element method that is based on the steady-state solution. Simple advection scheme was used for discretization of the fluid region. Method of lines and time integration was used for marching. The complexity of simulated physical phenomena was highly flexible and allowed to trade accuracy for computational time. In comparison with the very finely discretized model, highly comparable transients were obtained even for the thick accumulation wall

Improving the Accuracy and Scope of Control-Oriented Vapor Compression Cycle System Models

The benefits of applying advanced control techniques to vapor compression cycle systems are well know. The main advantages are improved performance and efficiency, the achievement of which brings both economic and environmental gains. One of the most significant hurdles to the practical application of advanced control techniques is the development of a dynamic system level model that is both accurate and mathematically tractable. Previous efforts in control-oriented modeling have produced a class of heat exchanger models known as moving-boundary models. When combined with mass flow device models, these moving-boundary models provide an excellent framework for both dynamic analysis and control design. This thesis contains the results of research carried out to increase both the accuracy and scope of these system level models. The improvements to the existing vapor compression cycle models are carried out through the application of various modeling techniques, some static and some dynamic, some data-based and some physics-based. Semiempirical static modeling techniques are used to increase the accuracy of both heat exchangers and mass flow devices over a wide range of operating conditions. Dynamic modeling techniques are used both to derive new component models that are essential to the simulation of very common vapor compression cycle systems and to improve the accuracy of the existing compressor model. A new heat exchanger model that accounts for the effects of moisture in the air is presented. All of these model improvements and additions are unified to create a simple but accurate system level model with a wide range of application. Extensive model validation results are presented, providing both qualitative and quantitative evaluation of the new models and model improvements.Air Conditioning and Refrigeration Project 17

Thermal Transients in District Heating Systems

Heat fluxes in a district heating pipeline systems need to be controlled on the scale from minutes to an hour to adjust to evolving demand. There are two principal ways to control the heat flux - keep temperature fixed but adjust velocity of the carrier (typically water) or keep the velocity flow steady but then adjust temperature at the heat producing source (heat plant). We study the latter scenario, commonly used for operations in Russia and Nordic countries, and analyze dynamics of the heat front as it propagates through the system. Steady velocity flows in the district heating pipelines are typically turbulent and incompressible. Changes in the heat, on either consumption or production sides, lead to slow transients which last from tens of minutes to hours. We classify relevant physical phenomena in a single pipe, e.g. turbulent spread of the turbulent front. We then explain how to describe dynamics of temperature and heat flux evolution over a network efficiently and illustrate the network solution on a simple example involving one producer and one consumer of heat connected by "hot" and "cold" pipes. We conclude the manuscript motivating future research directions.Comment: 31 pages, 7 figure

Connecting dispersion models and wall temperature prediction for laminar and turbulent flows in channels

In a former paper, Drouin et al. (2010) proposed a model for dispersion phenomena in heated channels that works for both laminar and turbulent regimes. This model, derived according to the double averaging procedure, leads to satisfactory predictions of mean temperature. In order to derive dispersion coefficients, the so called ‘‘closure problem’’ was solved, which gave us access to the temperature deviation at sub filter scale. We now propose to capitalize on this useful information in order to connect dispersion modeling to wall temperature prediction. As a first step, we use the temperature deviation modeling in order to connect wall to mean temperatures within the asymptotic limit of well established pipe flows. Since temperature in wall vicinity is mostly controlled by boundary conditions, it might evolve according to different time and length scales than averaged temperature. Hence, this asymptotic limit provides poor prediction of wall temperature when flow conditions encounter fast transients and stiff heat flux gradients. To overcome this limitation we derive a transport equation for temperature deviation. The resulting two-temperature model is then compared with fine scale simulations used as reference results. Wall temperature predictions are found to be in good agreement for various Prandtl and Reynolds numbers, from laminar to fully turbulent regimes and improvement with respect to classical models is noticeable

Optimal Control of Transient Flow in Natural Gas Networks

We outline a new control system model for the distributed dynamics of compressible gas flow through large-scale pipeline networks with time-varying injections, withdrawals, and control actions of compressors and regulators. The gas dynamics PDE equations over the pipelines, together with boundary conditions at junctions, are reduced using lumped elements to a sparse nonlinear ODE system expressed in vector-matrix form using graph theoretic notation. This system, which we call the reduced network flow (RNF) model, is a consistent discretization of the PDE equations for gas flow. The RNF forms the dynamic constraints for optimal control problems for pipeline systems with known time-varying withdrawals and injections and gas pressure limits throughout the network. The objectives include economic transient compression (ETC) and minimum load shedding (MLS), which involve minimizing compression costs or, if that is infeasible, minimizing the unfulfilled deliveries, respectively. These continuous functional optimization problems are approximated using the Legendre-Gauss-Lobatto (LGL) pseudospectral collocation scheme to yield a family of nonlinear programs, whose solutions approach the optima with finer discretization. Simulation and optimization of time-varying scenarios on an example natural gas transmission network demonstrate the gains in security and efficiency over methods that assume steady-state behavior

Computationally efficient stratified flow wet angle correlation for high resolution simulations

n high resolution two-phase pipe flow simulations, such as slug capturing simulation for liquid-gas pipe flow, explicit calculation of stratified flow wet angle has been proposed to improve computational speed of simulations. Most phenomenological and approximate models for obtaining reliable predictions for stratified flow wet angle employ iterative methods or contain long explicit equations which reduce computational efficiency of these models in high-resolution simulations. Therefore, the aim of this study is to adapt a simple mathematical model for predicting stratified flow wet angle to achieve computationally efficient high-resolution liquid-gas pipe flow simulations