A method, system, device and medium for drag correction of gas-liquid numerical simulation

Abstract

The present invention discloses a method, system, device, and medium for correcting drag in a gas-liquid numerical simulation, relating to the technical field of multiphase flow simulation. The method comprises: establishing a three-dimensional pipeline model; meshing the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model; establishing an initial computational fluid dynamics model based on the basic theory of the Euler-Euler two-fluid model; obtaining ambient pressure and determining a correction function based on the ambient pressure; correcting the initial computational fluid dynamics model based on the correction function to obtain a corrected computational fluid dynamics model; and performing an unsteady simulation of gas-liquid two-phase flow in a pipeline based on the corrected computational fluid dynamics model and the meshed three-dimensional pipeline model to obtain a gas holdup distribution in the pipeline. The present invention can correct the magnitude of interphase drag according to ambient pressure, thereby improving the accuracy of gas holdup prediction in high-pressure environments.

Description

Technical Field

[0001] The present invention relates to the technical field of multiphase flow simulation, and in particular to a method, system, equipment and medium for correcting drag in gas-liquid numerical simulation. Background Art

[0002] The oceans are rich in energy resources, including oil and natural gas, accounting for approximately 34% of total oil and gas resources. 70% of these resources are located in deepwater areas, presenting enormous development potential and significant extraction challenges. To address China's excessive dependence on foreign oil and gas, rapid and vigorous deepwater oil and gas exploration and development are crucial. Deepwater oil and gas extraction has become a crucial measure for safeguarding China's energy security.

[0003] Today, gas-liquid mixed transport technology has become the preferred option for deep-sea oil and gas transportation, significantly simplifying the extraction process, saving pipeline costs, and increasing oil well recovery rates. Its high efficiency, energy conservation, and environmental protection advantages have made it a promising option. Deep-sea oil and gas primarily consists of a gas-liquid mixture of natural gas and petroleum. Therefore, deep-sea oil and gas transportation is essentially a scientific problem involving the mixed transport of gas-liquid two-phase flows under high-pressure conditions.

[0004] Pressure changes significantly impact the dynamics of bubbles in two-phase flow, necessitating modifications to the gas-liquid two-phase numerical model under high-pressure conditions. However, current research is lacking, and existing numerical models suffer from inaccurate predictions of gas holdup under high-pressure conditions. Summary of the Invention

[0005] The purpose of the present invention is to provide a method, system, device and medium for correcting the drag force in gas-liquid numerical simulation, which can correct the magnitude of the interphase drag force according to the ambient pressure and improve the accuracy of gas holdup prediction under high pressure environment.

[0006] To achieve the above object, the present invention provides the following solutions:

[0007] A method for correcting drag in gas-liquid numerical simulation, comprising:

[0008] Establish a 3D pipeline model;

[0009] Using a polyhedron meshing method to mesh the three-dimensional pipeline model to obtain a meshed three-dimensional pipeline model;

[0010] Based on the basic theory of Euler-Euler two-fluid model, an initial computational fluid dynamics model was established;

[0011] Acquiring an ambient pressure, and determining a correction function according to the ambient pressure;

[0012] Correcting the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model;

[0013] Based on the modified computational fluid dynamics model and the meshed three-dimensional pipeline model, an unsteady simulation of gas-liquid two-phase flow in the pipeline is performed to obtain a gas holdup distribution in the pipeline; the gas holdup distribution represents the change in the volume fraction of the gas phase with radial position at a set position section of the three-dimensional pipeline model.

[0014] Optionally, 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;

[0015] The liquid phase continuity equation is specifically formulated as follows:

[0016]

[0017] The gas phase continuity equation is specifically formulated as follows:

[0018]

[0019] Among them, α 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; u l is the velocity vector of the liquid phase, u g 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 the time;

[0020] The liquid phase momentum equation is specifically formulated as follows:

[0021]

[0022] The gas phase momentum equation is specifically formulated as follows:

[0023]

[0024] Where 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 interphase forces between the liquid phase and the gas phase except the drag force, F gl is the sum of other interphase forces between the gas phase and the liquid phase except the drag force, and F lg =-F gl, the other interphase forces include lift, wall lubrication force, turbulent diffusion force and virtual mass force.

[0025] Optionally, the specific formula of the correction function is:

[0026]

[0027] Where P0 is the standard atmospheric pressure, P is the ambient pressure, and ξ is the correction function of the interphase momentum exchange coefficient.

[0028] Optionally, the revised computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a revised liquid phase momentum equation, and a revised gas phase momentum equation;

[0029] The modified liquid phase momentum equation is specifically formulated as follows:

[0030]

[0031] The modified gas phase momentum equation is specifically formulated as follows:

[0032]

[0033] Optionally, the three-dimensional pipeline model has a diameter of 0.15 meters and a total length of 6.6 meters.

[0034] Optionally, the total number of grids in the three-dimensional pipeline model after grid division is 535,372.

[0035] A gas-liquid numerical simulation drag correction system, comprising:

[0036] Pipeline model building module, used to build a three-dimensional pipeline model;

[0037] A model meshing module is used to mesh the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model;

[0038] Computational fluid dynamics model building module, used to build the initial computational fluid dynamics model based on the basic theory of Euler-Euler two-fluid model;

[0039] a correction function determination module, configured to obtain an ambient pressure and determine a correction function according to the ambient pressure;

[0040] A computational fluid dynamics model correction module, configured to correct the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model;

[0041] The gas-liquid numerical simulation module is used to perform unsteady simulation of gas-liquid two-phase flow in the pipe based on the modified computational fluid dynamics model and the meshed three-dimensional pipe model to obtain the gas holdup distribution in the pipe; the gas holdup distribution represents the change in the volume fraction of the gas phase with radial position at a set position section of the three-dimensional pipe model.

[0042] An electronic device includes 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 perform the above-mentioned gas-liquid numerical simulation drag correction method.

[0043] 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.

[0044] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0045] The gas-liquid numerical simulation drag correction method provided by the present invention obtains the ambient pressure, determines a correction function based on the ambient pressure, and then corrects the initial computational fluid dynamics model based on 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 holdup distribution in the pipe, and in particular, can effectively avoid the existing numerical model's erroneous prediction of the gas holdup in a high-pressure environment where the pressure intensity in the flow field is in the range of 0.5 MPa to 2 MPa. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0047] Figure 1 A flow chart of the drag correction method for gas-liquid numerical simulation provided by the present invention;

[0048] Figure 2 A specific flow chart of the gas-liquid numerical simulation drag correction method provided by the present invention;

[0049] Figure 3 A schematic structural diagram of a three-dimensional pipeline model provided by the present invention;

[0050] Figure 4 A schematic cross-sectional view of a three-dimensional pipeline model after meshing provided by the present invention;

[0051] Figure 5 A comparison diagram of gas holdup distribution at a pressure of 1.0 MPa provided by the present invention;

[0052] Figure 6 A comparison diagram of gas holdup distribution at a pressure of 2.0 MPa provided by the present invention;

[0053] Figure 7 This is a module diagram of the gas-liquid numerical simulation drag correction system provided by the present invention.

[0054] Explanation of symbols:

[0055] Pipeline model establishment module—1, model grid division module—2, computational fluid dynamics model establishment module—3, correction function determination module—4, computational fluid dynamics model correction module—5, gas-liquid numerical simulation module—6. DETAILED DESCRIPTION

[0056] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0057] Pressure changes significantly influence the dynamics of bubbles in two-phase flows, necessitating modifications to the gas-liquid two-phase numerical model under high-pressure conditions. Krishna first proposed that pressure influences two main aspects: on the one hand, an increase in pressure increases the density of the gas phase; on the other hand, an increase in pressure reduces the bubble's rise velocity, significantly affecting the bubble's motion characteristics in the main phase.

[0058] Therefore, establishing a modified model that considers the influence of bubbles under pressure and can uniformly describe the mechanism by which pressure affects gas-phase dynamic behavior will be of great significance for the development of numerical simulation technology for gas-liquid two-phase flows in high-pressure environments involved in deep-sea oil and gas development. This invention is a modified method that incorporates the influence of pressure on gas-liquid two-phase flow into the drag model.

[0059] Specifically, the purpose of the present invention is to provide a method, system, device and medium for correcting the drag force in gas-liquid numerical simulation, which can correct the magnitude of the interphase drag force according to the ambient pressure and improve the accuracy of gas holdup prediction under high-pressure environment.

[0060] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0061] Example 1

[0062] like Figure 1 and Figure 2 As shown, the present invention provides a method for correcting drag in a gas-liquid numerical simulation, comprising:

[0063] Step S1: Create a three-dimensional pipeline model.

[0064] In this embodiment, the diameter of the three-dimensional pipeline model is 0.15 meters and the total length is 6.6 meters. Figure 3 .

[0065] Step S2: meshing the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model.

[0066] In this embodiment, the cross section of the three-dimensional pipeline model after meshing is as follows: Figure 4 As shown in the figure, the polyhedron grid division method is adopted, the grid size is about 1.2 times the bubble diameter, and the total number of grids is 535372.

[0067] Step S3: Based on the basic theory of the Euler-Euler two-fluid model, an initial computational fluid dynamics (CFD) model is established.

[0068] Specifically, 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.

[0069] The liquid phase continuity equation is specifically formulated as follows:

[0070]

[0071] The gas phase continuity equation is specifically formulated as follows:

[0072]

[0073] Wherein, the subscripts l and g represent the liquid phase and the gas phase respectively, α 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; u l is the velocity vector of the liquid phase, u g is the velocity vector of the gas phase; S l is the liquid phase source term, S g is the gas phase source term, which is generally 0; t is time; ▽ represents the Hamiltonian operator.

[0074] The liquid phase momentum equation is specifically formulated as follows:

[0075]

[0076] The gas phase momentum equation is specifically formulated as follows:

[0077]

[0078] Where 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 sum of other interphase forces between the liquid phase and the gas phase except the drag force, F gl is the sum of other interphase forces between the gas phase and the liquid phase except the drag force, and F lg =-F gl The other interphase forces include lift, wall lubrication force, turbulent diffusion force and virtual mass force.

[0079] Specifically, there are many interaction forces between the gas phase and the liquid phase, including drag force, lift force, wall lubrication force, turbulent diffusion force and virtual mass force, among which the most important is the drag term β(u l -u g ), and β is the interphase momentum exchange coefficient, which needs to be corrected accordingly in the present invention.

[0080] Step S4: Acquire the ambient pressure, and determine a correction function according to the ambient pressure.

[0081] At different pressures, the gas holdup (i.e., the volume fraction of the gas phase) α in the simulation g After stabilization, α g It presents a certain functional relationship with the pressure P; according to 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.

[0082] The specific formula of the correction function is:

[0083]

[0084] Where P0 is the standard atmospheric pressure, P is the ambient pressure, and ξ is the correction function of the interphase momentum exchange coefficient.

[0085] Step S5: correcting the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model.

[0086] Specifically, the revised computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a revised liquid phase momentum equation, and a revised gas phase momentum equation.

[0087] The modified liquid phase momentum equation is specifically formulated as follows:

[0088]

[0089] The modified gas phase momentum equation is specifically formulated as follows:

[0090]

[0091] 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 pipeline is performed to obtain a gas holdup distribution in the pipeline; the gas holdup distribution represents how the volume fraction of the gas phase changes with radial position at a set position section of the three-dimensional pipeline model.

[0092] Specifically, an unsteady simulation of gas-liquid two-phase flow in a pipe is performed, and the corrected drag model is compiled into its interphase force model. The correction function is used to make a corrected prediction of the bubble drag under different pressure conditions, and finally the time-averaged value of the simulation is taken as the evaluation criterion. Figure 5 、 Figure 6 The average gas holdup distribution diagrams at pressures of 1.0 MPa and 2.0 MPa are given respectively. Figure 5 and Figure 6 It can be found 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.

[0093] Example 2

[0094] like Figure 7 As shown, in order to execute the corresponding method of the above embodiment 1 to achieve the corresponding functions and technical effects, a gas-liquid numerical simulation drag correction system is provided below, including:

[0095] The pipeline model building module 1 is used to build a three-dimensional pipeline model.

[0096] The model meshing module 2 is used to mesh the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model.

[0097] The computational fluid dynamics model establishment module 3 is used to establish an initial computational fluid dynamics model based on the basic theory of the Euler-Euler two-fluid model.

[0098] The correction function determination module 4 is configured to obtain the ambient pressure and determine the correction function according to the ambient pressure.

[0099] The computational fluid dynamics model correction module 5 is configured to correct the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model.

[0100] The gas-liquid numerical simulation module 6 is used to perform unsteady simulation of gas-liquid two-phase flow in the pipe based on the modified computational fluid dynamics model and the meshed three-dimensional pipe model to obtain the gas holdup distribution in the pipe; the gas holdup distribution represents the change in the volume fraction of the gas phase with radial position at a set position section of the three-dimensional pipe model.

[0101] Example 3

[0102] An embodiment of the present invention further provides an electronic device, comprising a memory and a processor, wherein the memory is configured to store a computer program, and the processor is configured to execute the computer program to enable the electronic device to execute the gas-liquid numerical simulation drag correction method of embodiment 1. The electronic device may be a server.

[0103] In addition, 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.

[0104] In summary, the present invention specifically discloses a pressure correction method for gas phase drag based on the Euler-Euler two-fluid model theory. During deep-sea oil and gas mixed transportation, the environment is high-pressure, and the gas phase exhibits different drag at different pressures, resulting in different phase holdup distributions. Therefore, the correction method of the present invention can be used to reasonably predict gas holdup. The drag correction method for gas-liquid numerical simulations under high-pressure conditions, provided by the present invention, is based on the fundamental theory of the Euler-Euler two-fluid model. Based on the relationship between the drag between the gas and liquid phases and the ambient pressure, a pressure-dependent correction relationship is derived through data fitting, enabling more reasonable prediction of the gas holdup distribution within a flow field at different pressures. Different correction coefficients are applied to different pressure ranges, resulting in different trends in the interphase drag at different pressures. This correction is designed for high-pressure environments and gas-liquid two-phase simulations with large pressure variation scales. The magnitude of the drag correction is determined based on the magnitude of the pressure in the flow field, effectively avoiding the inaccurate prediction of gas holdup by existing numerical models under high-pressure conditions.

[0105] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. Reference can be made to the common and similar parts between the various embodiments. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple, and the relevant parts can be referred to the method description.

[0106] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The above examples are only intended to help understand the method and core concept of the present invention. At the same time, those skilled in the art will find that the specific implementation methods and application scopes may vary based on the concept of the present invention. In summary, the contents of this specification should not be construed as limiting the present invention.

Claims

1. A method for drag correction in gas-liquid numerical simulation, characterized in that: include: Establish a 3D pipeline model; Using a polyhedron meshing method to mesh the three-dimensional pipeline model to obtain a meshed three-dimensional pipeline model; Based on the basic theory of Euler-Euler two-fluid model, an initial computational fluid dynamics model was established; 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; The liquid phase continuity equation is specifically formulated as follows: The gas phase continuity equation is specifically formulated as follows: Among them, α 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; u l is the velocity vector of the liquid phase, u g 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 the time; The liquid phase momentum equation is specifically formulated as follows: The gas phase momentum equation is specifically formulated as follows: Where 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 interphase forces between the liquid phase and the gas phase except the drag force, F gl is the sum of other interphase forces between the gas phase and the liquid phase except the drag force, and F lg =-F gl , the other interphase forces include lift force, wall lubrication force, turbulent diffusion force and virtual mass force; Obtain the ambient pressure and determine the correction function according to the ambient pressure; at different pressures, the volume fraction α of the gas phase in the simulation g After stabilization, α g It has a certain functional relationship with the pressure P. 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. The specific formula of the correction function is: Where P0 is the standard atmospheric pressure, P is the ambient pressure, and ξ is the correction function of the interphase momentum exchange coefficient; Correcting the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model; the corrected computational fluid dynamics model includes: a liquid phase continuity equation, a gas phase continuity equation, a corrected liquid phase momentum equation, and a corrected gas phase momentum equation; The modified liquid phase momentum equation is specifically formulated as follows: The modified gas phase momentum equation is specifically formulated as follows: Based on the modified computational fluid dynamics model and the meshed three-dimensional pipeline model, an unsteady simulation of gas-liquid two-phase flow in the pipeline is performed to obtain a gas holdup distribution in the pipeline; the gas holdup distribution represents the change in the volume fraction of the gas phase with radial position at a set position section of the three-dimensional pipeline model.

2. The drag correction method for gas-liquid numerical simulation according to claim 1, characterized in that: The diameter of the three-dimensional pipeline model is 0.15 meters and the total length is 6.6 meters.

3. The drag correction method for gas-liquid numerical simulation according to claim 2, characterized in that: The total number of grids in the three-dimensional pipeline model after grid division is 535,372.

4. A gas-liquid numerical simulation drag correction system, characterized in that: include: Pipeline model building module, used to build a three-dimensional pipeline model; A model meshing module is used to mesh the three-dimensional pipeline model using a polyhedron meshing method to obtain a meshed three-dimensional pipeline model; Computational fluid dynamics model building module, used to build the initial computational fluid dynamics model based on the basic theory of Euler-Euler two-fluid 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; The liquid phase continuity equation is specifically formulated as follows: The gas phase continuity equation is specifically formulated as follows: Among them, α 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; u l is the velocity vector of the liquid phase, u g 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 the time; The liquid phase momentum equation is specifically formulated as follows: The gas phase momentum equation is specifically formulated as follows: Where 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 interphase forces between the liquid phase and the gas phase except the drag force, F gl is the sum of other interphase forces between the gas phase and the liquid phase except the drag force, and F lg =-F gl , the other interphase forces include lift force, wall lubrication force, turbulent diffusion force and virtual mass force; The correction function determination module is used to obtain the ambient pressure and determine the correction function according to the ambient pressure; under different pressures, the volume fraction α of the gas phase in the simulation g After stabilization, α g It has a certain functional relationship with the pressure P. 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. The specific formula of the correction function is: Where P0 is the standard atmospheric pressure, P is the ambient pressure, and ξ is the correction function of the interphase momentum exchange coefficient; a computational fluid dynamics model correction module, configured to correct the initial computational fluid dynamics model according to the correction function to obtain a corrected computational fluid dynamics model; the corrected computational fluid dynamics model comprising: a liquid phase continuity equation, a gas phase continuity equation, a corrected liquid phase momentum equation, and a corrected gas phase momentum equation; The modified liquid phase momentum equation is specifically formulated as follows: The modified gas phase momentum equation is specifically formulated as follows: The gas-liquid numerical simulation module is used to perform unsteady simulation of gas-liquid two-phase flow in the pipe based on the modified computational fluid dynamics model and the meshed three-dimensional pipe model to obtain the gas holdup distribution in the pipe; the gas holdup distribution represents the change in the volume fraction of the gas phase with radial position at a set position section of the three-dimensional pipe model.

5. An electronic device, characterized in that: The 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 gas-liquid numerical simulation drag correction method according to any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that It stores a computer program, which, when executed by a processor, implements the gas-liquid numerical simulation drag correction method according to any one of claims 1 to 3.

Images