CN116484771B - Method and device for generating CFD grid of axial flow compressor - Google Patents
Method and device for generating CFD grid of axial flow compressor Download PDFInfo
- Publication number
- CN116484771B CN116484771B CN202310741561.6A CN202310741561A CN116484771B CN 116484771 B CN116484771 B CN 116484771B CN 202310741561 A CN202310741561 A CN 202310741561A CN 116484771 B CN116484771 B CN 116484771B
- Authority
- CN
- China
- Prior art keywords
- grid
- flow channel
- parameter
- circumferential
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 239000012530 fluid Substances 0.000 claims abstract description 105
- 239000003102 growth factor Substances 0.000 claims description 7
- 230000001172 regenerating effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 18
- 238000004364 calculation method Methods 0.000 description 5
- 239000007787 solid Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 2
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 description 1
- KLDZYURQCUYZBL-UHFFFAOYSA-N 2-[3-[(2-hydroxyphenyl)methylideneamino]propyliminomethyl]phenol Chemical compound OC1=CC=CC=C1C=NCCCN=CC1=CC=CC=C1O KLDZYURQCUYZBL-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- FEPMHVLSLDOMQC-UHFFFAOYSA-N virginiamycin-S1 Natural products CC1OC(=O)C(C=2C=CC=CC=2)NC(=O)C2CC(=O)CCN2C(=O)C(CC=2C=CC=CC=2)N(C)C(=O)C2CCCN2C(=O)C(CC)NC(=O)C1NC(=O)C1=NC=CC=C1O FEPMHVLSLDOMQC-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
The embodiment of the application discloses a method and a device for generating CFD grids of an axial-flow compressor; the method comprises the following steps: determining the fluid domain of a single flow channel in an axial flow compressor according to a geometric model of the axial flow compressor; selecting a grid dividing method corresponding to the fluid domain of the single flow channel; setting corresponding direction expanding parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division; and generating a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the flow direction parameter and the circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
Description
Technical Field
The embodiment of the application relates to the technical field of computational fluid dynamics (Computational Fluid Dynamics, CFD) grid generation, in particular to a method and a device for generating a CFD grid of an axial flow compressor.
Background
In fluid simulation calculation of an axial flow compressor, space discretization is needed to be realized based on grids, and the quality of the grids can directly influence the accuracy of overall calculation. Also, if the time spent by the iterative computation itself is not taken into account, typically more than half of the time is used to adjust and modify the grid. When the grid is adjusted, the grid generation parameters are required to be adjusted repeatedly, the grid is further regenerated, and then the quality of the grid is observed; if the adjusted grid quality does not meet the requirement, the grid generation parameters are needed to be adjusted again, and the above processes are repeated in sequence until the grid quality meets the requirement. Conceivably, depending on the scale of the number of grids, it takes generally tens of seconds to several minutes to generate one grid, and a lengthy grid generation process certainly slows down the overall work efficiency. Therefore, a method for accelerating the grid generation is needed.
Disclosure of Invention
In view of this, the embodiments of the present application are expected to provide a method and apparatus for generating CFD meshes of an axial flow compressor; the method can improve the CFD grid generation of the axial-flow compressor, and improve the CFD grid generation efficiency and the CFD grid quality.
The technical scheme of the embodiment of the application is realized as follows:
in a first aspect, an embodiment of the present application provides a method for generating a CFD grid of an axial flow compressor, where the method includes:
determining the fluid domain of a single flow channel in an axial flow compressor according to a geometric model of the axial flow compressor;
selecting a grid dividing method corresponding to the fluid domain of the single flow channel;
setting corresponding direction expanding parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division;
and generating a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the flow direction parameter and the circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
In a second aspect, an embodiment of the present application provides an apparatus for generating a CFD grid of an axial flow compressor, where the apparatus includes a determining portion, a selecting portion, a setting portion, and a grid dividing portion; wherein,,
the determining part is configured to determine the fluid domain of a single flow channel in the axial flow compressor according to a geometric model of the axial flow compressor;
the selecting part is configured to select a grid dividing method corresponding to the fluid domain of the single flow channel;
the setting part is configured to set corresponding spanwise parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division;
the grid dividing part is configured to generate a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the set flow direction parameter and the set circumferential parameter, and generate a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
The embodiment of the application provides a method and a device for generating CFD grids of an axial-flow compressor; determining the fluid domain of a single flow channel of the axial flow compressor according to the geometric model of the axial flow compressor, and selecting a corresponding grid dividing method; generating a preview grid corresponding to the fluid domain of the single flow channel based on the set spanwise parameter, the flow direction parameter and the circumferential parameter, so that a final 3D grid can be generated when the preview grid corresponding to the fluid domain of the single flow channel meets the set requirement; it can be understood that the 3D grid corresponding to the fluid domain of the whole flow channel of the axial flow compressor can be obtained by performing the periodical symmetric operation on the 3D grid of the fluid domain of the single flow channel. The CFD grid dividing method provided by the embodiment of the application can generate the preview grid through the set spanwise parameter, the set flow direction parameter and the set circumferential parameter; it can be appreciated that when the grid quality of the generated preview grid does not meet the set requirement, the grid quality can be adjusted by changing the grid division parameters until the preview grid meets the set requirement and then the final 3D grid is regenerated. In the CFD grid dividing method provided by the embodiment of the application, only the preview grid is needed to be generated when the grid is adjusted, and the final 3D grid is regenerated when the adjustment is basically finished, and the adjustment of the preview grid belongs to local grid adjustment, so that the number of grids needed to be generated is greatly reduced, the preview grid can be generated rapidly, and the generation rate of the whole grid is accelerated.
Drawings
Fig. 1 is a schematic flow chart of a method for generating CFD meshes of an axial flow compressor according to an embodiment of the present application;
fig. 2 is a schematic diagram of a geometric model of an axial compressor according to an embodiment of the present application;
FIG. 3 is a meridional view of the axial compressor shown in FIG. 2;
FIG. 4 is a schematic diagram of a full flow channel structure according to an embodiment of the present application;
FIG. 5 is a schematic diagram of an H-topology method according to an embodiment of the present application;
FIG. 6 is a schematic diagram of an OH topology method according to an embodiment of the present application;
FIG. 7 is a schematic diagram of a method for determining the shape of an OH-shaped topological single runner according to an embodiment of the present application;
FIG. 8 is a schematic view of the location of a first control point at the leading edge inlet provided by an embodiment of the present application;
FIG. 9 is a schematic illustration of the location of a second control point at the leading edge inlet provided by an embodiment of the present application;
FIG. 10 is a schematic illustration of the location of a third control point at the leading edge inlet provided by an embodiment of the present application;
FIG. 11 is a schematic illustration of a flow field of a single flow channel in an OH topology provided by an embodiment of the present application;
FIG. 12 is a schematic diagram of a multiple grid as a fine grid according to an embodiment of the present application;
FIG. 13 is a schematic diagram of a multiple grid as a medium grid according to an embodiment of the present application;
FIG. 14 is a schematic diagram of a multiple grid as a fine grid according to an embodiment of the present application;
FIG. 15 is a schematic view of an exponentially distributed spanwise grid provided by an embodiment of the present application;
FIG. 16 is a lower half of the spanwise grid of FIG. 15;
FIG. 17 is a schematic diagram of an OH sub-block according to an embodiment of the present application;
FIG. 18 is a schematic view of an O-block with a relative thickness of 0.15 according to an embodiment of the present application;
FIG. 19 is a schematic view of an O-block with a relative thickness of 0.5 according to an embodiment of the present application;
FIG. 20 is a schematic diagram of an O-block grid when the distribution coefficient of the O-block circumferential direction is 0.2 according to an embodiment of the present application;
FIG. 21 is a schematic diagram of an O-block grid when the distribution coefficient of the O-block circumferential direction is 0.8 according to an embodiment of the present application;
FIG. 22 is a schematic diagram of a preview grid corresponding to a fluid domain of a single flow channel according to an embodiment of the present application;
FIG. 23 is a schematic view of a 3D grid corresponding to a fluid domain of a single flow channel according to an embodiment of the present application;
fig. 24 is a schematic diagram of a device for generating CFD meshes of an axial flow compressor according to an embodiment of the present application;
fig. 25 is a schematic diagram of a specific hardware structure of a device for generating CFD meshes of an axial flow compressor according to an embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application.
Referring to fig. 1, a method for generating CFD grids of an axial flow compressor according to an embodiment of the present application is shown, where the method includes:
s101, determining a fluid domain of a single flow channel in an axial flow compressor according to a geometric model of the axial flow compressor;
s102, selecting a grid division method corresponding to the fluid domain of the single flow channel;
s103, setting corresponding spanwise parameters, flow direction parameters and circumferential parameters of the fluid domains of the single flow channels when grid division is carried out;
s104, generating a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the set flow direction parameter and the set circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
For the technical scheme shown in fig. 1, determining the fluid domain of a single flow channel of an axial flow compressor according to a geometric model of the axial flow compressor, and selecting a corresponding grid division method; generating a preview grid corresponding to the fluid domain of the single flow channel based on the set spanwise parameter, the flow direction parameter and the circumferential parameter, so that a final 3D grid can be generated when the preview grid corresponding to the fluid domain of the single flow channel meets the set requirement; it can be understood that the 3D grid corresponding to the fluid domain of the whole flow channel of the axial flow compressor can be obtained by performing the periodical symmetric operation on the 3D grid of the fluid domain of the single flow channel. The CFD grid dividing method provided by the embodiment of the application can generate the preview grid through the set spanwise parameter, the set flow direction parameter and the set circumferential parameter; it can be appreciated that when the grid quality of the generated preview grid does not meet the set requirement, the grid quality can be adjusted by changing the grid division parameters until the preview grid meets the set requirement and then the final 3D grid is regenerated. In the CFD grid dividing method provided by the embodiment of the application, only the preview grid is needed to be generated when the grid is adjusted, and the final 3D grid is regenerated when the adjustment is basically finished, and the adjustment of the preview grid belongs to local grid adjustment, so that the number of grids needed to be generated is greatly reduced, the preview grid can be generated rapidly, and the generation rate of the whole grid is accelerated.
For the solution shown in fig. 1, in some possible embodiments, the determining, according to a geometric model of an axial compressor, a fluid domain of a single flow channel in the axial compressor includes:
according to the geometric model of the axial flow compressor, determining an annular area surrounded by an inlet surface, an outlet surface, a hub surface, a casing surface and a blade surface of the axial flow compressor as a full flow channel;
obtaining the single flow channel based on the full flow channel division;
determining the shape of the single flow channel;
and obtaining the fluid domain corresponding to the single flow channel based on the shape topology of the single flow channel.
Specifically, referring to fig. 2, a geometric model of the axial flow compressor 2 is shown, and as shown in fig. 2, the geometric model of the axial flow compressor 2 includes a row of moving blades 21, a row of stationary blades 22, a hub 23 and a casing 24; it will be appreciated that the axial compressor 2 is generally composed of two major parts, the disk and blades connected to the rotary shaft of the axial compressor constitute the rotor of the axial compressor, and the casing, which is not rotated externally, and the blades connected to the casing constitute the stator of the axial compressor. The blades on the rotor are called blades 21 and the blades on the stator are called vanes 22. Each row of blades 21 and the immediately following row of vanes 22 constitute a stage of an axial compressor. It will be appreciated, therefore, that where the axial compressor 2 is comprised of multiple stages, fig. 2 shows a schematic diagram of only one stage of the axial compressor.
Note that, a solid arrow A1 in fig. 2 indicates a flow direction of the axial compressor 2; solid arrow B1 represents the circumferential direction of the axial compressor 2; the solid arrow C1 indicates the spanwise direction of the axial compressor 2.
As shown in fig. 3, which shows a meridian view of the axial compressor 2 shown in fig. 2, wherein the abscissa indicates the rotation axis direction Z in meters (m) and the ordinate indicates the spanwise radius R in meters (m). In fig. 3, a solid black line D1 indicates a blade row interface; 31 denotes a casing surface, 32 denotes an inlet surface, 33 denotes a hub surface, and 34 denotes an outlet surface. It will be appreciated that the inlet face 32 is generally defined by the starting position of the hub face 33 and the starting position of the casing face 31; the outlet face 34 is generally defined by the end position of the hub face 33 and the end position of the casing face 31; the blade row interface D1 divides the flow direction of the axial compressor into a plurality of "blade rows", so that each blade row can respectively consider the generated grid when specifically generating the grid, and the blade row interface D1 is generally located at the middle position of the two rows of blades.
It is generally desirable in practice to specify the inlet face 32, the outlet face 34, and the blade row interface D1. Of course, the above-mentioned process can be automatically completed by program, also can be manually set by user.
Based on the above, as shown in fig. 4, the region surrounded by the casing surface 31, the inlet surface 32, the hub surface 33, the outlet surface 34, and the blade row interface D1 can be referred to as a full flow path.
In addition, it should be noted that in the implementation process, it is generally assumed that the fluid area of the axial flow compressor has periodicity along the circumferential direction, so only the fluid area of a single flow channel needs to be calculated, and then the fluid area of the full flow channel can be obtained through periodic symmetry. Referring to fig. 5 and 6, two modes of full runner dividing single runner H-topology and OH-topology are shown, respectively.
For the above embodiments, in some examples, the determining the shape of the single flow channel includes:
the shape of the flow channel from the inlet throat to the outlet throat is determined by the shape of the blade;
the shape of the flow channel at the leading edge inlet and the trailing edge outlet is determined according to the set control point.
Optionally, the control point set at the leading edge inlet includes:
the first control point is positioned on an extension line of the throat flow direction, and the position coordinate of the first control point is determined by using a first parameter A in the following formula:
wherein, the value range of A is floating point number of [0,1 ];
the second control point is positioned at the left side of the front edge of the blade, and the position coordinates of the second control point are respectively determined along the two directions of the flow direction and the circumferential direction by using a second parameter B and a third parameter C in the following formula:
the third control point is located at the left side of the second control point and is used for controlling the tangential direction at the second control point, and the position coordinate of the third control point is determined by using a fourth parameter D in the following formula:
。
in the embodiment of the application, taking OH topology as an example, the shape of the flow channel from the inlet throat to the outlet throat is determined by the shape of the blade, and the shape of the flow channel of the non-blade segment is controlled by three control points, that is, the shape of the flow channel at the inlet of the front edge and the outlet of the tail edge can be controlled according to the three set control points respectively. Referring to FIG. 7, which illustrates an OH topology single-runner shape determination method, in an embodiment of the present application, the shape of a single runner is determined in an RTheta-M coordinate system, which, as understood, corresponds to expanding a cylindrical surface of a constant R coordinate of a cylindrical coordinate system (Z, R, theta) into a plane; the abscissa M corresponds to the Z axis of the cylindrical coordinates and is also the rotation axis direction. Wherein L1 in fig. 7 represents a flow path dividing line; l2 represents a flow passage dividing line with a control point; p represents the middle position of the flow passage dividing line; p (P) in Indicating the position of the inlet throat, P out Indicating the exit throat position.
It should be noted that, the above calculation formula only shows that the shape of the flow channel is controlled by 4 parameters at the front edge inlet; similarly, in the implementation process, the trailing edge outlet may also be determined by a set control point, and the calculation method is similar to that of the control point set at the leading edge inlet, so that the embodiment of the present application will not be described again.
Referring to fig. 8, the location of the first control point at the leading edge inlet is shown, where the first control point is shown as a black square (indicated by the dashed arrow in fig. 8).
Referring to fig. 9, the location of the second control point at the leading edge inlet is shown, where the second control point is shown as a black square (indicated by the dashed arrow in fig. 9).
Referring to fig. 10, the location of a third control point at the leading edge inlet is shown, where the third control point is shown as the leftmost black square in fig. 10 (indicated by the dashed arrow). Where α represents the flow channel inlet angle and β represents the blade inlet angle.
It should be noted that the above three control points are determined by four relative parameters based on the ratio, rather than directly specifying the absolute coordinates of the control points.
Finally, the flow field of a single flow channel of OH topology is shown in fig. 11.
For the technical solution described in fig. 1, in some possible embodiments, the selecting a mesh division method corresponding to the fluid domain of the single flow channel includes:
and a method of multiple meshing is adopted for the fluid domain of the single flow channel.
It will be appreciated that the multiple grid is a CFD acceleration convergence technique, and that there is no direct relationship between multiple grids and grid generation. Just to use multiple grids, the number of grids would be required to satisfy: the number of grids/(2 multiple grid levels) is a positive integer. Grid points = grid number +1.
Fig. 12, 13 and 14 show schematic diagrams of the multiple grids as a fine grid, a medium grid and a coarse grid, respectively. As can be seen from fig. 12, 13, 14, the number of grids from fine grids to coarse grids is 1 by 2 grids per coarsening in each direction, so the number of grids must be a positive integer multiple of (2 multiple grid levels).
In addition, the calculation flow of the multiple grids is as follows: firstly, calculating from a fine grid, then interpolating residual errors to a coarse grid through a limiting operator, after calculating the coarse grid, reconstructing back to the fine grid through interpolation of a continuation operator, thus finishing one round of correction of the fine grid, and then repeatedly cycling until convergence. Thus, a "generated grid" generates a fine grid of multiple grids.
For the technical solution described in fig. 1, in some possible embodiments, the setting of the spanwise parameter, the flow direction parameter, and the circumferential parameter corresponding to the fluid domain of the single flow channel when performing grid division includes:
setting corresponding direction-spreading network parameters when the fluid domain of the single flow channel is subjected to grid division;
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades;
and setting the flow direction parameter and the circumferential parameter.
For the foregoing embodiments, in some examples, the setting the corresponding spanwise parameter when the fluid domain of the single flow channel is gridded includes:
setting the corresponding spanwise grid width of the fluid domain of the single flow channel to be in an equal-ratio array when grid division is carried out, and calculating according to the following formula to obtain the spanwise parameters：
Wherein S represents a growth factor;representing a first layer mesh width;
wherein, when i takes on the value of N,
wherein L represents half of the height of the flow channel; n represents half the number of meshes; and is also provided with。
Referring to fig. 15, an exponentially distributed spanwise grid schematic is shown. The spanwise mesh was dense near the hub and the casing, and sparse near the average diameter. Fig. 16 shows the lower half of the spanwise grid of fig. 15, with the grid width increasing monotonically as seen in fig. 16.
It should be noted that, the exponential distribution of the grid is defined as: the ratio of the widths of two adjacent grids is equal to a constant, that is, the width of the spanwise grid is in an equal-ratio array.
Therefore, in the implementation process, half of the height of the flow channel is L, the growth factor S, half of the grid number is N, and the first layer of grid width. According to the sum formula of the equal ratio array:
when the value of i is N,
thus, L, S, N,And a relational expression between the two. Namely S and->Is a conversion formula of (a). And can calculate the spanwise grid distribution +.>，/>。
Furthermore, the gap grid is not considered in setting the spanwise parameters in the present embodiment.
For the above embodiments, in some examples, initializing the flow direction parameter and the circumferential parameter according to the spanwise parameter and the blade geometric model feature includes:
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades, wherein the initialization principle is as follows: making the average value of the spanwise mesh length equal to the average value of the flow direction and the circumferential mesh length;
wherein the flow direction parameters include: the number of H1 block flow direction nodes, the number of H4 block flow direction nodes, the number of H2 block circumferential nodes, the number of H3 block circumferential nodes, the number of O block circumferential nodes, the O block circumferential distribution coefficient, the O block circumferential starting node, the O block relative thickness, the number of O block normal nodes, the O block normal growth factor, the inlet block grid-pitch ratio, the inlet block node number, the outlet block grid-pitch ratio and the outlet block node number.
Referring to fig. 17, a schematic diagram of OH sub-blocks is shown, including an inlet block RU, an O block, an H1 block, an H2 block, an H3 block, an H4 block, and an outlet block CH. Thus, in the implementation process, the flow direction parameters can be obtained by setting the parameters of the inlet block RU, the O block, the H1 block, the H2 block, the H3 block, the H4 block and the outlet block CH.
On the other hand, the circumferential grid parameters are relatively independent, irrespective of the section height position. Therefore, after initialization, the circumferential grid parameters are irrelevant to the spanwise parameters and the geometric model characteristics of the blades, and can be freely set by a user.
It will be appreciated that after setting the spanwise, circumferential and flow direction parameters, a preview grid can be generated and visually displayed.
Therefore, for the technical solution described in fig. 1, in some possible embodiments, the generating, according to the set spanwise parameter, the flow direction parameter, and the circumferential parameter, the preview grid corresponding to the fluid domain of the single flow channel by using the selected grid dividing method, and generating, when the preview grid meets the set requirement, the 3D grid corresponding to the fluid domain of the single flow channel includes:
determining leaf height parameters;
acquiring the flow channel shape and the blade shape corresponding to the blade height;
inputting the relative thickness of the O blocks, and dividing the O blocks of the OH topology by using the relative thickness of the O blocks;
dividing the O block grid according to the number of O block ring nodes and the O block ring distribution coefficient;
setting other circumferential parameters, and drawing to obtain a preview grid corresponding to the fluid domain of the single flow channel;
when the preview grid meets the set requirement, regenerating a 3D grid corresponding to the fluid domain of the single flow channel by using the spanwise parameter, the circumferential parameter and the blade geometric model characteristic;
and solving a poisson equation to smooth the 3D grid corresponding to the fluid domain of the single flow channel.
It will be appreciated that in the implementation process, the following steps are required to generate the preview grid corresponding to the fluid domain of the single flow channel:
1) Inputting a set leaf height parameter, wherein the leaf height is in a value range of 0% to 100%. In particular embodiments, the blade root is typically 0% and the blade tip is typically 100%.
2) Obtaining a runner shape and a blade shape corresponding to the blade height;
3) The relative thickness of the O-block is entered, wherein the relative thickness of the O-block defaults to 0.15. Wherein fig. 18 and 19 show schematic views of the O-block with a relative thickness of 0.15 and 0.5, respectively. Meanwhile, dividing O blocks of an OH topology by using the relative thickness of the O blocks;
4) Dividing an O block grid according to the number of O block ring nodes and the O block ring distribution coefficient; wherein, fig. 20 and 21 show O-block grid diagrams when the O-block circumferential distribution coefficients are 0.2 and 0.8, respectively.
5) Setting other circumferential parameters, and drawing a preview grid corresponding to the fluid domain of the single flow channel, wherein the obtained preview grid corresponding to the fluid domain of the single flow channel is shown in fig. 22.
The preview grid corresponding to the fluid domain of the single flow channel generated above is not related to the 3D grid corresponding to the fluid domain of the single flow channel, and all the cross-sectional layers are recalculated when the 3D grid corresponding to the fluid domain of the single flow channel is generated.
In addition, in the preview grid interface, the leaf height positions of the sections may be switched, and only one section grid of the leaf height positions may be generated at a time. The preview grid generation speed is high, so that the grid generation waiting time in the grid adjustment stage can be reduced.
When the fluid engineer considers that the mesh quality of the preview mesh corresponding to the fluid domain of the single flow channel shown in fig. 22 meets the set requirement, the 3D mesh corresponding to the fluid domain of the single flow channel can be regenerated by using the spanwise parameter, the circumferential parameter and the blade geometric model feature, and the poisson equation smoothing process can be solved. Wherein,,
poisson's equation is to coordinate the physical domainTo calculate domain coordinates +.>The transformed elliptic equation is written as:
where P, Q, R is the source item used to control the grid distribution. UsingAs an argument, the equation is transformed and then solved.
The resulting 3D grid corresponding to the fluid domain of the single flow channel is shown in fig. 23. It can be understood that the 3D grid corresponding to the fluid domain of the single flow channel shown in fig. 23 is periodically symmetric, so that the 3D grid corresponding to the fluid domain of the full flow channel can be obtained.
In the process of generating the preview grid and the 3D grid corresponding to the fluid domain of the single flow channel, the complete grid generation flow is split into a plurality of sub-steps, each sub-step is visually displayed, and a fluid engineer can immediately see the modified influence range by only executing part of sub-steps after adjusting the grid generation parameters, so that the whole grid dividing process is direct, and the grid generation efficiency is improved. Secondly, in the embodiment of the application, the relative ratio parameter of the grid generation parameter is used, so that the preview grid quality can represent the 3D grid quality, only the local preview grid quality is needed to be checked in the grid division process, the grid division time is saved, and the global 3D grid is not needed to be generated when the grid is adjusted.
Based on the same inventive concept as the previous technical solution, referring to fig. 24, there is shown a composition of an apparatus 240 for generating CFD meshes of an axial flow compressor according to an embodiment of the present application, where the apparatus 240 includes a determining portion 2401, a selecting portion 2402, a setting portion 2403 and a mesh dividing portion 2404; wherein,,
the determining part 2401 is configured to determine a fluid domain of a single flow channel in an axial flow compressor according to a geometric model of the axial flow compressor;
the selecting portion 2402 is configured to select a grid division method corresponding to the fluid domain of the single flow channel;
the setting portion 2403 is configured to set a spanwise parameter, a flow direction parameter, and a circumferential parameter corresponding to the fluid domain of the single flow channel when meshing;
the meshing part 2404 is configured to generate a preview mesh corresponding to the fluid domain of the single flow channel by adopting the selected meshing method according to the set spanwise parameter, the set flow direction parameter and the set circumferential parameter, and generate a 3D mesh corresponding to the fluid domain of the single flow channel when the preview mesh meets the set requirement.
Optionally, the determining portion 2401 is configured to:
according to the geometric model of the axial flow compressor, determining an annular area surrounded by an inlet surface, an outlet surface, a hub surface, a casing surface and a blade surface of the axial flow compressor as a full flow channel;
obtaining the single flow channel based on the full flow channel division;
determining the shape of the single flow channel;
and obtaining the fluid domain corresponding to the single flow channel based on the shape topology of the single flow channel.
Optionally, the determining portion 2401 is configured to:
the shape of the flow channel from the inlet throat to the outlet throat is determined by the shape of the blade;
the shape of the flow channel at the leading edge inlet and the trailing edge outlet is determined according to the set control point.
Optionally, the determining portion 2401 is configured to:
the first control point is positioned on an extension line of the throat flow direction, and the position coordinate of the first control point is determined by using a first parameter A in the following formula:
wherein, the value range of A is floating point number of [0,1 ];
the second control point is positioned at the left side of the front edge of the blade, and the position coordinates of the second control point are respectively determined along the two directions of the flow direction and the circumferential direction by using a second parameter B and a third parameter C in the following formula:
the third control point is located at the left side of the second control point and is used for controlling the tangential direction at the second control point, and the position coordinate of the third control point is determined by using a fourth parameter D in the following formula:
。
optionally, the selecting portion 2402 is configured to:
and a method of multiple meshing is adopted for the fluid domain of the single flow channel.
Optionally, the setting portion 2403 is configured to:
setting corresponding direction-spreading network parameters when the fluid domain of the single flow channel is subjected to grid division;
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades;
and setting the flow direction parameter and the circumferential parameter.
Optionally, the setting portion 2403 is configured to:
setting the corresponding spanwise grid width of the fluid domain of the single flow channel to be in an equal-ratio array when grid division is carried out, and calculating according to the following formula to obtain the spanwise parameters：
Wherein S represents a growth factor;representing a first layer mesh width;
wherein, when i takes on the value of N,
wherein L represents half of the height of the flow channel; n represents half the number of meshes; and is also provided with。
Optionally, the setting portion 2403 is configured to:
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades, wherein the initialization principle is as follows: making the average value of the spanwise mesh length equal to the average value of the flow direction and the circumferential mesh length;
wherein the flow direction parameters include: the number of H1 block flow direction nodes, the number of H4 block flow direction nodes, the number of H2 block circumferential nodes, the number of H3 block circumferential nodes, the number of O block circumferential nodes, the O block circumferential distribution coefficient, the O block circumferential starting node, the O block relative thickness, the number of O block normal nodes, the O block normal growth factor, the inlet block grid-pitch ratio, the inlet block node number, the outlet block grid-pitch ratio and the outlet block node number.
Optionally, the meshing portion 2404 is configured to:
determining leaf height parameters;
acquiring the flow channel shape and the blade shape corresponding to the blade height;
inputting the relative thickness of the O blocks, and dividing the O blocks of the OH topology by using the relative thickness of the O blocks;
dividing the O block grid according to the number of O block ring nodes and the O block ring distribution coefficient;
setting other circumferential parameters, and drawing to obtain a preview grid corresponding to the fluid domain of the single flow channel;
when the preview grid meets the set requirement, regenerating a 3D grid corresponding to the fluid domain of the single flow channel by using the spanwise parameter, the circumferential parameter and the blade geometric model characteristic;
and solving a poisson equation to smooth the 3D grid corresponding to the fluid domain of the single flow channel.
It will be appreciated that in this embodiment, a "part" may be a part of a circuit, a part of a processor, a part of a program or software, etc., and of course may be a unit, or a module may be non-modular.
In addition, each component in the present embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional modules.
The integrated units, if implemented in the form of software functional modules, may be stored in a computer-readable storage medium, if not sold or used as separate products, and based on such understanding, the technical solution of the present embodiment may be embodied essentially or partly in the form of a software product, which is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) or processor to perform all or part of the steps of the method described in the present embodiment. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read Only Memory (ROM), a random access Memory (RAM, randomAccess Memory), a magnetic disk, an optical disk, or other various media capable of storing program codes.
The present embodiment thus provides a computer storage medium storing a program for generating an axial compressor CFD grid, where the program for generating an axial compressor CFD grid, when executed by at least one processor, implements the method steps for generating an axial compressor CFD grid in the above technical solution.
Referring to fig. 25, a specific hardware structure of a computing device 250 capable of implementing the apparatus 240 for generating an axial compressor CFD grid according to the foregoing apparatus 240 for generating an axial compressor CFD grid is shown, where the computing device 250 may be a wireless device, a mobile or cellular phone (including a so-called smart phone), a Personal Digital Assistant (PDA), a video game console (including a video display, a mobile video game device, a mobile video conference unit), a laptop computer, a desktop computer, a television set-top box, a tablet computing device, an electronic book reader, a fixed or mobile media player, and so on. The computing device 250 includes: a communication interface 2501, a memory 2502 and a processor 2503; the various components are coupled together by a bus system 2504. It is to be appreciated that the bus system 2504 is employed to enable connected communications between these components. The bus system 2504 includes a power bus, a control bus, and a status signal bus in addition to the data bus. For clarity of illustration, the various buses are labeled as bus system 2504 in fig. 25. Wherein,,
the communication interface 2501 is configured to receive and send signals during the process of receiving and sending information with other external network elements;
the memory 2502 for storing a computer program capable of running on the processor 2503;
the processor 2503 is configured to perform the following steps when executing the computer program:
determining the fluid domain of a single flow channel in an axial flow compressor according to a geometric model of the axial flow compressor;
selecting a grid dividing method corresponding to the fluid domain of the single flow channel;
setting corresponding direction expanding parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division;
and generating a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the flow direction parameter and the circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
It is to be appreciated that the memory 2502 in embodiments of the application can be either volatile memory or nonvolatile memory, or can include both volatile and nonvolatile memory. The nonvolatile Memory may be a Read-Only Memory (ROM), a Programmable ROM (PROM), an Erasable PROM (EPROM), an Electrically Erasable EPROM (EEPROM), or a flash Memory. The volatile memory may be random access memory (Random Access Memory, RAM) which acts as an external cache. By way of example, and not limitation, many forms of RAM are available, such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (Double Data Rate SDRAM), enhanced SDRAM (ESDRAM), synchronous DRAM (SLDRAM), and Direct RAM (DRRAM). The memory 2502 of the systems and methods described herein is intended to comprise, without being limited to, these and any other suitable types of memory.
While processor 2503 may be an integrated circuit chip with signal processing capabilities. In implementation, the steps of the methods described above may be performed by integrated logic circuitry in hardware or instructions in software in the processor 2503. The processor 2503 described above may be a general purpose processor, a digital signal processor (Digital Signal Processor, DSP), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a field programmable gate array (Field Programmable Gate Array, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present application may be embodied directly in the execution of a hardware decoding processor, or in the execution of a combination of hardware and software modules in a decoding processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in memory 2502, and processor 2503 reads information from memory 2502 and performs the steps of the method described above in conjunction with its hardware.
It is to be understood that the embodiments described herein may be implemented in hardware, software, firmware, middleware, microcode, or a combination thereof. For a hardware implementation, the processing units may be implemented within one or more application specific integrated circuits (Application Specific Integrated Circuits, ASIC), digital signal processors (Digital Signal Processing, DSP), digital signal processing devices (DSP devices, DSPD), programmable logic devices (Programmable Logic Device, PLD), field programmable gate arrays (Field-Programmable Gate Array, FPGA), general purpose processors, controllers, microcontrollers, microprocessors, other electronic units designed to perform the functions described herein, or a combination thereof.
For a software implementation, the techniques described herein may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes may be stored in a memory and executed by a processor. The memory may be implemented within the processor or external to the processor.
Specifically, the processor 2503 is further configured to execute the method steps of CFD grid generation of the axial flow compressor in the foregoing technical solution when running the computer program, which is not described herein.
It should be noted that: the technical schemes described in the embodiments of the present application may be arbitrarily combined without any collision.
The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (10)
1. A method of axial compressor CFD grid generation, the method comprising:
according to a geometric model of the axial flow compressor, determining a fluid domain of a single flow passage in the axial flow compressor through the shape of the single flow passage in the axial flow compressor based on an OH type topological method;
selecting a grid dividing method corresponding to the fluid domain of the single flow channel;
setting corresponding direction expanding parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division;
and generating a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the relative thickness of the O blocks in the flow direction parameter, the number of O block circumferential nodes, the O block circumferential distribution coefficient and the circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
2. The method of claim 1, wherein determining the fluid domain of a single flow channel in an axial flow compressor based on a geometric model of the axial flow compressor comprises:
according to the geometric model of the axial flow compressor, determining an annular area surrounded by an inlet surface, an outlet surface, a hub surface, a casing surface and a blade surface of the axial flow compressor as a full flow channel;
obtaining the single flow channel based on the full flow channel division;
determining the shape of the single flow channel;
and obtaining the fluid domain corresponding to the single flow channel based on the shape topology of the single flow channel.
3. The method of claim 2, wherein the determining the shape of the single flow channel comprises:
the shape of the flow channel from the inlet throat to the outlet throat is determined by the shape of the blade;
the shape of the flow channel at the leading edge inlet and the trailing edge outlet is determined according to the set control point.
4. A method according to claim 3, wherein the control point set at the leading edge inlet comprises:
the first control point is positioned on an extension line of the throat flow direction, and the position coordinate of the first control point is determined by using a first parameter A in the following formula:
wherein, the value range of A is floating point number of [0,1 ];
the second control point is positioned at the left side of the front edge of the blade, and the position coordinates of the second control point are respectively determined along the two directions of the flow direction and the circumferential direction by using a second parameter B and a third parameter C in the following formula:
the third control point is located at the left side of the second control point and is used for controlling the tangential direction at the second control point, and the position coordinate of the third control point is determined by using a fourth parameter D in the following formula:
。
5. the method according to claim 1, wherein the selecting the meshing method corresponding to the fluid domain of the single flow channel includes:
and a method of multiple meshing is adopted for the fluid domain of the single flow channel.
6. The method of claim 1, wherein the setting of the spanwise parameter, the flow direction parameter, and the circumferential parameter corresponding to the meshing of the fluid domains of the single flow channel comprises:
setting corresponding direction-spreading network parameters when the fluid domain of the single flow channel is subjected to grid division;
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades;
and setting the flow direction parameter and the circumferential parameter.
7. The method of claim 6, wherein the setting the corresponding spanwise parameter for meshing the fluid domain of the single flow channel comprises:
setting the corresponding spanwise grid width of the fluid domain of the single flow channel to be in an equal-ratio array when grid division is carried out, and calculating according to the following formula to obtain the spanwise parameters：
Wherein S represents a growth factor;representing a first layer mesh width;
wherein, when i takes on the value of N,
wherein L represents half of the height of the flow channel; n represents half the number of meshes; and is also provided with。
8. The method of claim 6, wherein initializing the flow direction parameters and the circumferential parameters based on the spanwise parameters and blade geometry model characteristics comprises:
initializing the flow direction parameters and the circumferential parameters according to the spanwise parameters and the geometric model characteristics of the blades, wherein the initialization principle is as follows: making the average value of the spanwise mesh length equal to the average value of the flow direction and the circumferential mesh length;
wherein the flow direction parameters include: the number of H1 block flow direction nodes, the number of H4 block flow direction nodes, the number of H2 block circumferential nodes, the number of H3 block circumferential nodes, the number of O block circumferential nodes, the O block circumferential distribution coefficient, the O block circumferential starting node, the O block relative thickness, the number of O block normal nodes, the O block normal growth factor, the inlet block grid-pitch ratio, the inlet block node number, the outlet block grid-pitch ratio and the outlet block node number.
9. The method according to claim 8, wherein the generating a preview grid corresponding to the fluid domain of the single flow channel by using the selected meshing method according to the set spanwise parameter, the flow direction parameter, and the circumferential parameter, and generating a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement includes:
determining leaf height parameters;
acquiring the flow channel shape and the blade shape corresponding to the blade height;
inputting the relative thickness of the O blocks, and dividing the O blocks of the OH topology by using the relative thickness of the O blocks;
dividing the O block grid according to the number of O block ring nodes and the O block ring distribution coefficient;
setting other circumferential parameters, and drawing to obtain a preview grid corresponding to the fluid domain of the single flow channel;
when the preview grid meets the set requirement, regenerating a 3D grid corresponding to the fluid domain of the single flow channel by using the spanwise parameter, the circumferential parameter and the blade geometric model characteristic;
and solving a poisson equation to smooth the 3D grid corresponding to the fluid domain of the single flow channel.
10. The CFD grid generating device for the axial flow compressor is characterized by comprising a determining part, a selecting part, a setting part and a grid dividing part; wherein,,
the determining part is configured to determine the fluid domain of a single flow passage in the axial flow compressor through the shape of the single flow passage in the axial flow compressor based on an OH type topology method according to a geometric model of the axial flow compressor;
the selecting part is configured to select a grid dividing method corresponding to the fluid domain of the single flow channel;
the setting part is configured to set corresponding spanwise parameters, flow direction parameters and circumferential parameters when the fluid domain of the single flow channel is subjected to grid division;
the grid dividing part is configured to generate a preview grid corresponding to the fluid domain of the single flow channel by adopting the selected grid dividing method according to the set spanwise parameter, the relative thickness of the O blocks in the flow direction parameter, the number of O block circumferential nodes, the O block circumferential distribution coefficient and the circumferential parameter, and generate a 3D grid corresponding to the fluid domain of the single flow channel when the preview grid meets the set requirement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310741561.6A CN116484771B (en) | 2023-06-21 | 2023-06-21 | Method and device for generating CFD grid of axial flow compressor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310741561.6A CN116484771B (en) | 2023-06-21 | 2023-06-21 | Method and device for generating CFD grid of axial flow compressor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116484771A CN116484771A (en) | 2023-07-25 |
| CN116484771B true CN116484771B (en) | 2023-08-25 |
Family
ID=87221823
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310741561.6A Active CN116484771B (en) | 2023-06-21 | 2023-06-21 | Method and device for generating CFD grid of axial flow compressor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116484771B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117669429B (en) * | 2024-01-31 | 2024-05-10 | 陕西空天信息技术有限公司 | Fluid simulation method and device for rotary machine and computer storage medium |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8061988B1 (en) * | 2007-11-19 | 2011-11-22 | Florida Turbine Technologies, Inc. | Process for conjugate CFD analysis of a turbine airfoil |
| KR20130143227A (en) * | 2012-06-21 | 2013-12-31 | 한국건설기술연구원 | Method for multi-site calibration on unitary watershed system |
| CN104166752A (en) * | 2014-07-04 | 2014-11-26 | 贵州大学 | Hydraulic torque converter full-flow-channel transient numerical simulation calculation method |
| KR101877167B1 (en) * | 2017-03-20 | 2018-08-07 | 인하대학교 산학협력단 | Axial compressor |
| WO2021246931A1 (en) * | 2020-06-03 | 2021-12-09 | Ingrid Cloud Ab | Node for providing a simulation model, system for illustrating estimated fluid movements around a structure, methods therefore and a computer program product |
| CN114492220A (en) * | 2021-11-23 | 2022-05-13 | 华北电力大学 | OpenFOAM platform-based method and system for predicting flow characteristics of ocean atmospheric boundary layer |
| CN114936484A (en) * | 2022-04-11 | 2022-08-23 | 中建八局第三建设有限公司 | Urban park ecological water system efficiency design simulation analysis method |
| WO2022257308A1 (en) * | 2021-06-10 | 2022-12-15 | 大连理工大学 | Joint simulation-based fuel assembly multi-subject structure design optimization method |
| CN115906718A (en) * | 2023-03-09 | 2023-04-04 | 陕西空天信息技术有限公司 | Rotating machinery CFD system |
| CN115935837A (en) * | 2022-11-10 | 2023-04-07 | 天津大唐国际盘山发电有限责任公司 | Numerical calculation method for simulating boiler tail flue gas flowing process |
| CN116244865A (en) * | 2023-05-11 | 2023-06-09 | 陕西空天信息技术有限公司 | Method and device for finite element modeling of axial flow impeller and computer storage medium |
| CN116244848A (en) * | 2023-02-10 | 2023-06-09 | 西安交通大学 | Semi-open centrifugal impeller parametric modeling generation method and grid topology method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6220797B2 (en) * | 2012-03-08 | 2017-10-25 | エンジン シミュレーション パートナーズ | Fluid dynamics system boundaries |
| US20140257765A1 (en) * | 2013-03-05 | 2014-09-11 | Livermore Software Technology Corporation | Numerical Simulation of FSI Using The Space-Time CE/SE Solver With A Moving Mesh For The Fluid Domain |
| EP3131025A1 (en) * | 2015-08-14 | 2017-02-15 | Siemens Aktiengesellschaft | Method for the prediction of surge in a gas compressor |
-
2023
- 2023-06-21 CN CN202310741561.6A patent/CN116484771B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8061988B1 (en) * | 2007-11-19 | 2011-11-22 | Florida Turbine Technologies, Inc. | Process for conjugate CFD analysis of a turbine airfoil |
| KR20130143227A (en) * | 2012-06-21 | 2013-12-31 | 한국건설기술연구원 | Method for multi-site calibration on unitary watershed system |
| CN104166752A (en) * | 2014-07-04 | 2014-11-26 | 贵州大学 | Hydraulic torque converter full-flow-channel transient numerical simulation calculation method |
| KR101877167B1 (en) * | 2017-03-20 | 2018-08-07 | 인하대학교 산학협력단 | Axial compressor |
| WO2021246931A1 (en) * | 2020-06-03 | 2021-12-09 | Ingrid Cloud Ab | Node for providing a simulation model, system for illustrating estimated fluid movements around a structure, methods therefore and a computer program product |
| WO2022257308A1 (en) * | 2021-06-10 | 2022-12-15 | 大连理工大学 | Joint simulation-based fuel assembly multi-subject structure design optimization method |
| CN114492220A (en) * | 2021-11-23 | 2022-05-13 | 华北电力大学 | OpenFOAM platform-based method and system for predicting flow characteristics of ocean atmospheric boundary layer |
| CN114936484A (en) * | 2022-04-11 | 2022-08-23 | 中建八局第三建设有限公司 | Urban park ecological water system efficiency design simulation analysis method |
| CN115935837A (en) * | 2022-11-10 | 2023-04-07 | 天津大唐国际盘山发电有限责任公司 | Numerical calculation method for simulating boiler tail flue gas flowing process |
| CN116244848A (en) * | 2023-02-10 | 2023-06-09 | 西安交通大学 | Semi-open centrifugal impeller parametric modeling generation method and grid topology method |
| CN115906718A (en) * | 2023-03-09 | 2023-04-04 | 陕西空天信息技术有限公司 | Rotating machinery CFD system |
| CN116244865A (en) * | 2023-05-11 | 2023-06-09 | 陕西空天信息技术有限公司 | Method and device for finite element modeling of axial flow impeller and computer storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 气膜冷却叶片热冲击分析网格生成技术;关鹏;艾延廷;石小江;包天南;;热能动力工程(10);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116484771A (en) | 2023-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116484771B (en) | Method and device for generating CFD grid of axial flow compressor | |
| CN115329696B (en) | Conservation type fixed wall boundary numerical simulation method and device based on non-body-attached grid | |
| CN116484772B (en) | Loss acquisition method, device, equipment and medium for through-flow design | |
| CN110914869B (en) | Method and system for automated generation of constraint curves in computer graphics | |
| CN116050031B (en) | Axial flow impeller blade design method and device, storage medium and electronic equipment | |
| CN116415377B (en) | She Panmo type generation method and device, electronic equipment and storage medium | |
| CN111176210A (en) | Numerical controller, machining path setting method, and computer-readable medium containing program | |
| CN114060102B (en) | Method and device for determining outlet metal angle of guide vane blade | |
| CN116502373B (en) | Runner calculation grid generation method and device | |
| CN117216910B (en) | Centrifugal impeller wheel disc model construction method and device and electronic equipment | |
| CN117171924A (en) | Blade design method, device, equipment and computer storage medium | |
| CN116244865B (en) | Method and device for finite element modeling of axial flow impeller and computer storage medium | |
| CN115563898A (en) | Solid wall boundary numerical simulation method based on non-body-attached grid and related equipment | |
| CN114840926B (en) | Method for generating three-dimensional finite element grid with complex tire patterns | |
| JP4981313B2 (en) | Three-dimensional shape processing apparatus, curved surface creation program, and curved surface creation method | |
| CN116305669B (en) | Method, device, equipment and medium for acquiring design parameters of axial flow impeller | |
| CN117669429B (en) | Fluid simulation method and device for rotary machine and computer storage medium | |
| CN117688699B (en) | Extension section runner calculation grid construction method and device | |
| CN113065292B (en) | Boundary model method for integrated turbine transition section unsteady flow calculation | |
| CN117315083B (en) | Electronic cam motion curve generation method and device | |
| CN108898656A (en) | A kind of free curve in three-dimensional modeling efficiently produces method | |
| EP4375868A1 (en) | Physics-aware meshing techniques for generating grids for objects | |
| CN114856930B (en) | Determination method and device for pneumatic parameters of double-wind-wheel wind turbine generator set and electronic equipment | |
| CN115719403B (en) | Method and device for constructing tool path, electronic equipment and storage medium | |
| CN117195593B (en) | Method, device, equipment and medium for acquiring gas flow parameters of blade |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |