CN117540507B - 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 PDF

Info

Publication number
CN117540507B
CN117540507B CN202410023043.5A CN202410023043A CN117540507B CN 117540507 B CN117540507 B CN 117540507B CN 202410023043 A CN202410023043 A CN 202410023043A CN 117540507 B CN117540507 B CN 117540507B
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.)
Active
Application number
CN202410023043.5A
Other languages
Chinese (zh)
Other versions
CN117540507A (en
Inventor
王子维
李彬
赵辉
曹杰
李白羽
张威龙
刘瑶
黄凯
吕尧兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202410023043.5A priority Critical patent/CN117540507B/en
Publication of CN117540507A publication Critical patent/CN117540507A/en
Application granted granted Critical
Publication of CN117540507B publication Critical patent/CN117540507B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Calculation method and device for wall surface distance of full-ring grid
Technical Field
The application belongs to the technical field of hydrodynamics, and particularly relates to a method and a device for calculating 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, the 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 the lower cost and the 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 of the application or the prior art solutions, the drawings which are used in the description of the embodiments or the prior art will be briefly described below, it being obvious that the drawings in the description below are only some of the embodiments described in the present application, and that other drawings can 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 application;
FIG. 2 is a diagram illustrating a single rotor full ring grid wall classification in accordance with an embodiment of the present application;
FIG. 3 is a schematic diagram of a space grid in accordance with an embodiment of the present application;
FIG. 4 is a schematic illustration of the inlet and outlet faces in accordance with 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 accordance with an embodiment of the present application;
FIG. 6 is a schematic view of a single-vane channel and full annular wall in accordance with 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
In order to make the objects, technical solutions and advantages of the present application more apparent, 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, embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
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 diagram illustrating classification of single-rotor full-ring grid walls according to an embodiment of the present application, wherein the single-rotor walls may be classified into hub walls, casing walls and blade walls. 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, assuming that the number of blades of a single rotor is N (N > 1), there is a space grid point P in a certain blade channel, whose coordinates in a cylindrical coordinate system are. Due to the periodicity of the single rotor space grid, there are spatial points/>, within a different blade channel. 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 existed due to the periodicity of the single rotor wall gridAt 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 channel The set of distance values from all wall points is U1:
---（6）
for any element PA i in set U, its corresponding wall point is a i,Ai in wall point set V. Let the coordinates of A i be . Due to the periodicity of the single rotor wall grid, there must be wall points. A j is in the set of wall points V.
Then the distance P1a j from P1 to a j is in set U1.
The method comprises the following steps:
;
;
since P1A j=PAi, for any element PA i in set U, there must be one element P1A j in set U1 equal to it. So there are:
---（7）
Similarly, for any element P1A k in set U1, its corresponding wall point is A k,Ak in set V of wall points. Let the coordinates of A k be . Due to the periodicity of the single rotor wall grid, there must be wall points. Thus, a l is in the set of wall points V.
Then the distance PA l of P from a l is in set U.
As a result of:
;
;
Because of P1A k=PAl, for any one element P1A k in set U1, there must be one element PA l in set U equal to it. 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 the full-ring grid is obtained by copying the single-channel grid along the circumferential direction, grid point numbers obtained by copying the space grid points in channels of different blades of the full-ring grid are { ncFullAnnulus 1,…,ncFullAnnulusN }, and N is the number of blades of the single rotor. Establishing a mapping relation between the full-ring space grid points and the corresponding single-channel space grid points:
(ncSinglePassage=SPToFAProjection(ncFullAnnulus1)
…
ncSinglePassage=SPToFAProjection(ncFullAnnulusN) —(10)
SPToFAProjection is a mapping function of the full-ring space grid points to the 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 surface distance of ncFullAnnulus # space grid points of the full-ring grid is as follows: WALLDISTANCE (NCFULLANNULUS).
The wall surface distance is calculated by the following formula:
ncSinglePassage=SPToFAProjection(ncFullAnnulus)；
wallDistance(ncFullAnnulus)=wallDistance(ncSinglePassage)---（11）
Wherein WALLDISTANCE (NCSINGLEPASSAGE) is the wall distance of the space grid point No. NCSINGLEPASSAGE 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 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 application. 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 exemplary embodiments according to 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 the 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 on … …," "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 "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 (6)

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;
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;
Wherein: the mapping relation between the pre-established single-channel inner space grid point and the full-ring space grid point is obtained by the following steps:
According to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
performing 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; the rotating copying of the space grid points in the single channel to generate the full-ring space grid points of the single rotor comprises the following steps:
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.
2. 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.
3. The method according to claim 2, wherein the method further comprises:
the number of each single-channel space grid in the full-ring space grid points is established.
4. 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;
The mapping module is used for 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;
wherein: the mapping module is used for:
According to a single rotor to be calculated, determining space grid points in a single channel on the single rotor;
Performing 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;
The mapping module is used for:
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.
5. 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-3.
6. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program which, when executed, implements the method for calculating the full-loop mesh wall surface distance according to any one of claims 1 to 3.
CN202410023043.5A 2024-01-08 2024-01-08 Calculation method and device for wall surface distance of full-ring grid Active CN117540507B (en)

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 CN117540507A (en) 2024-02-09
CN117540507B true 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)

* Cited by examiner, † Cited by third party
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

Patent Citations (10)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
Numerical Simulation and Evaluation of Velocity Fluctuations During Rotating Stall of a Centrifugal Pump;Andreas Lucius 等;Journal of Fluids Engineering;20110816;第133卷(第08期);第081102页 *
基于TBR模型的高水头混流式水轮机水力性能预测;孙龙刚 等;农业工程学报;20190408;第35卷(第07期);第62-69页 *
基于流热固耦合燃气轮机压气机叶轮强度分析;李广新 等;流体机械;20170430;第45卷(第04期);第28-32页 *
某膨胀循环发动机推力室冷却结构流场仿真分析;周伟;火箭推进;20150415;第41卷(第02期);第63-69页 *
湍流模拟壁面距离MPI/OpenMP混合并行计算方法;赵钟 等;空气动力学学报;20191215;第37卷(第06期);第883-892页 *
笛卡尔网格下精确高效的壁面距离计算方法;孟爽 等;空气动力学学报;20230428;第41卷(第07期);第93-101页 *
转子叶片气弹稳定性与强迫响应分析;韩乐 等;航空发动机;20210729;第47卷(第04期);第82-90页 *
转静干涉对转子叶片颤振特性的影响;杨慧 等;北京航空航天大学学报;20150401;第42卷(第02期);第258-264页 *
高水头混流式水轮机飞逸工况流动特性数值研究;麻全;中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑;20180215(第02期);第C037-210页 *

Also Published As

Publication number Publication date
CN117540507A (en) 2024-02-09

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
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
CN101403907A (en) Five-coordinate NC tool track fast generation method for scattered point cloud data
CN105184683A (en) Probability clustering method based on wind electric field operation data
Scharfenstein et al. Probabilistic CFD analysis of high pressure turbine blades considering real geometric effects
Frey et al. The discrete adjoint of a turbomachinery RANS solver
CN115774903A (en) Helicopter rotor disc load real-time generation method, device, equipment and medium
Rezvaya et al. Solving the hydrodynamical tasks using CFD programs
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
Pillai et al. Modelling wind turbine wakes at Middelgrunden wind farm
Ren et al. Investigation of compressor tip clearance flow based on the discontinuous Galerkin methods
Kaplan et al. Reverse engineering design of a hydraulic turbine runner
CN114841091A (en) Pneumatic optimization design method of centrifugal ventilator based on proxy model technology
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
CN111611685A (en) Actuating line method for simulating working flow field of axial flow exhaust fan of underground workshop of pumped storage power station
Rodriguez et al. A new inlet distortion and pressure loss based design of an intake system for stationary gas turbines
US7953563B2 (en) Determining effects of turbine blades on fluid motion
Chen et al. Scalable domain decomposition algorithms for simulation of flows passing full size wind turbine

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