CN109635405A - A kind of Multiphase Flow transient state calculation method and system based on space conservation - Google Patents
A kind of Multiphase Flow transient state calculation method and system based on space conservation Download PDFInfo
- Publication number
- CN109635405A CN109635405A CN201811479567.6A CN201811479567A CN109635405A CN 109635405 A CN109635405 A CN 109635405A CN 201811479567 A CN201811479567 A CN 201811479567A CN 109635405 A CN109635405 A CN 109635405A
- Authority
- CN
- China
- Prior art keywords
- fluid
- pipeline
- model
- conservation
- phase
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 57
- 238000004364 calculation method Methods 0.000 title claims abstract description 22
- 239000012530 fluid Substances 0.000 claims abstract description 259
- 238000004134 energy conservation Methods 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims description 43
- 230000004323 axial length Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 8
- 230000005501 phase interface Effects 0.000 claims description 4
- 238000009614 chemical analysis method Methods 0.000 claims 1
- 230000008602 contraction Effects 0.000 claims 1
- 238000003860 storage Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000004088 simulation Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 230000006854 communication Effects 0.000 description 5
- 230000006399 behavior Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of Multiphase Flow transient state calculation methods and system based on space conservation, comprising: obtains pipe parameter and fluid parameter；The transient data of heterogeneous fluid is obtained by hydrodynamic analysis method based on preset fluid and pipeline conservation model, continuity conservation model, conservation of momentum model and modeling of energy conservation according to the pipe parameter and fluid parameter, the present invention can be by increasing conservation model between conduit volume and fluid volume, the transient behaviour of incompressible and compressible fluid multi-fluid phase flowing under any flow velocity of analog.
Description
Technical field
The present invention relates to Multiphase Flow computing technique field more particularly to a kind of Multiphase Flow transient states based on space conservation
Calculation method and system.
Background technique
Multiphase Flow is widely present in chemical industry, nuclear industry and petroleum industry field.With automation science and technology high speed development,
Accurate efficient heterogeneous fluid pressure drop and the transient simulation prediction of phase fraction show important especially.It especially can in deep-sea oil industry
Accurately simulation and prediction stopping transportation, the complicated variable working condition problem such as restart and leak, simultaneously for highly effective and safe progress industrial production
It realizes green energy conservation consumption reduction, is of great significance.In addition, carrying out the transient multi-phase Numerical-Mode of high-precision, high efficiency, high reliability
It is quasi-, it can not only be produced with service industry, can also be technological design, improvement production efficiency, reduction operating cost etc., provide important
Technical support.
Traditional one-dimensional Multiphase Flow transient simulation constantly improves and changes on the basis of classical two-fluid model
Into, but various improved methods all have strong system suitability.For example, existing be suitable for incompressible fluid Multiphase Flow
Transient simulation method, when describing the Multiphase Flow problem of potential flows, it may appear that inclined to the simulation of pressure propagation problem
Difference；Analogy method suitable for low mach also will appear significant error during the solution of High Mach number problem.
Summary of the invention
The Multiphase Flow transient state calculation method based on space conservation that it is an object of the present invention to provide a kind of, the present invention
It is incompressible and compressible fluid multithread bulk flow under any flow velocity of analog by increase fluid and pipeline conservation model
Transient behaviour.The Multiphase Flow transient state computing system based on space conservation that it is another object of the present invention to provide a kind of.This
The further object of invention is to provide a kind of computer equipment, another purpose of the invention is that providing a kind of computer can
Read medium.
In order to reach the goals above, one aspect of the present invention discloses a kind of Multiphase Flow transient state calculating based on space conservation
Method, comprising:
Obtain pipe parameter and fluid parameter；
Preset fluid and pipeline conservation model, continuity conservation mould are based on according to the pipe parameter and fluid parameter
Type, conservation of momentum model and modeling of energy conservation obtain the transient data of heterogeneous fluid by hydrodynamic analysis method.
Preferably, the method also includes establishing the fluid to keep with pipeline conservation model, continuity conservation model, momentum
The step of permanent model and modeling of energy conservation.
Preferably, the fluid is established to specifically include with pipeline conservation model:
Determine the Conservation Relationship of the change rate of the long-pending expansion and shrinking percentage and fluid flow through conduit at any time of pipeline section；
The fluid and pipeline conservation model based on association phase fraction under heterogeneous fluid are established according to the Conservation Relationship.
Preferably, the Conservation Relationship are as follows:
Wherein, A is conduit cross-sectional area, and t is the time, and u is fluid velocity, and x is pipeline axial length.
Preferably, the fluid and pipeline conservation model are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, A is conduit cross-sectional area, and t is the time, and x is pipeline
Axial length.
Preferably, described that preset fluid and pipeline conservation model are based on according to the pipe parameter and fluid parameter, connected
Continuous property conservation model, conservation of momentum model and modeling of energy conservation obtain the wink of heterogeneous fluid by hydrodynamic analysis method
State data specifically include:
Single order upstreame scheme is selected to carry out pipeline model discrete；
It divides staggered-mesh and forms solving model；
According to the solving model and the fluid and pipeline conservation model, the continuity conservation model, the conservation of momentum
The transient data of heterogeneous fluid is calculated in model and modeling of energy conservation.
Preferably, the continuity conservation model are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, t is the time, and x is pipeline axial length.
Preferably, the conservation of momentum model are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, and k is in contact with i-th of fluid
The number of fluid phase, αiFor the phase fraction of i-th of fluid phase, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase
Degree, p is pipeline pressure, pikPressure between the phase interface for k fluid phase being mutually in contact with it for i-th of fluid, ГiIt is flowed for i-th
The frictional resistance of body phase and pipe wall, ГikIt is all i-th of fluids opposite with i-th of fluid is in contact k fluid mutually in pipeline
Axial effect resultant force in its preceding symbol is "+" if the force direction is contrary with fluid flow rate, if the active force side
It is "-" to then its preceding symbol identical as fluid flow rate direction, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is
Time, x are pipeline axial length.
Preferably, the modeling of energy conservation are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, p is pipeline pressure, EiIt is i-th of fluid compared to interior
Can, HiFor the ratio total enthalpy of i-th of fluid phase, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is the time, and x is
Pipeline axial length.
The Multiphase Flow transient state computing system based on space conservation that the invention also discloses a kind of, the system comprises models
Establish unit and model solution unit；
The model foundation unit is based on being associated with phase fraction between pipeline and the heterogeneous fluid for flowing through the pipeline for establishing
Fluid and pipeline conservation model, continuity conservation model, conservation of momentum model and modeling of energy conservation；
The model solution unit is used for according to the fluid and pipeline conservation model, the continuity conservation model, moves
It measures conservation model and modeling of energy conservation and the transient data of heterogeneous fluid is obtained by hydrodynamic analysis method.
The invention also discloses a kind of computer equipment, including memory, processor and storage are on a memory and can
The computer program run on a processor,
The processor realizes method as described above when executing described program.
The invention also discloses a kind of computer-readable mediums, are stored thereon with computer program,
The program realizes method as described above when being executed by processor.
The present invention can be achieved the Multiphase Flow transient state based on space conservation and calculate, i.e., by increasing the flowing of description fluid and pipe
Conservation Relationship between road is into conventional multiphase fluid model, incompressible under any flow velocity of analog and compressible fluid multithread
The Transient Flow of body phase.This method can overcome the calculating of conventional multiphase flow simulating be difficult to accurate description fluid and pipeline space it
Between relationship defect, it is advantageous that: incompressible under any flow velocity of analog and compressible fluid multi-fluid phase transient flow
It is dynamic, obtain high-precision, high efficiency, the transient multi-phase numerical simulation calculation of high reliability as a result, can be applied to chemical industry, nuclear industry,
The Multiphase Flow transient state of petroleum industry production calculates, and serves technological design in industry, production efficiency is improved and energy-saving
Etc..
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 shows a kind of stream of Multiphase Flow transient state one specific embodiment of calculation method based on space conservation of the present invention
Cheng Tu；
Fig. 2 shows a kind of one specific embodiments of Multiphase Flow transient state calculation method based on space conservation of the present invention to determine
The flow chart of conservation model between conduit volume and fluid volume；
Fig. 3 shows a kind of one specific embodiment pipeline section of Multiphase Flow transient state calculation method based on space conservation of the present invention
The schematic diagram of infinitesimal；
Fig. 4 shows a kind of one specific embodiment fluid of Multiphase Flow transient state calculation method based on space conservation of the present invention
The flow chart of dynamic analysis；
Fig. 5 shows a kind of one specific embodiment pipeline of Multiphase Flow transient state calculation method based on space conservation of the present invention
The schematic diagram of grid dividing；
Fig. 6 shows a kind of one specific embodiment classics of Multiphase Flow transient state calculation method based on space conservation of the present invention
The comparison diagram of original state, transient state and stable state in tap example；
Fig. 7 shows a kind of one specific embodiment fire hose of Multiphase Flow transient state calculation method based on space conservation of the present invention
One of the result figure of head example；
Fig. 8 shows a kind of one specific embodiment fire hose of Multiphase Flow transient state calculation method based on space conservation of the present invention
The two of the result figure of head example；
Fig. 9 shows a kind of one specific embodiment fire hose of Multiphase Flow transient state calculation method based on space conservation of the present invention
The three of the result figure of head example；
Figure 10 shows a kind of one specific embodiment of Multiphase Flow transient state computing system based on space conservation of the present invention
Structure chart；
Figure 11 shows a kind of Multiphase Flow transient state computing system another specific embodiment based on space conservation of the present invention
Structure chart；
Figure 12 shows the structural schematic diagram for being suitable for the computer equipment for being used to realize the embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
According to an aspect of the present invention, present embodiment discloses a kind of, and the Multiphase Flow transient state based on space conservation calculates
Method.As shown in Figure 1, in the present embodiment, the method 10 includes:
S100: pipe parameter and fluid parameter are obtained.Wherein, pipe parameter may include conduit cross-sectional area and duct length
Etc. parameters, fluid parameter may include the parameters such as the initial velocity of the number of fluid phase, the density of fluid and fluid.
S200: preset fluid is based on according to the pipe parameter and fluid parameter and is kept with pipeline conservation model, continuity
Permanent model, conservation of momentum model and modeling of energy conservation obtain the transient state number of heterogeneous fluid by hydrodynamic analysis method
According to.
The present invention by fluid in the duct flow mechanism intrinsic property study, in heterogeneous fluid model increase pipe
Conservation model between road volume and fluid volume is incompressible and compressible fluid multithread bulk flow under any flow velocity of analog
Transient behaviour improves conventional multiphase flow simulating and is difficult to the solution accuracy of Multiphase Flow transient data and high reliablity
The defect of relationship between accurate description fluid and pipeline space can be adapted in the wider array of Multiphase Flow calculating of range, alignment
The really Multiphase Flow problem in simulation chemical industry, nuclear industry and petroleum industry, is of great significance.
In a preferred embodiment, the Multiphase Flow transient state calculation method based on space conservation may also include described in foundation
The step of fluid and pipeline conservation model, continuity conservation model, conservation of momentum model and modeling of energy conservation.Specifically, such as
Shown in Fig. 2, the S100 is specific can include:
S110: determine that the conservation of the change rate of the long-pending expansion and shrinking percentage and fluid flow through conduit at any time of pipeline section closes
System.Preferably, for pipeline section infinitesimal as shown in Figure 3, the variation of conduit volume derived from the expansion at any time of pipeline section product and
Shrinking percentage, the expansion and shrinking percentage of pipeline section product at any time are as follows:
Wherein, A is conduit cross-sectional area, and Δ x is the length of pipeline infinitesimal, and t is the time.
When for only having a kind of fluid of fluid phase to flow through in pipeline section infinitesimal as shown in Figure 3, the fluid body of pipeline is flowed through
Long-pending variation is derived from the change rate of the inflow and outflow duct element, the change rate of fluid flow through conduit are as follows:
Wherein, A is conduit cross-sectional area, and Δ x is the length of pipeline infinitesimal, and u is fluid velocity.
Then the Conservation Relationship of the change rate of the long-pending expansion and shrinking percentage and fluid flow through conduit at any time of pipeline section can table
It is shown as:
Wherein, A is conduit cross-sectional area, and t is the time, and u is fluid velocity, and x is pipeline axial length.
S120: the fluid and pipeline conservation mould based on association phase fraction under heterogeneous fluid are established according to the Conservation Relationship
Type.It in the full space-time unique for considering pipeline, is filled always by least any fluid phase in pipeline, derives and establish association phase fraction
Fluid and pipeline conservation model, which may include two, and function of the one of them as fluid density, embodying can
The fluctuation property of compression fluid, wherein another function as fluid velocity, embodies the flow properties of incompressible fluid.
The present invention can improve two-fluid model, mutually all have preferable be applicable in compressible and incompressible fluid for multi-fluid phase
Property, multi-fluid phase simulation accuracy and high reliablity.
Preferably, the fluid and pipeline conservation model can are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, A is conduit cross-sectional area, and t is the time, and x is pipeline
Axial length.
In a preferred embodiment, the continuity conservation model of the multi-fluid are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, t is the time, and x is pipeline axial length.
The conservation of momentum model are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, and k is in contact with i-th of fluid
The number of fluid phase, αiFor the phase fraction of i-th of fluid phase, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase
Degree, p is pipeline pressure, pikPressure between the phase interface for k fluid phase being mutually in contact with it for i-th of fluid, ГiIt is flowed for i-th
The frictional resistance of body phase and pipe wall, ГikIt is all i-th of fluids opposite with i-th of fluid is in contact k fluid mutually in pipeline
Axial effect resultant force in its preceding symbol is "+" if the force direction is contrary with fluid flow rate, if the active force side
It is "-" to then its preceding symbol identical as fluid flow rate direction, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is
Time, x are pipeline axial length.
The modeling of energy conservation are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, p is pipeline pressure, EiIt is i-th of fluid compared to interior
Can, HiFor the ratio total enthalpy of i-th of fluid phase, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is the time, and x is
Pipeline axial length.
In a preferred embodiment, the transient state of heterogeneous fluid can be calculated by hydrodynamic method as shown in Figure 4
Data, the S300 are specific can include:
S310: selection single order upstreame scheme is discrete to pipeline model progress, as shown in Figure 5.
S320: it divides staggered-mesh and forms solving model.
S330: preset fluid is based on according to the pipe parameter and fluid parameter and is kept with pipeline conservation model, continuity
Permanent model, conservation of momentum model and modeling of energy conservation solve the solving model and obtain the transient data of heterogeneous fluid, such as
Pipeline pressure (the p in solving model at discrete point can be solved to obtainW、pP、pE), fluid density (αW、αP、αE) and fluid phase phase
Divide rate (ρW、ρE、ρP) and discrete subregion fluid velocity (uw、ue、uee)。
Original state, transition shape shown in Fig. 6 in the specific example of a tap, when fluid flows through faucet conduit
The comparison diagram of state and stable state.Multiphase Flow transient state calculation method through the invention obtains fluid as shown in figs. 7 to 9
Transient data as a result, can be seen that from such as Fig. 7-Fig. 9, using the Multiphase Flow transient state calculating side of the invention based on space conservation
Method and economics analysis solution data have the very high goodness of fit, indicate that analogy method of the invention can accurately and efficiently simulate multithread
Mutually transient behaviour in the duct.
According to another aspect of the present invention, the present embodiment also discloses based on a kind of Multiphase Flow transient state by space conservation
Calculation system.As shown in Figure 10, the system comprises parameter acquiring units 11 and model solution unit 12.
Wherein, the parameter acquiring unit 11 is for obtaining pipe parameter and fluid parameter.Wherein, pipe parameter may include
The parameters such as conduit cross-sectional area and duct length, fluid parameter may include the first of the number of fluid phase, the density of fluid and fluid
The parameters such as beginning speed.
The model solution unit 12 is used to be based on preset fluid and pipeline according to the pipe parameter and fluid parameter
Conservation model, continuity conservation model, conservation of momentum model and modeling of energy conservation are obtained by hydrodynamic analysis method
The transient data of heterogeneous fluid.
The present invention by fluid in the duct flow mechanism intrinsic property study, in heterogeneous fluid model increase pipe
Conservation model between road volume and fluid volume is incompressible and compressible fluid multithread bulk flow under any flow velocity of analog
Transient behaviour improves conventional multiphase flow simulating and is difficult to the solution accuracy of Multiphase Flow transient data and high reliablity
The defect of relationship between accurate description fluid and pipeline space can be adapted in the wider array of Multiphase Flow calculating of range, alignment
The really Multiphase Flow problem in simulation chemical industry, nuclear industry and petroleum industry, is of great significance.
In a preferred embodiment, as shown in figure 11, the system also includes model foundation unit 13, the model is built
Vertical unit can be used for establishing the fluid and pipeline conservation model, continuity conservation model, conservation of momentum model and the conservation of energy
Model.Preferably, model foundation unit 13 can determine the long-pending expansion and shrinking percentage and fluid flow through conduit at any time of pipeline section
Change rate Conservation Relationship, and according to the Conservation Relationship establish under heterogeneous fluid based on association phase fraction fluid and pipeline
Conservation model.It in the full space-time unique for considering pipeline, is filled always by least any fluid phase in pipeline, derives and establish association
The fluid and pipeline conservation model of phase fraction, which may include two, function of the one of them as fluid density, body
The fluctuation property of compressible fluid is showed, wherein another function as fluid velocity, embodies the stream of incompressible fluid
Dynamic attribute.The present invention can improve two-fluid model, for multi-fluid phase and compressible and incompressible fluid mutually all have compared with
Good applicability, multi-fluid phase simulation accuracy and high reliablity.
Preferably, for pipeline section infinitesimal as shown in Figure 3, the variation of conduit volume is derived from pipeline section product at any time swollen
Swollen and shrinking percentage, the expansion and shrinking percentage of pipeline section product at any time are as follows:
Wherein, A is conduit cross-sectional area, and Δ x is the length of pipeline infinitesimal, and t is the time.
When for only having a kind of fluid of fluid phase to flow through in pipeline section infinitesimal as shown in Figure 3, the fluid body of pipeline is flowed through
Long-pending variation is derived from the change rate of the inflow and outflow duct element, the change rate of fluid flow through conduit are as follows:
Wherein, A is conduit cross-sectional area, and Δ x is the length of pipeline infinitesimal, and u is fluid velocity.
Then the Conservation Relationship of the change rate of the long-pending expansion and shrinking percentage and fluid flow through conduit at any time of pipeline section can table
It is shown as:
Wherein, A is conduit cross-sectional area, and t is the time, and u is fluid velocity, and x is pipeline axial length.
Preferably, the fluid and pipeline conservation model can are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, A is conduit cross-sectional area, and t is the time, and x is pipeline
Axial length.
In a preferred embodiment, the continuity conservation model of the multi-fluid are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, t is the time, and x is pipeline axial length.
The conservation of momentum model are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, and k is in contact with i-th of fluid
The number of fluid phase, αiFor the phase fraction of i-th of fluid phase, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase
Degree, p is pipeline pressure, pikPressure between the phase interface for k fluid phase being mutually in contact with it for i-th of fluid, ГiIt is flowed for i-th
The frictional resistance of body phase and pipe wall, ГikIt is all i-th of fluids opposite with i-th of fluid is in contact k fluid mutually in pipeline
Axial effect resultant force in its preceding symbol is "+" if the force direction is contrary with fluid flow rate, if the active force side
It is "-" to then its preceding symbol identical as fluid flow rate direction, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is
Time, x are pipeline axial length.
The modeling of energy conservation are as follows:
Wherein, i is fluids within pipes phase number, and i=1,2 ..n, n are positive integer, αiFor the phase point of i-th of fluid phase
Rate, ρiFor the density of i-th of fluid phase, uiFor the speed of i-th of fluid phase, p is pipeline pressure, EiIt is i-th of fluid compared to interior
Can, HiFor the ratio total enthalpy of i-th of fluid phase, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is the time, and x is
Pipeline axial length.
In a preferred embodiment, the model solution unit 12 is further useful for selection single order upstreame scheme to pipe
Road carries out discrete, division staggered-mesh formation solving model, according to the pipe parameter and fluid parameter based on preset fluid
The solving model is solved with pipeline conservation model, continuity conservation model, conservation of momentum model and modeling of energy conservation to obtain
The transient data of heterogeneous fluid.Such as it can solve to obtain pipeline pressure (p in solving model at discrete pointW、pP、pE), fluid it is close
Spend (αW、αP、αE) and fluid phase phase fraction (ρW、ρE、ρP) and discrete subregion fluid velocity (uw、ue、uee)。
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer equipment, specifically, computer is set
It is standby for example can for personal computer, laptop computer, cellular phone, camera phone, smart phone, personal digital assistant,
Media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or these equipment
In any equipment combination.
Computer equipment specifically includes memory, processor and storage on a memory simultaneously in a typical example
The computer program that can be run on a processor is realized when the processor executes described program and is held as described above by client
Capable method, alternatively, the processor realizes the method executed as described above by server when executing described program.
Below with reference to Figure 12, it illustrates the structures for the computer equipment 600 for being suitable for being used to realize the embodiment of the present application to show
It is intended to.
As shown in figure 12, computer equipment 600 includes central processing unit (CPU) 601, can be read-only according to being stored in
Program in memory (ROM) 602 is loaded into random access storage device (RAM) from storage section 608) program in 603
And execute various work appropriate and processing.In RAM603, also it is stored with system 600 and operates required various program sum numbers
According to.CPU601, ROM602 and RAM603 are connected with each other by bus 604.Input/output (I/O) interface 605 is also connected to
Bus 604.
I/O interface 605 is connected to lower component: the importation 606 including keyboard, mouse etc.；It is penetrated including such as cathode
The output par, c 607 of spool (CRT), liquid crystal ultramagnifier (LCD) etc. and loudspeaker etc.；Storage section 608 including hard disk etc.；
And including such as LAN card, the communications portion 609 of the network interface card of modem etc..Communications portion 609 via such as because
The network of spy's net executes communication process.Driver 610 is also connected to I/O interface 606 as needed.Detachable media 611, such as
Disk, CD, magneto-optic disk, semiconductor memory etc. are mounted on as needed on driver 610, in order to read from thereon
Computer program be mounted as needed such as storage section 608.
Particularly, according to an embodiment of the invention, may be implemented as computer above with reference to the process of flow chart description
Software program.For example, the embodiment of the present invention includes a kind of computer program product comprising be tangibly embodied in machine readable
Computer program on medium, the computer program include the program code for method shown in execution flow chart.At this
In the embodiment of sample, which can be downloaded and installed from network by communications portion 609, and/or from removable
Medium 611 is unloaded to be mounted.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM),
Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each unit can be realized in the same or multiple software and or hardware when application.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method of element, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (12)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811479567.6A CN109635405B (en) | 2018-12-05 | 2018-12-05 | A method and system for transient calculation of multiphase flow based on space conservation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811479567.6A CN109635405B (en) | 2018-12-05 | 2018-12-05 | A method and system for transient calculation of multiphase flow based on space conservation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109635405A true CN109635405A (en) | 2019-04-16 |
| CN109635405B CN109635405B (en) | 2020-12-18 |
Family
ID=66071305
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811479567.6A Active CN109635405B (en) | 2018-12-05 | 2018-12-05 | A method and system for transient calculation of multiphase flow based on space conservation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109635405B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110059411A (en) * | 2019-04-19 | 2019-07-26 | 中国海洋石油集团有限公司 | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation |
| CN111400950A (en) * | 2020-03-17 | 2020-07-10 | 中国石油大学(北京) | Transient flow simulation method and device for hydrate slurry multiphase pipeline |
| CN113792449A (en) * | 2021-08-10 | 2021-12-14 | 中国船舶重工集团公司第七一九研究所 | Nuclear reactor monitoring method and system |
| CN116011352A (en) * | 2021-10-22 | 2023-04-25 | 中国人民解放军国防科技大学 | Method for unconditionally maintaining strong stability of hyperbolic conservation law of computational fluid dynamics |
| CN117669407A (en) * | 2023-10-30 | 2024-03-08 | 中国石油大学(北京) | Construction method and construction device for pipeline production multiphase flow thermodynamic model |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120118077A1 (en) * | 2010-08-24 | 2012-05-17 | Invensys Systems, Inc. | Multiphase metering system |
| CN102943620A (en) * | 2012-08-27 | 2013-02-27 | 中国石油大学(华东) | Pressure-controlled drilling method based on drilling annulus wellbore multi-phase flow computing |
| CN103149012A (en) * | 2013-03-08 | 2013-06-12 | 中国石油大学(北京) | Experimental device and method for simulating flow characteristic of coal bed gas gathering and transportation pipeline |
| CN104594834A (en) * | 2014-12-01 | 2015-05-06 | 中国石油大学(华东) | Method for monitoring drilling overflow condition of deepwater oil-based drilling fluid |
| CN104714049A (en) * | 2013-12-17 | 2015-06-17 | 国际商业机器公司 | Computer based fluid flow velocity estimation from concentrations of a reacting constituent for products and services |
| CN105387895A (en) * | 2015-10-16 | 2016-03-09 | 中国石油天然气股份有限公司规划总院 | Long-throat-and-neck venturi tube moisture flow measurement method based on flowing density |
-
2018
- 2018-12-05 CN CN201811479567.6A patent/CN109635405B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120118077A1 (en) * | 2010-08-24 | 2012-05-17 | Invensys Systems, Inc. | Multiphase metering system |
| CN102943620A (en) * | 2012-08-27 | 2013-02-27 | 中国石油大学(华东) | Pressure-controlled drilling method based on drilling annulus wellbore multi-phase flow computing |
| CN103149012A (en) * | 2013-03-08 | 2013-06-12 | 中国石油大学(北京) | Experimental device and method for simulating flow characteristic of coal bed gas gathering and transportation pipeline |
| CN104714049A (en) * | 2013-12-17 | 2015-06-17 | 国际商业机器公司 | Computer based fluid flow velocity estimation from concentrations of a reacting constituent for products and services |
| CN104594834A (en) * | 2014-12-01 | 2015-05-06 | 中国石油大学(华东) | Method for monitoring drilling overflow condition of deepwater oil-based drilling fluid |
| CN105387895A (en) * | 2015-10-16 | 2016-03-09 | 中国石油天然气股份有限公司规划总院 | Long-throat-and-neck venturi tube moisture flow measurement method based on flowing density |
Non-Patent Citations (3)
| Title |
|---|
| BERENDT,A等: "Transition of the electrohydrodynamic two-phase flow into the single-phase flow in a needle-to-plate negative corona discharge in the finite-volume chamber", 《JOURNAL OF ELECTROSTATICS》 * |
| K.MINAMI等: "Transient two-phase flow behavior in pipelines-experiment and modeling", 《INTERNATIONAL JOURNAL OF MULTIPHASE FLOW》 * |
| 宋尚飞 等: "基于LedaFlow的天然气凝析液管道瞬变过程模拟", 《油气储运》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110059411A (en) * | 2019-04-19 | 2019-07-26 | 中国海洋石油集团有限公司 | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation |
| CN111400950A (en) * | 2020-03-17 | 2020-07-10 | 中国石油大学(北京) | Transient flow simulation method and device for hydrate slurry multiphase pipeline |
| CN111400950B (en) * | 2020-03-17 | 2022-03-15 | 中国石油大学(北京) | Hydrate slurry multiphase pipeline transient flow simulation method and device |
| CN113792449A (en) * | 2021-08-10 | 2021-12-14 | 中国船舶重工集团公司第七一九研究所 | Nuclear reactor monitoring method and system |
| CN116011352A (en) * | 2021-10-22 | 2023-04-25 | 中国人民解放军国防科技大学 | Method for unconditionally maintaining strong stability of hyperbolic conservation law of computational fluid dynamics |
| CN117669407A (en) * | 2023-10-30 | 2024-03-08 | 中国石油大学(北京) | Construction method and construction device for pipeline production multiphase flow thermodynamic model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109635405B (en) | 2020-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109635405A (en) | A kind of Multiphase Flow transient state calculation method and system based on space conservation | |
| Benra et al. | A comparison of one‐way and two‐way coupling methods for numerical analysis of fluid‐structure interactions | |
| Simpson et al. | Evaluating the approximation of the affinity laws and improving the efficiency estimate for variable speed pumps | |
| Ramos et al. | Flat vs. curved rigid-lid LES computations of an open-channel confluence | |
| Dongari et al. | Effects of curvature on rarefied gas flows between rotating concentric cylinders | |
| Iliev et al. | Optimization of Francis turbines for variable speed operation using surrogate modeling approach | |
| CN107524921A (en) | Products Batch Transportation Pipeline batching interface tracking and device | |
| Stansby | Limitations of depth-averaged modeling for shallow wakes | |
| Meziou et al. | Low-dimensional modeling of transient two-phase flow in pipelines | |
| Venutelli | Analysis of dynamic wave model for unsteady flow in an open channel | |
| Zhang et al. | Modeling and simulation of pigging for a gas pipeline using a bypass pig | |
| CN107526914B (en) | Variable-watershed flow field calculation method of tilting-pad sliding bearing based on structured dynamic grid | |
| Yuan | Numerical investigation of wavy and spiral Taylor-Görtler vortices in medium spherical gaps | |
| Rosén et al. | Prediction of turbulent high Schmidt number mass transfer using a low Reynolds number k—ϵ turbulence model | |
| Zolfaghary Azizi et al. | Comparison of numerical methods for two-fluid model for gas–liquid transient flow regime and its application in slug modeling initiation | |
| CN107992642A (en) | Take the Forecasting Methodology and device of liquid critical flow velocity | |
| CN112949112A (en) | Rotor-sliding bearing system lubrication basin dynamic grid parallel computing method | |
| Furnes | Flow in Pelton turbines | |
| Patil et al. | Centrifugal pump head prediction using affinity laws modified for viscosity | |
| Li et al. | An improved smoothed particle hydrodynamics method for modeling multiphase flows | |
| Ooms et al. | Numerical study of laminar core‐annular flow in a torus and in a 90° pipe bend | |
| CN110232222B (en) | Deposition tube flow field analysis method and system | |
| Pochylý et al. | Rotational flow in centrifugal pump meridian using curvilinear coordinates | |
| Yang et al. | A direct numerical simulation method for solid-solid collision and coupling with fluid | |
| Zanotti et al. | A preconditioning mass matrix to avoid the ill-posed two-fluid model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information |
Inventor after: Shi Bohui Inventor after: Shi Guobin Inventor after: Kang Qi Inventor after: Gong Jing Inventor after: Wu Haihao Inventor after: Li Xiaoping Inventor before: Shi Bohui Inventor before: Shi Guobin Inventor before: Kang Qi Inventor before: Gong Jing Inventor before: Wu Haihao Inventor before: Li Xiaoping |
|
| CB03 | Change of inventor or designer information |