Analytical 'steady-state'-based derivation and clarification of the courant-friedrichs-lewy condition for pipe flow

被引:10
作者
Kowalczuk, Zdzislaw [1 ]
Tatara, Marek Sylwester [1 ]
机构
[1] Gdansk Univ Technol, Fac Elect Telecommun & Informat, Dept Robot & Decis Syst, Ul Narutowicza 11-12, PL-80233 Gdansk, Poland
关键词
Computational physics; Optimization; Pipeline flow process; Emulation; Courant number; FLUID-FLOW; STABILITY; MODEL; SCHEMES; 1D;
D O I
10.1016/j.jngse.2021.103953
中图分类号
TE [石油、天然气工业]; TK [能源与动力工程];
学科分类号
0807 ; 0820 ;
摘要
This article addresses the problem of choosing the optimal discretization grid for emulating fluid flow through a pipeline. The aggregated basic flow model is linearized near the operating point obtained from the steady state analytic solution of the differential equations under consideration. Based on this model, the relationship between the Courant number (mu) and the stability margin is examined. The numerically set coefficient mu(opt), ensuring the maximum margin of stability, is analyzed in terms of the physical and technological parameters of the flow. As a result of this analysis, a specific formula is obtained based on parameters describing the mechanics (geometry and physics) of the flow through the pipeline, which leads to the optimal value of the Courant number, separately for smooth and rough pumping conditions. A more detailed analysis of the distribution of the optimal mu coefficient in relation to the parameters of the pipeline flow mechanics shows four cases to consider when determining the coefficient mu(opt). Surprisingly, in three cases, the CFL condition is insufficient, which is expressed in the form of the proposed procedure for choosing the optimal value of mu. The final dichotomous model is derived from the Monte Carlo simulation results in which the effect of each parameter on the optimal Courant number is estimated and consolidated. Taking into account the recognized general laws of physics and using numerical methods and mathematical analysis, simple and useful analytical relationships describing the flow process are obtained. In addition, computer simulations are performed to verify the correctness of the proposed procedure, as well as a number of other considerations related to the modeling of fluid flow in transport pipelines.
引用
收藏
页数:16
相关论文
共 39 条
[1]  
Alghurabi A., 2020, J NAT GAS SCI ENG
[2]   On stability of staggered schemes [J].
Bauer, Amy L. ;
Loubere, Raphaeerl ;
Wendroff, Burton .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2008, 46 (02) :996-1011
[3]   LEAK DETECTION METHODS FOR PIPELINES [J].
BILLMANN, L ;
ISERMANN, R .
AUTOMATICA, 1987, 23 (03) :381-385
[4]   Simplified 1D Incompressible Two-Fluid Model with Artificial Diffusion for Slug Flow Capturing in Horizontal and Nearly Horizontal Pipes [J].
Bonzanini, Arianna ;
Picchi, Davide ;
Poesio, Pietro .
ENERGIES, 2017, 10 (09)
[5]  
Bridson R., 2015, FLUID SIMULATION COM
[6]  
Brogan W.L, 1991, MODERN CONTROL THEOR
[7]   A conservative overlap method for multi-block parallelization of compact finite-volume schemes [J].
Capuano, F. ;
Mastellone, A. ;
De Angelis, E. M. .
COMPUTERS & FLUIDS, 2017, 159 :327-337
[8]   An unconditionally stable 3-D ADI-MRTD method free of the CFL stability condition [J].
Chen, ZZ ;
Zhang, JZ .
IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS, 2001, 11 (08) :349-351
[9]   ON PARTIAL DIFFERENCE EQUATIONS OF MATHEMATICAL PHYSICS [J].
COURANT, R ;
FRIEDRICHS, K ;
LEWY, H .
IBM JOURNAL OF RESEARCH AND DEVELOPMENT, 1967, 11 (02) :215-+
[10]   Numerical method of characteristics for semilinear partial functional differential systems [J].
Czernous, W. .
JOURNAL OF NUMERICAL MATHEMATICS, 2008, 16 (01) :1-21