WO2021225173A1 - ガス流量解析方法 - Google Patents

ガス流量解析方法 Download PDF

Info

Publication number
WO2021225173A1
WO2021225173A1 PCT/JP2021/017599 JP2021017599W WO2021225173A1 WO 2021225173 A1 WO2021225173 A1 WO 2021225173A1 JP 2021017599 W JP2021017599 W JP 2021017599W WO 2021225173 A1 WO2021225173 A1 WO 2021225173A1
Authority
WO
WIPO (PCT)
Prior art keywords
gas
tank
flow rate
pressure
tanks
Prior art date
Application number
PCT/JP2021/017599
Other languages
English (en)
French (fr)
Inventor
清 菅井
Original Assignee
Sugai Kiyoshi
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sugai Kiyoshi filed Critical Sugai Kiyoshi
Publication of WO2021225173A1 publication Critical patent/WO2021225173A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]

Definitions

  • the present invention relates to a gas flow rate analysis method for determining the gas pressure of a plurality of tanks and the gas flow rate flowing through the gas pipeline. More specifically, the gas pressure and gas of a plurality of tanks subordinately connected by two or less gas pipelines.
  • the present invention relates to a gas flow analysis method for obtaining a gas flow rate flowing through a pipeline.
  • Non-Patent Document 1 simulates a boiler to which pipes, valves, and pumps are connected in an electric power plant, and modularizes and models a liquid flow rate system into a device having a mathematical formula having the same characteristics. The distribution is calculated.
  • This model has two fixed pressure points P 0 and P 9 and internal 1 to 8 branch pressure points, and the other branch pressure points P 1 , P 4 and P 8 are minute. Find the pressure fluctuation at each branch point when it changes.
  • the flow rate of the incompressible liquid flowing in the pipe connecting the tank 1 and the tank 2 is obtained based on the following Bernoulli's equation.
  • the pressure fluctuation of the entire boiler system is obtained by performing a matrix calculation using a Jacobian matrix and a convergence calculation of both pressure and flow rate by a computer.
  • Non-Patent Document 1 obtains the fluid pressure at each connection point of the modeled system and the liquid flow rate of the pipeline, and the pressure of the gas tank subordinately connected in the gas pipeline of the present application. It does not determine the flow rate flowing through the gas pipeline.
  • Non-Patent Document 3 is a paper explaining the calculation formula of the control valve, and the flow rate Q is expressed by the following formula using the Cv value.
  • the gas flowing in the pipe in the gas flow rate system is in the turbulent flow region, and the state in which friction occurs between the gas and the pipe wall is the same as that of the adjusting valve. Therefore, for example, the gas flowing in the pipe connecting the tank 1 and the tank 2
  • the flow rate can be calculated by applying the following equations (4) according to the equations (2) and (3).
  • Flow rate Q coefficient x SQRT ((P 1 ) 2- (P 2 ) 2 ) (4)
  • SQRT () represents the square root of the numerical value in (). (The same applies below.)
  • the present invention has been made in view of the conventional problems, and an object of the present invention is the gas pressure of a plurality of tanks subordinately connected by two or less gas pipelines and the gas flow rate flowing through the gas pipelines.
  • the present invention made to achieve the above object is a gas flow rate analysis method for obtaining the gas pressure of a plurality of tanks subordinately connected by two or less gas pipelines and the gas flow rate flowing through the gas pipeline by an analyzer. hand, The analyzer obtains the gas pressure P of the tank from the following [Equation 1].
  • the initial value of the gas flow coefficient that flows in the gas pipeline connecting N and M obtained in [Equation 2] is calculated from the initial flow rate between tank N and tank M and the initial value of the gas pressure between tank N and tank M.
  • Initial value setting stage to obtain and assign to variable memory as variable value
  • the analyzer sets a value DXM, which represents a change in the flow rate of gas flowing through each gas pipeline, for each predetermined calculation time ( ⁇ t).
  • DXM current gas flow rate flowing through the gas pipeline-calculated by the gas flow rate before ⁇ t flowing through the gas pipeline
  • Gas amount XM previous gas amount XM + change in gas flow rate
  • DXM x calculation time step ⁇ t To obtain the gas amount XM of the tank connecting each gas pipeline, Including a step of repeating the determination of gas pressures of a plurality of tanks a predetermined number of times according to [Equation 1].
  • Some of the gas pipelines are equipped with actuators, and the analysis device further includes a step of setting the opening degree of the actuators. It is characterized in that the flowing gas flow rate is obtained.
  • the actuator in multiple tanks that are subordinately connected by two or less gas pipelines, the actuator should be set at an arbitrary location in the gas pipeline in response to the operator's input, and accident analysis should be performed when the actuator is fully closed. It is characterized by.
  • the gas pressure and gas of the tank within a predetermined time when the actuator opening is changed by using a simple calculation formula. It is possible to know the fluctuation of the gas flow rate flowing through the pipeline.
  • the actuator can be set at an arbitrary location in the gas pipeline, and an accident analysis can be performed when the actuator opening is fully closed.
  • FIG. 1 It is a block diagram of the tank used for the gas flow rate analysis flow by this invention. This is the first half of the main routine used in the gas flow rate analysis flow according to the present invention. This is the latter half of the main routine used in the gas flow rate analysis flow according to the present invention. This is the first half of the subroutine used in the gas flow rate analysis flow according to the present invention. This is the latter half of the subroutine used in the gas flow rate analysis flow according to the present invention.
  • FIG. 1 is a configuration diagram of a tank for obtaining the gas pressure of a plurality of tanks subordinately connected by two or less gas pipelines according to the present invention and the gas flow rate flowing through the gas pipeline.
  • Fig. 1 there are two gas pipelines C and D from tanks 3 to 4, because the dual system is used as a failure countermeasure for important gas pipelines.
  • nitrogen gas is used as the gas.
  • nitrogen gas is simply referred to as gas.
  • the gas in the tank 1 flows into the tank 3 through the gas pipeline A.
  • the gas in the tank 2 flows into the tank 3 through the gas pipeline B.
  • the gas in the tank 3 flows into the tank 4 through the gas pipeline C and the gas pipeline D, which are dual systems.
  • the flow rate of the gas in the gas pipeline C and the gas pipeline D is controlled by the opening degree of the actuator C and the actuator D.
  • DXM current gas flow rate flowing through the gas pipeline-gas flow rate before ⁇ t flowing through the gas pipeline
  • gas pressure P in the tank can be calculated by the following formula according to [Equation 1].
  • Gas pressure P (gas amount XM x gas-specific constant RC x gas temperature T) / tank volume V
  • the processing flow consists of the main routine shown in FIG. 2 and the subroutine shown in FIG.
  • the main routine is a repetitive calculation, and the calculation is performed every ⁇ t in the calculation time interval and repeated a predetermined number of times.
  • the number of repetitions was set to, for example, 70.
  • step S1 FLAG1 and FLAG2 are set to 1, the number of repeated operations I is set to 0, and the operation time t is set to 0.
  • the initial values of the parameters set in step S3 include the gas pressures P1 to P4 of the tanks 1 to 4, the tank volumes V1 to V4 of the tanks 1 to 4, and the initial flow rate WA1 of the gas pipelines A to D. , WB1, WC1, WD1, each of the gas temperatures T1 to T4 of the tanks 1 to 4, the gas-specific constant RC, the opening degree AKA of each of the actuators C and D, the opening degree AKB, and the calculation time step ⁇ t.
  • the initial flow rates WA1, WB1, WC1, and WD1 of the gas pipelines A to D are determined by referring to the calculation formula (2), the measured values of similar systems in the past, or the measured values of the operating system. And.
  • WAA1 WA1
  • WBB1 WB1
  • WDD WC1 + WD1.
  • the initial values of the respective gas amounts XM1 to XM3 of the tanks 1 to 3 set in step S4 are the parameters of [Equation 1], the gas pressures P1 to P3 set in step S3, the tank volumes V1 to V3, and the gas-specific constants. Obtained using the initial values of RC and gas temperatures T1 to T3. Further, the initial values of the flow coefficients AG and BG of the gas pipelines A and B are obtained by using the initial values of the initial flow rates WA1 and WB1 and the gas pressures P1 to P3 set in step S3 in the parameter of [Equation 2]. These initial values are recorded in a predetermined memory for data processing after the calculation.
  • Flow coefficient AG Pipeline A flow rate (WA1) / SQRT (pressure (P1) 2 -pressure (P3) 2 )
  • Flow coefficient BG Pipeline B flow rate (WB1) / SQRT (pressure (P2) 2 -pressure (P3) 2 )
  • WAA1 AG (flow coefficient) x SQRT ((pressure P1) 2- (pressure P3) 2 )
  • WBB1 BG (flow coefficient) x SQRT ((pressure P2) 2- (pressure P3) 2 )
  • step S7 the gas outflow amounts WAA1 and WBB1 obtained in step S6 are used to obtain the gas outflow change amounts DXM1 and DXM2 of the tanks 1 and 2.
  • Gas outflow change DXM1 current gas outflow (WAA1) -previous gas outflow (WAA1)
  • Gas outflow change DXM2 current gas outflow (WBB1) -previous gas outflow (WBB1)
  • WAA1 current gas outflow
  • WBB1 current gas outflow
  • WBB1 current gas outflow
  • the previous gas outflow amount (WAA1) and the previous gas outflow amount (WBB1) are read from the memory area stored under the same name at the time of the previous calculation and used. At the time of the first calculation, the value stored as the initial value is read out and used.
  • the current gas outflow amount (WAA1) and the current gas outflow amount (WBB1) are stored in the previous memory outflow amount (WAA1) and the previous gas outflow amount (WBB1), which are the above memory areas.
  • the amount of change in gas outflow is the amount of change in gas amount in the tank.
  • step S8 the subroutine is executed in step S8.
  • the subroutine shown in FIG. 3 will be described later.
  • step S10 the gas amounts XM1 to XM3 of the tanks 1 to 3 are obtained for each calculation time step ⁇ t by using the gas amount changes DXM1 to DXM3 in the tanks 1 to 3 obtained in steps S7 and S9.
  • Gas amount of tank 1 XM1 previous gas amount XM1 + DXM1 ⁇ calculation time step ⁇ t
  • Gas amount of tank 2 XM2 previous gas amount XM2 + DXM2 ⁇ calculation time step ⁇ t
  • Gas amount of tank 3 XM3 previous gas amount XM3 + DXM3 ⁇ calculation time step ⁇ t
  • step S11 the gas pressures P1 to P3 of each tank are obtained from the result obtained in step S10 by using [Equation 1].
  • Gas pressure P1 (XM1 x gas-specific RC x gas temperature T1) / tank volume V1
  • Gas pressure P2 (XM2 x gas-specific RC x gas temperature T2) / tank volume V2
  • Gas pressure P3 (XM3 x gas-specific RC x gas temperature T3) / tank volume V3
  • step S12 the gas pressures P1 to P3 obtained by the above formula are processed according to the conditions.
  • P3 ⁇ P1, P3 P1
  • P3 P2
  • step S13 the gas outflow amount WAA1, WBB1, the sum WCC of the gas outflow amount WAA1 and WBB1, the gas outflow change amount DXM1 to DXM3, the gas amounts XM1 to XM3 of each of the tanks 1 to 3, and the gas obtained by the above calculation.
  • the pressures P1 to P3 are stored in a predetermined memory area in association with the calculation time t for data processing such as graph display after the calculation.
  • This subroutine is to obtain the sum WDD of the inflow amount of the gas flowing from the tank 3 to the tank 4 in the gas pipelines C and D in which the actuators C and D are installed, and the pressure of the tank 4.
  • the initial values of the gas pressure P4, the tank volume V4, the gas-specific constant RC, and the gas temperature T4 set in step S3 of the main routine are used. 1] is used to obtain. Then, the initial values of the flow coefficient CG and DG are set using the initial values of the initial flow rates WC1, WD1 and the gas pressures P3 and P4 that flow from the tank 3 to the tank 4 via the gas pipeline C and the gas pipeline D [Equation 2]. ] To obtain. These initial values are recorded in a predetermined memory for data processing after the calculation.
  • Gas amount of tank 4 XM4 pressure P4 x tank volume V4 / (constant RC x gas temperature T4)
  • Flow coefficient CG WC1 / (0.8 ⁇ SQRT ((pressure P3) 2- (pressure P4) 2 ))
  • Flow coefficient DG WD1 / (0.8 ⁇ SQRT ((pressure P3) 2- (pressure P4) 2 ))
  • step S17 the number of calculations I is updated.
  • step S18 If the preset number of calculations is 3 in step S18, the process proceeds to step S19 to change the opening degrees (opening AKA, opening AKB) of the actuators C and D to the changed values of 0.795. ..
  • step S20 when the number of calculations is other than 3, the opening degrees (opening AKA, opening AKB) of the actuators C and D are both set to normal values of 0.8.
  • the number of calculations for changing the opening degree (opening AKA, opening degree AKB) of the actuators C and D and the actuator opening degree can be set to any number.
  • step S21 based on the flow coefficient CG and DG obtained by the initial value calculation and the set opening degrees (opening C and opening D) of the actuators C and D, from the tank 3 according to [Equation 2].
  • the gas inflow amounts WCC1 and WDD1 flowing into the tank 4 are obtained.
  • Gas inflow WCC1 CG x AKA x SQRT ((pressure P3) 2- (pressure P4) 2 )
  • Gas inflow WDD1 DG x AKB x SQRT ((pressure P3) 2- (pressure P4) 2 )
  • step S22 the sum WDD of the gas inflow amount flowing from the tank 3 to the tank 4 is obtained by the gas inflow amount WCC1 and WDD1 obtained in step S21.
  • Gas inflow WDD WCC1 + WDD1 WDD is returned to the main routine as a return value.
  • step S23 the gas inflow change amount DXM4 in the tank 4 is obtained based on the gas inflow amount WDD obtained in step S22 and the previous gas inflow amount WDD.
  • Gas inflow change amount DXM4 current gas inflow amount WDD-previous gas inflow amount WDD
  • the previous gas inflow amount WDD is read from the memory area stored under the same name at the time of the previous calculation and used. At the time of the first calculation, the value stored as the initial value is read out and used. Then, when the calculation here is completed, the gas inflow amount WDD this time is stored in the previous gas inflow amount WDD, which is the memory area.
  • the amount of change in gas inflow is the amount of change in gas amount.
  • step S24 the gas amount XM4 of the tank 4 is obtained based on the gas inflow change amount DXM4 obtained in step S23.
  • Gas amount of tank 4 XM4 previous gas amount XM4 + DXM4 ⁇ calculation time step ⁇ t
  • step S25 the gas pressure P4 of the tank 4 can be obtained from the gas amount XM4 obtained in step S24 using [Equation 1].
  • Gas pressure P4 (XM4 x gas-specific constant RC x gas temperature T4) / tank volume V4
  • step S26 the gas pressures P1 to P3 obtained by the above formula are processed according to the conditions.
  • P4 ⁇ 0 P4 0
  • P4 ⁇ P3 P4 P3
  • step S27 the gas inflow amount WCC1, the gas inflow amount WDD1, the gas inflow amount WDD, the gas amount XM4 of the tank 4, and the gas pressure P4 obtained by the above calculation are for data processing such as graph display after the calculation. , It is stored in a predetermined memory area in relation to the calculation time t.
  • the data obtained from the simulation can be used as design data when determining the size of the tank volume, pump capacity, pipe diameter, etc. of a system having multiple tanks that are subordinately connected by two or less gas pipelines. .. It can also be used to determine the ratings of flowmeters, pressure gauges, etc. installed in gas pipelines.
  • the present invention is not limited to the above embodiment, and the gas of two or less systems is not limited to the technical scope of the present invention. It can be implemented in various ways for multiple tanks that are subordinately connected in a pipeline.
  • the present invention can be applied to the field of gas flow rate analysis for obtaining the gas pressure of a plurality of tanks subordinately connected by two or less gas pipelines and the gas flow rate flowing through the gas pipeline.

Abstract

２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるガス流量解析方法を提供する。　本発明の複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるガス流量解析方法は、タンクのガス圧力Ｐを、　ＰＶ＝ＸＭ×ＲＣ×Ｔ　（Ｐ；圧力、Ｖ；タンク容積、ＸＭ；ガス量、ＲＣ；ガス固有の定数、Ｔ；ガス温度） から求め、　タンク１，２間を結ぶガス管路に流れる流量を、　流量＝流量係数×ＳＱＲＴ（（タンク１のガス圧力）２－（タンク２のガス圧力）２） から求める段階と、　予め定められた演算時間（Δｔ）ごとに、ガス管路に流れるガス流量の変化から、前記ガス管路を結ぶタンクのガス量ＸＭを求め、前記複数タンクのガス圧力を求める段階と、ガス管路のいくつかがアクチュエータを備え、アクチュエータの開度を設定する段階を含むことを特徴とする。

Description

ガス流量解析方法
　本発明は、複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるガス流量解析方法に関し、より詳しくは、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるガス流量解析方法に関する。
　一般に、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めることは、流体方程式に関連し計算機を用いたとしても容易でない。
　非特許文献１は、電力プラントにおいて、配管、弁、ポンプが接続されたボイラのシミュレーションを行うものであり、同じ特性の数式を持つ機器にモジュール化し、モデル化された液体流量系の圧力と流量分布を求めるものである。
　このモデルは、Ｐ０及びＰ９の２つの固定圧力点と内部の１から８の分岐圧力点を有していて、それ以外の分岐圧力点であるＰ１、Ｐ４およびＰ８が微小に変化した場合の各分岐点の圧力変動を求める。
　非特許文献１の圧力変動の計算は、例えばタンク１とタンク２を結ぶ配管の中を流れる非圧縮性液体に対する流量を下記ベルヌーイの式を基に求める。
　　　流量Ｑ　∝　ＳＱＲＴ（Ｐ１－Ｐ２）　（１）
　
　ボイラ系全体の圧力変動は、計算機によりヤコビアン行列を用いた行列演算と圧力と流量双方の収束演算を行って求める。
　しかし、非特許文献１は、モデル化されたシステムの各接続点における流体圧力と管路の液体流量を求めるものであって、本願が目的とするガス管路で従属接続されるガスタンクの圧力とガス管路に流れる流量を求めるものではない。
　近年は、流体の解析ツールとして各種のシミュレーションソフトが世に販売されていてガスを扱うこともできるが、高価でかつ実際の運用場所に合わせて使用するには、ガスを含む流体全般に関する幅広い知識とシミュレーションソフトをこなす専門技術を必要とする。
　しかるに、ガス流量を扱うに当たり参考になる式として、プロセス設備の配管の途中に挿入される調整弁の容量（Ｃｖ値）を求める計算式がある。
　調整弁は液体、気体などの流体を通過させたり遮断したりする目的で使用される。
　非特許文献２は、気体のＣｖ値を求める式として、手計算での確認が行える下記の米国ＦＣＩ（Ｆｌｕｉｄ　Ｃｏｎｔｒｏｌｓ　Ｉｎｓｔｉｔｕｔｅ　Ｉｎｃ．）が定めた式（２）を挙げている。
　（別途、計算式には計算精度が高いＩＥＣ（国際電気標準会議）規格があるが、この式の計算はコンピュータを用いて行う必要がある。）
Figure JPOXMLDOC01-appb-I000001
　　　但し　
　　　　　　Ｑ　：　最大流量（ｍ３／ｈｒ）　（１５．６℃，１０１．３ｋＰａ・Ａ）
　　　　　　Ｇ　：　比重　（空気＝１）
　　　　　Ｔｆ　：　ガス温度（℃）
　　　　　Ｐ１　：　１次側絶対圧力（ｋｇｆ／ｃｍ２　ａｂｓ）
　　　　　Ｐ２　：　２次側絶対圧力（ｋｇｆ／ｃｍ２　ａｂｓ）
　　　　　ΔＰ　：　Ｐ１－Ｐ２（ｋｇｆ／ｃｍ２）
　
　一方、非特許文献３は、調節弁の計算式について説明した論文であり、流量ＱをＣｖ値を用い下式で表している。
Figure JPOXMLDOC01-appb-I000002
　上記論文では、管内を流れる弁入口および弁出口の間の気体に対して、エネルギー保存則を適用し、流れが乱流領域にあり、弁通過時に気体と管壁との間で摩擦が生ずるとして、摩擦を伴う流れに対する一般式を適用し、（３）式を求めている。
　（２）式と（３）式が、係数２８７と２９０の微小な差異を除き一致することは明らかである。
　（２）式で、条件「ΔＰ　＜　Ｐ１／２　　の場合」を付しているのは、気体（ガス）においては、ΔＰ（＝Ｐ１－Ｐ２）＞　Ｐ１／２　の条件では弁の縮流部の流速が音速に達して流量が飽和し、それ以上いくら差圧を大きく与えても流量が増加しないことを考慮してなされたためであり、調整弁に対してのみ適用される。
　ガス流量システムにおける管内を流れるガスは乱流領域にあり、ガスと管壁との間で摩擦が生ずる状態は調整弁と同一であるので、例えばタンク１とタンク２を結ぶ配管の中を流れるガス流量の計算は（２）、（３）式に依り、下記（４）式を適用して行える。
　　　流量Ｑ＝係数×ＳＱＲＴ（（Ｐ１）２－（Ｐ２）２）　（４）
　
　但し、ＳＱＲＴ（　）は（　）内の数値の平方根を表す。（以下同様。）
　なお、タンク内のガスはＰＶ＝ＮＲＴの公式を満たす。（特許文献１）
特開平３－２０９０９７号公報
プラントシミュレーション　長崎大学大学院生産科学研究科　松井信生　ＷＷＷ．ｅｅｅ．ｎａｇａｓａｋｉ－ｕ．ａｃ．ｊｐ＞ＰＤＦ＞ＳｉｍｕｌａｔｉｏｎＭａｔｓｕｉ１ アズビル株式会社　製品仕様書のダウンロード　コントロールバルブ（調整弁）／操作端　ポートサイズとＣｖ値計算式　ＳＰ１－ＣＶ２０１４－Ｓｅｃ１０－１　１．３　気体のＣｖ値計算式 調整弁の計算式について　浅野　勇　昭和３６年９月　計測第９号　５５３～５６０頁
　本発明は、従来の問題点に鑑みてなされたものであって、本発明の目的は、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を、式（４）を用い簡易な数値計算で求めるガス流量解析方法を提供することにある。
　上記目的を達成するためになされた本発明は、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を解析装置によって求めるガス流量解析方法であって、
　前記解析装置が、タンクのガス圧力Ｐを、下記［数１］から求め、
［数１］
　ＰＶ＝ＸＭ×ＲＣ×Ｔ
　（ここで　Ｐ；ガス圧力、Ｖ；タンク容積、ＸＭ；ガス量、ＲＣ；ガス固有の定数、Ｔ；ガス温度）
　タンクＮ，Ｍ間を結ぶガス管路に流れるガス流量を、下記［数２］から求める段階と、
［数２］
　流量＝流量係数×ＳＱＲＴ（（タンクＮのガス圧力）２－（タンクＭのガス圧力）２）
 
［数１］で求められる前記複数タンクの中の各タンクのガス量ＸＭの初期値を、各タンクのガス圧力、タンク容積、ガス固有の定数及びガス温度の初期値から求め、
［数２］で求められるＮとＭを結ぶガス管路に流れるガス流量係数の初期値を、タンクＮとタンクＭの間の初期流量、及び、タンクＮとタンクＭのガス圧力の初期値から求めて、変数値として変数メモリに代入する初期値設定段階と、
　前記解析装置が、予め定められた演算時間（Δｔ）ごとに、各ガス管路に流れるガス流量の変化を表す値ＤＸＭを、
　ＤＸＭ＝ガス管路に流れる今回のガス流量－ガス管路に流れるΔｔ前のガス流量
により求め、
　ガス量ＸＭ＝前回のガス量ＸＭ＋ガス流量の変化ＤＸＭ×演算時間刻みΔｔ
により、各ガス管路を結ぶタンクのガス量ＸＭを求め、
［数１］に従い、複数タンクのガス圧力を求めることを所定回数繰り返す段階と、を含み、
　ガス管路のいくつかがアクチュエータを備え、前記解析装置がアクチュエータの開度を設定する段階を更に含む、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めることを特徴とする。
　また、２系統以下のガス管路で従属接続される複数タンクにおいて、オペレータの入力に対応してアクチュエータをガス管路の任意の場所に設定し、アクチュエータを全閉した場合の事故解析を行うことを特徴とする。
　本発明によれば、２系統以下のガス管路で従属接続される複数タンクにおいて、簡易な計算式を用いて、アクチュエータ開度を変更した場合の、所定時間内のタンクのガス圧力及び、ガス管路に流れるガス流量の変動を知ることができる。
　また、アクチュエータをガス管路の任意の場所に設定し、アクチュエータ開度を全閉した場合の事故解析を行うことができる。
本発明によるガス流量解析フローに用いられるタンクの構成図である。 本発明によるガス流量解析フローに用いられるメインルーチンの前半部である。 本発明によるガス流量解析フローに用いられるメインルーチンの後半部である。 本発明によるガス流量解析フローに用いられるサブルーチンの前半部である。 本発明によるガス流量解析フローに用いられるサブルーチンの後半部である。
　以下、本発明の実施形態について、ガス管路の所定の箇所で、所定の時間にアクチュエータ開度を変更した場合に、各タンクのガス圧力及び、ガス管路に流れるガス流量がどのように変動するかを調べる解析方法について説明する。
　図１は、本発明による２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるタンクの構成図である。
　　図１で、タンク３から４へのガス管路がＣ、Ｄの２系統あるのは、重要なガス管路に対しては故障対策として２重系が用いられるからである
　ガスの種類として、窒素ガス、アルゴンガス等のガスを用いる構成を形成することができる。
　本発明の実施例は、ガスとして窒素ガスを用いる。以下窒素ガスを単にガスと呼ぶ。
　図１によれば、タンク１にあるガスがガス管路Ａを介してタンク３に流入する。同様に、タンク２にあるガスがガス管路Ｂを介してタンク３に流入する。更に、タンク３のガスは２重系であるガス管路Ｃ、ガス管路Ｄを介してタンク４に流入する。
　本発明では、ガス管路Ｃ，ガス管路ＤのガスをアクチュエータＣ、アクチュエータＤの開度により流量制御する。
　解析対象となるガスがタンク間を流れる場合の、各タンクにおけるガス圧力及び、ガス管路に流れるガス流量を、演算時間Δｔごとに求める。
　演算時間を累積していくことにより、所定時間内の各タンクのガス圧力及び、ガス管路に流れるガス流量の変動を解析することができる。
　各タンクにおけるガス圧力及び、ガス管路に流れるガス流量を求めるに当たり、代表的な計算式を始めに説明する。
　一般にガスは［数１］を満たす。
　［数１］
　ＰＶ＝ＸＭ×ＲＣ×Ｔ
　ここで、
　Ｐ；ガス圧力、Ｖ；タンク容積、ＸＭ；ガス量、ＲＣ；ガス固有の定数、Ｔ；ガス温度である。
　一方、タンク間に流れるガス流量を、例えば、タンク１と３の間のガス管路Ａに流れる流量ＷＡ（例えば毎秒当たり）を、前述した式（４）に倣い下記［数２］で求める。
　［数２］
　流量ＷＡ＝ＡＧ（流量係数）×ＳＱＲＴ（（タンク１の圧力）２－（タンク３の圧力）２）
　演算時間刻みΔｔごとのガス流量の変化を表す値であるＤＸＭ（例えば毎秒当たり）は，下記［数３］で求められる。
　［数３］
　ＤＸＭ＝ガス管路に流れる今回のガス流量－ガス管路に流れるΔｔ前のガス流量
　上記ＤＸＭは、タンク間を結ぶガス流量の変化なので、タンク内のガス量ＸＭの変化をもたらす。
　従って、演算時間刻みΔｔごとのタンク内のガス量ＸＭは下記［数４］で求められる。
　［数４］
　ガス量ＸＭ＝前回のガス量ＸＭ＋ガス流量の変化ＤＸＭ×演算時間刻みΔｔ
　タンク内のガス量ＸＭが分れば、タンク内のガス圧力Ｐは、［数１］に従い、下記の式で求められる。
　ガス圧力Ｐ＝（ガス量ＸＭ×ガス固有の定数ＲＣ×ガス温度Ｔ）／タンク容積Ｖ
　以下、ガス流量解析方法についてフローチャートを基に説明する。
処理フローは図２に示すメインルーチンと図３に示すサブルーチンからなる。
　始めに図２に示すメインルーチンについて説明する。ステップごとに符号Ｓ１～Ｓ１３を付けた。以下すべてステップＳ１～ステップＳ１３とする。
　メインルーチンは繰り返し演算になっており、演算時間刻みΔｔごとに演算を行い所定回数繰り返す。実施例では、繰り返し回数を例えば７０とした。
　ステップＳ１では、ＦＬＡＧ１、ＦＬＡＧ２を１に、繰り返し演算回数　Ｉを０に、演算時間ｔを０に設定する。
　ステップＳ２の分岐では、ＦＬＡＧ１＝１であるので所定のパラメータの初期値を設定する初期値設定のステップＳ３に進む。
　ステップＳ３で設定するパラメータの初期値には、タンク１～４のそれぞれのガス圧力Ｐ１～Ｐ４、タンク１～４のそれぞれのタンク容積Ｖ１～Ｖ４、ガス管路Ａ～Ｄのそれぞれの初期流量ＷＡ１、ＷＢ１、ＷＣ１、ＷＤ１、タンク１～４のそれぞれのガス温度Ｔ１～Ｔ４、ガス固有の定数ＲＣ、アクチェータＣ，Ｄそれぞれの開度ＡＫＡ、開度ＡＫＢ、演算時間刻みΔｔがある。
　ここで、ガス管路Ａ～Ｄそれぞれの初期流量ＷＡ１、ＷＢ１、ＷＣ１、ＷＤ１は、計算式（２）、過去の類似システムの実測値、或いは稼働中のシステムの実測値を参考にして定めるものとする。
　また、後の演算で用いられる初回演算時の前回ガス流出入量は、前回ガス流出量（ＷＡＡ１）＝ＷＡ１、前回ガス流出量（ＷＢＢ１）＝ＷＢ１、前回ガス流入量ＷＤＤ＝ＷＣ１＋ＷＤ１として、所定のメモリー領域に記憶する。
　ステップＳ４で設定するタンク１～３のそれぞれのガス量ＸＭ１～ＸＭ３の初期値は［数１］のパラメータに、ステップＳ３で設定したガス圧力Ｐ１～Ｐ３、タンク容積Ｖ１～Ｖ３、ガス固有の定数ＲＣ，ガス温度Ｔ１～Ｔ３の各初期値を用いて求める。
　また、ガス管路Ａ、Ｂの流量係数ＡＧ、ＢＧの初期値は、［数２］のパラメータにステップＳ３で設定した初期流量ＷＡ１，ＷＢ１及びガス圧力Ｐ１～Ｐ３の初期値を用いて求める。これらの初期値は、演算後のデータ処理のため所定のメモリに記録される。
　
　タンク１のガス量ＸＭ１＝圧力Ｐ１×タンク容積Ｖ１／（定数ＲＣ×ガス温度Ｔ１）
　タンク２のガス量ＸＭ２＝圧力Ｐ２×タンク容積Ｖ２／（定数ＲＣ×ガス温度Ｔ２）
　タンク３のガス量ＸＭ３＝圧力Ｐ３×タンク容積Ｖ３／（定数ＲＣ×ガス温度Ｔ３）
　
　流量係数ＡＧ＝管路Ａ流量（ＷＡ１）／ＳＱＲＴ（圧力（Ｐ１）２－圧力（Ｐ３）２）
　流量係数ＢＧ＝管路Ｂ流量（ＷＢ１）／ＳＱＲＴ（圧力（Ｐ２）２－圧力（Ｐ３）２）
　
　初期値の設定が終了すると、ステップＳ５でＦＬＡＧ１＝０とする。
　ステップＳ６では、演算時間ｔ＝ｔ＋演算時間刻みΔｔ
とした後、ステップＳ４の初期値計算で求めた流量係数ＡＧ、ＢＧを用い、タンク１から３及び、タンク２から３へのガス管路Ａ，Ｂのガス流出量ＷＡＡ１，ＷＢＢ１（例えば毎秒当たり）を［数２］を用いて求める。
　
　ＷＡＡ１＝ＡＧ（流量係数）×ＳＱＲＴ（（圧力Ｐ１）２－（圧力Ｐ３）２）
　ＷＢＢ１＝ＢＧ（流量係数）×ＳＱＲＴ（（圧力Ｐ２）２－（圧力Ｐ３）２）
　
　ステップＳ７では、ステップＳ６で得られたガス流出量ＷＡＡ１、ＷＢＢ１を用い、タンク１，２のガス流出変化量ＤＸＭ１、ＤＸＭ２を求める。
　
　ガス流出変化量ＤＸＭ１＝今回ガス流出量（ＷＡＡ１）－前回ガス流出量（ＷＡＡ１）
　ガス流出変化量ＤＸＭ２＝今回ガス流出量（ＷＢＢ１）－前回ガス流出量（ＷＢＢ１）
　
　ここで、前回ガス流出量（ＷＡＡ１）、前回ガス流出量（ＷＢＢ１）は、前回演算時に同名称で記憶されたメモリ領域から読み出して用いる。初回演算時は、初期値として記憶された値を読みだして用いる。
　そして、ここでの演算が終わると、今回ガス流出量（ＷＡＡ１）、今回ガス流出量（ＷＢＢ１）は、上記メモリ領域である前回ガス流出量（ＷＡＡ１）、前回ガス流出量（ＷＢＢ１）に記憶される。
　
　ガス流出変化量はとりも直さずタンクのガス量変化量である。
　
　続いて、ステップＳ８でサブルーチンを実行する。図３で示すサブルーチンについては後述する。
　ステップＳ９では、タンク３のガス量の変化量ＤＸＭ３を、タンク１、２からタンク３へ流出するガス流出量ＷＡＡ１とＷＢＢ１の和ＷＣＣから、サブルーチンの戻り値でありタンク３からタンク４へ流入するガス管路Ｃ，Ｄのガス流入量の和ＷＤＤを引いて求める。
　ガス量の変化量ＤＸＭ３＝ＷＣＣ－ＷＤＤ
　　　　　　　　ＷＣＣ＝ＷＡＡ１＋ＷＢＢ１
　
　ステップＳ１０では、ステップＳ７，ステップＳ９で求められたタンク１～３内のガス量変化量ＤＸＭ１～ＤＸＭ３を用い、演算時間刻みΔｔごとのタンク１～３のガス量ＸＭ１～ＸＭ３を求める。
　
　タンク１のガス量ＸＭ１＝前回ガス量ＸＭ１＋ＤＸＭ１×演算時間刻みΔｔ
　タンク２のガス量ＸＭ２＝前回ガス量ＸＭ２＋ＤＸＭ２×演算時間刻みΔｔ
　タンク３のガス量ＸＭ３＝前回ガス量ＸＭ３＋ＤＸＭ３×演算時間刻みΔｔ
　
　ステップＳ１１において、ステップＳ１０で得られた結果から［数１］を用いて、各タンクのガス圧力Ｐ１～Ｐ３を求める。
　
　ガス圧力Ｐ１＝（ＸＭ１×ガス固有のＲＣ×ガス温度Ｔ１）／タンク容積Ｖ１
　ガス圧力Ｐ２＝（ＸＭ２×ガス固有のＲＣ×ガス温度Ｔ２）／タンク容積Ｖ２
　ガス圧力Ｐ３＝（ＸＭ３×ガス固有のＲＣ×ガス温度Ｔ３）／タンク容積Ｖ３
　
　ステップＳ１２では、上記計算式で得られたガス圧力Ｐ１～Ｐ３に対し、条件に応じた処理を行う。
　
　Ｐ１～Ｐ３が≦０の場合、Ｐ１～Ｐ３＝０　　　　　　　　　
　Ｐ３≧Ｐ１の場合、Ｐ３＝Ｐ１
　Ｐ３≧Ｐ２の場合、Ｐ３＝Ｐ２
　
　ステップＳ１３で、以上の演算で得られた、ガス流出量ＷＡＡ１，ＷＢＢ１、ガス流出量ＷＡＡ１とＷＢＢ１の和ＷＣＣ、ガス流出変化量ＤＸＭ１～ＤＸＭ３、タンク１～３それぞれのガス量ＸＭ１～ＸＭ３、ガス圧力Ｐ１～３は、演算後のグラフ表示等のデータ処理のため、演算時間ｔに関連づけられて所定のメモリ領域に記憶される。
　最後に演算が終了したか否かの判断を行い、Ｉが７０に達していなければ元のステップＳ２へ戻る。
　次に、サブルーチンについて図３のフローチャートを元に説明する。
このサブルーチンは、アクチュエータＣ、Ｄが設置されているガス管路Ｃ、Ｄの、タンク３からタンク４へ流入するガスの流入量の和ＷＤＤ、及び、タンク４の圧力を求めることにある。
　最初の分岐ステップＳ１４では、メインルーチンのステップＳ１において設定したＦＬＡＧ２＝１なので次のステップＳ１５に進む。
　ステップＳ１５で設定するタンク４のガス量ＸＭ４の初期値は、メインルーチンのステップＳ３で設定したガス圧力Ｐ４、タンク容積Ｖ４、ガス固有の定数ＲＣ，ガス温度Ｔ４の各初期値を用い、［数１］を用いて求める。
　ついで流量係数ＣＧ、ＤＧの初期値を、ガス管路Ｃ及びガス管路Ｄを介しタンク３からタンク４に流入する初期流量ＷＣ１、ＷＤ１及びガス圧力Ｐ３，Ｐ４の初期値を用いて［数２］から求める。これらの初期値は、演算後のデータ処理のため所定のメモリに記録される。
　
　タンク４のガス量ＸＭ４＝圧力Ｐ４×タンク容積Ｖ４／（定数ＲＣ×ガス温度Ｔ４）
　
　流量係数ＣＧ＝ＷＣ１／（０．８×ＳＱＲＴ（（圧力Ｐ３）２－（圧力Ｐ４）２））
　流量係数ＤＧ＝ＷＤ１／（０．８×ＳＱＲＴ（（圧力Ｐ３）２－（圧力Ｐ４）２））
　
　ステップＳ１６で、初期値計算を終え、ＦＬＡＧ２＝０とする。
　ステップＳ１７で、計算回数　Ｉを更新する。
　ステップＳ１８で、あらかじめ設定された計算回数が３である場合には、ステップＳ１９に進みアクチュエータＣ、Ｄの開度（開度ＡＫＡ，開度ＡＫＢ）をいずれも変更値の０．７９５に変更する。
　ステップＳ２０で、計算回数が３以外の場合にはアクチュエータＣ、Ｄの開度（開度ＡＫＡ，開度ＡＫＢ）をいずれも通常値０．８に設定する。
　なお、アクチュエータＣ、Ｄの開度（開度ＡＫＡ，開度ＡＫＢ）を変更する計算回数およびアクチュエータ開度は任意の数に設定可能である。
　ステップＳ２１で、初期値計算で得られた流量係数ＣＧ、ＤＧ及び、設定されたアクチュエータＣ、Ｄの開度（開度Ｃ，開度Ｄ）を元に、［数２］に従い、タンク３からタンク４へ流入するガス流入量ＷＣＣ１、ＷＤＤ１を求める。
　
　ガス流入量ＷＣＣ１＝ＣＧ×ＡＫＡ×ＳＱＲＴ（（圧力Ｐ３）２－（圧力Ｐ４）２）
　ガス流入量ＷＤＤ１＝ＤＧ×ＡＫＢ×ＳＱＲＴ（（圧力Ｐ３）２－（圧力Ｐ４）２）
　
　ステップＳ２２で、ステップＳ２１で求めたガス流入量ＷＣＣ１、ＷＤＤ１によりタンク３からタンク４へ流入するガス流入量の和ＷＤＤを求める。
　
　ガス流入量ＷＤＤ＝ＷＣＣ１＋ＷＤＤ１
　
　なお、ＷＤＤは戻り値としてメインルーチンへ返される。
　
　ステップＳ２３で、ステップＳ２２で得られたガス流入量ＷＤＤと前回のガス流入量ＷＤＤを元に、タンク４内のガス流入変化量ＤＸＭ４を求める。
　
　ガス流入変化量ＤＸＭ４＝今回ガス流入量ＷＤＤ－前回ガス流入量ＷＤＤ
　
　ここで、前回ガス流入量ＷＤＤは、前回演算時に同名称で記憶されたメモリ領域から読み出して用いる。初回演算時は、初期値として記憶された値を読みだして用いる。
　そして、ここでの演算が終わると、今回ガス流入量ＷＤＤは、上記メモリ領域である前回ガス流入量ＷＤＤに記憶される。
　
　ガス流入変化量はとりも直さずガス量変化量である。
　
　ステップＳ２４で、ステップＳ２３で得られたガス流入変化量ＤＸＭ４を元に、タンク４のガス量ＸＭ４を求める。
　
　タンク４のガス量ＸＭ４＝前回ガス量ＸＭ４＋ＤＸＭ４×演算時間刻みΔｔ
　
　ステップＳ２５で、ステップＳ２４で得られたガス量ＸＭ４から、［数１］を用いてタンク４のガス圧力Ｐ４が求まる。
　
　ガス圧力Ｐ４＝（ＸＭ４×ガス固有の定数ＲＣ×ガス温度Ｔ４）／タンク容積Ｖ４
　
　ステップＳ２６で、上記計算式で得られたガス圧力Ｐ１～Ｐ３に対し、条件に応じた処理を行う。　
　
　Ｐ４≦０の場合　　Ｐ４＝０
　Ｐ４≧Ｐ３の場合　Ｐ４＝Ｐ３
　
　ステップＳ２７で、以上の演算で得られた、ガス流入量ＷＣＣ１、ガス流入量ＷＤＤ１、ガス流入量ＷＤＤ、タンク４のガス量ＸＭ４、ガス圧力Ｐ４は、演算後のグラフ表示等のデータ処理のため、演算時間ｔに関連ずけられ所定のメモリ領域に記憶される。
　上記解析方法は、２系統以下のガス管路で従属接続される複数タンクにおいて、ガス管路の所定の箇所に設定されたアクチュエータの開度を変更した場合に、各タンクのガス圧力及び、ガス管路に流れるガス流量がどのように変動するかを調べたものである。
　シミュレーションから得られたデータは、２系統以下のガス管路で従属接続される複数タンクを有するシステムのタンク容積、ポンプ容量、配管径等の大きさを決める場合の設計データとして役立たせることができる。
また、ガス管路に装備する流量計、圧力計等の定格を定めるために用いることができる。
　また、本解析方法を用いれば、アクチュエータをガス管路の任意の場所に設定し、設定された管路を全閉した場合の２系統以下のガス管路で従属接続される複数タンクを有するシステムの事故解析を行うことが出来る。
解析で得られるデータは、事故対策に役立たせることができる。
　以上、本発明の実施形態について図面を参照しながら説明したが、本発明は、上記の実施形態に限られるものではなく、本発明の技術的範囲から逸脱しない範囲内で、２系統以下のガス管路で従属接続される複数タンクについて多様に実施することができる。
　本発明は、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めるガス流量解析の分野に適用できる。

Claims (2)

  1. 　２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を解析装置によって求めるガス流量解析方法であって、
    　前記解析装置が、タンクのガス圧力Ｐを、下記［数１］から求め、
    ［数１］
    　ＰＶ＝ＸＭ×ＲＣ×Ｔ
    　（ここで　Ｐ；ガス圧力、Ｖ；タンク容積、ＸＭ；ガス量、ＲＣ；ガス固有の定数、Ｔ；ガス温度）
     
    　タンクＮ，Ｍ間を結ぶガス管路に流れるガス流量を、下記［数２］から求める段階と、
    ［数２］
    　流量＝流量係数×ＳＱＲＴ（（タンクＮのガス圧力）２－（タンクＭのガス圧力）２）　
     
    　前記［数１］で求められる前記複数タンクの中の各タンクのガス量ＸＭの初期値を、各タンクのガス圧力、タンク容積、ガス固有の定数及びガス温度の初期値から求め、
    　前記［数２］で求められるタンクＮとタンクＭを結ぶ各ガス管路に流れるガス流量係数の初期値を、タンクＮとタンクＭの間の初期流量、及び、タンクＮとタンクＭのガス圧力の初期値から求めて、変数値としてメモリーに記憶する初期値設定段階と、
    　前記解析装置が、予め定められた演算時間（Δｔ）ごとに、各ガス管路に流れるガス流量の変化を表す値ＤＸＭを
    　ＤＸＭ＝ガス管路に流れる今回のガス流量－ガス管路に流れるΔｔ前のガス流量
    により求め、
    　ガス量ＸＭ＝前回のガス量ＸＭ＋ガス流量の変化ＤＸＭ×演算時間刻みΔｔ
    により、各ガス管路を結ぶタンクのガス量ＸＭを求め、
    ［数１］に従い、複数タンクのガス圧力を求めることを所定回数繰り返す段階と、を含み、
    　ガス管路のいくつかがアクチュエータを備え、前記解析装置がアクチュエータの開度を設定する段階を更に含む、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めることを特徴とするガス流量解析方法。
     
  2. 　前記アクチュエータの開度を設定する段階において、オペレータの入力に対応して前記アクチュエータを前記ガス管路の任意の場所に設定し、前記アクチュエータを全閉した場合の事故解析を行う、２系統以下のガス管路で従属接続される複数タンクのガス圧力及び、ガス管路に流れるガス流量を求めることを特徴とする請求項１に記載のガス流量解析方法。
PCT/JP2021/017599 2020-05-08 2021-05-07 ガス流量解析方法 WO2021225173A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020-082785 2020-05-08
JP2020082785A JP6908758B1 (ja) 2020-05-08 2020-05-08 ガス流量解析方法

Publications (1)

Publication Number Publication Date
WO2021225173A1 true WO2021225173A1 (ja) 2021-11-11

Family

ID=76967251

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/017599 WO2021225173A1 (ja) 2020-05-08 2021-05-07 ガス流量解析方法

Country Status (2)

Country Link
JP (1) JP6908758B1 (ja)
WO (1) WO2021225173A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03209097A (ja) * 1989-10-02 1991-09-12 Robert M Pierson 高圧のガスを圧力容器に迅速に充填する方法と装置
WO2011155539A1 (ja) * 2010-06-09 2011-12-15 三菱重工業株式会社 数値解析装置、要素生成プログラム、及び数値解析方法
JP2017041124A (ja) * 2015-08-20 2017-02-23 株式会社東芝 シミュレーション装置およびシミュレーション方法
WO2019031371A1 (ja) * 2017-08-07 2019-02-14 日本電気株式会社 状態分析装置、状態分析方法及び記録媒体

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03209097A (ja) * 1989-10-02 1991-09-12 Robert M Pierson 高圧のガスを圧力容器に迅速に充填する方法と装置
WO2011155539A1 (ja) * 2010-06-09 2011-12-15 三菱重工業株式会社 数値解析装置、要素生成プログラム、及び数値解析方法
JP2017041124A (ja) * 2015-08-20 2017-02-23 株式会社東芝 シミュレーション装置およびシミュレーション方法
WO2019031371A1 (ja) * 2017-08-07 2019-02-14 日本電気株式会社 状態分析装置、状態分析方法及び記録媒体

Also Published As

Publication number Publication date
JP6908758B1 (ja) 2021-07-28
JP2021177350A (ja) 2021-11-11

Similar Documents

Publication Publication Date Title
US9864383B2 (en) Method and system for determining characteristic parameters of a hydraulic network
CA2466932C (en) Flow measurement module and method
Behrooz et al. Modeling and state estimation for gas transmission networks
Wu et al. An empirical discharge coefficient model for orifice flow
Grace et al. Experimental parametric equation for the prediction of valve coefficient (Cv) for choke valve trims
Billmann et al. Leak detection methods for pipelines
Meniconi et al. Pressure reducing valve characterization for pipe system management
Diener et al. Sizing of throttling device for gas/liquid two‐phase flow part 2: Control valves, orifices, and nozzles
WO2021225173A1 (ja) ガス流量解析方法
Dempster et al. Prediction of the flow and force characteristics of safety relief valves
Jablonská et al. Flow characteristics of control valve for different strokes
Ozdamar et al. An Experimental and Numerical Study on Pressure Drop Coefficient of Ball Valves
Kovar et al. Simulation modelling of water-supply network as mechatronic system
Muftah CFD Modeling of elbow and orifice meters
Reader-Harris et al. Orifice Discharge Coefficient
CN114415748B (zh) 考虑空化影响的套筒调节阀阻力系数估计方法及估计系统
Rasaienia Accuracy Improvement in Sizing Control Valves Using Neural Network
Lončar et al. Numerical models for analysis of hydraulic transients
Likitha et al. CFD Analysis on Different Shapes of Concentric Orifice Plate for Turbulent Flow
Kang et al. A Study on the Comparison of Loss Coefficient by 1-inch Ball and Glove Valve Opening Ratio
Benaissa et al. MODIFICATION OF THE GEOMETRY OF THE FLOW CONDITIONER
Xu et al. Process variability reduction in flow control and wedge flowmeter performance evaluation
Sandberg et al. Venturi flowmeter performance installed downstream of the branch of a tee junction
Yusof et al. Validation of Flow Impact to Detect the Energy Loss in Ball Valve
Nabil et al. Experimental and numerical effects of flow hydraulics and pipe geometry on leakage behavior of laboratory water network distribution systems

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21800572

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21800572

Country of ref document: EP

Kind code of ref document: A1