WO2024148671A1 - 一种气液数值模拟曳力修正方法、系统、设备及介质 - Google Patents

一种气液数值模拟曳力修正方法、系统、设备及介质 Download PDF

Info

Publication number
WO2024148671A1
WO2024148671A1 PCT/CN2023/080761 CN2023080761W WO2024148671A1 WO 2024148671 A1 WO2024148671 A1 WO 2024148671A1 CN 2023080761 W CN2023080761 W CN 2023080761W WO 2024148671 A1 WO2024148671 A1 WO 2024148671A1
Authority
WO
WIPO (PCT)
Prior art keywords
gas
phase
model
liquid
fluid dynamics
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.)
Ceased
Application number
PCT/CN2023/080761
Other languages
English (en)
French (fr)
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.)
Jiangsu University
Original Assignee
Jiangsu University
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 Jiangsu University filed Critical Jiangsu University
Priority to GB2313400.0A priority Critical patent/GB2619201B/en
Publication of WO2024148671A1 publication Critical patent/WO2024148671A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Definitions

  • the present invention relates to the technical field of multiphase flow simulation, and in particular to a method, system, equipment and medium for gas-liquid numerical simulation drag correction.
  • the ocean contains abundant energy resources such as oil and natural gas, accounting for about 34% of the total oil and gas resources, and 70% of them are in deep sea areas, with huge development potential and difficulty in exploitation.
  • energy resources such as oil and natural gas
  • it is necessary to vigorously carry out China's deep-sea oil and gas exploration and development as soon as possible.
  • Deep-sea oil and gas exploitation has become an important measure to ensure China's energy security.
  • deep-sea oil and gas transportation gas-liquid mixed transportation technology has gradually become the preferred solution for deep-sea oil and gas transportation because it greatly simplifies the mining process, saves pipeline laying costs, and improves oil well recovery. It has the advantages of high efficiency, energy saving, and environmental protection.
  • the main components of deep-sea oil and gas are a gas-liquid mixture of natural gas and petroleum. Therefore, deep-sea oil and gas transportation is essentially a scientific problem of gas-liquid two-phase flow mixed transportation under high-pressure conditions.
  • the purpose of the present invention is to provide a method, system, device and medium for correcting the drag in gas-liquid numerical simulation, which can correct the magnitude of the interphase drag according to the ambient pressure and improve the accuracy of predicting the gas holdup under high pressure environment.
  • the present invention provides the following solutions:
  • a method for correcting drag in gas-liquid numerical simulation comprising:
  • the three-dimensional pipeline model is meshed by using a polyhedral meshing method to obtain a mesh
  • the gas holdup distribution represents the change of the volume fraction of the gas phase with the radial position at the set position section of the three-dimensional pipeline model.
  • the initial computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a liquid phase momentum equation and a gas phase momentum equation;
  • liquid phase continuity equation is specifically formulated as follows:
  • the gas phase continuity equation is specifically formulated as follows:
  • ⁇ l is the volume fraction of the liquid phase
  • ⁇ g is the volume fraction of the gas phase
  • ⁇ l is the density of the liquid phase
  • ⁇ g is the density of the gas phase
  • ul is the velocity vector of the liquid phase
  • ug is the velocity vector of the gas phase
  • S l is the liquid phase source term
  • S g is the gas phase source term
  • t is time;
  • liquid phase momentum equation is specifically formulated as follows:
  • the gas phase momentum equation is specifically formulated as follows:
  • g is the gravitational acceleration vector
  • ⁇ l is the liquid phase stress tensor
  • ⁇ g is the gas phase stress tensor
  • ⁇ is the phase momentum exchange coefficient
  • F lg is the sum of other phase forces between the liquid phase and the gas phase except the drag force
  • F gl is the sum of other phase forces between the gas phase and the liquid phase except the drag force
  • F lg -F gl
  • Other interphase forces include lift, wall lubrication force, turbulent diffusion force and virtual mass force.
  • P0 is the standard atmospheric pressure
  • P is the ambient pressure
  • ⁇ is the correction function of the phase momentum exchange coefficient
  • the modified computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a modified liquid phase momentum equation and a modified gas phase momentum equation;
  • the modified liquid phase momentum equation has the following specific formula:
  • the modified gas phase momentum equation is specifically formulated as follows:
  • the three-dimensional pipeline model has a diameter of 0.15 meters and a total length of 6.6 meters.
  • the total number of grids in the three-dimensional pipeline model after grid division is 535372.
  • a gas-liquid numerical simulation drag correction system comprising:
  • Pipeline model building module used to build a three-dimensional pipeline model
  • a model meshing module used to mesh the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model
  • a computational fluid dynamics model building module is used to build an initial computational fluid dynamics model based on the basic theory of the Euler-Euler two-fluid model
  • a correction function determination module used for obtaining the ambient pressure and determining the correction function according to the ambient pressure
  • a computational fluid dynamics model correction module used to correct the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model
  • the gas-liquid numerical simulation module is used to perform unsteady simulation of gas-liquid two-phase flow in the pipe according to the modified computational fluid dynamics model and the meshed three-dimensional pipeline model to obtain the gas content distribution in the pipe; the gas content distribution represents the change of the volume fraction of the gas phase with the radial position at the set position section of the three-dimensional pipeline model.
  • An electronic device comprises a memory and a processor, wherein the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the above-mentioned gas-liquid numerical simulation drag correction method.
  • a computer-readable storage medium stores a computer program, which implements the above-mentioned gas-liquid numerical simulation drag correction method when executed by a processor.
  • the present invention discloses the following technical effects:
  • the gas-liquid numerical simulation drag correction method obtained by the present invention obtains the ambient pressure, determines the correction function according to the ambient pressure, and then corrects the initial computational fluid dynamics model according to the correction function. It can correct the magnitude of the phase drag under different pressure conditions, so that the corrected computational fluid dynamics model can more accurately predict the gas content distribution in the pipe, and especially can effectively avoid the existing numerical model's erroneous prediction of the gas content in a high-pressure environment where the pressure intensity in the flow field is in the range of 0.5Mpa to 2Mpa.
  • FIG1 is a flow chart of a method for correcting drag in gas-liquid numerical simulation provided by the present invention
  • FIG2 is a specific flow chart of the gas-liquid numerical simulation drag correction method provided by the present invention.
  • FIG3 is a schematic structural diagram of a three-dimensional pipeline model provided by the present invention.
  • FIG4 is a cross-sectional schematic diagram of a three-dimensional pipeline model after meshing provided by the present invention.
  • FIG5 is a comparison diagram of gas holdup distribution at a pressure of 1.0 MPa provided by the present invention.
  • FIG6 is a comparison diagram of gas holdup distribution at a pressure of 2.0 MPa provided by the present invention.
  • FIG. 7 is a module diagram of the gas-liquid numerical simulation drag correction system provided by the present invention.
  • Pipeline model building module 1, model grid division module—2, computational fluid dynamics model building module—3, correction function determination module—4, computational fluid dynamics model correction module—5, gas-liquid numerical simulation module—6.
  • the change in pressure has a great influence on the dynamic model of bubbles in two-phase flow, and the gas-liquid two-phase numerical model needs to be modified under high pressure.
  • the influence of pressure is mainly in two aspects: on the one hand, the increase in pressure increases the density of the gas phase, and on the other hand, the increase in pressure reduces the rising speed of bubbles, which has a great influence on the movement characteristics of bubbles in the main phase.
  • the present invention is a correction method that introduces the influence of pressure on gas-liquid two-phase flow into the drag model.
  • the purpose of the present invention is to provide a gas-liquid numerical simulation drag correction method, system, equipment and medium, which can correct the magnitude of interphase drag according to the ambient pressure and improve the accuracy of gas content prediction under high pressure environment.
  • the present invention provides a method for correcting drag in a gas-liquid numerical simulation, comprising:
  • Step S1 Establish a three-dimensional pipeline model.
  • the diameter of the three-dimensional pipeline model is 0.15 meters and the total length is 6.6 meters.
  • the structure thereof is shown in FIG3 .
  • Step S2 meshing the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model.
  • the cross section of the three-dimensional pipeline model after meshing is shown in FIG4 .
  • a polyhedral meshing method is adopted, the mesh size is about 1.2 times of the bubble diameter, and the total number of meshes is 535,372.
  • Step S3 Based on the basic theory of the Euler-Euler two-fluid model, establish an initial computational fluid dynamics (CFD) model.
  • CFD computational fluid dynamics
  • the initial computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a liquid phase momentum equation and a gas phase momentum equation.
  • liquid phase continuity equation is specifically formulated as follows:
  • the gas phase continuity equation is specifically formulated as follows:
  • ⁇ l is the volume fraction of the liquid phase
  • ⁇ g is the volume fraction of the gas phase
  • ⁇ l is the density of the liquid phase
  • ⁇ g is the density of the gas phase
  • ul is the velocity vector of the liquid phase
  • ug is the velocity vector of the gas phase
  • S l is the liquid phase source term
  • S g is the gas phase source term, and they are generally 0
  • t is time
  • liquid phase momentum equation is specifically formulated as follows:
  • the gas phase momentum equation is specifically formulated as follows:
  • g is the gravitational acceleration vector
  • ⁇ l is the liquid phase stress tensor
  • ⁇ g is the gas phase stress tensor
  • ⁇ is the phase momentum exchange coefficient, which determines the magnitude of the drag force
  • F lg is the relative velocity of the liquid phase to the gas phase in addition to the drag force.
  • the other interphase forces include lift, wall lubrication force, turbulent diffusion force, virtual mass force, etc.
  • Step S4 Acquire the ambient pressure, and determine a correction function according to the ambient pressure.
  • the gas holdup i.e., the volume fraction of the gas phase
  • ⁇ g and the pressure P show a certain functional relationship; based on this relationship, a fitting algorithm is used to fit the pressure correction function ⁇ .
  • the specific correction method divides the influence of pressure into two intervals. When the ambient pressure is less than 1MPa, the correction factor adopts a linear formula; when the ambient pressure is greater than 1MPa, the correction factor adopts a power function formula.
  • P0 is the standard atmospheric pressure
  • P is the ambient pressure
  • ⁇ is the correction function of the phase momentum exchange coefficient
  • Step S5 Correcting the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model.
  • the modified computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a modified liquid phase momentum equation and a modified gas phase momentum equation.
  • the modified liquid phase momentum equation has the following specific formula:
  • the modified gas phase momentum equation is specifically formulated as follows:
  • Step S6 Based on the modified computational fluid dynamics model and the meshed three-dimensional pipeline model, an unsteady simulation of the gas-liquid two-phase flow in the pipe is performed to obtain the gas content distribution in the pipe; the gas content distribution represents the change of the volume fraction of the gas phase with the radial position at the set position section of the three-dimensional pipeline model.
  • Figures 5 and 6 show the average gas holdup distribution at pressures of 1.0 MPa and 2.0 MPa, respectively. It can be found from Figures 5 and 6 that the gas holdup distribution predicted by the correction method proposed in the present invention is more consistent with the experimentally measured gas holdup distribution, that is, the present invention can more accurately predict the gas holdup distribution under high pressure conditions in the pipe.
  • a gas-liquid numerical simulation drag correction system including:
  • the pipeline model building module 1 is used to build a three-dimensional pipeline model.
  • the model meshing module 2 is used to mesh the three-dimensional pipeline model using a polyhedral meshing method to obtain a meshed three-dimensional pipeline model.
  • the computational fluid dynamics model building module 3 is used to build an initial computational fluid dynamics model based on the basic theory of the Euler-Euler two-fluid model.
  • the correction function determination module 4 is used to obtain the ambient pressure and determine the correction function according to the ambient pressure.
  • the computational fluid dynamics model correction module 5 is used to correct the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model.
  • the gas-liquid numerical simulation module 6 is used to perform unsteady simulation of the gas-liquid two-phase flow in the pipe according to the modified computational fluid dynamics model and the three-dimensional pipeline model after meshing, and obtain the gas content distribution in the pipe; the gas content distribution represents the change of the volume fraction of the gas phase with the radial position at the set position section of the three-dimensional pipeline model.
  • An embodiment of the present invention further provides an electronic device, comprising a memory and a processor, wherein the memory is used to store a computer program, and the processor is used to run the computer program to enable the electronic device to execute the gas-liquid numerical simulation drag correction method in embodiment 1.
  • the electronic device may be a server.
  • the present invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the gas-liquid numerical simulation drag correction method in the first embodiment.
  • the present invention specifically discloses a pressure correction method for gas phase drag based on the Euler-Euler two-fluid model theory.
  • the environment is a high-pressure environment, and the gas phase has different drags under different pressures, resulting in different phase content distributions.
  • the correction method of the present invention can be used to reasonably predict the gas content.
  • the gas-liquid numerical simulation drag correction method provided by the present invention is suitable for high-pressure conditions. Based on the basic theory of the Euler-Euler two-fluid model, according to the relationship between the drag between the gas and liquid phases and the environmental pressure, the correction relationship affected by pressure is given by data fitting, which can more reasonably predict the gas content distribution in the flow basin under different pressures.
  • phase drag shows different trends under different pressures.
  • This correction is designed for high-pressure environments and gas-liquid two-phase simulations with a large pressure change scale.
  • the magnitude of the correction drag is determined according to the intensity of the pressure in the flow field, which can effectively avoid the existing numerical model from making incorrect predictions of the gas content under high-pressure environments.
  • each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments.
  • the same or similar parts between the embodiments can be referred to each other.
  • the description is relatively simple, and the relevant parts can be referred to the method part.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开一种气液数值模拟曳力修正方法、系统、设备及介质，涉及多相流模拟技术领域，该方法包括：建立三维管道模型；采用多面体网格划分方式，对三维管道模型进行网格划分，得到网格划分后的三维管道模型；基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型；获取环境压力，并根据环境压力确定修正函数；根据修正函数对初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型；根据修正后的计算流体动力学模型和网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布。本发明能够根据环境压力修正相间曳力的大小，提高对高压环境下气含率预测的准确性。

Description

一种气液数值模拟曳力修正方法、系统、设备及介质 技术领域
本发明涉及多相流模拟技术领域，特别是涉及一种气液数值模拟曳力修正方法、系统、设备及介质。
背景技术
海洋中蕴藏着丰富的石油、天然气等能源资源，约占油气资源总量的34％，而其中70％又都处于深海区域，开发潜力以及开采难度都十分巨大。为了解决中国油气对外依存度过高的问题，需要尽快大力开展中国的深海油气勘探与开发，深海油气开采已成为保障中国能源安全的重要举措。
现如今，在深海油气输运中，由于气液混输技术大幅简化了开采工艺流程，节省了管道铺设成本、提高了油井采收率，具有高效、节能、环保的优点，已逐渐成为深海油气输运的首选方案。深海油气的主要成分为天然气和石油的气液混合物，因此深海油气输运本质上是一种高压条件下气液两相流混输的科学问题。
压力的变化对于两相流中气泡的动力学模型有着很大的影响，在高压环境下需要对气液两相数值模型进行修正。然而现有技术尚且缺乏相关研究，现有数值模型在对高压环境下的气含率预测中，存在准确率不高的问题。
发明内容
本发明的目的是提供一种气液数值模拟曳力修正方法、系统、设备及介质，能够根据环境压力修正相间曳力的大小，提高对高压环境下气含率预测的准确性。
为实现上述目的，本发明提供了如下方案：
一种气液数值模拟曳力修正方法，包括：
建立三维管道模型；
采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网 格划分后的三维管道模型；
基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型；
获取环境压力，并根据所述环境压力确定修正函数；
根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型；
根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
可选地，所述初始计算流体动力学模型包括：液相连续性方程、气相连续性方程、液相动量方程和气相动量方程；
所述液相连续性方程，具体公式为：
所述气相连续性方程，具体公式为：
其中，αl为液相的体积分数，αg为气相的体积分数；ρl为液相的密度，ρg为气相的密度；ul为液相的速度矢量，ug为气相的速度矢量；Sl为液相源项，Sg为气相源项；t为时间；
所述液相动量方程，具体公式为：
所述气相动量方程，具体公式为：
其中，g为重力加速度矢量；τl为液相应力张量，τg为气相应力张量；β为相间动量交换系数；Flg为除曳力之外，液相相对于气相的其他相间力之和，Fgl为除曳力之外，气相相对于液相的其他相间力之和，且Flg＝-Fgl，所述其 他相间力包括升力、壁面润滑力、湍流扩散力和虚拟质量力。
可选地，所述修正函数的具体公式为：
其中，P0为标准大气压，P为环境压力，ξ为相间动量交换系数的修正函数。
可选地，所述修正后的计算流体动力学模型包括：液相连续性方程、气相连续性方程、修正后的液相动量方程和修正后的气相动量方程；
所述修正后的液相动量方程，具体公式为：
所述修正后的气相动量方程，具体公式为：
可选地，所述三维管道模型的管径为0.15米，总长为6.6米。
可选地，所述网格划分后的三维管道模型中的总网格数为535372个。
一种气液数值模拟曳力修正系统，包括：
管道模型建立模块，用于建立三维管道模型；
模型网格划分模块，用于采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网格划分后的三维管道模型；
计算流体动力学模型建立模块，用于基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型；
修正函数确定模块，用于获取环境压力，并根据所述环境压力确定修正函数；
计算流体动力学模型修正模块，用于根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型；
气液数值模拟模块，用于根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
一种电子设备，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行上述的气液数值模拟曳力修正方法。
一种计算机可读存储介质，其存储有计算机程序，所述计算机程序被处理器执行时实现上述的气液数值模拟曳力修正方法。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供的气液数值模拟曳力修正方法，通过获取环境压力，并根据环境压力确定修正函数，进而根据修正函数对初始计算流体动力学模型进行修正，能够对不同压力条件下的相间曳力的大小进行修正，从而使修正后的计算流体动力学模型能够更加精确地预测管内的气含率分布，特别是能够有效避免现有数值模型对于流场中压力的强度在0.5Mpa～2Mpa范围内的高压环境下的气含率的错误预测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的气液数值模拟曳力修正方法的流程图；
图2为本发明提供的气液数值模拟曳力修正方法的具体流程图；
图3为本发明提供的三维管道模型的结构示意图；
图4为本发明提供的网格划分后的三维管道模型的截面示意图；
图5为本发明提供的1.0MPa的压力下气含率分布对比图；
图6为本发明提供的2.0MPa的压力下气含率分布对比图；
图7为本发明提供的气液数值模拟曳力修正系统的模块图。
符号说明：
管道模型建立模块—1，模型网格划分模块—2，计算流体动力学模型建
立模块—3，修正函数确定模块—4，计算流体动力学模型修正模块—5，气液数值模拟模块—6。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
压力的变化对于两相流中气泡的动力学模型有着很大的影响，在高压环境下需要对气液两相数值模型进行修正。Krishna首先提出压力的影响主要在两个方面：一方面压力的增加增大了气相密度，另一方面压力的增加减小了气泡上升速度，对气泡在主相中的运动特征产生较大影响。
因此，建立考虑气泡在压力影响下的修正模型并且能够统一描述压力对气相动力学行为的影响机理，这对于深海油气开发中涉及的高压环境气液两相数值模拟技术的发展将有重要的意义。本发明是一种将压力对气液两相流的影响引入到曳力模型中的修正方法。
具体地，本发明的目的是提供一种气液数值模拟曳力修正方法、系统、设备及介质，能够根据环境压力修正相间曳力的大小，提高对高压环境下气含率预测的准确性。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
如图1及图2所示，本发明提供一种气液数值模拟曳力修正方法，包括：
步骤S1：建立三维管道模型。
在本实施例中，所述三维管道模型的管径为0.15米，总长为6.6米，其结构参见图3。
步骤S2：采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网格划分后的三维管道模型。
在本实施例中，网格划分后的三维管道模型的截面如图4所示，采用多面体网格划分方式，网格大小约为气泡直径的1.2倍，总网格数为535372。
步骤S3：基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学(Computational Fluid Dynamics，CFD)模型。
具体地，所述初始计算流体动力学模型包括：液相连续性方程、气相连续性方程、液相动量方程和气相动量方程。
所述液相连续性方程，具体公式为：
所述气相连续性方程，具体公式为：
其中，下标l和g分别表示液相和气相，αl为液相的体积分数，αg为气相的体积分数；ρl为液相的密度，ρg为气相的密度；ul为液相的速度矢量，ug为气相的速度矢量；Sl为液相源项，Sg为气相源项，一般均为0；t为时间；表示哈密顿算子。
所述液相动量方程，具体公式为：
所述气相动量方程，具体公式为：
其中，g为重力加速度矢量；τl为液相应力张量，τg为气相应力张量；β为相间动量交换系数，决定着曳力的大小；Flg为除曳力之外，液相相对于气 相的其他相间力之和，Fgl为除曳力之外，气相相对于液相的其他相间力之和，且Flg＝-Fgl，所述其他相间力包括升力、壁面润滑力、湍流扩散力和虚拟质量力等。
具体地，气相与液相存在的相互作用力有多种，包括曳力、升力、壁面润滑力、湍流扩散力和虚拟质量力等，其中最主要的是曳力项β(ul-ug)，而β为作为相间动量交换系数，在本发明中，需要进行相应修正。
步骤S4：获取环境压力，并根据所述环境压力确定修正函数。
在不同压力下，模拟中气含率(即气相的体积分数)αg趋于稳定后，αg与压力P呈现一定的函数关系；根据其关系，采用拟合算法拟合出压力的修正函数ξ，具体的修正方法将压力的影响分为两个区间，当环境压力小于1MPa时，修正因子采用线性式；当环境压力大于1MPa时，修正因子采用幂函数式。
所述修正函数的具体公式为：
其中，P0为标准大气压，P为环境压力，ξ为相间动量交换系数的修正函数。
步骤S5：根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型。
具体地，所述修正后的计算流体动力学模型包括：液相连续性方程、气相连续性方程、修正后的液相动量方程和修正后的气相动量方程。
所述修正后的液相动量方程，具体公式为：
所述修正后的气相动量方程，具体公式为：
步骤S6：根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
具体地，进行管内气液两相流的非定常模拟，将修正过后的曳力模型编译到其相间力模型中，利用修正函数对不同压力条件下的气泡曳力进行修正预测，最终取模拟的时均值作为评价标准。图5、图6分别给出了压力在1.0MPa和2.0MPa下的平均气含率分布图，由图5及图6可以发现，采用本发明所提出的修正方法所预测得到的气含率分布，与实验测得的气含率分布情况更加吻合，即本发明可以更加精确的预测管内高压条件下的气含率分布。
实施例二
如图7所示，为了执行上述实施例一对应的方法，以实现相应的功能和技术效果，下面提供一种气液数值模拟曳力修正系统，包括：
管道模型建立模块1，用于建立三维管道模型。
模型网格划分模块2，用于采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网格划分后的三维管道模型。
计算流体动力学模型建立模块3，用于基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型。
修正函数确定模块4，用于获取环境压力，并根据所述环境压力确定修正函数。
计算流体动力学模型修正模块5，用于根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型。
气液数值模拟模块6，用于根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
实施例三
本发明实施例还提供一种电子设备，包括存储器和处理器，该存储器用于存储计算机程序，该处理器用于运行计算机程序以使电子设备执行实施例一中的气液数值模拟曳力修正方法。所述电子设备可以是服务器。
另外，本发明还提供一种计算机可读存储介质，其存储有计算机程序，该计算机程序被处理器执行时实现实施例一中的气液数值模拟曳力修正方法。
综上所述，本发明具体公开了一种基于欧拉-欧拉双流体模型理论气相曳力的压力修正方法，在深海油气混输过程中，所处的环境即为高压环境，并且，不同的压力下气相具有不同的曳力，导致不同的相含率分布，则可利用本发明的修正方法对气含率进行合理预测。本发明提供的适用于高压条件下的气液数值模拟曳力修正方法，基于欧拉-欧拉双流体模型基本理论，根据气液相间曳力随环境压力的变化关系，通过数据拟合方式给出了受压力影响的修正关系，可以更加合理的预测不同压力下的流域内气含率分布。针对不同压力范围采取不同的修正系数，在不同的压力下相间曳力呈现不同的趋势，这一修正是针对高压环境以及具有较大的压力变化尺度的气液两相模拟而设计的，根据流场中压力的强度大小来确定修正曳力的大小，可以有效避免现有数值模型对于高压环境下对气含率的错误预测。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (9)

  1. 一种气液数值模拟曳力修正方法，其特征在于，包括：
    建立三维管道模型；
    采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网格划分后的三维管道模型；
    基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型；
    获取环境压力，并根据所述环境压力确定修正函数；
    根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型；
    根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
  2. 根据权利要求1所述的气液数值模拟曳力修正方法，其特征在于，所述初始计算流体动力学模型包括：液相连续性方程、气相连续性方程、液相动量方程和气相动量方程；
    所述液相连续性方程，具体公式为：
    所述气相连续性方程，具体公式为：
    其中，αl为液相的体积分数，αg为气相的体积分数；ρl为液相的密度，ρg为气相的密度；ul为液相的速度矢量，ug为气相的速度矢量；Sl为液相源项，Sg为气相源项；t为时间；
    所述液相动量方程，具体公式为：
    所述气相动量方程，具体公式为：
    其中，g为重力加速度矢量；τl为液相应力张量，τg为气相应力张量；β为相间动量交换系数；Flg为除曳力之外，液相相对于气相的其他相间力之和，Fgl为除曳力之外，气相相对于液相的其他相间力之和，且Flg＝-Fgl，所述其他相间力包括升力、壁面润滑力、湍流扩散力和虚拟质量力。
  3. 根据权利要求2所述的气液数值模拟曳力修正方法，其特征在于，所述修正函数的具体公式为：
    其中，P0为标准大气压，P为环境压力，ξ为相间动量交换系数的修正函数。
  4. 根据权利要求3所述的气液数值模拟曳力修正方法，其特征在于，所述修正后的计算流体动力学模型包括：液相连续性方程、气相连续性方程、修正后的液相动量方程和修正后的气相动量方程；
    所述修正后的液相动量方程，具体公式为：
    所述修正后的气相动量方程，具体公式为：
  5. 根据权利要求1所述的气液数值模拟曳力修正方法，其特征在于，所述三维管道模型的管径为0.15米，总长为6.6米。
  6. 根据权利要求5所述的气液数值模拟曳力修正方法，其特征在于，所述网格划分后的三维管道模型中的总网格数为535372个。
  7. 一种气液数值模拟曳力修正系统，其特征在于，包括：
    管道模型建立模块，用于建立三维管道模型；
    模型网格划分模块，用于采用多面体网格划分方式，对所述三维管道模型进行网格划分，得到网格划分后的三维管道模型；
    计算流体动力学模型建立模块，用于基于欧拉-欧拉双流体模型基本理论，建立初始计算流体动力学模型；
    修正函数确定模块，用于获取环境压力，并根据所述环境压力确定修正函数；
    计算流体动力学模型修正模块，用于根据所述修正函数对所述初始计算流体动力学模型进行修正，得到修正后的计算流体动力学模型；
    气液数值模拟模块，用于根据所述修正后的计算流体动力学模型和所述网格划分后的三维管道模型，进行管内气液两相流的非定常模拟，得到管内的气含率分布；所述气含率分布表示在所述三维管道模型的设定位置截面处，气相的体积分数随径向位置的变化情况。
  8. 一种电子设备，其特征在于，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行如权利要求1至6中任意一项所述的气液数值模拟曳力修正方法。
  9. 一种计算机可读存储介质，其特征在于，其存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至6中任意一项所述的气液数值模拟曳力修正方法。
PCT/CN2023/080761 2023-01-13 2023-03-10 一种气液数值模拟曳力修正方法、系统、设备及介质 Ceased WO2024148671A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
GB2313400.0A GB2619201B (en) 2023-01-13 2023-03-10 A gas-liquid numerical simulation drag force correction method and system, device, and medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202310041937.2A CN116306342B (zh) 2023-01-13 2023-01-13 一种气液数值模拟曳力修正方法、系统、设备及介质
CN202310041937.2 2023-01-13

Publications (1)

Publication Number Publication Date
WO2024148671A1 true WO2024148671A1 (zh) 2024-07-18

Family

ID=86836753

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/080761 Ceased WO2024148671A1 (zh) 2023-01-13 2023-03-10 一种气液数值模拟曳力修正方法、系统、设备及介质

Country Status (2)

Country Link
CN (1) CN116306342B (zh)
WO (1) WO2024148671A1 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119044866A (zh) * 2024-09-24 2024-11-29 无锡学院 一种大气电场监测结果的修订方法
CN119249971A (zh) * 2024-12-04 2025-01-03 天目山实验室 一种激光空化气泡的数值仿真方法、装置、设备及介质
CN119249972A (zh) * 2024-12-04 2025-01-03 天目山实验室 一种气液冰三相流动分析方法、装置、设备及介质
CN119323165A (zh) * 2024-12-18 2025-01-17 浙江理工大学 一种多尺度颗粒流动计算方法
CN119623355A (zh) * 2025-02-12 2025-03-14 哈尔滨工业大学（威海） 基于虚拟流体与弹性体耦合作用的旋流器结构优化方法
CN119849375A (zh) * 2025-02-21 2025-04-18 哈尔滨工程大学 一种适应于核能系统一体化汽水分离装置性能分析的耦合求解方法、程序、设备及存储介质
CN120337829A (zh) * 2025-06-20 2025-07-18 浙江理工大学 基于多参考系旋转流道内固液两相数值模拟方法
CN120524869A (zh) * 2025-07-23 2025-08-22 上正阀门集团有限公司 一种可优化防堵塞调节截止阀结构的流场分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120053908A1 (en) * 2010-08-27 2012-03-01 Doyub Kim Method for simulating dispersed bubble flow
CN106682348A (zh) * 2017-01-09 2017-05-17 福州大学 采用低雷诺数湍流模型计算筛板萃取塔液液流场的方法
CN107132156A (zh) * 2017-05-05 2017-09-05 西安石油大学 一种颗粒密度和粒径动态变化流化床的模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115310384A (zh) * 2022-08-24 2022-11-08 清华大学 一种基于阻力系数修正的固液两相流数值计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120053908A1 (en) * 2010-08-27 2012-03-01 Doyub Kim Method for simulating dispersed bubble flow
CN106682348A (zh) * 2017-01-09 2017-05-17 福州大学 采用低雷诺数湍流模型计算筛板萃取塔液液流场的方法
CN107132156A (zh) * 2017-05-05 2017-09-05 西安石油大学 一种颗粒密度和粒径动态变化流化床的模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI GUANG, YANG XIAOGANG, DAI GANCE: "CFD simulation of gas-liquid flow in bubble column", CIESC JOURNAL, HUAXUE GONGYE CHUBANSHE, CN, vol. 59, no. 8, 15 August 2008 (2008-08-15), CN , pages 1958 - 1965, XP093190214, ISSN: 0438-1157 *
LI MENG; LI XIANG-YANG; WANG HONG-ZHI; XIE YONG-BING; CAO HONG-BIN: "Numerical Simulation of Gas-Liquid Two-phase Flow in a Bubble Column with Various Drag Models", THE CHINESE JOURNAL OF PROCESS ENGINEERING, KEXUE CHUBANSHE, BEIJING, CN, vol. 15, no. 02, 15 April 2015 (2015-04-15), CN , pages 181 - 189, XP009556039, ISSN: 1009-606X *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119044866A (zh) * 2024-09-24 2024-11-29 无锡学院 一种大气电场监测结果的修订方法
CN119249971A (zh) * 2024-12-04 2025-01-03 天目山实验室 一种激光空化气泡的数值仿真方法、装置、设备及介质
CN119249972A (zh) * 2024-12-04 2025-01-03 天目山实验室 一种气液冰三相流动分析方法、装置、设备及介质
CN119323165A (zh) * 2024-12-18 2025-01-17 浙江理工大学 一种多尺度颗粒流动计算方法
CN119623355A (zh) * 2025-02-12 2025-03-14 哈尔滨工业大学（威海） 基于虚拟流体与弹性体耦合作用的旋流器结构优化方法
CN119849375A (zh) * 2025-02-21 2025-04-18 哈尔滨工程大学 一种适应于核能系统一体化汽水分离装置性能分析的耦合求解方法、程序、设备及存储介质
CN120337829A (zh) * 2025-06-20 2025-07-18 浙江理工大学 基于多参考系旋转流道内固液两相数值模拟方法
CN120524869A (zh) * 2025-07-23 2025-08-22 上正阀门集团有限公司 一种可优化防堵塞调节截止阀结构的流场分析方法

Also Published As

Publication number Publication date
CN116306342B (zh) 2025-09-30
CN116306342A (zh) 2023-06-23

Similar Documents

Publication Publication Date Title
WO2024148671A1 (zh) 一种气液数值模拟曳力修正方法、系统、设备及介质
Aizinger et al. A discontinuous Galerkin method for two-dimensional flow and transport in shallow water
Zhao et al. Numerical modeling of local scour below a piggyback pipeline in currents
CN105912753B (zh) 基于强度折减法的海底边坡三维稳定性分析方法
Gong et al. Application and prospects of multi-phase pipeline simulation technology in empowering the intelligent oil and gas fields
Ortega et al. On the dynamic response of flexible risers caused by internal slug flow
CN105955928A (zh) 基于cfd预报船舶阻力的计算方法
US20240311535A1 (en) Hydrodynamic free-surface lattice boltzmann simulation method and system, and storage medium
CN105138731A (zh) 一种水合物分解引起海底斜坡不稳定性评价系统及方法
Zhang et al. Progressive bridge collapse analysis under both scour and floods by coupling simulation in structural and hydraulic fields. Part Ⅰ: Numerical solver
Seo et al. Numerical analyses on the formation, propagation, and deformation of landslide tsunami using LS-DYNA and NWT
CN113111608B (zh) 一种新型局部湍流脉动生成方法
Soydan et al. An improved direct forcing immersed boundary method with integrated mooring algorithm for floating offshore wind turbines
CN119227580B (zh) 一种采矿规划方法、系统、存储介质和设备
CN120951848A (zh) 预测多孔介质流动与变形的模拟方法、装置、设备及介质
Hou Various remeshing arrangements for two-dimensional finite element crack closure analysis
Pal et al. An overview of the numerical approaches to water hammer modelling: the ongoing quest for practi‑cal and accurate numerical approaches. Water. 2021; 13 (1): 1597
Caron et al. Sensitivity analysis of finite volume simulations of a breaking dam problem
CN118313309A (zh) 一种飞行器的数值仿真方法、系统、存储介质和电子设备
GB2619201A (en) A gas-liquid numerical simulation drag force correction method and system
CN111241682A (zh) 管道截面流动预测分析方法和装置
Zhang et al. Numerical analysis on the effects of a submerged bottom-mounted barrier in the head-on collision of two solitary waves
CN116482758B (zh) 流固边界耦合介质弹性波正演横波奇异性处理方法及装置
Amouzgar et al. A GPU-accelerated shallow flow model for tsunami simulations
Nandi et al. A solution verification study for URANS simulations of flow over a 5: 1 rectangular cylinder using grid convergence index and least squares procedures

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 202313400

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20230310

WWE Wipo information: entry into national phase

Ref document number: 2313400.0

Country of ref document: GB

WWP Wipo information: published in national office

Ref document number: 2313400.0

Country of ref document: GB

WWG Wipo information: grant in national office

Ref document number: 2313400.0

Country of ref document: GB

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

Ref document number: 23915437

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 23915437

Country of ref document: EP

Kind code of ref document: A1