CN117540507A - Calculation method and device for wall surface distance of full-ring grid - Google Patents
Calculation method and device for wall surface distance of full-ring grid Download PDFInfo
- Publication number
- CN117540507A CN117540507A CN202410023043.5A CN202410023043A CN117540507A CN 117540507 A CN117540507 A CN 117540507A CN 202410023043 A CN202410023043 A CN 202410023043A CN 117540507 A CN117540507 A CN 117540507A
- Authority
- CN
- China
- Prior art keywords
- space grid
- distance
- wall
- full
- point
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 28
- 238000013507 mapping Methods 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims description 46
- 238000004590 computer program Methods 0.000 claims description 12
- 230000003362 replicative effect Effects 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
The application discloses a calculation method and a calculation device for a full-ring grid wall surface distance, wherein a first space grid point in a single channel of a single rotor is obtained; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
Description
Technical Field
The application belongs to the technical field of hydrodynamics, and particularly relates to a calculation method and a calculation device for a full-ring grid wall surface distance.
Background
Impeller machines are widely used in the energy power industry, such as aircraft engines, gas turbines, and the like. The flow field of the impeller machine is analyzed, and the flow field can be used for evaluating the working characteristics of the impeller machine and improving the design of the impeller machine. Compared to the aerodynamic test of the impeller machinery, three-dimensional simulation of the flow field of the impeller machinery based on CFD (Computational Fluid Dynamics ) has been largely applied to the analysis of the flow field of the impeller machinery due to its lower cost and faster speed.
Single rotors are used as the basic component of impeller machines, and extensive numerical simulation studies have been conducted in academia and industry. In simulating certain working conditions, such as near stall working conditions, unstable working conditions and the like, the flow field structure occupies a plurality of blade channels, and simulation is generally required to be performed based on a full ring network lattice.
The flow state of the flow field is classified into laminar flow and turbulent flow, and the real flow in the nature is basically a turbulent flow problem. Currently, when turbulence is simulated based on CFD, the simulation is generally performed based on a turbulence model. Common turbulence models include SA models, SST models, etc., all of which require calculation of wall distances, i.e., the closest distance of a spatial point to the wall.
The current methods for calculating the wall distance are mainly of two types: one is to solve partial differential equations, and the other is to calculate directly by geometric methods.
The method for solving the wall surface distance based on the partial differential equation is difficult to apply to complex shapes. The method mainly aims at solving the partial differential equation on the complex appearance, has poor stability, and obtains an approximate solution, and the accuracy of the wall surface distance directly influences the stability, the convergence and the accuracy of the calculation result of the CFD solver.
In the geometric method for directly calculating the wall distance, the first is a direct search method, and the second is a hierarchical search method.
For the direct search method, this is achieved by calculating the distance between the spatial point and each wall point and taking the minimum value. Let the number of spatial points be m, the number of wall points be n, the time complexity of such a method be O (m x n).
For the hierarchical search method, the mainstream algorithm is the box method. For the box method, wall points are generally classified into different boxes based on a Cartesian coordinate range. When the wall surface distance of the space point is calculated, a plurality of boxes closest to the space point are found first, and then the wall surface points in the boxes are traversed, so that the shortest distance is calculated.
When the existing calculation method is adopted to calculate the wall distance of the single-rotor full-ring network grid, the calculation efficiency is low.
Disclosure of Invention
The invention aims to provide a method, a device, a terminal device and a storage medium for calculating the wall surface distance of a full-ring grid, so as to solve the defects in the prior art.
In a first aspect, an embodiment of the present invention provides a method for calculating a wall distance of a full-ring grid, where the method includes:
acquiring a first space grid point in a single channel of a single rotor;
calculating the distance from the first space grid point to each wall point of the single rotor;
determining the minimum distance among the distances of the wall points as a target wall distance;
and mapping the target wall distance to the full-ring space grid point according to a pre-established mapping relation between the single-channel space grid point and the full-ring space grid point.
Optionally, the mapping relationship between the pre-established single-channel inner space grid point and the full-channel inner space grid point is obtained by the following manner:
according to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel to generate full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
Optionally, the rotationally replicating the spatial grid points in the single channel to generate full-loop spatial grid points of the single rotor includes:
determining a number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first and second spatial grid points differ by one or more blade periods, and the first and second blade channels are different blade channels;
determining a first wall point and calculating a first distance between the first wall point and the first spatial grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and the second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between a first space grid point and a second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
Optionally, the mapping relationship between the single-channel inner space grid point and the full-ring space grid point is one-to-many, and the mapping relationship between the full-ring space grid point and the single-channel inner space grid point is one-to-many.
Optionally, the method further comprises:
the number of each single-channel space grid in the full-ring space grid points is established.
In a second aspect, an embodiment of the present invention provides a computing device for a wall distance of a full-ring grid, where the device includes:
the acquisition module is used for acquiring a first space grid point in a single channel of the single rotor;
a calculation module for calculating the distance from the first space grid point to each wall point of the single rotor;
the determining module is used for determining the minimum distance in the distances of the wall points as a target wall distance;
and the mapping module is used for mapping the target wall distance to the full-ring space grid point according to the pre-established mapping relation between the single-channel space grid point and the full-ring space grid point.
Optionally, the mapping module is configured to:
according to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel to generate full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
Optionally, the mapping module is configured to:
determining a number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first and second spatial grid points differ by one or more blade periods, and the first and second blade channels are different blade channels;
determining a first wall point and calculating a first distance between the first wall point and the first spatial grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and the second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between a first space grid point and a second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
Optionally, the mapping relationship between the single-channel inner space grid point and the full-ring space grid point is one-to-many, and the mapping relationship between the full-ring space grid point and the single-channel inner space grid point is one-to-many.
Optionally, the mapping module is further configured to:
the number of each single-channel space grid in the full-ring space grid points is established.
In a third aspect, an embodiment of the present invention provides a terminal device, including: at least one processor and memory;
the memory stores a computer program; the at least one processor executes the computer program stored in the memory to implement the method for calculating the full-loop mesh wall surface distance provided in the first aspect.
In a fourth aspect, an embodiment of the present invention provides a computer readable storage medium, where a computer program is stored, where the computer program is executed to implement the method for calculating the wall surface distance of the full-ring grid provided in the first aspect.
The embodiment of the invention has the following advantages:
the embodiment of the invention provides a calculation method, a calculation device, terminal equipment and a storage medium for the wall surface distance of a full-ring grid, wherein a first space grid point in a single channel of a single rotor is obtained; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
Drawings
In order to more clearly illustrate the embodiments or prior art solutions of the present application, the drawings that are required for the description of the embodiments or prior art will be briefly described below, it being apparent that the drawings in the following description are only some of the embodiments described in the present application, and that other drawings may be obtained according to these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a flowchart of a method for calculating a full-loop mesh wall distance according to an embodiment of the present disclosure;
FIG. 2 is a diagram illustrating classification of single rotor full ring grid walls according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a space grid in an embodiment of the present application;
FIG. 4 is a schematic illustration of an inlet face and an outlet face in an embodiment of the present application;
FIG. 5 is a schematic diagram of a single-bladed channel spatial grid and a full-annular spatial grid in an embodiment of the present application;
FIG. 6 is a schematic view of a single-vane channel and full annular wall in an embodiment of the present application;
FIG. 7 is a block diagram of an embodiment of a computing device for full ring grid wall distance in accordance with the present invention;
fig. 8 is a schematic structural view of a terminal device of the present invention.
Detailed Description
For the purposes, technical solutions and advantages of the present application, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
An embodiment of the invention provides a calculation method of a full-ring grid wall surface distance, which is used for calculating the wall surface distance of a rotor. The execution body of the embodiment is a computing device of the full-ring grid wall surface distance, and the computing device is arranged on a terminal device, for example, the terminal device at least comprises a computer terminal and the like.
Referring to fig. 1, a step flow chart of an embodiment of a method for calculating a full-loop mesh wall surface distance according to the present invention is shown, and the method specifically may include the following steps:
s101, acquiring a first space grid point in a single channel of a single rotor;
specifically, the embodiment of the invention calculates the wall surface distance for a single rotor, wherein the single rotor has a plurality of single channels, each single channel has an independent space region, and a plurality of space grid points, namely a first space grid point, exist in the space region.
S102, calculating the distance from the first space grid point to each wall point of the single rotor;
specifically, the single rotor comprises a plurality of blades, the terminal equipment acquires wall points on the single blade, and the terminal equipment calculates the distances between the first space grid point and the plurality of wall points according to the space coordinates of the first space grid point and the space coordinates of the wall points, so that the plurality of wall distances can be obtained.
The wall distance of a space grid point (typically, the cell center of a space grid cell) is the minimum distance to all wall grid points (hereinafter, referred to as wall points).
S103, determining the minimum distance among the distances of all the wall points as a target wall distance;
specifically, the terminal device selects a minimum value among the wall distances among the plurality of wall distances, and determines the wall distance minimum value as the target wall distance.
S104, mapping the target wall distance to the full-ring space grid point according to the pre-established mapping relation between the single-channel space grid point and the full-ring space grid point.
Specifically, according to the geometric periodic characteristics of the single rotor, a mapping relation between space grid points in a single channel and full-ring space grid points is established, only a single-channel target wall surface distance is calculated, and then the single-channel target wall surface distance is mapped to the full-ring space grid points, namely, a certain single-rotor full-ring grid is taken as an example object; firstly, generating a single-channel grid of the single rotor; generating a full ring grid of a single rotor through rotary replication, and recording the mapping relation between a single-channel grid unit and the full ring grid unit in the process; extracting wall points of the full-ring grid; calculating the minimum distance from the single-channel grid to the wall points of the full-ring grid; mapping the wall distance of the single-channel space grid to the full-ring space grid; in this way, the calculation efficiency is improved.
According to the calculation method of the full-ring grid wall surface distance, the first space grid points in a single channel of a single rotor are obtained; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
The method for calculating the wall surface distance of the full-ring grid provided by the embodiment of the invention is further described in a further embodiment of the invention.
Fig. 2 is a schematic diagram illustrating classification of single-rotor full-ring grid wall surfaces according to an embodiment of the present application, where the single-rotor wall surfaces may be classified into hub wall surfaces, casing wall surfaces, and blade wall surfaces. A schematic of the spatial grid is shown in fig. 3. Each single channel of the single rotor is a distributed area of the spatial grid through the enclosed space enveloped by the inlet face and the outlet face, as shown in fig. 4.
The wall distance of a space grid point (typically, the cell center of a space grid cell) is the minimum distance to all wall grid points (hereinafter, referred to as wall points).
The full ring mesh is replicated circumferentially from a single-blade channel mesh, as shown in FIG. 5. The relationship between the single channel mesh and the full ring mesh is shown in fig. 6.
Optionally, the pre-established mapping relationship between the space grid points in the single channel and the space grid points in the full-channel is obtained by the following way:
according to the single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel, generating full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
Optionally, rotationally replicating the spatial grid points within the single channel to generate full-ring spatial grid points of the single rotor, including:
determining the number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first spatial grid point and the second spatial grid point differ by one or more vane cycles, and the first vane passage and the second vane passage are different vane passages;
determining a first wall point and calculating a first distance between the first wall point and a first space grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and a second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between the first space grid point and the second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
Optionally, the mapping relationship between the single-channel inner space grid point and the full-ring space grid point is one-to-many, and the mapping relationship between the full-ring space grid point and the single-channel inner space grid point is one-to-many.
Specifically, it is assumed that the number of blades of a single rotor is N (N>1) Within a vane passage there are spatial grid points P whose coordinates in a cylindrical coordinate system are. Due to the periodicity of the single rotor space grid, there are space points +.>. I.e., P1 and P are in different blade passes, but the coordinates differ by one or more blade cycles.
For any wall pointThe distance between the P point and the A point is PA, and the wall point is necessarily present due to the periodicity of the single rotor wall grid>At this time, the distance between P1 and A1 is P1A1;
wherein:
since PA and P1A1 are both greater than 0, there are:
---（3）
let the set of all wall points be V:
---（4）
wherein m is the number of wall points.
With space grid points in a certain vane passageThe set U of distance values from all wall points is:
---（5）
space grid point for another blade channelThe set of distance values from all wall points is U1:
---（6）
for any one element PA in set U i The corresponding wall point is A i ，A i In the set of wall points V. Let A i Is the coordinates of (a). Due to the periodicity of the single rotor wall grid, there must be wall points +.>。A j In the set of wall points V.
Then P1 and A j Distance P1A of (2) j In set U1.
The method comprises the following steps:
;
;
due to P1A j =PA i So for any one element PA in set U i Necessarily there is an element P1A equal to it j In set U1. So there are:
---（7）
similarly, for any one element P1A in the set U1 k The corresponding wall point is A k ，A k In the set of wall points V. Let A k Is the coordinates of (a). Due to the periodicity of the single rotor wall grid, there must be wall points. Thus, A l In the set of wall points V.
Then P and A l Distance PA of (2) l In set U.
As a result of:
;
;
due to P1A k =PA l So for any one element P1A in the set U1 k Necessarily there is an element PA equal to it l In set U. So there are:
---（8）
so that:
---（9）
the wall distance of the spatial grid point P is the minimum value of all the distance values in the set U, and the wall distance of the spatial grid point P1 is the minimum value of all the distance values in the set U1. Therefore, P and P1 are equal to the wall distance.
Therefore, only the wall distances of all the space grid points of one blade channel are calculated, and the wall distances of the space grid points of other blade channels can be obtained.
Step 2: and establishing a mapping relation between the single-channel grid cell and the full-ring grid cell.
Because the full-ring grid is formed by copying a single-channel grid along the circumferential direction, the mapping relation between the two space grid points can be established in the copying process. Considering that a certain space grid point in the single-channel grid can be mapped to the space grid point at the same position of each channel of the whole channel, the mapping relation between the space grid point in the single channel and the space grid point in the whole channel is one-to-many, and the mapping relation between the space grid point in the whole channel and the space grid point in the single channel is one-to-many.
Let ncSinglePassage be the number of a certain spatial grid point of the single channel. After obtaining full-ring grids based on single-channel grids along circumferential direction, in channels of different blades of the full-ring grids, grid point numbers obtained by copying the space grid points are { ncFullAnnulus }, respectively 1 ,…,ncFullAnnulus N N is the number of blades of a single rotor. Establishing a mapping relation between the full-ring space grid points and the corresponding single-channel space grid points:
(ncSinglePassage=SPToFAProjection(ncFullAnnulus 1 )
…
ncSinglePassage=SPToFAProjection(ncFullAnnulus N ) —(10)
SPToFAProjections are mapping functions of full-loop space grid points to corresponding single-channel space grid points. Step 3: and calculating to obtain the wall surface distance of the space grid point in the single channel.
As shown in fig. 6, the wall surface grid is a full-ring grid, and based on the existing general wall surface distance algorithm, the wall surface distance of the space grid point in the single channel can be obtained.
Step 4: the wall distances of the space grid points in the single channel are mapped to the full-ring space grid points.
The wall distance of the ncfullane number space grid point of the full-loop grid is as follows: wallDistance (ncFullAnnulus).
The wall surface distance is calculated by the following formula:
ncSinglePassage=SPToFAProjection(ncFullAnnulus)；
wallDistance(ncFullAnnulus)=wallDistance(ncSinglePassage)---（11）
wallDistance (ncSinglePassage) is the wall distance of the space grid point of the ncSinglePassage number of the single-channel grid.
Optionally, the method further comprises:
the number of each single-channel space grid in the full-ring space grid points is established.
The embodiment of the invention provides a wall distance generation method for a single-rotor full-ring grid, which fully considers geometric periodic characteristics of an impeller mechanical space grid and a wall grid, only calculates the wall distance of the grid of one blade channel, and expands the wall distance of the space grid of all blade channels to obtain the wall distance.
It should be noted that, for simplicity of description, the method embodiments are shown as a series of acts, but it should be understood by those skilled in the art that the embodiments are not limited by the order of acts, as some steps may occur in other orders or concurrently in accordance with the embodiments. Further, those skilled in the art will appreciate that the embodiments described in the specification are presently preferred embodiments, and that the acts are not necessarily required by the embodiments of the invention.
According to the calculation method of the full-ring grid wall surface distance, the first space grid points in a single channel of a single rotor are obtained; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
Another embodiment of the present invention provides a computing device for a full-ring grid wall surface distance, configured to execute the computing method for a full-ring grid wall surface distance provided in the foregoing embodiment.
Referring to fig. 7, there is shown a block diagram of an embodiment of a full-loop mesh wall distance computing device of the present invention, which may include the following modules: an acquisition module 701, a calculation module 702, a determination module 703 and a mapping module 704, wherein:
the acquiring module 701 is configured to acquire a first spatial grid point in a single channel of a single rotor;
the calculating module 702 is configured to calculate a distance from the first spatial grid point to each wall point of the single rotor;
the determining module 703 is configured to determine a minimum distance among the distances of the wall points as a target wall distance;
the mapping module 704 is configured to map the target wall distance to the full-ring space grid point according to a pre-established mapping relationship between the single-channel space grid point and the full-ring space grid point.
The computing device for the full-ring grid wall surface distance provided by the embodiment of the invention obtains the first space grid point in the single channel of the single rotor; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
A further embodiment of the present invention further provides a computing device for the wall distance of the full-ring grid provided in the above embodiment.
Optionally, the mapping module is configured to:
according to the single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel, generating full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
Optionally, the mapping module is configured to:
determining the number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first spatial grid point and the second spatial grid point differ by one or more vane cycles, and the first vane passage and the second vane passage are different vane passages;
determining a first wall point and calculating a first distance between the first wall point and a first space grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and a second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between the first space grid point and the second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
Optionally, the mapping relationship between the single-channel inner space grid point and the full-ring space grid point is one-to-many, and the mapping relationship between the full-ring space grid point and the single-channel inner space grid point is one-to-many.
Optionally, the mapping module is further configured to:
the number of each single-channel space grid in the full-ring space grid points is established.
For the device embodiments, since they are substantially similar to the method embodiments, the description is relatively simple, and reference is made to the description of the method embodiments for relevant points.
The computing device for the full-ring grid wall surface distance provided by the embodiment of the invention obtains the first space grid point in the single channel of the single rotor; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
An embodiment of the present invention provides a terminal device, configured to execute the method for calculating the wall surface distance of the full-ring mesh provided in the above embodiment.
Fig. 8 is a schematic structural view of a terminal device of the present invention, as shown in fig. 8, the terminal device comprising: at least one processor 801 and memory 802;
the memory stores a computer program; the at least one processor executes the computer program stored in the memory to implement the method for calculating the wall distance of the full-loop grid provided in the above embodiment.
The terminal equipment provided by the embodiment obtains the first space grid point in the single channel of the single rotor; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
Still another embodiment of the present application provides a computer readable storage medium storing a computer program, where the computer program when executed implements the method for calculating the wall surface distance of the full-ring grid provided in any one of the above embodiments.
According to the computer-readable storage medium of the present embodiment, by acquiring a first spatial grid point within a single channel of a single rotor; calculating the distance from the first space grid point to each wall point of the single rotor; determining the minimum distance among the distances of all the wall points as a target wall distance; according to the pre-established mapping relation between the space grid points in the single channel and the full-ring space grid points, mapping the target wall distance to the full-ring space grid points, according to the geometric periodic characteristics of the space grid and the wall grid of the impeller machinery, only calculating the wall distance of the grid of one blade channel, and then expanding to obtain the wall distances of the space grid of all the blade channels, so that the calculation efficiency of the wall distances can be improved.
It should be noted that the foregoing detailed description is exemplary and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments in accordance with the present application. As used herein, the singular is intended to include the plural unless the context clearly indicates otherwise. Furthermore, it will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, steps, operations, devices, components, and/or groups thereof.
It should be noted that the terms "first," "second," and the like in the description and claims of the present application and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are capable of operation in sequences other than those illustrated or otherwise described herein.
Furthermore, the terms "comprise" and "have," as well as any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those elements but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways, such as rotated 90 degrees or at other orientations, and the spatially relative descriptors used herein interpreted accordingly.
In the above detailed description, reference is made to the accompanying drawings, which form a part hereof. In the drawings, like numerals typically identify like components unless context indicates otherwise. The illustrated embodiments described in the detailed description, drawings, and claims are not meant to be limiting. Other embodiments may be utilized, and other changes may be made, without departing from the spirit or scope of the subject matter presented herein.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method for calculating a full-loop mesh wall distance, the method comprising:
acquiring a first space grid point in a single channel of a single rotor;
calculating the distance from the first space grid point to each wall point of the single rotor;
determining the minimum distance among the distances of the wall points as a target wall distance;
and mapping the target wall distance to the full-ring space grid point according to a pre-established mapping relation between the single-channel space grid point and the full-ring space grid point.
2. The method according to claim 1, wherein the pre-established mapping relationship between the space grid points in the single channel and the space grid points in the full channel is obtained by:
according to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel to generate full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
3. The method of claim 2, wherein rotationally replicating the spatial grid points within the single channel to generate full-loop spatial grid points for a single rotor comprises:
determining a number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first and second spatial grid points differ by one or more blade periods, and the first and second blade channels are different blade channels;
determining a first wall point and calculating a first distance between the first wall point and the first spatial grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and the second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between a first space grid point and a second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
4. The method of claim 1, wherein the single-channel inner space grid point and the full-channel inner space grid point are mapped one-to-many, and the full-channel inner space grid point and the single-channel inner space grid point are mapped one-to-one.
5. The method according to claim 4, wherein the method further comprises:
the number of each single-channel space grid in the full-ring space grid points is established.
6. A computing device for full-loop mesh wall distance, the device comprising:
the acquisition module is used for acquiring a first space grid point in a single channel of the single rotor;
a calculation module for calculating the distance from the first space grid point to each wall point of the single rotor;
the determining module is used for determining the minimum distance in the distances of the wall points as a target wall distance;
and the mapping module is used for mapping the target wall distance to the full-ring space grid point according to the pre-established mapping relation between the single-channel space grid point and the full-ring space grid point.
7. The apparatus of claim 6, wherein the mapping module is configured to:
according to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
and carrying out rotary copying on the space grid points in the single channel to generate full-ring space grid points of the single rotor, and recording the mapping relation between the space grid points and the full-ring space grid points in the single channel.
8. The apparatus of claim 7, wherein the mapping module is configured to:
determining a number of blades of the single rotor;
a first space grid point is arranged in the first blade channel;
a second space grid point is arranged in the second blade channel; the coordinates of the first and second spatial grid points differ by one or more blade periods, and the first and second blade channels are different blade channels;
determining a first wall point and calculating a first distance between the first wall point and the first spatial grid point;
determining a second wall point according to the periodicity of the single rotor wall geometry, and calculating a second distance between the second wall point and the second space grid point;
determining a first distance minimum value as a target first distance according to the plurality of first distances;
determining a second distance minimum value as a target second distance according to the plurality of second distances;
under the condition that the first distance of the target is the same as the second distance of the target, establishing a corresponding relation between a first space grid point and a second space grid point;
and determining the full-ring space grid point of the single rotor according to the corresponding relation between the first space grid point and the second space grid point.
9. A terminal device, comprising: at least one processor and memory;
the memory stores a computer program; the at least one processor executes the computer program stored by the memory to implement the method of computing full-loop mesh wall distance of any one of claims 1-5.
10. A computer readable storage medium, wherein a computer program is stored in the computer readable storage medium, said computer program when executed implementing the method for calculating the full-loop mesh wall surface distance according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410023043.5A CN117540507B (en) | 2024-01-08 | 2024-01-08 | Calculation method and device for wall surface distance of full-ring grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410023043.5A CN117540507B (en) | 2024-01-08 | 2024-01-08 | Calculation method and device for wall surface distance of full-ring grid |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117540507A true CN117540507A (en) | 2024-02-09 |
| CN117540507B CN117540507B (en) | 2024-06-04 |
Family
ID=89782646
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410023043.5A Active CN117540507B (en) | 2024-01-08 | 2024-01-08 | Calculation method and device for wall surface distance of full-ring grid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117540507B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101882177A (en) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | Aeroelastic stability fluid-structure interaction prediction method of turbo-machine changed interblade phase angles |
| CN108733922A (en) * | 2018-05-20 | 2018-11-02 | 中国人民解放军海军工程大学 | Vacuole computational methods under a kind of PODDED PROPULSOR deflection operating mode |
| CN112100937A (en) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | Efficient overlapped grid parallel assembly method based on wall surface distance |
| CN112613126A (en) * | 2020-12-29 | 2021-04-06 | 中国航天空气动力技术研究院 | Wall surface distance calculation method and device applied to aircraft appearance |
| CN113515902A (en) * | 2021-06-23 | 2021-10-19 | 中国人民解放军国防科技大学 | Method and device for initializing full-loop unsteady simulation of inner flow of gas compressor |
| CN114756974A (en) * | 2022-06-13 | 2022-07-15 | 中国空气动力研究与发展中心计算空气动力研究所 | Wall distance calculation method considering object plane normal information |
| CN114943167A (en) * | 2022-05-31 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, system, medium and equipment for calculating wall surface distance of structural grid |
| CN115357849A (en) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Method and device for calculating wall surface distance under Cartesian grid |
| CN115470511A (en) * | 2022-11-15 | 2022-12-13 | 中国空气动力研究与发展中心计算空气动力研究所 | Cartesian grid generation method and device suitable for multi-component model |
| CN115859536A (en) * | 2023-03-02 | 2023-03-28 | 中国航发四川燃气涡轮研究院 | Asynchronous vibration frequency locking numerical simulation method for rotor blade of gas compressor |
-
2024
- 2024-01-08 CN CN202410023043.5A patent/CN117540507B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101882177A (en) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | Aeroelastic stability fluid-structure interaction prediction method of turbo-machine changed interblade phase angles |
| CN108733922A (en) * | 2018-05-20 | 2018-11-02 | 中国人民解放军海军工程大学 | Vacuole computational methods under a kind of PODDED PROPULSOR deflection operating mode |
| CN112100937A (en) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | Efficient overlapped grid parallel assembly method based on wall surface distance |
| CN112613126A (en) * | 2020-12-29 | 2021-04-06 | 中国航天空气动力技术研究院 | Wall surface distance calculation method and device applied to aircraft appearance |
| CN113515902A (en) * | 2021-06-23 | 2021-10-19 | 中国人民解放军国防科技大学 | Method and device for initializing full-loop unsteady simulation of inner flow of gas compressor |
| CN114943167A (en) * | 2022-05-31 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, system, medium and equipment for calculating wall surface distance of structural grid |
| CN114756974A (en) * | 2022-06-13 | 2022-07-15 | 中国空气动力研究与发展中心计算空气动力研究所 | Wall distance calculation method considering object plane normal information |
| CN115357849A (en) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Method and device for calculating wall surface distance under Cartesian grid |
| CN115470511A (en) * | 2022-11-15 | 2022-12-13 | 中国空气动力研究与发展中心计算空气动力研究所 | Cartesian grid generation method and device suitable for multi-component model |
| CN115859536A (en) * | 2023-03-02 | 2023-03-28 | 中国航发四川燃气涡轮研究院 | Asynchronous vibration frequency locking numerical simulation method for rotor blade of gas compressor |
Non-Patent Citations (9)
| Title |
|---|
| ANDREAS LUCIUS 等: "Numerical Simulation and Evaluation of Velocity Fluctuations During Rotating Stall of a Centrifugal Pump", JOURNAL OF FLUIDS ENGINEERING, vol. 133, no. 08, 16 August 2011 (2011-08-16), pages 081102 * |
| 周伟: "某膨胀循环发动机推力室冷却结构流场仿真分析", 火箭推进, vol. 41, no. 02, 15 April 2015 (2015-04-15), pages 63 - 69 * |
| 孙龙刚 等: "基于TBR模型的高水头混流式水轮机水力性能预测", 农业工程学报, vol. 35, no. 07, 8 April 2019 (2019-04-08), pages 62 - 69 * |
| 孟爽 等: "笛卡尔网格下精确高效的壁面距离计算方法", 空气动力学学报, vol. 41, no. 07, 28 April 2023 (2023-04-28), pages 93 - 101 * |
| 李广新 等: "基于流热固耦合燃气轮机压气机叶轮强度分析", 流体机械, vol. 45, no. 04, 30 April 2017 (2017-04-30), pages 28 - 32 * |
| 杨慧 等: "转静干涉对转子叶片颤振特性的影响", 北京航空航天大学学报, vol. 42, no. 02, 1 April 2015 (2015-04-01), pages 258 - 264 * |
| 赵钟 等: "湍流模拟壁面距离MPI/OpenMP混合并行计算方法", 空气动力学学报, vol. 37, no. 06, 15 December 2019 (2019-12-15), pages 883 - 892 * |
| 韩乐 等: "转子叶片气弹稳定性与强迫响应分析", 航空发动机, vol. 47, no. 04, 29 July 2021 (2021-07-29), pages 82 - 90 * |
| 麻全: "高水头混流式水轮机飞逸工况流动特性数值研究", 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑, no. 02, 15 February 2018 (2018-02-15), pages 037 - 210 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117540507B (en) | 2024-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Raciti Castelli et al. | Modeling strategy and numerical validation for a Darrieus vertical axis micro-wind turbine | |
| Potsdam et al. | Unstructured mesh CFD aerodynamic analysis of the NREL Phase VI rotor | |
| Van Zante et al. | Recommendations for achieving accurate numerical simulation of tip clearance flows in transonic compressor rotors | |
| Kersken et al. | Time-linearized and time-accurate 3D RANS methods for aeroelastic analysis in turbomachinery | |
| CN102938003B (en) | Method for predicting aeroelasticity stability numerical value of turbomachinery with error frequency included | |
| CN109753716B (en) | Nuclear/thermal power turboset fluid excitation numerical calculation method and system based on flow field simulation | |
| Doubrawa et al. | Optimization-based calibration of FAST. Farm parameters against large-eddy simulations | |
| CN105184683A (en) | Probability clustering method based on wind electric field operation data | |
| CN115774903A (en) | Helicopter rotor disc load real-time generation method, device, equipment and medium | |
| Rezvaya et al. | Solving the hydrodynamical tasks using CFD programs | |
| Digraskar | Simulations of flow over wind turbines | |
| Turgut et al. | Factors influencing computational predictability of aerodynamic losses in a turbine nozzle guide vane flow | |
| CN117540507B (en) | Calculation method and device for wall surface distance of full-ring grid | |
| Young et al. | Accounting for eccentricity in compressor performance predictions | |
| Jin | Numerical simulation of wind turbine wakes based on actuator line method in NEK5000 | |
| Sergio et al. | Optimization methodology assessment for the inlet velocity profile of a hydraulic turbine draft tube: part II—performance evaluation of draft tube model | |
| Pillai et al. | Modelling wind turbine wakes at Middelgrunden wind farm | |
| Schuff et al. | Coupled Mode Flutter Analysis of Turbomachinery Blades Using an Adaptation of the p–k Method | |
| Kaplan et al. | Reverse engineering design of a hydraulic turbine runner | |
| Wilkinson et al. | Approximating urban wind interference | |
| Engels-Putzka et al. | Forced response sensitivity analysis using an adjoint harmonic balance solver | |
| Zieße et al. | Validation of a Synthetic-Eddy Method for Modelling Incoming Wakes in Scale-Resolving Simulations of Turbomachinery Cascades | |
| Marta et al. | Blade shape Optimization using a RANS Discrete Adjoint solver | |
| CN117454721B (en) | Wind power plant wake superposition effect evaluation method and medium based on digital simulation experiment | |
| Rodriguez et al. | A new inlet distortion and pressure loss based design of an intake system for stationary gas turbines |
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 |