论文标题
非理想气体的稳态网络流动方程的数值解
Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas
论文作者
论文摘要
我们为非理想气体制定了一个稳态网络流问题,该问题将注入率和网络中的淋巴结压力与管道流动相关联。对于这个问题,我们介绍并证明了关于符合单调性特性的广泛类别的非理想压力密度关系的广义解决方案唯一性的定理。此外,我们开发了一种用于稳态问题的数值解的牛顿 - 拉夫森算法,该方程式的系统非二键化使其成为可能。已开发的算法已在基准实例上进行了广泛的测试,并证明可以稳健地收敛到广义溶液。先前的结果表明,由于对初始猜测具有极大的敏感性,因此很难通过牛顿 - 拉夫森方法来解决理想气体的稳态网络流动方程。相比之下,我们发现稳态问题的非限制化是牛顿 - 拉夫森方法稳健收敛的关键。我们根据解决方案的独特性确定标准,在该解决方案的独特性下,非线性求解器发现的非物理广义溶液的存在意味着物理溶液不存在,即问题的不可行。最后,我们比较了基于理想和非理想状态方程的压力和流解决方案,以证明需要在实践中应用后者。本文开发的求解器是开源的,可用于学术和研究社区以及行业。
We formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.