CN116502373A

CN116502373A - Runner calculation grid generation method and device

Info

Publication number: CN116502373A
Application number: CN202310786762.8A
Authority: CN
Inventors: ***; 赵文强; 刘涛; 魏征; 郝帅
Original assignee: Shaanxi Aerospace Information Technology Co ltd
Current assignee: Shaanxi Aerospace Information Technology Co ltd
Priority date: 2023-06-30
Filing date: 2023-06-30
Publication date: 2023-07-28
Anticipated expiration: 2043-06-30
Also published as: CN116502373B

Abstract

The disclosure provides a method and a device for generating a runner calculation grid, relates to the technical field of computers, and can be applied to a scene of performance analysis of a designed axial flow impeller machine. The runner calculation grid generation method comprises the following steps: obtaining impeller design parameters of an axial flow impeller machine; determining a meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters; determining calculation grid characteristic data based on a meridian structure diagram; and determining a calculation grid corresponding to the single runner through calculation grid characteristic data. The technical scheme of the embodiment of the disclosure can improve the generation efficiency of the flow channel calculation grid of the axial flow impeller machine and improve the accuracy and stability of the generated calculation grid.

Description

Runner calculation grid generation method and device

Technical Field

The disclosure relates to the technical field of computers, in particular to a runner calculation grid generation method and a runner calculation grid generation device.

Background

The impeller machinery through-flow calculation is a numerical solution technology for analyzing the quasi-three-dimensional attribute of the air flow, in particular to an axial flow impeller machinery with complex flow channel structures such as an axial flow compressor, and the radial direction and the spanwise direction distribution condition of the air flow parameters in the flow channel can be observed through the through-flow calculation, and the realization of all functions is based on the generated flow channel calculation grid.

At present, the related axial flow impeller machine through-flow grid computing method mainly gives grid division size according to the geometric model of the axial flow impeller machine, but in the scheme, grid setting is needed again after the geometric model of the axial flow impeller machine is optimized and improved, so that the workload is large, and the generating efficiency of the flow channel computing grid is low.

It should be noted that the information disclosed in the above background section is only for enhancing understanding of the background of the present disclosure and thus may include information that does not constitute prior art known to those of ordinary skill in the art.

Disclosure of Invention

An object of an embodiment of the present disclosure is to provide a method for generating a runner computation grid, a device for generating a runner computation grid, an electronic device, and a computer readable storage medium, so as to effectively improve the efficiency of generating a runner computation grid of an axial flow impeller machine, reduce the workload required for generating a runner computation grid, and reduce the labor cost.

Other features and advantages of the present disclosure will be apparent from the following detailed description, or may be learned in part by the practice of the disclosure.

According to a first aspect of an embodiment of the present disclosure, there is provided a runner computation grid generating method, including:

Obtaining impeller design parameters of an axial flow impeller machine;

determining a meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters;

determining calculation grid characteristic data based on the meridian structure diagram;

and determining the computing grid corresponding to the single flow channel through the computing grid characteristic data.

According to a second aspect of the embodiments of the present disclosure, there is provided a flow channel calculation grid generating apparatus, including:

the impeller design parameter acquisition module is used for acquiring impeller design parameters of the axial flow impeller machine;

the meridian structure diagram determining module is used for determining a meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters;

the grid characteristic data determining module is used for determining and calculating grid characteristic data based on the meridian structure diagram;

and the runner calculation grid generation module is used for determining the calculation grid corresponding to the single runner through the calculation grid characteristic data.

According to a third aspect of embodiments of the present disclosure, there is provided an electronic device, comprising: a processor; and a memory having stored thereon computer readable instructions which when executed by the processor implement the runner computation grid generation method of any of the above.

According to a fourth aspect of embodiments of the present disclosure, there is provided a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements a runner computation grid generation method according to any one of the above.

The technical scheme provided by the embodiment of the disclosure can comprise the following beneficial effects:

according to the flow channel calculation grid generation method in the example embodiment of the disclosure, the impeller design parameters of the axial flow impeller machine can be obtained, then the meridian structure diagram of the axial flow impeller machine on a single flow channel can be determined according to the impeller design parameters, grid characteristic data corresponding to the calculation grid can be determined based on the meridian structure diagram, and further the calculation grid corresponding to the single flow channel can be determined through the calculation grid characteristic data. On the one hand, a meridian direction structure diagram can be obtained according to the impeller design parameters, and then a calculation grid corresponding to a single flow channel of the axial-flow impeller machine can be directly generated according to the meridian direction structure diagram, and the generation of the flow channel calculation grid can be automatically realized only by inputting the impeller design parameters in the generation process, so that the generation efficiency of the calculation grid is effectively improved, the manual participation is reduced, and the labor cost is reduced; on the other hand, after the design optimization or improvement of the axial flow impeller machine, the calculation grids are not required to be re-divided according to the geometric model of the axial flow impeller machine, and the re-generation of the calculation grids can be completed by only adjusting the impeller design parameters of the axial flow impeller machine, so that the generation efficiency of the calculation grids is further improved, the workload required by the calculation grid division is reduced, the time required by the calculation grid division is shortened, and the design efficiency of the axial flow impeller machine is further improved; on the other hand, the calculation grids are generated through the meridian structure diagram, the standardization and standardization of the dividing process of the calculation grids are realized, the problem that the divided calculation grids have large difference due to different experiences of different designers can be effectively avoided, the accuracy and stability of the generated calculation grids are effectively improved, and the robustness of the design analysis result of the axial-flow impeller machine is ensured.

It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.

Drawings

The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description, serve to explain the principles of the disclosure. It will be apparent to those of ordinary skill in the art that the drawings in the following description are merely examples of the disclosure and that other drawings may be derived from them without undue effort.

Fig. 1 is a schematic diagram of a system architecture of an exemplary application environment to which a runner computation grid generation method and apparatus of an embodiment of the present disclosure may be applied.

Fig. 2 schematically illustrates a flow diagram of a runner computation grid generation method according to some embodiments of the present disclosure.

Fig. 3 schematically illustrates a flow diagram for generating a meridional block diagram according to some embodiments of the present disclosure.

Fig. 4 schematically illustrates a structural schematic of sections of a meridional block diagram according to some embodiments of the present disclosure.

Fig. 5 schematically illustrates a schematic diagram of segment sizes according to some embodiments of the present disclosure.

Fig. 6 schematically illustrates a schematic diagram of constructing a meridian curve in accordance with some embodiments of the present disclosure.

Fig. 7 schematically illustrates a flow diagram for generating a computational grid according to some embodiments of the present disclosure.

FIG. 8 schematically illustrates a flow diagram for generating a computational grid according to further embodiments of the present disclosure.

Fig. 9 schematically illustrates a structural schematic of a set grid node orientation according to some embodiments of the present disclosure.

Fig. 10 schematically illustrates a schematic view of a tangential thickness of a airfoil according to some embodiments of the disclosure.

FIG. 11 schematically illustrates a structural schematic of a first blade side grid node coordinate point set in accordance with some embodiments of the present disclosure.

FIG. 12 schematically illustrates a structural schematic of a second blade side grid node coordinate point set in accordance with some embodiments of the present disclosure.

Fig. 13 schematically illustrates a structural schematic of a computational grid according to some embodiments of the present disclosure.

Fig. 14 schematically illustrates a structural schematic of a grid-densified computational grid according to some embodiments of the present disclosure.

Fig. 15 schematically illustrates a schematic diagram of a runner computation grid generating device according to some embodiments of the present disclosure.

Fig. 16 schematically illustrates a structural schematic diagram of a computer system of an electronic device according to some embodiments of the present disclosure.

Fig. 17 schematically illustrates a schematic diagram of a computer-readable storage medium according to some embodiments of the present disclosure.

In the drawings, the same or corresponding reference numerals indicate the same or corresponding parts.

Detailed Description

Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the present specification. Rather, they are merely examples of apparatus and methods consistent with some aspects of the present description as detailed in the accompanying claims.

The terminology used in the description presented herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the description. As used in this specification and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.

It should be understood that although the terms first, second, third, etc. may be used in this specification to describe various information, these information should not be limited to these terms. These terms are only used to distinguish one type of information from another. For example, the first information may also be referred to as second information, and similarly, the second information may also be referred to as first information, without departing from the scope of the present description. The word "if" as used herein may be interpreted as "at … …" or "at … …" or "responsive to a determination", depending on the context.

Example embodiments will now be described more fully with reference to the accompanying drawings. However, the exemplary embodiments may be embodied in many forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of the example embodiments to those skilled in the art.

Furthermore, the described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a thorough understanding of embodiments of the disclosure. One skilled in the relevant art will recognize, however, that the disclosed aspects may be practiced without one or more of the specific details, or with other methods, components, devices, steps, etc. In other instances, well-known methods, devices, implementations, or operations are not shown or described in detail to avoid obscuring aspects of the disclosure.

Moreover, the drawings are only schematic illustrations and are not necessarily drawn to scale. The block diagrams depicted in the figures are merely functional entities and do not necessarily correspond to physically separate entities. That is, the functional entities may be implemented in software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor devices and/or microcontroller devices.

As shown in fig. 1, the system architecture 100 may include one or more of a desktop computer 101, a portable computer 102, a smart phone 103, and other terminal devices, a network 104, and a server 105. The network 104 is the medium used to provide communication links between the terminal devices and the server 105. The network 104 may include various connection types, such as wired, wireless communication links, or fiber optic cables, among others. The terminal device may be a variety of electronic devices having data processing capabilities with a display screen for providing a user with an interface for inputting impeller design parameters and a computing grid for displaying an axial flow impeller machine, including but not limited to desktop computers, portable computers, smartphones, and the like as described above. It should be understood that the number of terminal devices, networks and servers in fig. 1 is merely illustrative. There may be any number of terminal devices, networks, and servers, as desired for implementation. For example, the server 105 may be a server cluster formed by a plurality of servers.

The method for generating the runner computation grid provided by the embodiment of the disclosure is generally executed by a terminal device, and accordingly, the runner computation grid generating device is generally arranged in the terminal device. However, it is easily understood by those skilled in the art that the method for generating the runner computation grid provided in the embodiment of the present disclosure may also be performed by the server 105, and accordingly, the runner computation grid generating device may also be provided in the server 105, which is not particularly limited in the present exemplary embodiment.

Further, it should be understood that the runner computation grid generation method of the embodiments of the present disclosure may be configured as a software module. In some implementations, the runner computation grid generation scheme of the present disclosure may be deployed separately to enable direct generation of a single runner computation grid of an axial flow turbomachine by inputting impeller design parameters of the axial flow turbomachine. In other implementation scenarios, the flow channel computing grid generation scheme of the present disclosure may be deployed in other software, as a functional module of the software, for example, in design analysis software of an axial flow impeller machine, and the application mode of the flow channel computing grid generation method is not particularly limited in the present disclosure.

Computational grid generation techniques have evolved into an important component of Computational Fluid Dynamics (CFD) and computational heat transfer chemistry (NHT). In computational fluid dynamics, the geometry of discrete points distributed in a flow field according to a certain rule is called computational grid, and the process of generating these nodes is called computational grid generation. The computational grid generation is a tie connecting a geometric model and a numerical algorithm, the geometric model can be solved numerically only when being divided into computational grids with certain standards, in general, the denser the computational grid division is, the more accurate the result is, but the larger the required computational load is, and the more time is consumed.

In the related art, when the grid is calculated by the axial flow impeller machine, the grid dividing size is mainly set according to the geometric model of the axial flow impeller machine, but in the scheme, grid setting is needed again after the geometric model of the axial flow impeller machine is optimized and improved, the workload is large, and the generation efficiency of the runner calculation grid is low.

Based on one or more problems in the related art, in the present exemplary embodiment, a flow channel calculation grid generation method is provided first, and a detailed description will be given below taking a terminal device to execute the method as an example. Fig. 2 schematically illustrates a flow diagram of a runner computation grid generation method according to some embodiments of the present disclosure. Referring to fig. 2, the flow channel calculation grid generation method may include the steps of:

Step S210, obtaining impeller design parameters of an axial flow type impeller machine;

step S220, determining a meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters;

step S230, determining calculation grid characteristic data based on the meridian structure diagram;

step S240, determining a computing grid corresponding to the single flow channel according to the computing grid feature data.

According to the flow channel calculation grid generation method in the embodiment, on one hand, a meridian direction structure diagram can be obtained according to the impeller design parameters, and further calculation grids corresponding to single flow channels of the axial flow impeller machine can be directly generated according to the meridian direction structure diagram, the generation of the flow channel calculation grids can be automatically realized only by inputting the impeller design parameters in the generation process, the generation efficiency of the calculation grids is effectively improved, the manual participation is reduced, and the labor cost is reduced; on the other hand, after the design optimization or improvement of the axial flow impeller machine, the calculation grids are not required to be re-divided according to the geometric model of the axial flow impeller machine, and the re-generation of the calculation grids can be completed by only adjusting the impeller design parameters of the axial flow impeller machine, so that the generation efficiency of the calculation grids is further improved, the workload required by the calculation grid division is reduced, the time required by the calculation grid division is shortened, and the design efficiency of the axial flow impeller machine is further improved; on the other hand, the calculation grids are generated through the meridian structure diagram, the standardization and standardization of the dividing process of the calculation grids are realized, the problem that the divided calculation grids have large difference due to different experiences of different designers can be effectively avoided, the accuracy and stability of the generated calculation grids are effectively improved, and the robustness of the design analysis result of the axial-flow impeller machine is ensured.

Next, a flow channel calculation grid generation method in the present exemplary embodiment will be further described.

In step S210, an impeller design parameter of the axial flow type impeller machine is obtained.

In an exemplary embodiment of the present disclosure, an axial flow type impeller machine refers to an impeller machine that flows in an axial direction through a fluid medium, for example, the axial flow type impeller machine may be a one-stage axial flow compressor or a multi-stage axial flow compressor, and of course, the axial flow type impeller machine may also be an axial flow water turbine, an axial flow compressor, an axial flow mixing pump, an axial flow water pump, or the like, and the exemplary embodiment is not particularly limited as to the type of the axial flow type impeller machine.

The impeller design parameter refers to geometric dimension data required in the design process of the axial flow impeller machine, for example, the impeller design parameter may be length data of each part of the axial flow impeller machine, such as a casing molded line length, a hub molded line length, a blade section length, a stator blade section length, a blade plate length, etc., and may also be height data of each part, such as a runner inlet height, a runner outlet height, a blade inlet and outlet height, a stator blade inlet and outlet height, etc., and of course, the impeller design parameter may also be a connection angle of each part surface connection, etc., and the type and number of the impeller design parameter are not limited in particular in this example embodiment.

The input impeller design parameters can be obtained through a pre-provided input interface, the impeller design parameters can also be obtained through directly analyzing the geometric model of the designed axial-flow impeller machine, and of course, the axial-flow impeller machine in a real scene can also be scanned through three-dimensional imaging equipment, so that the corresponding impeller design parameters are obtained.

In step S220, a radial structure of the axial flow impeller machine on the single flow channel is determined according to the impeller design parameters.

In an example embodiment of the present disclosure, the meridian structure diagram refers to a feature structure diagram of an axial flow type impeller machine in the meridian direction, for example, the meridian structure diagram may be a partial structure diagram of the axial flow type impeller machine in the meridian direction only including an external molded line, or may be an integral structure diagram of the axial flow type impeller machine in the meridian direction, and the content of the meridian structure diagram is not particularly limited in the example embodiment.

Because the axial flow impeller machine has symmetry in blade design distribution, and further the axial flow impeller machine has symmetry in flow channel distribution, the integral flow channel calculation grid of the axial flow impeller machine can be obtained only by obtaining the calculation grid of the axial flow impeller machine on a single flow channel, and for convenience in calculation and explanation, the embodiment of the disclosure takes the calculation grid for generating the single flow channel as an example.

In step S230, calculated grid feature data is determined based on the meridian structure map.

In an exemplary embodiment of the present disclosure, the calculation grid feature data refers to feature data describing the positions of grid nodes required for generating a calculation grid, for example, the calculation grid feature data may be calculation grid node coordinates determined based on a meridian structure diagram, or may be a spatial position relationship between grid nodes, and of course, the calculation grid feature data may be other types of feature data describing the positions of grid nodes, which is not particularly limited in this exemplary embodiment.

For example, the number of grid nodes to be set on the meridian structure chart can be determined according to the blade size data on the meridian structure chart, and then grid node setting is performed on the meridian structure chart through the determined number of grid nodes, so that calculated grid characteristic data such as calculated grid node coordinates, spatial position relations among grid nodes and the like can be obtained; of course, the determination of the grid node positions may also be implemented by other means on a meridian structure map, which is not particularly limited in this example embodiment.

In step S240, a computing grid corresponding to the single flow channel is determined according to the computing grid characteristic data.

In an example embodiment of the present disclosure, after the calculation grid feature data is obtained, a calculation grid node position in a three-dimensional space may be determined according to characteristics of the axial flow type impeller machine and the calculation grid feature data, and then a calculation grid of a single flow channel corresponding to the axial flow type impeller machine may be drawn through the calculation grid node position in the three-dimensional space, and further a global calculation grid corresponding to the axial flow type impeller machine may be generated according to the calculation grid of the single flow channel.

The cuboid computing grid can be generated through computing grid feature data, the hexahedral computing grid can also be generated through computing grid feature data, of course, any shape computing grid such as prismatic grid, pyramid grid and the like can also be generated through computing grid feature data, and specifically, the self-defining setting can be performed according to the design characteristics of the axial flow impeller machine, and the shape of the generated computing grid is not particularly limited in the embodiment.

The radial structure diagram is obtained through the impeller design parameters, and then the calculation grid corresponding to the single flow channel of the axial flow impeller machine can be directly generated according to the radial structure diagram, the generation of the flow channel calculation grid can be automatically realized only by inputting the impeller design parameters in the generation process, the generation efficiency of the calculation grid is effectively improved, the manual participation is reduced, and the labor cost is reduced; meanwhile, after design optimization or improvement of the axial flow impeller machine, the calculation grids are not required to be divided again according to the geometric model of the axial flow impeller machine, and the calculation grids can be regenerated by adjusting impeller design parameters of the axial flow impeller machine, so that the generation efficiency of the calculation grids is further improved, the workload required by calculation grid division is reduced, the time required by calculation grid division is shortened, and the design efficiency of the axial flow impeller machine is further improved.

Next, the steps of step S110 to step S140 will be described in detail.

In an exemplary embodiment of the present disclosure, the impeller design parameters may include blade parameters and non-blade parameters, where the blade parameters refer to data related to the blade structure on the axial flow impeller machine, and for example, the blade parameters may include the height of the blade, the length of the blade, the thickness of the blade edge plate, the number of blades, the spacing of the blades, etc., and of course, the required blade parameters are not the same as the required blade parameters in different application scenarios, and the exemplary embodiment does not make a special limitation on the type of the blade parameters. Taking an axial-flow compressor as an example, the blade parameters may be blade parameters related to a movable blade or blade parameters related to a stationary blade.

The non-blade parameters refer to related data which do not belong to the blade structure on the axial-flow impeller machine, for example, the non-blade parameters can comprise a casing molded line size, a hub molded line size, a runner inlet height, a runner outlet height and the like, and similarly, the required non-blade parameters are different according to different application scenes, and the type of the blade parameters is not particularly limited in the embodiment.

The radial structure diagram of the axial flow type impeller machine on a single flow channel can be determined according to the impeller design parameters through the steps in fig. 3, and referring to fig. 3, the method specifically can include:

Step S310, determining a blade section of the axial flow impeller machine on a meridian plane of a single flow channel according to the blade parameters;

step S320, determining a non-blade section of the axial flow impeller machine on a meridian plane of a single flow channel according to the non-blade parameters;

step S330, determining a meridian structure diagram of the axial flow impeller machine on the single flow channel through the blade section and the non-blade section.

The blade section refers to an area of the axial flow type impeller machine, which belongs to a blade structure on a meridian plane of a single flow channel, and the blade section can comprise a movable blade section and a stationary blade section, and can be drawn according to blade parameters.

The non-blade section refers to an area of the axial flow impeller machine, which does not belong to a blade structure, on the meridian plane of the single flow channel, and the non-blade section on the meridian plane of the single flow channel can be drawn according to non-blade parameters.

After the blade section and the non-blade section are obtained, the blade section and the non-blade section can be spliced according to the position relation between the blade section and the non-blade section, so that a meridian structure diagram of the axial flow impeller machine on a single flow channel is obtained.

Referring to fig. 4, a meridional structure diagram of an axial flow type impeller machine on a single flow channel may be determined according to impeller design parameters, wherein the impeller design parameters may include blade parameters and non-blade parameters, wherein the blade parameters may include a blade length 401, a blade height 402, a vane length 403, a vane height 404, etc. of the axial flow type impeller machine in a meridional direction, and the non-blade parameters may include a casing profile length 405, a hub profile length 406, a flow channel inlet height 407, a flow channel outlet height 408, etc. of the axial flow type impeller machine in a meridional direction.

Specifically, the vane segments of the axial flow turbomachine on the meridian plane of the single flow path may be determined based on vane parameters, for example, the vane segments 410 and vane segments 412 may be obtained, the non-vane segments of the axial flow turbomachine on the meridian plane of the single flow path may be determined based on non-vane parameters, for example, the flow path inlet segment 409, the cavity segments 411 between the vanes and vanes, and the flow path outlet segment 413 may be obtained. After all the blade sections and the non-blade sections are obtained, the radial structure diagram of the axial flow type impeller machine on a single flow channel can be obtained by splicing according to the position relation between the blade sections and the non-blade sections.

Each section of the meridian structure diagram is determined through the impeller design parameters, and then the meridian structure diagram of the axial-flow impeller machine on a single flow channel is obtained through splicing of each section, so that division of each structure on the meridian structure diagram is realized, drawing of calculation grids of the axial-flow impeller machine can be realized according to the characteristics of different sections, and the accuracy of the drawn calculation grids is ensured.

In an example embodiment of the present disclosure, the computing grid feature data may include meridian grid node coordinates, where meridian grid node coordinates refer to coordinates corresponding to grid nodes of a meridian up computing grid, and determining the computing grid feature data based on a meridian structure map may be achieved by:

the total meridian node quantity can be determined according to the impeller design parameters of each section in the meridian structure diagram; the total meridian node quantity is distributed, and the zone meridian node quantity corresponding to each zone is obtained; constructing a meridian curve corresponding to a meridian structure diagram; and carrying out equidistant interpolation on the meridian curve according to the number of the regional meridian nodes to obtain the meridian grid node coordinates.

The total meridian node number refers to the total number of grid nodes which should be set on the meridian structure diagram, and can be calculated according to the impeller design parameters of each section in the meridian structure diagram. The number of zone meridian nodes refers to the number of grid nodes which should be set in each zone on a meridian structure diagram, and the number of zone meridian nodes corresponding to each zone can be obtained by distributing the total number of meridian nodes.

The method can construct a meridian curve corresponding to the meridian structure diagram, and perform equidistant interpolation on the meridian curve according to the calculated number of the regional meridian nodes, and can obtain the meridian grid node coordinates or the position relationship among the meridian grid nodes by establishing a coordinate system corresponding to the meridian structure diagram.

Alternatively, the average size corresponding to each section can be calculated according to the impeller design parameters of each section in the meridian structure diagram, then the preset number of spanwise grid nodes can be obtained, and the total meridian node number can be determined according to the average size and the number of spanwise grid nodes.

Referring to fig. 5, taking a bucket segment as an example, any one of the radial block diagrams may include a segment inlet height 501, a segment outlet height 502, a segment casing length 503, and a segment hub length 504, and the average size may include an average height and an average length, for example, the average height may be calculated by the relation (1), and the average length may be calculated according to the relation (2):

H(n)=[h(n1)+h(n2)]/2 （1）

L(n)=[l(n1)+l(n2)]/2 （2）

where H (n) may represent the average height of the nth section, H (n 1) may represent the section inlet height 501 of the nth section, H (n 2) may represent the section outlet height 502 of the nth section, l (n) may represent the average length of the nth section, l (n 1) may represent the section casing length 503 of the nth section, and l (n 2) may represent the section hub length 504 of the nth section.

After obtaining the average size corresponding to each segment, the total meridian node number may be determined according to the average size and the number of spanwise grid nodes, for example, the total meridian node number may be calculated by the relationships (3) to (5):

J = K×(L/ H)×f （3）

L=∑L(n) （4）

H=∑H(n)/n （5）

wherein J can represent the total meridian node number, and the result is taken up by an integer; k can represent the number of the spanwise grid nodes and can be preset, for example, the value of K can be 28, and of course, the value of K can be self-defined according to the actual application situation, and the embodiment is not limited to this; f can represent a scale factor, generally a constant greater than 1, and the total meridian node quantity reaches a working condition expected value through the set f, and can be specifically set in a self-defined manner according to an actual application scene, wherein the value of f is not particularly limited; l may represent the average length of the full segment, L (n) may represent the average length of the nth segment, H may represent the average height of the full segment, H (n) may represent the average height of the nth segment, n may represent the number of segments in the meridional structure map, and Σ may represent the summation operation.

Optionally, the total meridian node number can be allocated by the following steps to obtain the zone meridian node number corresponding to each zone:

The estimated meridian node number corresponding to each segment may be determined based on the average size and the total meridian node number; and acquiring preset multiple grid coefficients, and determining the number of zone meridian nodes corresponding to each zone according to the estimated number of the meridian nodes and the multiple grid coefficients.

The estimated meridian node number refers to an estimated value obtained by distributing the total meridian node number according to the average size, the regional meridian node number refers to the actual meridian grid node number obtained by adjusting the estimated meridian node number of each region according to the requirement, namely the meridian grid node number finally used for calculating grid division.

For example, the allocation of the total meridian node number based on the average size may be implemented by the relation (6) to obtain the estimated meridian node number corresponding to each segment:

J(n)= (J-1)×[L(n) / L] （6）

where J (n) may represent the estimated meridian node number for the nth zone, and the result is an upward integer, J may represent the total meridian node number, L (n) may represent the average length of the nth zone, and L may represent the average length of the full zone.

Further, accurate adjustment of the estimated meridian node number according to the multiple grid coefficients can be achieved through a relation (7) to obtain the segment meridian node number corresponding to each segment:

J(nλ)= [J(n) / K(λ)]×K(λ)+1 （7）

Where J (nλ) may represent the number of zone meridian nodes of the nth zone, J (n) may represent the number of estimated meridian nodes of the nth zone, K (λ) may represent multiple grid coefficients, and K (λ) generally takes the power of λ of 2, where the range of λ values may be 1,2,3, and of course, other forms of multiple grid coefficients may be set, and the present exemplary embodiment is not limited thereto.

Optionally, after the number of the section meridian nodes is obtained, a meridian curve corresponding to the meridian structure diagram can be constructed, and then equidistant interpolation can be performed on the meridian curve according to the number of the section meridian nodes, so as to obtain the meridian grid node coordinates.

The meridian curve can be obtained by constructing quasi-orthogonal lines in each section of the meridian structure diagram, equidistant interpolation points on the quasi-orthogonal lines and sequentially connecting the equidistant interpolation points according to the direction of the flow channel; of course, the meridian curve may be obtained by sequentially connecting the non-uniform interpolation points according to the distance from the boundary layer (i.e., the hub molded line and the casing molded line) and the flow channel direction, which is not particularly limited in this exemplary embodiment.

Referring to fig. 6, in the meridian structure chart, a quasi-orthogonal line 601 in each section may be determined, and N points, that is, equidistant interpolation points 602, may be equidistantly interpolated on the quasi-orthogonal line 601 in each section, where N may be a positive integer greater than 1, and generally may be custom set according to an application situation, or may be randomly generated according to a random generation range, and the value of N is not limited in this example embodiment.

The spline curves can be used for sequentially connecting N points to form N meridian curves (the N meridian curves can comprise hub molded lines and casing molded lines), equidistant interpolation can be carried out on the meridian curves according to the number of the meridian nodes of the sections, the meridian grid nodes corresponding to the sections can be obtained, and then the meridian grid node coordinates (z, r) corresponding to the sections can be obtained.

Obtaining sector meridian nodes corresponding to each sector through impeller design parameters of each sector in a meridian structure diagram, further carrying out equidistant interpolation on meridian curves according to the number of the sector meridian nodes to obtain meridian grid node coordinates, effectively ensuring that the obtained meridian grid nodes accord with the design of the axial-flow impeller machine, ensuring that the divided calculation grids match with flow channels of the axial-flow impeller machine, and improving the accuracy of the generated flow channel calculation grids; and the standardization and standardization of the generation of the calculation grid nodes are realized, the accuracy of the determined grid nodes is ensured, and meanwhile, the generation efficiency of the calculation grid is improved.

In an example embodiment of the present disclosure, determining a computing grid corresponding to a single runner according to computing grid feature data may be implemented through the steps in fig. 7, and referring to fig. 7, may specifically include:

step S710, determining polar angle coordinate values according to the meridian grid node coordinates;

step S720, determining grid node column coordinates through the meridian grid node coordinates and the polar angle coordinate values;

step S730, generating a calculation grid corresponding to the single flow channel based on the grid node cylindrical coordinates.

The polar angle coordinate value refers to a coordinate value corresponding to the meridian grid node in the theta direction in the cylindrical coordinate system, namely, the theta coordinate value, the grid node cylindrical coordinate refers to a coordinate value corresponding to the grid node when the grid node is converted into the cylindrical coordinate system, the grid node cylindrical coordinate can be constructed through the meridian grid node coordinate and the polar angle coordinate value, and further after the grid node cylindrical coordinate is obtained, the corresponding grid node can be connected based on the grid node cylindrical coordinate, so that the calculation grid corresponding to the axial-flow impeller machinery single flow channel is obtained.

Alternatively, the polar coordinate values may include a blade segment polar coordinate value and a non-blade segment polar coordinate value, where the blade segment polar coordinate value refers to a polar coordinate value corresponding to a meridian grid node coordinate in the blade segment, and the non-blade segment polar coordinate value refers to a polar coordinate value corresponding to a meridian grid node coordinate in the non-blade segment.

The coordinate system conversion can be carried out on the noon grid node coordinates, the space grid node coordinates are determined, and the polar angle coordinate values of the blade segments are calculated according to the space grid node coordinates. The space grid node coordinates refer to grid node coordinates obtained by converting meridian grid node coordinates into a cartesian coordinate system, and of course, the space grid node coordinates may also be grid node coordinates obtained by converting meridian grid node coordinates into other types of space coordinate systems, and the type of the space coordinate system is not particularly limited in this example embodiment.

For example, by converting the meridian grid node coordinates (z, r) into a cartesian coordinate system, the spatial grid node coordinates (x, y, z) of the blade section on the pressure surface and the suction surface can be obtained, and then the blade segment polar angle coordinate value can be calculated by the spatial grid node coordinates (x, y, z) of the meridian grid node in the cartesian coordinate system, for example, the blade segment polar angle coordinate value can be calculated by the relation (8) and the relation (9):

θ(p) = acrtan[y(p)/x(p)] （8）

θ(s)= acrtan[y(s)/x(s)] （9）

wherein θ (p) may represent θ coordinate values of grid nodes on the blade section pressure surface in a cylindrical coordinate system, θ(s) may represent θ coordinate values of grid nodes on the blade section suction surface in a cylindrical coordinate system, x (p) and y (p) may represent x and y coordinate values of grid nodes on the blade section pressure surface in a Cartesian coordinate system, respectively, and x(s) and y(s) may represent x and y coordinate values of grid nodes on the blade section suction surface in a Cartesian coordinate system, respectively.

Alternatively, the polar angle coordinate value of the non-blade segment corresponding to each grid node on the non-blade segment can be determined by setting a constant polar angle coordinate value (θ coordinate value), that is, the polar angle coordinate value of the non-blade segment can be determined by the positional relationship of the non-blade segment and the blade segment in the meridian structure diagram and the polar angle coordinate value of the blade segment. For example, the θ coordinate value of the non-blade section preceding the blade section in the flow path direction may be set as the θ coordinate value of the blade section inlet, and the θ coordinate value of the non-blade section following the blade section in the flow path direction may be set as the θ coordinate value of the blade section outlet. Of course, the polar coordinate values corresponding to each grid node on the non-blade segment may also be determined by other methods, and the method of calculating the polar coordinate values of the non-blade segment according to the present exemplary embodiment is not particularly limited.

After obtaining the polar angle coordinate values of the blade segments and the polar angle coordinate values of the non-blade segments, grid node cylindrical coordinates (z, r, theta) can be constructed according to meridian grid node coordinates (z, r) and the polar angle coordinate values theta, and then a calculation grid corresponding to a single flow channel can be generated through the grid node cylindrical coordinates (z, r, theta).

In an example embodiment of the present disclosure, the polar angle coordinate values may include a first polar angle coordinate value and a second polar angle coordinate value, wherein the first polar angle coordinate value may refer to a polar angle coordinate value corresponding to a grid node on a pressure side in the polar angle coordinate values of the blade segment, and the second polar angle coordinate value may refer to a polar angle coordinate value corresponding to a grid node on a suction side in the polar angle coordinate values of the blade segment. It should be understood that "first" and "second" herein are used only to distinguish between the polar coordinate values corresponding to the grid nodes on the pressure side and the polar coordinate values corresponding to the grid nodes on the suction side, and are not limited in any particular way, the first polar coordinate value may also represent the polar coordinate value corresponding to the grid nodes on the suction side, and the second polar coordinate value may also represent the polar coordinate value corresponding to the grid nodes on the pressure side, and should not be limited in any way to this exemplary embodiment.

The generation of the calculation grid corresponding to the single flow channel based on the grid node cylindrical coordinates can be realized through the steps in fig. 8, and referring to fig. 8, the method specifically may include:

step S810, setting a grid node direction, and determining a first blade side grid node coordinate point set of the grid node column coordinates in the grid node direction;

Step S820, determining the tangential thickness of the blade profile according to the first polar angular coordinate value and the second polar angular coordinate value;

step S830, converting the first blade side grid node coordinate point set based on the tangential thickness of the blade profile, to obtain a second blade side grid node coordinate point set;

step S840, connecting coordinate points between the first blade side grid node coordinate point set and the second blade side grid node coordinate point set, and generating a calculation grid corresponding to the single flow channel.

The grid node direction refers to a unit vector preset to describe a positional direction relationship between grid nodes, for example, a certain angular point position in a single-channel fluid domain of an axial-flow impeller machine may be used as a starting point of the grid node direction, and unit vectors in different directions are specified by the starting point to obtain the grid node direction corresponding to the single-channel fluid domain. By determining the directions of the grid nodes, the position relation among the grid nodes can be described, and the subsequent generation of the calculation grid is facilitated.

Referring to fig. 9, for a single-channel fluid field 901 corresponding to an axial-flow type impeller machine, a certain corner point in the single-channel fluid field 901 may be used as an initial origin of a grid node direction, where i may represent a grid node order value distributed along a circumferential direction, j may represent a grid node order value distributed along a meridian direction, k may represent a grid node order value distributed along a spanwise direction, and initial positions of i, j and k are all 0, and further, positions of any grid node in the grid may be calculated by i, j and k, for example, positions of any grid node may be represented as v [ i ] [ j ] [ k ].

The first blade side grid node coordinate point set refers to a coordinate point set corresponding to grid nodes on the pressure surface side of the runner in the grid node direction, and the second blade side grid node coordinate point set refers to a coordinate point set corresponding to grid nodes on the suction surface side of the runner in the grid node direction. It should be understood that "first" and "second" herein are used only to distinguish the coordinate point set corresponding to the grid node on the pressure side of the flow path from the coordinate point set corresponding to the grid node on the suction side of the flow path, and are not intended to have any special meaning, and should not be construed as limiting the present exemplary embodiment.

Referring to fig. 10, for a blade 1001, a blade profile tangential thickness 1002 may be determined by a first polar angular coordinate value and a second polar angular coordinate value, for example, the blade profile tangential thickness 1002 may be determined by the relation (10):

dT= [(θ(p)×r- θ(s)×r] （10）

where dT may represent the tangential thickness 1002 of the airfoil of the blade 1001, θ (p) may represent a first polar coordinate value, i.e., a polar coordinate value corresponding to the grid node on the pressure side of the blade segment polar coordinate values, θ(s) may represent a second polar coordinate value, i.e., a polar coordinate value corresponding to the grid node on the suction side of the blade segment polar coordinate values, and r may represent an r coordinate value in the grid node cylindrical coordinates (z, r, θ).

After the grid node column coordinates corresponding to the pressure surface side of the blade section are obtained, the grid node column coordinates corresponding to the pressure surface side of the non-blade section can be obtained by extension, and the grid node column coordinates corresponding to the pressure surface side of the blade section and the non-blade section are converted to the grid node direction, so that a first blade side grid node coordinate point set v [0] [ j ] [ k ] is obtained.

The first blade side grid node coordinate point set v [0] [ j ] [ k ] can be converted based on the blade tangential thickness dT to obtain the second blade side grid node coordinate point set v [1] [ j ] [ k ], and in the conversion process, θ (p) in the first blade side grid node coordinate point set is mainly converted into θ (s ') in the second blade side grid node coordinate point set in the same grid node direction, for example, conversion between θ (p) and θ (s') can be achieved through a relational expression (11):

θ(s')=θ(p)+r×(2π / b) – dT （11）

wherein θ (s') may represent a polar coordinate value in the second blade side grid node coordinate point set, that is, a coordinate point set corresponding to the grid node on the suction side of the flow channel, θ (p) may represent a polar coordinate value in the first blade side grid node coordinate point set, r may represent an r coordinate value in grid node cylindrical coordinates (z, r, θ), b may represent the number of blades corresponding to the axial flow impeller machine, and dT may represent a tangential thickness of the blade profile of the blade.

Referring to fig. 11, 12 and 13, a first set of blade-side grid node coordinates 1102 in the flow field 1101 may be converted by the relation (11), that is, v [0] [ j ] [ k ] is converted according to the tangential thickness of the blade profile, to obtain a second set of blade-side grid node coordinates 1202 in the flow field 1201, that is, v [1] [ j ] [ k ], and may further pass through the first set of blade-side grid node coordinates 1102 and the second set of blade-side grid node coordinates 1202.

After the first blade side grid node coordinate point set 1102 and the second blade side grid node coordinate point set 1202 are obtained, coordinate points corresponding to the first blade side grid node coordinate point set 1102 and the second blade side grid node coordinate point set 1202 can be connected, and a calculation grid 1301 corresponding to a single flow channel of the axial flow type impeller machine can be obtained.

In an example embodiment of the present disclosure, the grid densification process may be performed on the computing grid 1301 obtained in fig. 13, and since the grid densification needs to be performed near the hub and the casing in the practical application process, the purpose of this is that, because a boundary layer exists near the wall surface in the flow channel, by encrypting the grid nodes near the wall surface, the hydrodynamic characteristics of the fluid flow may be reflected more truly.

The computational grid may be grid-densified by:

the method comprises the steps that a preset growth factor and the number of the spanwise grid nodes can be obtained, and the grid scale of each spanwise grid layer is determined according to the length of a spanwise curve of a calculated grid, the growth factor and the number of the spanwise grid nodes; and interpolating the spanwise curve through the grid scale to obtain spanwise interpolation point coordinates, and further generating a grid densely processed calculation grid based on the spanwise interpolation point coordinates.

For example, the grid scale of the first spanwise grid layer may be calculated by relation (12):

S=δ ₁ (1+a+a ² +···+a ^k/2 )×2 （12）

wherein S can represent the length, delta of the spanwise curve corresponding to any spanwise curve ₁ The grid scale of the first spanwise grid layer can be represented, a can represent a preset growth factor, generally a constant larger than 1, the growth factor can be set in a self-defined mode according to an actual application scene, and the value of a is not limited in particular; k may represent a spanwise mesh node.

The grid scale delta of the first spanwise grid layer of any spanwise curve can be determined according to the relation (12) ₁ Thereafter, the grid dimension delta of the grid layer can be expanded through the first layer ₁ And growth factor a can calculate any spanwise lattice The grid scale delta of the layers can interpolate the spanwise curve through the grid scale delta of any spanwise grid layer to obtain spanwise interpolation point coordinates, and further, a grid densely processed calculation grid can be generated based on the spanwise interpolation point coordinates, and as shown in fig. 14, the grid densely processed calculation grid 1401 can be obtained through the grid densely processing of the calculation grid 1301 in fig. 13, and the grid densely processed calculation grid 1401 is the calculation grid finally used for axial flow impeller mechanical analysis and research.

It should be noted that although the steps of the methods of the present disclosure are illustrated in the accompanying drawings in a particular order, this does not require or imply that the steps must be performed in that particular order or that all of the illustrated steps be performed in order to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step to perform, and/or one step decomposed into multiple steps to perform, etc.

In addition, in the present exemplary embodiment, a flow channel calculation grid generating apparatus is also provided. Referring to fig. 15, the flow path calculation grid generating apparatus 1500 includes: the impeller design parameter acquisition module 1510, the meridional structure map determination module 1520, the grid feature data determination module 1530, and the runner calculation grid generation module 1540. Wherein:

The impeller design parameter acquisition module 1510 is configured to acquire impeller design parameters of the axial flow impeller machine;

the meridian structure diagram determining module 1520 is configured to determine a meridian structure diagram of the axial flow type impeller machine on a single flow channel according to the impeller design parameter;

the grid feature data determination module 1530 is configured to determine calculated grid feature data based on the meridian structure map;

the runner computation grid generating module 1540 is configured to determine a computation grid corresponding to the single runner according to the computation grid feature data.

In one exemplary embodiment of the present disclosure, based on the foregoing, the impeller design parameters include a vane parameter and a non-vane parameter; the meridian-structure-map determination module 1520 is configured to:

determining a blade section of the axial flow impeller machine on a meridian plane of a single flow channel according to the blade parameters;

determining a non-blade section of the axial flow impeller machine on a meridian plane of a single flow channel according to the non-blade parameter;

determining a meridional structure of the axial flow turbomachine on the single flow channel through the blade sections and the non-blade sections.

In one exemplary embodiment of the present disclosure, calculating grid feature data includes meridian grid node coordinates based on the foregoing approach; the mesh feature data determination module 1530 is configured to:

Determining the total meridian node quantity according to the impeller design parameters of each section in the meridian structure diagram;

the total meridian node quantity is distributed to obtain the zone meridian node quantity corresponding to each zone;

constructing a meridian curve corresponding to the meridian structure diagram;

and carrying out equidistant interpolation on the meridian curve according to the number of the meridian nodes of the section to obtain the meridian grid node coordinates.

In one exemplary embodiment of the present disclosure, based on the foregoing scheme, the mesh feature data determination module 1530 is configured to:

calculating the average size corresponding to each section according to the impeller design parameters of each section in the meridian structure diagram;

and acquiring the preset number of the spanwise grid nodes, and determining the total number of meridian nodes according to the average size and the number of the spanwise grid nodes.

determining the estimated meridian node number corresponding to each section based on the average size and the total meridian node number;

and acquiring a preset multiple grid coefficient, and determining the number of zone meridian nodes corresponding to each zone according to the estimated number of meridian nodes and the multiple grid coefficient.

In one exemplary embodiment of the present disclosure, based on the foregoing scheme, the runner computation grid generation module 1540 is configured to:

determining polar angle coordinate values according to the meridian grid node coordinates;

determining grid node column coordinates through the meridian grid node coordinates and the polar angle coordinate values;

and generating a calculation grid corresponding to the single flow channel based on the grid node cylindrical coordinates.

In one exemplary embodiment of the present disclosure, based on the foregoing aspects, the polar angle coordinate values include blade segment polar angle coordinate values and non-blade segment polar angle coordinate values, the runner computation grid generation module 1540 is configured to:

performing coordinate system conversion on the meridian grid node coordinates to determine space grid node coordinates;

calculating polar angle coordinate values of the leaf segments according to the space grid node coordinates; and

and determining the polar angle coordinate value of the non-blade segment through the position relation of the non-blade segment and the blade segment in the meridian structural diagram and the polar angle coordinate value of the blade segment.

In one exemplary embodiment of the present disclosure, based on the foregoing, the polar angular coordinate values include a first polar angular coordinate value and a second polar angular coordinate value; the runner computation grid generation module 1540 is configured to:

Setting a grid node direction, and determining a first blade side grid node coordinate point set of the grid node column coordinates in the grid node direction;

determining the tangential thickness of the blade profile according to the first polar angle coordinate value and the second polar angle coordinate value;

converting the first blade side grid node coordinate point set based on the tangential thickness of the blade profile to obtain a second blade side grid node coordinate point set;

and connecting coordinate points between the first blade side grid node coordinate point set and the second blade side grid node coordinate point set to generate a calculation grid corresponding to the single flow channel.

In an exemplary embodiment of the present disclosure, based on the foregoing scheme, the runner computing grid generating device 1500 further includes a computing grid encryption module that may be used to:

performing grid densification treatment on the computing grid;

the grid densification processing for the computing grid comprises the following steps:

acquiring a preset growth factor and the number of the spanwise grid nodes;

determining the grid scale of each spanwise grid layer according to the length of the spanwise curve of the calculated grid, the growth factors and the number of nodes of the spanwise grid;

Interpolating the spanwise curve through the grid scale to obtain spanwise interpolation point coordinates;

and generating a calculation grid after grid densification processing based on the spanwise interpolation point coordinates.

The specific details of each module of the above-mentioned medium-channel computing grid generating device are already described in detail in the corresponding channel computing grid generating method, so that they will not be described in detail here.

It should be noted that although several modules or units of the runner computation grid generating device are mentioned in the above detailed description, this division is not mandatory. Indeed, the features and functionality of two or more modules or units described above may be embodied in one module or unit in accordance with embodiments of the present disclosure. Conversely, the features and functions of one module or unit described above may be further divided into a plurality of modules or units to be embodied.

In addition, in an exemplary embodiment of the present disclosure, an electronic device capable of implementing the above-described runner computation grid generation method is also provided.

Those skilled in the art will appreciate that the various aspects of the present disclosure may be implemented as a system, method, or program product. Accordingly, various aspects of the disclosure may be embodied in the following forms, namely: an entirely hardware embodiment, an entirely software embodiment (including firmware, micro-code, etc.) or an embodiment combining hardware and software aspects may be referred to herein as a "circuit," module "or" system.

An electronic device 1600 according to such an embodiment of the present disclosure is described below with reference to fig. 16. The electronic device 1600 shown in fig. 16 is merely an example and should not be construed as limiting the functionality and scope of use of embodiments of the present disclosure.

As shown in fig. 16, the electronic device 1600 is embodied in the form of a general purpose computing device. The components of the electronic device 1600 may include, but are not limited to: the at least one processing unit 1610, the at least one memory unit 1620, a bus 1630 connecting the different system components (including the memory unit 1620 and the processing unit 1610), and a display unit 1640.

Wherein the storage unit stores program code that is executable by the processing unit 1610 such that the processing unit 1610 performs steps according to various exemplary embodiments of the present disclosure described in the above-described "exemplary methods" section of the present specification. For example, the processing unit 1610 may perform step S210 shown in fig. 2 to obtain impeller design parameters of an axial flow type impeller machine; step S220, determining a meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters; step S230, determining calculation grid characteristic data based on the meridian structure diagram; step S240, determining a computing grid corresponding to the single flow channel according to the computing grid feature data.

The memory unit 1620 may include readable media in the form of volatile memory units, such as Random Access Memory (RAM) 1621 and/or cache memory 1622, and may further include Read Only Memory (ROM) 1623.

The storage unit 1620 may also include a program/utility 1624 having a set (at least one) of program modules 1625, such program modules 1625 including, but not limited to: an operating system, one or more application programs, other program modules, and program data, each or some combination of which may include an implementation of a network environment.

Bus 1630 may be a local bus representing one or more of several types of bus structures including a memory unit bus or memory unit controller, a peripheral bus, an accelerated graphics port, a processing unit, or using any of a variety of bus architectures.

Electronic device 1600 may also communicate with one or more external devices 1670 (e.g., keyboard, pointing device, bluetooth device, etc.), with one or more devices that enable a user to interact with electronic device 1600, and/or with any device (e.g., router, modem, etc.) that enables electronic device 1600 to communicate with one or more other computing devices. Such communication may occur through an input/output (I/O) interface 1650. Also, electronic device 1600 can communicate with one or more networks such as a Local Area Network (LAN), a Wide Area Network (WAN) and/or a public network, such as the Internet, through network adapter 1660. As shown, network adapter 1660 communicates with other modules of electronic device 1600 over bus 1630. It should be appreciated that although not shown, other hardware and/or software modules may be used in connection with electronic device 1600, including, but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, data backup storage systems, and the like.

From the above description of embodiments, those skilled in the art will readily appreciate that the example embodiments described herein may be implemented in software, or in combination with the necessary hardware. Thus, the technical solution according to the embodiments of the present disclosure may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (may be a CD-ROM, a U-disk, a mobile hard disk, etc.) or on a network, and includes several instructions to cause a computing device (may be a personal computer, a server, a terminal device, or a network device, etc.) to perform the method according to the embodiments of the present disclosure.

In an exemplary embodiment of the present disclosure, a computer-readable storage medium having stored thereon a program product capable of implementing the method described above in the present specification is also provided. In some possible embodiments, the various aspects of the present disclosure may also be implemented in the form of a program product comprising program code for causing a terminal device to carry out the steps according to the various exemplary embodiments of the disclosure as described in the "exemplary methods" section of this specification, when the program product is run on the terminal device.

Referring to fig. 17, a program product 1700 for implementing the above-described runner computation grid generation method according to an embodiment of the present disclosure is described, which may employ a portable compact disc read-only memory (CD-ROM) and include program code, and may be run on a terminal device, such as a personal computer. However, the program product of the present disclosure is not limited thereto, and in this document, a readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.

The program product may employ any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. The readable storage medium can be, for example, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples (a non-exhaustive list) of the readable storage medium would include the following: an electrical connection having one or more wires, a portable disk, a hard disk, random Access Memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

The computer readable signal medium may include a data signal propagated in baseband or as part of a carrier wave with readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A readable signal medium may also be any readable medium that is not a readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.

Program code embodied on a readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

Program code for carrying out operations of the present disclosure may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computing device, partly on the user's device, as a stand-alone software package, partly on the user's computing device, partly on a remote computing device, or entirely on the remote computing device or server. In the case of remote computing devices, the remote computing device may be connected to the user computing device through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computing device (e.g., connected via the Internet using an Internet service provider).

Furthermore, the above-described figures are only schematic illustrations of processes included in the method according to the exemplary embodiments of the present disclosure, and are not intended to be limiting. It will be readily appreciated that the processes shown in the above figures do not indicate or limit the temporal order of these processes. In addition, it is also readily understood that these processes may be performed synchronously or asynchronously, for example, among a plurality of modules.

From the above description of embodiments, those skilled in the art will readily appreciate that the example embodiments described herein may be implemented in software, or in combination with the necessary hardware. Thus, the technical solution according to the embodiments of the present disclosure may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (may be a CD-ROM, a U-disk, a mobile hard disk, etc.) or on a network, and includes several instructions to cause a computing device (may be a personal computer, a server, a touch terminal, or a network device, etc.) to perform the method according to the embodiments of the present disclosure.

Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This application is intended to cover any adaptations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.

It is to be understood that the present disclosure is not limited to the precise arrangements and instrumentalities shown in the drawings, and that various modifications and changes may be effected without departing from the scope thereof. The scope of the present disclosure is limited only by the appended claims.

Claims

1. A runner computing grid generation method, comprising:

obtaining impeller design parameters of an axial flow impeller machine;

2. The flow channel computing grid generation method of claim 1, wherein the impeller design parameters include blade parameters and non-blade parameters; the determining the meridian structure diagram of the axial flow impeller machine on a single flow channel according to the impeller design parameters comprises the following steps:

3. The runner computation grid generation method of claim 2, wherein the computation grid feature data includes meridian grid node coordinates; the determining the calculated grid characteristic data based on the meridian structure diagram comprises the following steps:

constructing a meridian curve corresponding to the meridian structure diagram;

4. The method of generating a grid for runner computation according to claim 3, wherein said determining the total meridian node number based on the impeller design parameters of each segment in said meridian structure map comprises

5. The method for generating the runner computation grid according to claim 4, wherein the allocating the total meridian node number to obtain the segment meridian node number corresponding to each segment includes:

6. The flow channel computing grid generation method according to claim 3, wherein the determining the computing grid corresponding to the single flow channel by the computing grid feature data includes:

7. The method of claim 6, wherein the polar values include blade segment polar values and non-blade segment polar values, wherein determining polar values from the meridional grid node coordinates comprises:

8. The flow channel calculation grid generation method according to claim 6 or 7, wherein the polar angle coordinate values include a first polar angle coordinate value and a second polar angle coordinate value; the generating the calculation grid corresponding to the single flow channel based on the grid node cylindrical coordinates comprises the following steps:

9. The flow channel computing grid generation method according to claim 1, characterized in that the method further comprises:

performing grid densification treatment on the computing grid;

acquiring a preset growth factor and the number of the spanwise grid nodes;

10. A runner computation mesh generation apparatus, comprising: