CN109635405B - Multiphase flow transient calculation method and system based on space conservation - Google Patents

Multiphase flow transient calculation method and system based on space conservation Download PDF

Info

Publication number
CN109635405B
CN109635405B CN201811479567.6A CN201811479567A CN109635405B CN 109635405 B CN109635405 B CN 109635405B CN 201811479567 A CN201811479567 A CN 201811479567A CN 109635405 B CN109635405 B CN 109635405B
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.)
Active
Application number
CN201811479567.6A
Other languages
Chinese (zh)
Other versions
CN109635405A (en
Inventor
史博会
石国赟
康琪
宫敬
吴海浩
李晓平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201811479567.6A priority Critical patent/CN109635405B/en
Publication of CN109635405A publication Critical patent/CN109635405A/en
Application granted granted Critical
Publication of CN109635405B publication Critical patent/CN109635405B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 invention provides a multiphase flow transient calculation method and a system based on space conservation, which comprises the following steps: acquiring pipeline parameters and fluid parameters; according to the pipeline parameters and the fluid parameters, transient data of the multiphase fluid are obtained through a fluid dynamics analysis method based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model.

Description

Multiphase flow transient calculation method and system based on space conservation
Technical Field
The invention relates to the technical field of multi-phase flow calculation, in particular to a multi-phase flow transient calculation method and system based on space conservation.
Background
Multiphase flow is widely found in the chemical, nuclear and petroleum industries. With the rapid development of automation technology, accurate and efficient transient simulation prediction of multiphase fluid pressure drop and phase fraction is particularly important. Especially, the method can accurately simulate and predict the problems of complex variable working conditions such as stop transportation, restart, leakage and the like in the deep sea petroleum industry, and has important significance for efficiently and safely carrying out industrial production and realizing green energy conservation and consumption reduction. In addition, the transient multiphase numerical simulation with high precision, high efficiency and high reliability is developed, so that the method not only can serve industrial production, but also can provide important technical support for process design, improvement of production efficiency, reduction of operation cost and the like.
The traditional one-dimensional multiphase flow transient simulation is continuously improved and improved on the basis of a classical two-fluid model, but various improved methods have strong system applicability. For example, when describing the problem of the multiphase flow of the incompressible flow, the current method is suitable for the transient simulation method of the multiphase flow of the incompressible fluid, and the simulation deviation of the pressure propagation problem occurs; the simulation method suitable for low mach numbers also has significant errors in the process of solving the problem of high mach numbers.
Disclosure of Invention
The invention aims to provide a multiphase flow transient calculation method based on space conservation. Another object of the present invention is to provide a multiphase flow transient computing system based on space conservation. It is a further object of the present invention to provide a computer apparatus and it is a further object of the present invention to provide a computer readable medium.
In order to achieve the above object, the present invention discloses a multiphase flow transient calculation method based on space conservation, which includes:
acquiring pipeline parameters and fluid parameters;
and obtaining transient data of the multiphase fluid by a fluid dynamics analysis method based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model according to the pipeline parameters and the fluid parameters.
Preferably, the method further comprises the step of establishing the fluid and pipe conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model.
Preferably, the establishing of the fluid and pipeline conservation model specifically comprises:
determining the conservation relation between the expansion and contraction rate of the pipeline sectional area along with time and the change rate of the fluid flowing through the pipeline;
and establishing a fluid and pipeline conservation model based on the associated phase fraction under the multiphase fluid according to the conservation relation.
Preferably, the conservation relation is:
Figure BDA0001893038690000021
wherein A is the cross-sectional area of the pipeline, t is time, u is the fluid velocity, and x is the axial length of the pipeline.
Preferably, the fluid and pipe conservation model is:
Figure BDA0001893038690000022
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiThe velocity of the ith fluid phase, A is the cross-sectional area of the conduit, t is the time, and x is the axial length of the conduit.
Preferably, the obtaining of the transient data of the multiphase fluid by the fluid dynamics analysis method based on the preset fluid and pipeline conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model according to the pipeline parameters and the fluid parameters specifically includes:
selecting a first-order windward format to disperse the pipeline model;
dividing the staggered grids to form a solution model;
and calculating to obtain transient data of the multiphase fluid according to the solving model, the fluid and pipeline conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model.
Preferably, the continuity conservation model is:
Figure BDA0001893038690000023
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, t is the time, and x is the axial length of the pipe.
Preferably, the momentum conservation model is:
Figure BDA0001893038690000031
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, k is the number of fluid phases in contact with the ith fluid, and αiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, pikIs the pressure between the phase interfaces of the ith fluid phase and the k fluid phases with which it is in contact, fiIs the friction between the ith fluid phase and the pipe wall, rikThe acting force of all k fluids in contact with the ith fluid relative to the ith fluid phase in the axial direction of the pipeline is marked as "+" if the acting force direction is opposite to the flow velocity direction of the fluids, the former mark is marked as "-" if the acting force direction is the same as the flow velocity direction of the fluids, theta is the inclination angle between the pipeline and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipeline.
Preferably, the energy conservation model is:
Figure BDA0001893038690000032
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, EiIs the i-th fluid phase internal energy, HiIs the specific total enthalpy of the ith fluid phase, theta is the angle of inclination between the pipe and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipe.
The invention also discloses a multiphase flow transient computing system based on space conservation, which comprises a model establishing unit and a model solving unit;
the model establishing unit is used for establishing a fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model based on the correlation phase fraction between the pipeline and the multiphase fluid flowing through the pipeline;
the model solving unit is used for obtaining transient data of the multiphase fluid through a fluid dynamics analysis method according to the fluid and pipeline conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model.
The invention also discloses a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor,
the processor, when executing the program, implements the method as described above.
The invention also discloses a computer-readable medium, having stored thereon a computer program,
which when executed by a processor implements the method as described above.
The invention can realize the multiphase flow transient calculation based on space conservation, namely, the transient flow of incompressible and compressible fluid and multi-fluid phase under any flow rate can be simulated by adding the conservation relation between the described fluid flow and the pipeline into the traditional multiphase fluid model. The method can overcome the defect that the traditional multiphase flow simulation calculation is difficult to accurately describe the relationship between the fluid and the pipeline space, and has the advantages that: the method can simulate the transient flow of incompressible and compressible fluid multi-fluid phases at any flow rate, obtain the transient multi-phase numerical simulation calculation result with high precision, high efficiency and high reliability, can be applied to the multi-phase flow transient calculation in the production of chemical industry, nuclear industry and petroleum industry, and serves the aspects of process design, production efficiency improvement, energy conservation, consumption reduction and the like in the industry.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of an embodiment of a multiphase flow transient calculation method based on space conservation according to the present invention;
FIG. 2 is a flow chart of a method for multiphase flow transient calculation based on space conservation according to an embodiment of the present invention for determining a conservation model between a pipe volume and a fluid volume;
FIG. 3 is a schematic diagram of pipe section micro-elements according to an embodiment of the present invention;
FIG. 4 is a flow chart of a fluid dynamics analysis of an embodiment of a method for multiphase flow transient computation based on space conservation in accordance with the present invention;
FIG. 5 is a schematic diagram illustrating pipe meshing according to an embodiment of the present invention;
FIG. 6 is a diagram showing a comparison of initial state, transient state and steady state in a classical faucet algorithm according to an embodiment of the present invention;
FIG. 7 is a graph showing the results of an example faucet algorithm according to an embodiment of the present invention;
FIG. 8 is a second graph showing the results of an example faucet based on a space conservation-based multiphase flow transient calculation method of an embodiment of the present invention;
FIG. 9 is a third graph showing the results of an example faucet calculation of an embodiment of a multiphase flow transient calculation method based on space conservation according to the present invention;
FIG. 10 is a block diagram illustrating one embodiment of a space conservation based multiphase flow transient computing system in accordance with the present invention;
FIG. 11 is a block diagram illustrating another embodiment of a multiphase flow transient computing system based on space conservation in accordance with the present invention;
FIG. 12 shows a schematic block diagram of a computer device suitable for use in implementing embodiments of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
According to one aspect of the invention, the embodiment discloses a multiphase flow transient calculation method based on space conservation. As shown in fig. 1, in this embodiment, the method 10 includes:
s100: pipeline parameters and fluid parameters are obtained. The pipe parameters may include pipe cross-sectional area, pipe length and other parameters, and the fluid parameters may include the number of fluid phases, the density of the fluid, the initial velocity of the fluid and other parameters.
S200: and obtaining transient data of the multiphase fluid by a fluid dynamics analysis method based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model according to the pipeline parameters and the fluid parameters.
According to the invention, through the intrinsic characteristic research of the flowing mechanism of the fluid in the pipeline, the conservation model between the pipeline volume and the fluid volume is added in the multiphase fluid model, the transient state of the flow of incompressible and compressible fluid and multiple fluids at any flow rate can be simulated, the solving accuracy and reliability of the multiphase flow transient data are high, the defect that the relation between the fluid and the pipeline space is difficult to accurately describe in the traditional multiphase flow simulation is overcome, the method can be applied to the multiphase flow calculation with a wider range, and has important significance for accurately simulating the multiphase flow problem in the chemical industry, the nuclear industry and the petroleum industry.
In a preferred embodiment, the multiphase flow transient calculation method based on space conservation may further include the step of establishing a conservation model of the fluid and the pipeline, a conservation model of continuity, a conservation model of momentum, and a conservation model of energy. Specifically, as shown in fig. 2, the S100 may specifically include:
s110: the conservation of the rate of change of fluid flow through the pipe is determined as a function of the expansion and contraction of the pipe cross-sectional area over time. Preferably, for a pipe segment infinitesimal as shown in fig. 3, the change in pipe volume results from the expansion and contraction of the pipe cross-sectional area over time as follows:
Figure BDA0001893038690000051
wherein A is the cross-sectional area of the pipeline, Deltax is the length of the micro-element of the pipeline, and t is time.
For a fluid flow of only one fluid phase in a pipe segment infinitesimal as shown in fig. 3, the change in volume of the fluid flowing through the pipe results from the rate of change of flow into and out of the pipe element:
Figure BDA0001893038690000061
wherein A is the cross-sectional area of the pipeline, Deltax is the length of the micro-element of the pipeline, and u is the fluid velocity.
The conservation of the rate of change of fluid flow through the conduit with respect to the rate of expansion and contraction of the cross-sectional area of the conduit over time can be expressed as:
Figure BDA0001893038690000062
wherein A is the cross-sectional area of the pipeline, t is time, u is the fluid velocity, and x is the axial length of the pipeline.
S120: and establishing a fluid and pipeline conservation model based on the associated phase fraction under the multiphase fluid according to the conservation relation. Taking into account that the pipe is always filled with at least any one fluid phase in the pipe in the full space-time range of the pipe, a conservation model of the fluid and the pipe is derived for establishing the related phase fraction, wherein the conservation model can comprise two terms, wherein one term represents the fluctuation property of the compressible fluid as a function of the density of the fluid, and the other term represents the flow property of the incompressible fluid as a function of the velocity of the fluid. The invention can improve a double-fluid model, has better applicability to both multi-fluid phases and compressible and incompressible fluid phases, and has high simulation accuracy and reliability of the multi-fluid phases.
Preferably, the fluid and pipe conservation model may be:
Figure BDA0001893038690000063
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiThe velocity of the ith fluid phase, A is the cross-sectional area of the conduit, t is the time, and x is the axial length of the conduit.
In a preferred embodiment, the continuous conservation model for multiple fluids is:
Figure BDA0001893038690000064
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, t is the time, and x is the axial length of the pipe.
The momentum conservation model is as follows:
Figure BDA0001893038690000071
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, k is the number of fluid phases in contact with the ith fluid, and αiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, pikIs the pressure between the phase interfaces of the ith fluid phase and the k fluid phases with which it is in contact, fiIs the friction between the ith fluid phase and the pipe wall, rikFor all k fluids in contact with the ith fluid relative to the ith fluid phaseThe resultant force in the axial direction of the pipeline is marked as "+" before the direction of the acting force is opposite to the flow velocity direction of the fluid, if the direction of the acting force is the same as the flow velocity direction of the fluid, the sign is "-", theta is the inclination angle between the pipeline and the horizontal line, g is the gravity acceleration, t is the time, and x is the axial length of the pipeline.
The energy conservation model is as follows:
Figure BDA0001893038690000072
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, EiIs the i-th fluid phase internal energy, HiIs the specific total enthalpy of the ith fluid phase, theta is the angle of inclination between the pipe and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipe.
In a preferred embodiment, the transient data of the multiphase fluid may be calculated by a fluid dynamics method as shown in fig. 4, and the S300 may specifically include:
s310: the first order windward format is selected to discretize the pipe model as shown in fig. 5.
S320: and dividing the staggered grids to form a solution model.
S330: solving the solution model based on the preset fluid and pipeline conservation model, continuity conservation model, momentum conservation model and energy conservation model according to the pipeline parameters and the fluid parameters to obtain transient data of the multiphase fluid, such as pipeline pressure (p) at discrete points in the solution modelW、pP、pE) Fluid density (alpha)W、αP、αE) And phase fraction of fluid phase (p)W、ρE、ρP) And the fluid velocity (u) of the discrete sub-regionsw、ue、uee)。
FIG. 6 is a graph showing a comparison of initial, transient and steady state flow of fluid through a faucet tube in a specific embodiment of the faucet. The transient data results of the fluid shown in fig. 7-9 are obtained by the multiphase flow transient calculation method, and as can be seen from fig. 7-9, the multiphase flow transient calculation method based on space conservation and theoretical analytic solution data have high goodness of fit, which indicates that the simulation method can accurately and efficiently simulate the transient state of a multi-flow phase in a pipeline.
According to another aspect of the present invention, the present embodiment further discloses a multiphase flow transient computing system based on space conservation. As shown in fig. 10, the system includes a parameter acquisition unit 11 and a model solving unit 12.
The parameter obtaining unit 11 is configured to obtain a pipeline parameter and a fluid parameter. The pipe parameters may include pipe cross-sectional area, pipe length and other parameters, and the fluid parameters may include the number of fluid phases, the density of the fluid, the initial velocity of the fluid and other parameters.
The model solving unit 12 is configured to obtain transient data of the multiphase fluid by a fluid dynamics analysis method based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model according to the pipeline parameters and the fluid parameters.
According to the invention, through the intrinsic characteristic research of the flowing mechanism of the fluid in the pipeline, the conservation model between the pipeline volume and the fluid volume is added in the multiphase fluid model, the transient state of the flow of incompressible and compressible fluid and multiple fluids at any flow rate can be simulated, the solving accuracy and reliability of the multiphase flow transient data are high, the defect that the relation between the fluid and the pipeline space is difficult to accurately describe in the traditional multiphase flow simulation is overcome, the method can be applied to the multiphase flow calculation with a wider range, and has important significance for accurately simulating the multiphase flow problem in the chemical industry, the nuclear industry and the petroleum industry.
In a preferred embodiment, as shown in fig. 11, the system further comprises a model establishing unit 13, which can be used to establish the fluid and pipe conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model. Preferably, the model establishing unit 13 may determine a conservation relation between the expansion and contraction rate of the pipe cross-sectional area with time and the change rate of the fluid flowing through the pipe, and establish a conservation model of the fluid and the pipe based on the associated phase fraction in the multiphase fluid according to the conservation relation. Taking into account that the pipe is always filled with at least any one fluid phase in the pipe in the full space-time range of the pipe, a conservation model of the fluid and the pipe is derived for establishing the related phase fraction, wherein the conservation model can comprise two terms, wherein one term represents the fluctuation property of the compressible fluid as a function of the density of the fluid, and the other term represents the flow property of the incompressible fluid as a function of the velocity of the fluid. The invention can improve a double-fluid model, has better applicability to both multi-fluid phases and compressible and incompressible fluid phases, and has high simulation accuracy and reliability of the multi-fluid phases.
Preferably, for a pipe segment infinitesimal as shown in fig. 3, the change in pipe volume results from the expansion and contraction of the pipe cross-sectional area over time as follows:
Figure BDA0001893038690000081
wherein A is the cross-sectional area of the pipeline, Deltax is the length of the micro-element of the pipeline, and t is time.
For a fluid flow of only one fluid phase in a pipe segment infinitesimal as shown in fig. 3, the change in volume of the fluid flowing through the pipe results from the rate of change of flow into and out of the pipe element:
Figure BDA0001893038690000091
wherein A is the cross-sectional area of the pipeline, Deltax is the length of the micro-element of the pipeline, and u is the fluid velocity.
The conservation of the rate of change of fluid flow through the conduit with respect to the rate of expansion and contraction of the cross-sectional area of the conduit over time can be expressed as:
Figure BDA0001893038690000092
wherein A is the cross-sectional area of the pipeline, t is time, u is the fluid velocity, and x is the axial length of the pipeline.
Preferably, the fluid and pipe conservation model may be:
Figure BDA0001893038690000093
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiThe velocity of the ith fluid phase, A is the cross-sectional area of the conduit, t is the time, and x is the axial length of the conduit.
In a preferred embodiment, the continuous conservation model for multiple fluids is:
Figure BDA0001893038690000094
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, t is the time, and x is the axial length of the pipe.
The momentum conservation model is as follows:
Figure BDA0001893038690000095
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, k is the number of fluid phases in contact with the ith fluid, and αiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, pikK for the ith fluid phase in contact therewithPressure between the phase boundaries of the fluid phases, riIs the friction between the ith fluid phase and the pipe wall, rikThe acting force of all k fluids in contact with the ith fluid relative to the ith fluid phase in the axial direction of the pipeline is marked as "+" if the acting force direction is opposite to the flow velocity direction of the fluids, the former mark is marked as "-" if the acting force direction is the same as the flow velocity direction of the fluids, theta is the inclination angle between the pipeline and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipeline.
The energy conservation model is as follows:
Figure BDA0001893038690000101
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, EiIs the i-th fluid phase internal energy, HiIs the specific total enthalpy of the ith fluid phase, theta is the angle of inclination between the pipe and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipe.
In a preferred embodiment, the model solving unit 12 is further configured to select a first-order windward format to discretize the pipeline, divide the staggered grids to form a solving model, and solve the solving model based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model according to the pipeline parameters and the fluid parameters to obtain transient data of the multiphase fluid. For example, the pipeline pressure (p) at discrete points in the solution model may be solvedW、pP、pE) Fluid density (alpha)W、αP、αE) And phase fraction of fluid phase (p)W、ρE、ρP) And the fluid velocity (u) of the discrete sub-regionsw、ue、uee)。
The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. A typical implementation device is a computer device, which may be, for example, a personal computer, a laptop computer, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
In a typical example, the computer device specifically comprises a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method performed by the client as described above when executing the program, or the processor implementing the method performed by the server as described above when executing the program.
Referring now to FIG. 12, shown is a schematic block diagram of a computer device 600 suitable for use in implementing embodiments of the present application.
As shown in fig. 12, the computer apparatus 600 includes a Central Processing Unit (CPU)601 which can perform various appropriate works and processes according to a program stored in a Read Only Memory (ROM)602 or a program loaded from a storage section 608 into a Random Access Memory (RAM)) 603. In the RAM603, various programs and data necessary for the operation of the system 600 are also stored. The CPU601, ROM602, and RAM603 are connected to each other via a bus 604. An input/output (I/O) interface 605 is also connected to bus 604.
The following components are connected to the I/O interface 605: an input portion 606 including a keyboard, a mouse, and the like; an output section 607 including a Cathode Ray Tube (CRT), a liquid crystal feedback (LCD), and the like, and a speaker and the like; a storage section 608 including a hard disk and the like; and a communication section 609 including a network interface card such as a LAN card, a modem, or the like. The communication section 609 performs communication processing via a network such as the internet. The driver 610 is also connected to the I/O interface 606 as needed. A removable medium 611 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 610 as necessary, so that a computer program read out therefrom is mounted as necessary on the storage section 608.
In particular, according to an embodiment of the present invention, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, embodiments of the invention include a computer program product comprising a computer program tangibly embodied on a machine-readable medium, the computer program comprising program code for performing the method illustrated in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 609, and/or installed from the removable medium 611.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims (11)

1. A multiphase flow transient calculation method based on space conservation is characterized by comprising the following steps:
acquiring pipeline parameters and fluid parameters;
and obtaining transient data of the multiphase fluid by a fluid dynamics analysis method according to the pipeline parameters and the fluid parameters based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model, wherein the fluid and pipeline conservation model is a conservation model of the fluid and pipeline based on the associated phase fraction under the multiphase fluid by determining the conservation relation between the expansion and contraction rate of the cross-sectional area of the pipeline along with time and the fluid volume change rate of the fluid flowing through the pipeline and establishing the obtained multiphase fluid according to the conservation relation.
2. The method of claim 1, further comprising the step of establishing the fluid and conduit conservation model, the continuity conservation model, the momentum conservation model, and the energy conservation model.
3. The method of claim 2, wherein the conservation relation is:
Figure FDA0002744968070000011
wherein A is the cross-sectional area of the pipeline, t is time, u is the fluid velocity, and x is the axial length of the pipeline.
4. The method of any one of claims 1-3, wherein the fluid and conduit conservation model is:
Figure FDA0002744968070000012
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiThe velocity of the ith fluid phase, A is the cross-sectional area of the conduit, t is the time, and x is the axial length of the conduit.
5. The method of claim 1, wherein the obtaining transient data of the multiphase fluid by a fluid dynamics analysis method according to the pipeline parameters and the fluid parameters based on a preset fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model specifically comprises:
selecting a first-order windward format to disperse the pipeline model;
dividing the staggered grids to form a solution model;
and calculating to obtain transient data of the multiphase fluid according to the solving model, the fluid and pipeline conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model.
6. The method of claim 1, wherein the continuity conservation model is:
Figure FDA0002744968070000021
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, t is the time, and x is the axial length of the pipe.
7. The method of claim 1, wherein the momentum conservation model is:
Figure FDA0002744968070000022
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, k is the number of fluid phases in contact with the ith fluid, and αiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, pikIs the pressure between the phase interfaces of the ith fluid phase and the k fluid phases with which it is in contact, fiIs the friction between the ith fluid phase and the pipe wall, rikFor all k fluids in contact with the ith fluid phase in the axial direction of the pipe relative to the ith fluid phaseAnd applying resultant force, wherein the former sign is "+" if the direction of the acting force is opposite to the direction of the flow velocity of the fluid, the former sign is "-" if the direction of the acting force is the same as the direction of the flow velocity of the fluid, theta is the inclination angle between the pipeline and the horizontal line, g is the gravity acceleration, t is the time, and x is the axial length of the pipeline.
8. The method of claim 1, wherein the energy conservation model is:
Figure FDA0002744968070000023
wherein i is the number of fluid phases in the pipeline, i is 1,2,. n, n is a positive integer, alphaiIs the phase fraction, p, of the ith fluid phaseiIs the density of the ith fluid phase, uiIs the velocity of the ith fluid phase, p is the line pressure, EiIs the i-th fluid phase internal energy, HiIs the specific total enthalpy of the ith fluid phase, theta is the angle of inclination between the pipe and the horizontal line, g is the gravitational acceleration, t is the time, and x is the axial length of the pipe.
9. A multi-phase flow transient computing system based on space conservation is characterized by comprising a model establishing unit and a model solving unit;
the model establishing unit is used for establishing a fluid and pipeline conservation model, a continuity conservation model, a momentum conservation model and an energy conservation model based on the correlation phase fraction between the pipeline and the multiphase fluid flowing through the pipeline;
the model solving unit is used for obtaining transient data of the multiphase fluid through a fluid dynamics analysis method according to the fluid and pipeline conservation model, the continuity conservation model, the momentum conservation model and the energy conservation model, wherein the fluid and pipeline conservation model is a conservation relation which is obtained by determining the expansion and contraction rate of the cross-sectional area of the pipeline along with time and the volume change rate of the fluid flowing through the pipeline, and the conservation model of the fluid and the pipeline based on the associated phase fraction under the multiphase fluid is established according to the conservation relation.
10. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor,
the processor, when executing the program, implements the method of any of claims 1-8.
11. A computer-readable medium, having stored thereon a computer program,
the program when executed by a processor implementing the method according to any one of claims 1-8.
CN201811479567.6A 2018-12-05 2018-12-05 Multiphase flow transient calculation method and system based on space conservation Active CN109635405B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811479567.6A CN109635405B (en) 2018-12-05 2018-12-05 Multiphase flow transient calculation method and system 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 Multiphase flow transient calculation method and system based on space conservation

Publications (2)

Publication Number Publication Date
CN109635405A CN109635405A (en) 2019-04-16
CN109635405B true 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 Multiphase flow transient calculation method and system based on space conservation

Country Status (1)

Country Link
CN (1) CN109635405B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
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
CN111400950B (en) * 2020-03-17 2022-03-15 中国石油大学(北京) Hydrate slurry multiphase pipeline transient flow simulation method and device
CN113792449B (en) * 2021-08-10 2023-12-19 中国船舶重工集团公司第七一九研究所 Nuclear reactor monitoring method and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8621937B2 (en) * 2010-08-24 2014-01-07 Invensys Systems, Inc. Multiphase metering system
CN102943620B (en) * 2012-08-27 2013-08-28 中国石油大学(华东) Pressure-controlled drilling method based on drilling annulus wellbore multi-phase flow computing
CN103149012B (en) * 2013-03-08 2015-07-01 中国石油大学(北京) Experimental device and method for simulating flow characteristic of coal bed gas gathering and transportation pipeline
CN104594834B (en) * 2014-12-01 2016-12-21 中国石油大学(华东) A kind of monitoring method of deep water oil base drilling fluid drilling well flooded conditions

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Also Published As

Publication number Publication date
CN109635405A (en) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109635405B (en) Multiphase flow transient calculation method and system based on space conservation
Wilson et al. Comprehensive approach to verification and validation of CFD simulations—Part 2: Application for RANS simulation of a cargo/container ship
Vidakovic Statistical modeling by wavelets
Fleckenstein et al. Simplified modeling of the influence of surfactants on the rise of bubbles in VOF-simulations
Ramos et al. Flat vs. curved rigid-lid LES computations of an open-channel confluence
Vítkovský et al. Numerical error in weighting function-based unsteady friction models for pipe transients
Meziou et al. Low-dimensional modeling of transient two-phase flow in pipelines
Fabricius et al. A comparison of the roughness regimes in hydrodynamic lubrication
Youssef et al. Application of statistical analysis techniques to pipeline on-bottom stability analysis
Schmandt et al. Loss coefficients for periodically unsteady flows in conduit components: illustrated for laminar flow in a circular duct and a 90 degree bend
CN109826599B (en) Method and system for determining gas injection quantity in thickened oil recovery
Song et al. A surrogate model for shallow water equations solvers with deep learning
Korzilius et al. Modeling liquid slugs accelerating in inclined conduits
Zhang et al. Modeling and simulation of pigging for a gas pipeline using a bypass pig
Hadj-Taı¨ eb et al. Validation of hyperbolic model for water-hammer in deformable pipes
Dai et al. Accuracy and reliability of piecewise-constant method in studying the responses of nonlinear dynamic systems
Liou Numerical properties of the discrete gas cavity model for transients
CN110232222B (en) Sediment tube flow field analysis method and system
CN108510111A (en) A kind of method and device of prediction annular space supercharging
Chuang et al. An analytical solution of groundwater flow in wedge-shaped aquifers with estuarine boundary conditions
Pawar et al. Spatiotemporal Down-scaling for Multiphase Flow in Porous Media Using Implicit Hypernetwork Neural Representation
Jing et al. Flow behavior of two-dimensional wet foam: Effect of foam quality
Ba geri et al. Evaluating the effect of surface tension on two-phase flow in horizontal pipes
Zhang et al. Modeling of slug dissipation and generation in gas-liquid hilly-terrain pipe flow
Paternina-Verona et al. Two-dimensional simulation of emptying manoeuvres in water pipelines with admitted air

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