Hydraulic and thermal coupling liquid ammonia pipeline flow characteristic simulation method and device

Abstract

The invention discloses a hydraulic and thermal coupling liquid ammonia pipeline flow characteristic simulation method and a hydraulic and thermal coupling liquid ammonia pipeline flow characteristic simulation device. The method comprises the steps that on the basis of a liquid ammonia pipeline physical model, a liquid ammonia pipeline flow unsteady-state hydraulic and thermal coupling model is established and comprises an energy conservation equation, a momentum conservation equation, a continuity equation and a soil heat conduction equation; according to the current grid density, carrying out two-dimensional finite element grid division and one-dimensional pipe section discretization processing on the coupling model; based on the current initial condition and the boundary condition, solving the discrete coupling model to obtain flow field data of each node; judging whether the flow field data of each node converges or not; if not, the initial condition and the boundary condition are modified, and / or the grid density is modified; and solving again to obtain the flow field data of each node. According to the method, the flow characteristics of the liquid ammonia pipeline are accurately calculated, and reliable fluid mechanics information is provided for design and safe and stable operation of the liquid ammonia long-distance pipeline.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ammonia storage and transportation, and particularly relates to a method and device for simulating flow characteristics of a liquid ammonia pipeline coupled with hydraulic and thermal forces. BACKGROUND

[0002] Ammonia is widely concerned as an efficient and economical energy storage and hydrogen carrier in the context of green and low-carbon energy transformation. Due to the leakage risk of liquid ammonia, long-distance transportation pipelines are often laid underground. Under the coupling effect of the surrounding soil environment, the average parameter calculation method used for conventional oil and gas pipelines is no longer applicable to the description of the flow characteristics of the liquid ammonia pipeline. At the same time, the high volume expansion and temperature sensitivity of ammonia make the pipeline always in a hydraulic and thermal alternating state, which further leads to the difficulty in predicting the distribution of the liquid ammonia physical field along the pipeline. In view of the complexity of the liquid ammonia pipeline flow process, the research on the flow characteristics of the liquid ammonia pipeline is relatively scarce. SUMMARY

[0003] In order to at least partially solve the technical problems existing in the prior art, the present application is made by the inventors, and through specific embodiments, a method and device for simulating flow characteristics of a liquid ammonia pipeline coupled with hydraulic and thermal forces are provided, which realizes accurate calculation of the flow characteristics of the liquid ammonia pipeline and provides reliable fluid mechanics information for the design and safe and stable operation of the long-distance liquid ammonia pipeline.

[0004] In the first aspect, an embodiment of the present application provides a method for simulating flow characteristics of a liquid ammonia pipeline coupled with hydraulic and thermal forces, comprising:

[0005] On the basis of a physical model of the liquid ammonia pipeline, a non-steady-state hydraulic and thermal coupling model of the liquid ammonia pipeline is established, including an energy conservation equation, a momentum conservation equation, a continuity equation, and a soil heat conduction equation;

[0006] According to the current grid density, two-dimensional finite element grid division and one-dimensional pipe segment discretization processing are performed on the coupling model;

[0007] Based on the current initial conditions and boundary conditions, the discretized coupling model is solved to obtain flow field data of each node;

[0008] It is determined whether the flow field data of each node converges;

[0009] If not, the current initial conditions and boundary conditions are modified, and / or the current grid density is modified; and the flow field data of each node is re-solved.

[0010] In the second aspect, an embodiment of the present application provides a device for simulating flow characteristics of a liquid ammonia pipeline coupled with hydraulic and thermal forces, comprising:

[0011] The coupling model establishing module is configured to establish a non-steady-state hydraulic-thermal coupling model of the liquid ammonia pipeline flow on the basis of the physical model of the liquid ammonia pipeline, including an energy conservation equation, a momentum conservation equation, a continuity equation and a soil heat conduction equation.

[0012] The discretization module is configured to perform two-dimensional finite element mesh division and one-dimensional pipe segment discretization processing on the coupling model according to the current grid density.

[0013] The flow field data solving module is configured to solve the discretized coupling model to obtain flow field data of each node based on the current initial condition and the boundary condition.

[0014] The judgment module is configured to judge whether the flow field data of each node converges.

[0015] The correction module is configured to modify the current initial condition and the boundary condition and / or modify the current grid density if the judgment module judges that the flow field data of each node does not converge.

[0016] In a third aspect, an embodiment of the present application provides a computer storage medium, which stores computer executable instructions, and the computer executable instructions are executed by a processor to implement the simulation method for the hydraulic-thermal coupling liquid ammonia pipeline flow characteristics.

[0017] In a fourth aspect, an embodiment of the present application provides a server, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor implements the simulation method for the hydraulic-thermal coupling liquid ammonia pipeline flow characteristics when executing the program.

[0018] The simulation method for the hydraulic-thermal coupling liquid ammonia pipeline flow characteristics provided by the embodiment of the present application has at least the following beneficial effects:

[0019] The simulation method for the hydraulic-thermal coupling liquid ammonia pipeline flow characteristics provided by the embodiment of the present application has at least the following beneficial effects:

[0020] Additional features and advantages of the present application will be set forth in the description that follows, and in part will be apparent from the description, or can be learned by practice of the application. The objectives and other advantages of the present application will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.

[0021] The technical solutions of the present application are described in further detail below with reference to the accompanying drawings and examples. BRIEF DESCRIPTION OF DRAWINGS

[0022] The accompanying drawings are included to provide a further understanding of the present application and are incorporated in and constitute a part of the specification, illustrate embodiments of the present application and, together with the description, serve to explain the present application. In the drawings:

[0023] Figure 1 Flow chart of the method for simulating flow characteristics of the hydromechanically and thermally coupled liquid ammonia pipeline in the embodiment of the present application;

[0024] Figure 2 Structural schematic diagram of the physical model of the liquid ammonia pipeline in the embodiment of the present application;

[0025] Figure 3 Pressure and temperature relative error diagrams of the model at different monitoring points in the embodiment of the present application;

[0026] Figure 4 Structural schematic diagram of the device for simulating flow characteristics of the hydromechanically and thermally coupled liquid ammonia pipeline in the embodiment of the present application. DETAILED DESCRIPTION

[0027] Exemplary embodiments of the present disclosure will be described more fully hereinafter with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it is to be understood that the present disclosure is not to be limited to the embodiments shown in the drawings, which are provided for the purpose of explanation only. Rather, the present disclosure is to cover all modifications, equivalents, and alternatives falling within the scope of the present disclosure. The description of the exemplary embodiments is intended to be illustrative, and not to limit the scope of the present disclosure. As such, the following detailed description is not restrictively to the accompanying drawings.

[0028] It is to be understood that the terminology used in the present application is solely for the purpose of describing particular embodiments of the present application and is not intended to be limiting of the present application. Additionally, for the purposes of the present application, ranges of values are to be understood as being inclusive of the values that are recited within the range. Each smaller range that falls within a broader range is to be understood to be specifically disclosed. The upper and lower limits of these smaller ranges are independently combinable with each other as well as the upper and lower limits of the broader range, to form further smaller ranges within the broader ranges.

[0029] Unless otherwise indicated, all technical and scientific terms used herein have the same meaning as those commonly understood by one of ordinary skill in the art to which this application belongs. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present application, the preferred methods and materials are described. All publications mentioned in this specification are herein incorporated by reference to disclose and describe the methods and / or materials in connection with which the publications are cited. The citation of any reference is not an admission that it is prior art with respect to the present application.

[0030] The embodiment of the present application provides a water and heat coupled liquid ammonia pipeline flow characteristic simulation method and device, and accurately calculates the liquid ammonia pipeline flow characteristic.

[0031] Embodiments

[0032] The embodiment of the present application provides a water and heat coupled liquid ammonia pipeline flow characteristic simulation method, and a flow thereof is shown in Figure 1 The method comprises the following steps:

[0033] Step S11: on the basis of a liquid ammonia pipeline flow physical model, a liquid ammonia pipeline flow unsteady water and heat coupled model is established.

[0034] First, a liquid ammonia pipeline physical model is determined, including pipeline geometric parameters and upstream boundary related parameters. Referring to Figure 2 A structure diagram of the liquid ammonia pipeline physical model is shown.

[0035] On this basis, a liquid ammonia pipeline flow unsteady water and heat coupled model is established, including:

[0036] 1. Energy conservation equation

[0037]

[0038]

[0039] In formula (1)-(3), h is field enthalpy, J / kg; m is fluid mass, kg; U is flow rate, m / s; M is mole fraction; subscript i represents the i-th fluid, and subscripts t and z represent time and space, respectively; S is enthalpy source, J; Q wall is heat flow through the pipe wall, J; p is fluid pressure, kPa; T l is fluid temperature, K; φ is fugacity coefficient, C p,l is fluid constant-pressure specific heat capacity, J / (kg·℃); C p0 , C p1 , C p2 , C p3 and C p4 are different specific heat coefficients, which are respectively 4.238, -0.004215, 2.041x10-6 , -2.126 x 10 -9 and 2.61 x 10 -11 ; R is the ideal gas constant, taken as 8.314 J / (mol K).

[0040] 2. Momentum conservation equation

[0041]

[0042] In formula (4), H(%) is the cross-sectional liquid holdup, p l (kg / m) is the fluid density, v(m / s) is the flow rate, τ(N / m) is the shear force between the fluid and the pipe wall, g(m / s 2 ) is the acceleration of gravity, θ(°) is the pipe inclination angle, and x is the pipe axial direction.

[0043] 3. Continuity equation

[0044]

[0045] 4. Soil heat conduction equation

[0046]

[0047] In formula (6), p s (kg / m) is the soil density, c p,s (J / (kg·℃)) is the soil constant-pressure specific heat capacity, Ts(℃) is the soil temperature, λs(W(m·℃)) is the soil thermal conductivity, and y is the pipe radial direction.

[0048] Step S12: Based on the liquid ammonia pipe flow physical model, the two-dimensional finite element grid division and one-dimensional pipe segment discretization processing are performed on the coupled model according to the set grid density.

[0049] Further, the two-dimensional finite element grid division can be performed on the coupled model based on the Delaunay triangulation method. More precise grid elements are used in the vicinity of the pipe wall, so as to accurately predict the soil temperature field. The Delaunay triangulation method follows the “minimum angle maximum” and “empty circumscribed circle” criteria for division for any given planar point set.

[0050] The two-dimensional grid division is performed in the pipe radial direction (cross section), the one-dimensional grid division is performed in the pipe axial direction, and finally the solution is performed under the comprehensive discrete results of the two. Through the grid division, the continuous calculation domain is divided into discrete grids, and then the discrete solution of the basic equation in the calculation domain in step S13 is realized.

[0051] Step S13: Based on the current initial condition and boundary condition, the discrete coupled model is solved to obtain the flow field data of each node.

[0052] The coupled model after discretization is solved by finite difference method.

[0053] Further, the energy conservation equation and the soil heat conduction equation after discretization are combined to obtain a pipe temperature field control equation:

[0054]

[0055] J q ·n=-k(T s -T amb ) (11)

[0056]

[0057] wherein, A l is the cross-sectional area of the fluid, A s is the cross-sectional area of the soil, A wall is the cross-sectional area of the pipe wall, p wall is the density of the pipe wall, c p,wall is the specific heat capacity of the pipe wall at constant pressure, T wall is the pipe wall temperature, Q s is the heat flow through the soil, and d is an operator representing the change, J q (W / m 2 ) is the heat flux at the boundary of the finite element, n is the unit normal vector at the boundary of the finite element, and represents the internal region of the finite element, k is the heat transfer coefficient, T amb is the initial ambient temperature (℃), R s is the pipe burial depth, R w is the outer diameter of the pipe wall, K s is the soil thermal conductivity, the superscript n1 represents the nth time step, the subscript j represents the jth node, and the subscript wall represents the pipe wall.

[0058] The heat transfer coefficient k (W / ㎡·℃) and the initial ambient temperature Tamb (℃) are user-defined. Equations (10) and (11) are used to solve J q ·n in equation (10), and equation (12) is used to iteratively solve the ambient and pipe wall temperature vectors at the next time step for the node.

[0059] Through the solution of the equation, the flow field data of each node is obtained, including the pressure, flow rate, and temperature field distribution data of each node.

[0060] Step S14: Determine whether the flow field data of each node converges.

[0061] The convergence criterion is that the flow field data tends to be smooth and no longer changes.

[0062] Multiple monitoring points can be set, and when the fluctuation curves of the temperature, pressure and flow rate of the monitoring points tend to be flat and no longer change, it is determined that the flow field data of each node converges.

[0063] If the step S14 determines yes, the prediction result can be output; if the step S14 determines no, the step S15 is executed, and / or the step S16 is executed.

[0064] The step S15: modifying the current initial condition and boundary condition.

[0065] The initial condition and boundary condition of the coupling model are set as a mass flow node at the inlet of the pipeline and a pressure outlet node at the outlet of the pipeline.

[0066] For example, the initial condition and boundary condition can be set as follows: taking the pipeline material as low carbon steel, the thermal conductivity coefficient is 46.4 W / (m·℃), the pipe length is 10 km, the inlet temperature is 20℃, the initial ambient temperature is 15℃, the inlet flow rate is 0.9 m / s, and the outlet pressure is 1.2 MPa.

[0067] The initial condition or the boundary condition can be modified based on the tunicate swarm optimization algorithm (the process of searching for the best food source position by the tunicate swarm is analogous to parameter optimization).

[0068] The step S16: modifying the current grid density.

[0069] If only the step S15 is executed and the step S16 is not executed, after the step S15, the step S13 is executed; if only the step S16 is executed and the step S15 is not executed, after the step S16, the step S12 is executed; if both the step S15 and the step S16 are executed, after the execution, the step S12 is executed. The flow field data of each node is obtained by re-solution.

[0070] The simulation method for the flow characteristics of the liquid ammonia pipeline provided by the embodiment of the present application can more accurately calculate the distribution law of the physical field in the flow process of the pipeline transported liquid ammonia, the numerical simulation method used is more practical, the simulation accuracy can be greatly improved, reliable fluid mechanics information is provided for the design and safe and stable operation of the long-distance pipeline for transporting liquid ammonia, and the method has important significance in the field of ammonia storage and transportation technology.

[0071] In some embodiments, the error of the final output prediction result and the experimental result can also be compared and analyzed, and the initial conditions and boundary conditions of the model and / or the grid density can be further modified.

[0072] Referring to Figure 3 As shown in the figure, three monitoring points are selected, and comparison results of comparing the simulation results with the actual operation data are shown. The results show that the average relative error of the model in temperature prediction of different monitoring points is 3.6%, the maximum relative error is 4.3%, the average relative error of pressure prediction is 1.97%, and the maximum relative error is 2.5%, which meets the engineering actual calculation accuracy requirement, so the model has high accuracy.

[0073] Based on the inventive concept of the application, the embodiment of the application further provides a water and heat coupled liquid ammonia pipeline flow characteristic simulation device. Figure 4 As shown in the figure, the device comprises:

[0074] The coupling model establishing module 41 is configured to establish a liquid ammonia pipeline flow non-steady-state water and heat coupled model on the basis of the liquid ammonia pipeline physical model, and the model comprises an energy conservation equation, a momentum conservation equation, a continuity equation and a soil heat conduction equation.

[0075] The discretization module 42 is configured to perform two-dimensional finite element grid division and one-dimensional pipeline segment discretization processing on the coupling model according to the current grid density.

[0076] The flow field data solving module 43 is configured to solve the discretized coupling model based on the current initial conditions and boundary conditions to obtain flow field data of each node.

[0077] The judgment module 44 is configured to judge whether the flow field data of each node converges.

[0078] The correction module 45 is configured to modify the current initial conditions and boundary conditions and / or modify the current grid density if the judgment module 44 judges that the flow field data of each node does not converge, and the flow field data solving module 43 is further configured to re-solve to obtain the flow field data of each node.

[0079] As to the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be described in detail here.

[0080] Based on the inventive concept of the application, the embodiment of the application further provides a computer storage medium, wherein the computer storage medium stores computer executable instructions, and the computer executable instructions are executed by a processor to implement the water and heat coupled liquid ammonia pipeline flow characteristic simulation method.

[0081] Based on the inventive concept of the present invention, the present invention also provides a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the above-mentioned method for simulating the flow characteristics of a liquid ammonia pipeline by hydraulic-thermal coupling.

[0082] Unless otherwise specifically stated, terms such as processing, calculation, operation, determination, display, etc., may refer to the actions and / or processes of one or more processing or computing systems or similar devices that represent the manipulation and conversion of data representing physical (e.g., electronic) quantities within the registers or memory of the processing system into other data similarly representing physical quantities within the memory, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals can be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description can be represented by voltage, current, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0083] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process may be rearranged without departing from the scope of this disclosure. The appended method claims provide elements of various steps in an exemplary order and are not intended to limit the scope to the specific order or hierarchy described.

[0084] In the detailed description above, various features are combined together in a single embodiment to simplify this disclosure. This approach to disclosure should not be construed as reflecting an intention that embodiments of the claimed subject matter require more features than those stated in each claim. Rather, as reflected in the appended claims, the invention is presented with fewer features than all of the features in a single disclosed embodiment. Therefore, the appended claims are hereby clearly incorporated into the detailed description, wherein each claim stands alone as a preferred embodiment of the invention.

[0085] Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein can be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans can implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.

[0086] The steps of a method or algorithm described in connection with the embodiments disclosed herein can be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium can be integral to the processor. The processor and the storage medium can reside in an ASIC. The ASIC can reside in a user terminal. In the alternative, the processor and the storage medium can reside as discrete components in a user terminal.

[0087] For a software implementation, the techniques described herein can be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes can be stored in memory units and executed by processors. The memory unit can be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.

[0088] The above description includes one or more examples of the embodiments. Of course, not all possible combinations of components or method steps described above can be claimed as embodiments. One of ordinary skill in the art can recognize that modifications and variations of the embodiments described herein are possible and are within the scope of the present disclosure. It is therefore intended that the embodiments described herein be considered in all respects as only illustrative and not restrictive. That is, the embodiments described herein are intended to cover all alternatives, modifications and equivalents. In addition, while the term "comprise" has been used in the detailed description and claims, the word "comprise" and variations such as "comprises" or "comprising," when used in the specification and / or claims, can mean "consist of or "consist essentially of. Furthermore, the use of the term "or," as used in the detailed description and claims, is to be interpreted as a non-exclusive "or," unless explicitly indicated to the contrary.

Claims

1. A method for simulating the flow characteristics of liquid ammonia pipelines through hydraulic-thermal coupling, characterized in that, include: Based on the physical model of the liquid ammonia pipeline, an unsteady hydraulic-thermal coupling model of liquid ammonia pipeline flow is established, including the energy conservation equation, momentum conservation equation, continuity equation and soil thermal conductivity equation. Based on the current mesh density, the coupled model is subjected to two-dimensional finite element mesh generation and one-dimensional pipe segment discretization. Based on the current initial and boundary conditions, the discretized coupled model is solved to obtain the flow field data for each node; Determine whether the flow field data at each node has converged; If not, modify the current initial and boundary conditions, and / or modify the current mesh density; The flow field data for each node is obtained by resolving the problem.

2. The method as described in claim 1, characterized in that, The energy conservation equation is as follows: In equations (1)-(3), h is the field enthalpy, m is the fluid mass, U is the flow velocity, M is the mole fraction, the subscript i represents the i-th fluid, and the subscripts t and z represent time and space, respectively. S is the enthalpy source, Q wall Let p be the heat flow through the pipe wall, and T be the fluid pressure. l Where φ is the fluid temperature, φ is the fugacity coefficient, and C is the fluid temperature. p,l C is the specific heat capacity at constant pressure of the fluid. p0 C p1 C p2 C p3 and C p4 These represent different specific heat coefficients, where R is the ideal gas constant.

3. The method as described in claim 2, characterized in that, The momentum conservation equation is: In equation (4), H is the cross-sectional liquid holdup, ρ l Let ρ be the fluid density, v be the flow velocity, τ be the shear force between the fluid and the pipe wall, g be the gravitational acceleration, θ be the pipe inclination angle, and x be the pipe axial direction.

4. The method as described in claim 3, characterized in that, The continuity equation is:

5. The method as described in claim 4, characterized in that, The soil thermal conductivity equation is as follows: In equation (6), ρ s For soil density, c p,s λs is the soil's specific heat capacity at constant pressure, Ts is the soil temperature, λs is the soil's thermal conductivity, and y is the pipe's radial direction.

6. The method as described in claim 4, characterized in that, The solution to the discretized coupled model includes: The coupled model after simultaneous discretization is solved using the finite difference method.

7. The method as described in claim 6, characterized in that, The coupled model after simultaneous discretization includes: By simultaneously solving the discretized energy conservation equation and the soil heat conduction equation, the governing equation for the pipe temperature field is obtained as follows: J q =-λ s ▽T s (10) J q ·n=-k(T s -T amb ) (11) Among them, A l Let A be the cross-sectional area of ​​the fluid. s Let A be the cross-sectional area of ​​the soil. wall ρ is the cross-sectional area of ​​the pipe wall. wall c is the pipe wall density. p,wall T is the isobaric specific heat capacity of the pipe wall. wall Q is the pipe wall temperature. s For the heat flow through the soil, δ is an operator representing the change, J q Let T be the heat flux at the finite element boundary, k be the heat transfer coefficient, n be the unit normal vector at the finite element boundary, Ω represent the interior region of the finite element, and T be the heat flux at the finite element boundary. amb Let R be the initial value of the ambient temperature. s For the pipeline burial depth, R w K is the outer diameter of the pipe wall. s The soil thermal conductivity is represented by the superscript n1, which represents the nth time step, the subscript j, which represents the jth node, and the subscript wall, which represents the pipe wall.

8. The method as described in claim 1, characterized in that, The coupled model is meshed using two-dimensional finite element methods, including: The coupled model is meshed using a two-dimensional finite element method based on the Delaunay triangulation method.

9. The method according to any one of claims 1 to 8, characterized in that, The obtained flow field data for each node includes: The pressure, flow velocity, and temperature field distribution data of each node are obtained.

10. A hydraulically and thermally coupled device for simulating the flow characteristics of liquid ammonia pipelines, characterized in that, include: The coupling model building module is used to build an unsteady hydraulic-thermal coupling model of liquid ammonia pipe flow based on the physical model of the liquid ammonia pipeline, including the energy conservation equation, momentum conservation equation, continuity equation and soil heat conduction equation. The discretization module is used to perform two-dimensional finite element mesh generation and one-dimensional pipe segment discretization processing on the coupled model according to the current mesh density. The flow field data solving module is used to solve the discretized coupled model based on the current initial and boundary conditions to obtain the flow field data of each node. The judgment module is used to determine whether the flow field data of each node has converged; The correction module is used to modify the current initial conditions and boundary conditions, and / or modify the current mesh density if the judgment module determines that it is not true; the process data solving module is also used to re-solve to obtain the flow field data of each node.

11. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the hydraulic-thermal coupling liquid ammonia pipeline flow characteristic simulation method according to any one of claims 1 to 9.

12. A server, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the method for simulating the flow characteristics of a liquid ammonia pipeline by hydraulic-thermal coupling as described in any one of claims 1 to 9.

Images