An analysis method for T-type gas-liquid two-phase pipe flow pressure fluctuation attenuator

Abstract

The present application relates to the technical field of fluid mechanics, and discloses an analysis method for a T-shaped gas-liquid two-phase pipeline pressure pulsation attenuator, which comprises the following steps: S1: establishing a T-shaped gas-liquid two-phase pipeline geometric model and a corresponding calculation model; S2: performing geometric model verification, analyzing the flow pattern and flow field characteristics of the T-shaped gas-liquid two-phase pipeline flow by respectively changing the gas-liquid velocity and the gas phase inflow; S3: simulating the high-frequency noise response of the pressure generated by the T-shaped gas-liquid two-phase pipeline; S4: simulating H-shaped attenuators of various sizes, and comparatively analyzing the attenuation effect of a specific pressure frequency band under different sizes of the resonant cavity; and S5: obtaining the physical structure parameters of the attenuator for better attenuating the pressure pulsation of the T-shaped gas-liquid two-phase pipeline flow by comparatively analyzing the attenuation effects of different specifications of the H-shaped attenuator. The present application considers various influencing factors of the pipeline flow characteristics in the T-shaped gas-liquid two-phase pipeline, simulates the pipeline flow in a real scene, and analyzes the H-shaped attenuator with a better attenuation effect.

Description

Technical Field

[0001] The present invention relates to the technical field of fluid mechanics, and in particular to an analysis method for a T-type gas-liquid two-phase pipe flow pressure pulsation attenuator. Background Art

[0002] Gas-liquid two-phase flow systems are widely used in oil pipelines, chemical plants, ship pipelines, nuclear power plants, and other applications requiring gas-liquid exchange, mixing, and separation. Numerous reports indicate that gas-liquid two-phase flow in pipelines can cause significant pipeline vibration, wear, failure, and noise.

[0003] Gas-liquid two-phase flow can cause severe pipe vibration and generate high noise levels. In ship operations, the vibration and noise caused by fluid transmission in pipelines is a significant factor directly impacting the ship's noise performance, and controlling pipeline vibration and noise is a crucial component of ship vibration and noise reduction engineering. After hazardous gases are treated in the combustion chamber of a ship's generator, the remaining gas is a mixture of carbon dioxide and oxygen, with carbon dioxide comprising 90%. This gas is then mixed into seawater through a T-shaped pipe. The dissolution of the carbon dioxide creates a complex gas-liquid two-phase flow, generating significant noise pollution.

[0004] There are few studies on the pressure pulsation and attenuation analysis of T-type gas-liquid two-phase pipe flow. The main focus of T-type pipeline research is on the two-phase separation phenomenon, and little attention is paid to the mixing of pure liquid and pure gas phases. Most of the related research focuses on the two-phase mixing in microchannel T-tubes or reactors. There are few studies on the gas-liquid two-phase mixing of T-type pipes with pipe diameters of about 50mm to 3mm. Therefore, there are even fewer pressure pulsation attenuation devices for this type of T-type gas-liquid two-phase pipe flow. H-type attenuator, such as Figure 2 As shown in the figure, as the main research object, it mainly uses the oil column resonance principle to form anti-phase waves to attenuate pulsation in a specific frequency band. It has a simple structure, good attenuation effect, and is widely used. In addition to achieving the purpose of noise reduction by avoiding the generation of flow types such as intermittent flow, stirring flow, and wavy flow, the existing technology has not yet achieved the noise control of two-phase flow by changing the internal specific structure of the H-type attenuator, so as to install an attenuator to control the flow field pulsation in the pipe. Summary of the Invention

[0005] In order to overcome or alleviate one or more of the above technical problems, the purpose of the present invention is to provide an analysis method for a T-type gas-liquid two-phase pipe flow pressure pulsation attenuator, which analyzes the corresponding H-type attenuator for suppressing high-frequency frequency response based on the characteristics of the simulated T-type gas-liquid two-phase pipe flow, and compares the parameters to select the optimal parameters so that the pressure pulse is more effectively suppressed, thereby greatly reducing the noise generated by the mixing of the two-phase flow.

[0006] The application provides the following technical solutions:

[0007] An analysis method for a T-shaped gas-liquid two-phase pipe flow pressure fluctuation attenuator, comprising the following steps:

[0008] S1: establishing a T-shaped gas-liquid two-phase pipe geometry model and a corresponding calculation model thereof;

[0009] S2: verifying the geometry model, setting physical parameters and boundary conditions of the T-shaped gas-liquid two-phase pipe, and analyzing flow pattern and flow field characteristics of the T-shaped gas-liquid two-phase pipe flow by respectively changing gas-liquid velocity and gas phase inflow;

[0010] S3: taking a pressure or flow velocity fluctuation signal in M sequence form as an excitation source of a gas-liquid inlet of the T-shaped gas-liquid two-phase pipe to simulate a pressure high-frequency noise response generated by the T-shaped gas-liquid two-phase pipe;

[0011] S4: simulating H-shaped attenuators of various sizes, and comparing and analyzing attenuation effects of specific pressure bands under different sizes of resonant cavities;

[0012] S5: obtaining physical structure parameters of an attenuator for better attenuating T-shaped gas-liquid two-phase pipe flow pressure fluctuation through comparison and analysis of attenuation effects of different specifications of H-shaped attenuators.

[0013] According to some embodiments, the T-shaped gas-liquid two-phase pipe geometry model in the S1 step is simulated by a simulation software, and the simulation steps are as follows:

[0014] S11: calculating geometry construction, drawing by using three-dimensional model drawing software to obtain a geometry model;

[0015] S12: calculating a mesh division scheme, importing the obtained geometry model into the simulation software to perform mesh processing of the model;

[0016] S13: calculating a required model and selecting parameters, for two-phase flow simulation, selecting a solver of the three-dimensional model drawing software and a corresponding turbulence model, a hybrid k-ε model, an interphase interaction turbulence model, and a flow field solving algorithm including PISO, PIMPLE and SIMPLE algorithms;

[0017] S14: setting physical parameters of a flow field, the physical parameters including interphase mass and heat transfer physical parameters required for two-phase flow;

[0018] S15: specifying boundary conditions and initial conditions, and the influence of changes in initial fields and boundary parameters on the flow field is an analysis focus of simulation;

[0019] S16: specifying a discrete format, and the discrete object is each term in a control equation;

[0020] S17: Control the calculated parameters, the purpose is to ensure the calculation accuracy, data sampling accuracy and sampling length, and complete the numerical simulation calculation of the flow field at the fastest speed.

[0021] According to some embodiments, in step S4, the calculation formula of the resonance angular frequency of the H-type attenuator is:

[0022]

[0023] In the formula, ω is the resonance angular frequency, a is the sound speed of the two-phase flow, gas or liquid, D1 is the neck diameter of the attenuator, L1 is the neck length, and V2 is the resonance cavity volume.

[0024] According to some embodiments, in step S1, the calculation model includes a gas-liquid two-phase flow control equation, a gas-liquid two-phase flow turbulence model, a boundary layer and wall function, and a solubility model, and the gas-liquid two-phase flow control equation includes an average continuity equation and an average momentum conservation equation, and an average energy conservation equation.

[0025] Compared with the prior art, the present application has the following beneficial effects:

[0026] The present application compares and analyzes the attenuator by simulating the pipe flow in the T-shaped gas-liquid two-phase pipeline. In the simulation process, various influencing factors of the pipe flow characteristics in the T-shaped gas-liquid two-phase pipeline are fully considered. In order to more accurately describe the frequency response of the two-phase flow system under different excitations and obtain more responses under different frequency components, the M sequence signal is added to the inlet boundary. The grid independence analysis is performed on the numerical model, the calculation efficiency and the accuracy of the results are balanced, and the most suitable grid division scheme is selected for subsequent calculation. In order to verify the accuracy of the flow pattern prediction of the calculation model, the gas-liquid two-phase flow in the T-shaped gas-liquid two-phase pipeline is also simulated and analyzed.

[0027] On the basis of grid division, different gas-liquid inlet velocity boundary conditions are selected, and the flow pattern change of the T-shaped pipeline, the outlet pressure, the outlet gas-liquid velocity fluctuation characteristics, the amplitude change and the change of the outlet solubility efficiency are analyzed. Then the variation law of the flow characteristics under different gas phase inflow angles and different environmental pressures is investigated. In this way, the pressure fluctuation is reduced and the CO2 solubility efficiency is improved.

[0028] Unstable pressure and velocity excitation signals, such as M sequence pseudo-random signals, are applied to the inlet of the T-shaped pipeline to simulate the frequency response of the pipeline. Several different pulsation attenuators and different arrangement modes are compared and analyzed to ensure that only liquid or gas is in the attenuator cavity (i.e. the attenuator only attenuates the pulsation of gas or liquid). The insertion loss and transmission loss of each attenuator are obtained through simulation to evaluate the attenuation effect of the attenuator on the pressure pulsation noise. BRIEF DESCRIPTION OF DRAWINGS

[0029] Figure 1 A schematic diagram of a T-shaped gas-liquid two-phase pipeline provided by an embodiment of the present application.

[0030] Figure 2 A schematic diagram of an H-shaped attenuator mentioned in the background of the present application.

[0031] Figure 3 A flow chart of a geometric two-dimensional simulation method in step S1 of a design method of a T-shaped gas-liquid two-phase pipe flow pressure pulsation attenuator provided by an embodiment of the present application.

[0032] Figure 4 A size diagram of a T-shaped gas-liquid two-phase pipeline provided by an embodiment of the present application.

[0033] Figure 5 Table 1 of main flow boundary types and initial values provided by an embodiment of the present application.

[0034] Figure 6 A size diagram of a T-shaped gas-liquid two-phase pipeline with a gas phase inlet angle provided by an embodiment of the present application.

[0035] Figure 7 Table 2 of gas phase inlet angle and selection of reduced velocity provided by an embodiment of the present application.

[0036] Figure 8 A schematic diagram of T-U and T-D installation methods provided by an embodiment of the present application.

[0037] Figure 9 A size identification diagram of a resonant cavity of an H-shaped attenuator provided by an embodiment of the present application.

[0038] Figure 10 Table 3 of attenuator geometric sizes and resonant frequencies provided by an embodiment of the present application. DETAILED DESCRIPTION

[0039] The present application is described in detail below with reference to the embodiments and drawings, but it should be understood that the embodiments and drawings are only used to exemplarily describe the present application, and cannot constitute any limitation on the protection scope of the present application. All reasonable transformations and combinations within the scope of the inventive concept of the present application fall within the protection scope of the present application.

[0040] The present application is further described below with reference to the drawings.

[0041] Embodiment 1

[0042] This embodiment uses CO2-water as a gas-liquid two-phase, and describes the analysis process of a T-shaped gas-liquid two-phase pipe flow pressure pulsation attenuator.

[0043] S1: Establishing a T-shaped gas-liquid two-phase pipeline geometry model and its corresponding calculation model;

[0044] As Figure 1 , the structure of the T-shaped gas-liquid two-phase pipeline is shown in the figure. The left end is the liquid phase inlet, the upper end is the gas phase inlet, and the right end is the gas-liquid mixed outlet. The change of gas-liquid inlet flow velocity affects the simulation of T-shaped gas-liquid two-phase pipeline in terms of flow pattern, pressure fluctuation, and solubility.

[0045] Therefore, the geometry two-dimensional simulation of T-shaped gas-liquid two-phase pipeline is carried out by third-party software, such as Figure 3 The simulation steps are as follows:

[0046] (1) Calculation of the geometry construction, mainly using third-party three-dimensional model drawing software to draw; the third-party three-dimensional model drawing software is OpenFOAM;

[0047] (2) Calculation of the meshing scheme, which can choose to import the geometry model of the third-party software into the simulation software such as Fluent for mesh processing of the fitted model. If it is a simple two-dimensional model, step 1 can be skipped, and the blockMesh tool can be used to process the mesh of the flow field;

[0048] (3) Calculation of the required model and parameter selection. For two-phase flow simulation, OpenFOAM provides solvers such as twoPhaseEulerFoam, interFoam, reactingTwoPhaseEulerFoam, and corresponding turbulence models such as mixed k-ε model, interphase interaction turbulence model, etc. The flow field solving algorithm includes PISO, PIMPLE, SIMPLE algorithm;

[0049] (4) Setting of the physical parameters of the flow field. For two-phase flow, the interphase mass, heat transfer, and other physical parameters need to be specified;

[0050] (5) Specification of boundary conditions and initial conditions. The change of initial field and boundary parameters has a great influence on the flow field, which is the research focus of simulation;

[0051] (6) Specification of the discrete format, mainly for each term in the control equation;

[0052] (7) Parameters for controlling calculation, the main purpose of which is to ensure the calculation accuracy, data sampling accuracy, and sampling length, and to complete the numerical simulation calculation of the flow field at the fastest speed.

[0053] The calculation model includes gas-liquid two-phase flow control equation, gas-liquid two-phase flow turbulence model, boundary layer and wall function, and solubility model. The gas-liquid two-phase flow control equation includes average continuity equation, average momentum conservation equation, and average energy conservation equation, etc.

[0054] The T-shaped gas-liquid two-phase pipeline two-dimensional model as shown in Figure 4 is obtained through the S1 step, considering that the CO2-water two-phase flow system is compressible, and the reactingTwoPhaseEulerFoam solver is used for solving. As the standard two-phase solver of OpenFOAM, it is based on the Euler-Euler model and can process each equivalent as a continuous phase. In order to more accurately describe the turbulent parameters under different phase fractions, the hybrid k-ε model is selected as the turbulent model.

[0055] S2: Perform geometric model verification, set the physical parameters and boundary conditions of the T-shaped gas-liquid two-phase pipeline; analyze the flow pattern and flow field characteristics of the T-shaped gas-liquid two-phase pipeline flow by changing the gas-liquid velocity and gas phase inflow respectively;

[0056] In this embodiment, the gas phase is a mixture of gaseous water and CO2, and the liquid phase is a mixture of CO2 dissolved in water. In the initial state, the gas phase branch pipe section will be full of gas, and the horizontal pipe section will be full of liquid.

[0057] The boundary condition type and value of the main parameters are shown in Table 1 Figure 5 , in which P is the absolute pressure, T is the absolute temperature, U is the velocity vector, and its subscripts x, y, and z represent the three components of the velocity in the Cartesian coordinate system, mainly referring to the coordinate system in Figure 4 , and ω represents the mass fraction of the component, and the subscripts g / l indicate that the phase of the flow field is gas / liquid, and the subscripts c / w indicate that the component is CO2 / water, i.e. ωc,g represents the mass fraction of component CO2 in the liquid phase.

[0058] In the table, the unit with only numbers represents a fixed value (fixedValue), which is a Dirichlet boundary condition, i.e. the value of a certain field at a certain boundary is a constant, and the zero gradient (zeroGradient) is a Neumann boundary condition. The inletOutlet is a special outlet boundary condition, which is different from the zero gradient condition, which can set the backflow of the flow field to 0 or other fixed values. The remaining boundary conditions such as fixed flow (fixedFluxPressure), special pressure for relative pressure, and turbulent boundary conditions such as k wall function (kqRWallFunction) and epsilon wall function (epsilonWallFunction) ensure the closure of the turbulent conservation equation.

[0059] For the physical parameters, the calculation model is downward in the direction of gravity. In addition, the reactingTwoPhaseEulerFoam solver also needs to specify the phase parameter file and the thermal physical parameter file of the two-phase flow system. Among them, the phase parameter file specifies that the average diameter of the gas bubble is 0.1 mm.

[0060] The increase of gas superficial velocity will lead to the increase of cross-sectional void fraction in the horizontal main pipe, and thus more easily form the plug flow with continuous gas-liquid alternation. When the gas superficial velocity increases to a certain extent, the flow pattern in the pipe will evolve from plug flow to stratified flow. The increase of liquid superficial velocity will decrease the cross-sectional void fraction in the main pipe, and more easily form bubbles and plugs. The increase of gas superficial velocity will lead to the increase of plug size, and thus decrease the frequency of pressure fluctuation in the pipe. The larger gas momentum increases the amplitude of pressure fluctuation in the developing stage, but because it is more likely to develop into stratified flow, the amplitude of pressure fluctuation in the stable stage will decrease. Meanwhile, the increase of gas superficial velocity will significantly increase the CO2 dissolution in water, and the increase of liquid superficial velocity will promote the formation of shorter plugs, and thus slightly increase the frequency of pressure fluctuation in the pipe and slightly decrease the CO2 dissolution in water.

[0061] The effects of different gas inlet angles θ on the flow characteristics of T-type pipe are as follows:

[0062] As Figure 6 , this embodiment selects five groups of gas branch pipe inclination angles θ, i.e. 45°, 67°, 90°, 112° and 135°, to change the gas inlet direction. The selection of gas inlet angle and superficial velocity is shown in Table 2 in Figure 7 . Since the change of gas superficial velocity will greatly affect the change of pipe parameters, this embodiment also sets different gas-liquid velocities, and changes the environmental pressure to 1 MPa to simulate the pressure environment of seawater.

[0063] 1. When the environmental pressure is increased to 1 MPa and the gas-liquid mixing pipe section is short, the flow pattern in the T-type pipe will be more stable, and at a stable gas-liquid inlet velocity, the flow pattern in the pipe will be stable in stratified flow. A lower gas superficial velocity is not conducive to resisting the impact of liquid flow, and thus in the developing stage of the flow pattern, the gas flow will break into smaller bubbles. When the gas-liquid superficial velocity is kept constant, the developing stage of the flow pattern will produce more bubbles when the gas inlet angle is 90°, i.e. the horizontal component of the gas inlet velocity is zero.

[0064] 2. At a large gas superficial velocity, the closer θ is to 90°, the larger the amplitude of pressure fluctuation will be. With the increase of gas superficial velocity, the amplitude of pressure fluctuation will also increase significantly. The closer θ is to 90°, the larger the frequency of pressure fluctuation will be. When θ increases, the horizontal component of the gas inlet velocity will also increase, and the gas and liquid will flow to the outlet without being fully mixed, and the dissolution of gas in water will decrease.

[0065] S3: The pressure or flow rate fluctuation signal in the form of M sequence is applied to the T-type gas-liquid two-phase pipeline as the excitation source of the gas-liquid inlet of the T-type gas-liquid two-phase pipeline to simulate the high-frequency noise response generated by the T-type gas-liquid two-phase pipeline;

[0066] The M sequence is generated by a shift register, and a series of discrete functions f(t) changing with time is obtained, that is, a series of discrete functions changing with time is obtained, and the M sequence excitation signal of the corresponding flow field can be obtained by multiplying the amplitude of the inlet flow field.

[0067] As a special pseudo-static response signal, the M sequence signal can give the system an approximate static excitation, so as to observe the stability, vibration, resonance, attenuation characteristics and the like of the system, and the advantage is that a large number of frequency components can be obtained in one period, so that the simulation time is shortened. When a sine sweep signal is used as an input signal to excite the system, only a few frequencies can be identified each time.

[0068] The M sequence is generated by a generator, which is usually composed of the following components:

[0069] (1) Shift register: This is the core part of the M sequence generation, which is usually a finite state binary register, for example, an 8-bit register. It can store a series of continuous bit sequences {ak} = {a1, a2, a3,..., an}. At each clock cycle, the content of the register will move one bit to the left.

[0070] (2) Feedback taps: These are the positions of the signals extracted from some bit positions of the shift register. For an 8-bit register, there can be multiple feedback taps, such as C1, C2,..., Ck, where Ci represents whether to extract a signal from the ith bit (Ci = 1 means extraction, and Ci = 0 means no extraction). The bit positions selected by these taps will participate in the next modulo-2 addition operation.

[0071] (3) Modulo-2 adder: 2 addition is a special binary addition in which there is no carry operation, only XOR operation. The extracted feedback tap bits are operated through the modulo-2 adder, and the result is placed in the leftmost side of the shift register.

[0072] S4: Simulate various sizes of H-type attenuators, and compare and analyze the attenuation effect of a specific pressure frequency band under the size of the resonant cavity;

[0073] As shown in Figure 8 , the H-type attenuator has a resonant cavity which can be installed on the same side or on the opposite side of the gas phase inlet, as shown in Figure 8 .a is the same side, that is, the "T-U" installation mode; Figure 8 .b is the opposite side, that is, the "T-D" installation mode.

[0074] According to the lumped parameter method, the pulsation attenuator can be regarded as an inductance-capacitance series circuit in the circuit, wherein the neck of the attenuator is the inductance of the attenuated current, and the cavity is the capacitance of the attenuated voltage. When the frequency of the external pressure pulsation is consistent with the resonance frequency of the muffler, a better pulsation attenuation effect can be achieved. Through this approximation, the calculation formula of the resonance angular frequency of the attenuator can be finally obtained as follows:

[0075]

[0076] In the formula, ω is the resonance angular frequency, a is the sound speed of the two-phase flow, gas or liquid, D1 is the diameter of the neck of the attenuator, L1 is the length of the neck, and V2 is the volume of the resonance cavity. The specific schematic diagram is as shown in Figure 9 .

[0077] The sound speed used in the above calculation formula is the sound speed in the pure liquid or pure gas medium, so as to simulate the pure liquid or pure gas state in the cavity, and the resonance frequency of the attenuator is designed for the main frequency of the pressure pulsation under different inlet excitation modes. The attenuator size under the corresponding gas-liquid medium and excitation is shown in Table 3-attenuator geometric size and its resonance frequency. Figure 10

[0078] S5: Through the analysis and comparison of the attenuation effects of different specifications of H-type attenuators, the physical structure parameters of the attenuator for better attenuating the pressure pulsation of T-type gas-liquid two-phase pipe flow are obtained.

[0079] In the "T-U" installation mode, that is, the attenuator is placed upright above the T-type pipe, at this time, the cavity of the attenuator is almost filled with gas, and the attenuation effects of the attenuators at the characteristic frequencies of the T-type pipe are similar, and the attenuation amplitude is only about 13 dB.

[0080] In the "T-D" installation mode, that is, the attenuator is placed upside down below the T-type pipe, at this time, the cavity of the attenuator is almost filled with liquid. In order to improve the calculation efficiency, the reduced gas speed is reduced to Jg=0.161 m / s, and the attenuator arrangement distance is extended. Under this working condition, the T-type pipe has multiple characteristic frequencies at high frequencies; the attenuators 1 and 3 have good attenuation effects at the 3rd, 4th, 5th and 6th frequencies of the T-type pipe, that is, near 211.11, 436.11, 722.22 and 1016.67 Hz, wherein the attenuation effect of the attenuator 1 near the 4th frequency is the best, close to the designed frequency, and the attenuation amplitude is 27 dB; the attenuator 2 can better suppress the pressure pulsation at the 5th frequency, and the attenuation amplitude is 29 dB. The resonance frequency of the attenuator 2 and 4 filled with water is relatively high, and the attenuation effect near the four frequencies is general.

[0081] ​The above examples are only preferred embodiments of the present application, and the protection scope of the present application is not limited to the above examples. Any technical scheme falling within the concept of the present application belongs to the protection scope of the present application. It should be pointed out that improvements and refinements made by those of ordinary skill in the art without departing from the principles of the present application should also be considered as falling within the protection scope of the present application.

Claims

1. An analysis method for a T-type gas-liquid two-phase pipe flow pressure pulsation attenuator, characterized by: The following steps are involved: S1: Establish a T-type gas-liquid two-phase pipeline geometric model and its corresponding calculation model; S2: Verify the geometric model and set the physical parameters and boundary conditions of the T-shaped gas-liquid two-phase pipeline; analyze the flow pattern and flow field characteristics of the T-shaped gas-liquid two-phase pipeline by changing the gas-liquid velocity and gas inflow; S3: A pressure or flow rate pulsation signal in the form of an M sequence is applied to the T-shaped gas-liquid two-phase pipeline as an excitation source at the gas-liquid inlet of the T-shaped gas-liquid two-phase pipeline to simulate the high-frequency pressure noise response generated by the T-shaped gas-liquid two-phase pipeline; S4: Simulate H-type attenuators of various sizes and compare and analyze the attenuation effect of specific pressure frequency bands under different resonant cavity sizes; S5: By analyzing and comparing the attenuation effects of H-type attenuators of different specifications, the physical structural parameters of the attenuator that is best for attenuating pressure pulsations in T-type gas-liquid two-phase pipe flows are obtained; The T-shaped gas-liquid two-phase pipeline geometric model described in step S1 is subjected to a two-dimensional geometric simulation using simulation software. The simulation steps are as follows: S11: geometric construction of the calculation field, using 3D model drawing software to draw and obtain a geometric model; S12: The calculated meshing scheme is imported into the simulation software to perform mesh processing to fit the model; S13: Calculate the required model and parameter selection. For two-phase flow simulation, select the solver of the 3D model rendering software and the corresponding turbulence model, hybrid k-ε model, interphase interaction turbulence model, and the flow field solution algorithm including PISO, PIMPLE, and SIMPLE algorithms; S14: Setting the physical parameters of the flow field, including the physical parameters of mass and heat transfer between phases required for two-phase flow; S15: Specifying boundary conditions and initial conditions, and the impact of changes in initial field and boundary parameters on the flow field are the key points of simulation analysis; S16: Specification of discrete format, where the discrete objects are the terms in the control equation; S17: Control calculation parameters to ensure calculation accuracy, data sampling accuracy and sampling length while completing the numerical simulation calculation of the flow field at the fastest speed; The calculation model in step S1 includes the gas-liquid two-phase flow control equation, the gas-liquid two-phase flow turbulence model, the boundary layer and wall function and the solubility model. The gas-liquid two-phase flow control equation includes the average continuity equation, the average momentum conservation equation and the average energy conservation equation.

2. The analysis method for a T-type gas-liquid two-phase pipe flow pressure pulsation attenuator according to claim 1, characterized in that: In step S4, the calculation formula for the resonant angular frequency of the H-type attenuator is: Where ω is the resonant angular frequency, a is the speed of sound of the two-phase flow, gas, or liquid, D1 is the attenuator neck diameter, L1 is the neck length, and V2 is the resonant cavity volume.

Images