CN111626005B

CN111626005B - Flow-solid-thermal coupling simulation method and device based on virtual hole grid self-identification technology

Info

Publication number: CN111626005B
Application number: CN202010461830.XA
Authority: CN
Inventors: 高浩卜; 陈曦; 王珂
Original assignee: China Aero Engine Research Institute
Current assignee: China Aero Engine Research Institute
Priority date: 2020-05-27
Filing date: 2020-05-27
Publication date: 2022-12-06
Anticipated expiration: 2040-05-27
Also published as: CN111626005A

Abstract

In order to solve the problem of low simulation calculation precision in the prior art, the invention provides a fluid-solid-heat coupling simulation method and device based on a virtual hole grid self-identification technology, and the simulation precision of a source item method and the automation degree of engineering design application are improved. The method comprises the following steps: generating virtual cooling hole information of a non-porous flame tube wall mesh model according to the hole arrangement requirement of the cooling holes; determining an add source item grid based on the virtual cooling hole information; completing a fluid-solid-thermal coupling simulation based on the grid of added source items; the invention also discloses a corresponding device, a computer readable storage medium and computer equipment, which improve the accuracy of determining the added source item grid and improve the precision of the fluid-solid-thermal coupling simulation.

Description

Flow-solid-thermal coupling simulation method and device based on virtual hole grid self-identification technology

Technical Field

The disclosure relates to the field of simulation of flow combustion of an aircraft engine combustion chamber, in particular to a flow-solid-heat coupling simulation method and device based on a virtual hole grid self-identification technology.

Background

In an aircraft engine combustion system, the wall surface of a flame tube is influenced by high-temperature flame and gas in a combustion chamber for a long time, and is easy to damage. Divergent cooling is a mainstream method of cooling the wall surface of the flame tube. The cooling mode is dispersed and a large amount of apertures are distributed on the wall surface of the whole flame tube wall in a certain angle full-covering manner, on the one hand, wall cooling is carried out through a cooling airflow passing through a hole channel, on the other hand, a cooling airflow layer is formed near the inner wall of the flame tube, and certain isolation effect is played on the wall surface of a combustion chamber and high-temperature gas inside the combustion chamber.

In the simulation calculation of the combustion chamber of the aero-engine, the complexity of a calculation grid model of the combustion chamber is greatly increased due to a large number of small hole structures on the wall surface of the flame tube, the number of grids is greatly increased, and the calculation time cost is increased. Therefore, experts and scholars at home and abroad propose a method for adding source items such as mass, momentum, energy and the like to a non-porous flame tube wall grid model, and virtual holes are constructed on a wall surface grid to simulate the flow effect of an inlet and an outlet of a cooling hole and the cooling effect of airflow flowing through a hole channel. The existing source item method mainly comprises the following steps of determining a grid for adding source items: a direct calling method, wherein when a non-porous grid model is constructed, grids of an inlet surface and an outlet surface of a virtual hole are independently established, and the grid surface of the inlet or the outlet is directly called for adding source items; the geometric position relation exhaustion method exhausts all possible conditions of the geometric position relation between the grids and the hole-shaped lines, and on a non-hole grid model, the position relation between the grids and the hole-shaped lines is judged by traversing each grid, so that whether the grid is the grid needing to be added with the source items is determined.

Although the direct calling method can completely express the geometric information of the cooling hole shape line, because the area of the inlet and outlet surfaces of the built holes is smaller, a grid with smaller scale is required to be constructed, and the number of the grids is not obviously reduced; although the geometric position relation exhaustion method can realize a great reduction in the number of meshes, the geometric position relation exhaustion method may cause the omission of the meshes for adding source items to easily occur in the judgment of the geometric position relation between the polygonal surface meshes and the complex hole-shaped lines, which causes a problem of low simulation calculation accuracy.

Disclosure of Invention

In order to solve the problem of low simulation calculation precision in the prior art, the invention provides a fluid-solid-thermal coupling simulation method and device based on a virtual hole grid self-identification technology, and the simulation precision is improved.

In a first aspect of the disclosure, a fluid-solid-thermal coupling simulation method based on a virtual hole grid self-identification technology includes: generating virtual cooling hole information of a flame tube wall mesh model with a non-porous structure according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement; determining an addition source item mesh based on the virtual cooling hole information, the addition source item mesh comprising an entrance face mesh of virtual cooling holes, an exit face mesh of virtual cooling holes, and a solid domain heat sink region mesh; completing a fluid solid thermal coupling simulation based on the grid of added source items.

Optionally, generating virtual cooling hole information of the non-porous flame tube wall mesh model according to the hole arrangement requirement of the cooling holes includes: acquiring hole arrangement requirements and hole basic information of cooling holes, wherein the hole basic information comprises: hole center position coordinate information, hole direction vector information and hole cross section shape information of a first hole in each row of holes; if the hole arrangement requirements comprise the requirements of avoiding the main combustion holes and the mixing holes, generating hole center coordinate information and expansion line point cloud coordinate information of each cooling hole according to the hole basic information and the hole arrangement requirements of the cooling holes;

judging whether the cooling holes meet a first preset condition and a second preset condition, and generating the virtual cooling hole information based on hole basic information of the cooling holes meeting the first preset condition and the second preset condition;

wherein, judging whether the cooling hole meets the first preset condition and the second preset condition comprises: if the hole center point and the point cloud of the expansion line of the cooling hole are both in the solid domain grid of the flame tube wall, judging that the cooling hole meets a first preset condition, otherwise, judging that the cooling hole does not meet the first preset condition; when the side wall surface mesh of each main combustion hole or each blending hole is orthographically projected to the point cloud plane of the cooling hole, if a side wall surface mesh projection graph exists on the point cloud plane of the cooling hole, the side wall surface mesh projection graph is inside an expansion point cloud line of the cooling hole or is intersected with the expansion point cloud line of the cooling hole, and the projection length of the side wall surface mesh projection does not exceed half of the hole depth of the cooling hole, the cooling hole is judged not to meet a second preset condition, and otherwise, the cooling hole is judged to meet the second preset condition; and the side wall surface grid projection graph is a graph formed by orthographic projection of the side wall surface grid of the main combustion hole or the mixing hole to a point cloud plane of the cooling hole.

Optionally, the generating the virtual cooling hole information based on the hole basic information of the cooling hole that meets the first preset condition and meets the second preset condition includes: if the hole arrangement requirement comprises a requirement that the intersecting holes and the adjacent holes are replaced by the vertical holes, generating the vertical holes based on the cooling holes which meet the first preset condition and do not meet the second preset condition; generating a vertical hole center point and an expanded line point cloud according to the information of the vertical hole, judging whether the vertical hole meets a first preset condition and a second preset condition, if so, replacing the original cooling hole with the vertical hole, otherwise, deleting the cooling hole, not replacing the vertical hole, and generating hole information of all the cooling holes based on all the original cooling holes and the vertical holes which meet the first preset condition and the second preset condition.

Optionally, determining to add the source item mesh based on the virtual cooling hole information includes:

orthographically projecting each entrance wall surface grid to a point cloud plane of a virtual cooling hole, and if an entrance wall surface grid projection graph is inside a point cloud line of the virtual cooling hole or is intersected with the point cloud line of the virtual cooling hole, and the projection length of the entrance wall surface grid projection is not more than half of the depth of the virtual cooling hole, determining the entrance wall surface grid as the entrance surface grid of the virtual cooling hole; orthographically projecting each outlet wall surface grid to a point cloud plane of a virtual cooling hole, and if the outlet wall surface grid projection graph is inside or intersected with a point cloud line of the virtual cooling hole and the projection length of the outlet wall surface grid projection does not exceed half of the hole depth of the virtual cooling hole, determining the outlet wall surface grid as the outlet surface grid of the virtual cooling hole; and judging each grid in the solid area of the wall of the flame tube, and if the central point of the grid is positioned in the virtual hole, determining the grid as the solid area heat sink area grid of the virtual cooling hole.

Optionally, calculating the effective flow area of the virtual cooling hole based on the added source item grid;

calculating the effective flow area of the virtual cooling hole based on the grid of added source items includes: projecting the added source item mesh onto a point cloud plane of a corresponding virtual cooling hole to form an added source item mesh projection graph, sequentially connecting mesh internal vertexes of the added source item mesh projection graph, a set intersection point and a shape line point cloud point to form an overlapped polygon, and calculating an effective flow area according to the overlapped polygon, wherein the set intersection point is an intersection point of a mesh edge of the added source item mesh projection graph and a point cloud shape line; completing a fluid solid thermal coupling simulation based on the grid of added source items comprises: performing a fluid-solid-thermal coupling simulation based on adding a source term grid and the effective flow area.

Optionally, the method further includes: determining an inlet data acquisition region and an outlet data acquisition region of a virtual cooling hole based on the virtual cooling hole information;

acquiring flow-thermodynamic data for an inlet data acquisition region and an outlet data acquisition region, the flow-thermodynamic data including static pressure, velocity, and temperature;

judging the gas flow direction of the virtual cooling hole according to flow-thermodynamic data, and calculating and obtaining the source item surface density of a fluid domain grid control equation, wherein the control equation comprises at least one of a continuity equation, a momentum equation and an energy equation; completing the flow-solid-thermal coupling simulation based on adding the source item grid and the effective flow area comprises the following steps: calculating a first source term of a fluid domain adding source term grid based on the gas flow direction of the virtual cooling holes, the source term areal density and the grid effective flow area, adding the first source term to the fluid domain grid control equation, calculating a second source term of a solid domain heat sink area grid energy equation according to the virtual hole gas hole flow and the temperature, and adding the second source term to the energy equation of the solid domain heat sink area grid.

Optionally, calculating an addition source term of the fluid domain addition source term grid based on the virtual cooling hole gas flow direction, the source term areal density, and the grid effective flow area includes: calculating the size of a third source item of a current calculation step of adding a source item grid according to the gas flow direction of the virtual cooling hole, the area density of the source item and the effective circulation area of the grid, acquiring the size of a fourth source item of a previous calculation step, and calculating the first source item through a relaxation factor method according to the size of the third source item and the size of the fourth source item; calculating the source term of the solid domain heat sink area grid energy equation according to the flow and the temperature of the virtual holes and the air holes comprises the following steps: and calculating the size of a fifth source item of the current calculation step of the solid region heat sink region grid according to the flow and the temperature of the virtual holes and the air holes, acquiring the size of a sixth source item of the previous calculation step, and calculating the second source item by a relaxation factor method according to the size of the fifth source item and the size of the sixth source item.

In a second aspect of the present disclosure, a fluid-solid-thermal coupling simulation apparatus based on a virtual hole grid self-identification technology includes:

the virtual cooling hole information confirmation module is used for generating virtual cooling hole information of the flame cylinder wall mesh model with the non-porous structure according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement;

an adding source item grid determining module, configured to determine an adding source item grid based on the virtual cooling hole information, where the adding source item grid includes an inlet face grid of a virtual cooling hole, an outlet face grid of a virtual cooling hole, and a solid region heat sink region grid;

a computation module to complete a fluid-solid-thermal coupling simulation based on the grid of added source items.

In a third aspect of the disclosure, a computer readable storage medium, having stored thereon a computer program which, when executed by a processor, carries out the steps of the method of any one of the first aspects of the disclosure.

In a fourth aspect of the present disclosure, a computer device comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the method of any one of the first aspect of the present disclosure when executing the computer program.

The invention is improved on the basis of the existing source item method fluid-solid-thermal coupling simulation method. According to the technical scheme of the embodiment of the disclosure, hole center coordinate information, hole direction vector information and hole cross section shape information in virtual cooling hole information are used as bases, a unified and general calculation mode is adopted for various geometric position conditions, and the added source item grids and the effective flow area of the added source item grids are determined. In addition, at the cooling hole in-process of arranging, to some design principles that aeroengine combustion chamber flame tube cooling hole was arranged, compare the mode that current simple array was arranged and was added manual adjustment, realized automatic identification and the function of arranging, improve the degree of automation that engineering design used.

Drawings

The accompanying drawings, which are included to provide a further understanding of the disclosure and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the disclosure and together with the description serve to explain the principles of the disclosure.

Fig. 1 is a flow chart of a fluid-solid-thermal coupling simulation method based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure;

FIG. 2 is a schematic view of a cooling hole of an embodiment of the present disclosure;

FIG. 3 is a schematic view of an expanded line point cloud calculation of a cooling hole of an embodiment of the present disclosure;

FIG. 4 is a schematic illustration of a rotational arrangement of cooling holes according to an embodiment of the present disclosure;

FIG. 5 is a schematic drawing arrangement of cooling holes of an embodiment of the present disclosure;

FIG. 6 is a flow chart of virtual cooling hole information generation in accordance with an embodiment of the present disclosure;

FIG. 7 is a side wall surface mesh projection schematic view of a primary fire hole of an embodiment of the present disclosure;

FIG. 8 is a schematic diagram of a grid projection computation of an embodiment of the present disclosure;

FIG. 9 is a schematic diagram of one effect of automatically identifying and deleting original cooling holes according to an embodiment of the disclosure;

FIG. 10 is another virtual cooling hole information generation flow diagram of an embodiment of the present disclosure;

FIG. 11 is a schematic diagram of one effect of a vertical hole in place of an original cooling hole according to embodiments of the present disclosure;

FIG. 12 is a schematic diagram of an effective flow area calculation for adding a grid of source items according to an embodiment of the present disclosure;

FIG. 13 is a schematic representation of the effect of characterizing pore cross-sectional area of an embodiment of the present disclosure compared to the prior art;

FIG. 14 is a schematic view of a data acquisition region and a solid-domain heat sink region of an embodiment of the present disclosure;

FIG. 15 is a block diagram of a fluid solid heat coupling simulation device based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure;

FIG. 16 is another flow diagram of a fluid solid heat coupling simulation method based on a virtual hole grid self-identification technique in accordance with an embodiment of the present disclosure;

FIG. 17 is a diagram of one substep of a fluid solid thermal coupling simulation method based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure;

FIG. 18 is a diagram of another substep of a fluid solid thermal coupling simulation method based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure;

FIG. 19 is a diagram of another substep of a fluid solid thermal coupling simulation method based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure;

fig. 20 is a diagram of another substep of a fluid solid thermal coupling simulation method based on a virtual hole grid self-identification technique according to an embodiment of the present disclosure.

Detailed Description

The present disclosure will be described in further detail with reference to the drawings and embodiments. It is to be understood that the specific embodiments described herein are for purposes of illustration only and are not to be construed as limitations of the present disclosure. It should be further noted that, for the convenience of description, only the portions relevant to the present disclosure are shown in the drawings.

It should be noted that the embodiments and features of the embodiments in the present disclosure may be combined with each other without conflict. The present disclosure will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.

Here, basic information of the wells, well information, and well arrangement appearing in the embodiments disclosed in the present application are explained.

Basic information of the well: well information for the first well in each row of wells.

Pore information: hole center position coordinates, hole direction vectors, and hole cross-section shape information.

Hole arrangement: in the first case: a simple rotating array or a stretched array. In the second case: after the array, the main combustion holes and the mixing holes are automatically identified and avoided, and the interference or similar holes are automatically removed. In a third case: in general, the cooling holes are all inclined holes, when the inclined holes interfere or are close to the main combustion holes or the mixing holes, the vertical holes are used for replacing the inclined holes, whether the vertical holes interfere or are close to each other is judged, if the vertical holes do not interfere or are close to each other, the vertical holes are used for automatically replacing the inclined holes, and the inclined holes are not directly removed as in the second case.

According to the basic information of the holes, the holes are arranged according to the hole arrangement, and the hole information of all the holes can be generated.

Example 1:

referring to fig. 1 to fig. 3, the embodiment discloses a fluid-solid-thermal coupling simulation method based on a virtual pore grid self-identification technology, which includes:

y1: generating virtual cooling hole information of a flame tube wall mesh model with a non-porous structure according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement;

y2: determining an addition source item mesh based on the virtual cooling hole information, the addition source item mesh comprising an inlet face mesh of virtual cooling holes, an outlet face mesh of virtual cooling holes, and a solid domain heat sink region mesh;

y3: completing a fluid solid thermal coupling simulation based on the grid of added source items.

The invention is improved on the basis of the existing source term method fluid-solid-thermal coupling simulation (calculation) method. According to the technical scheme of the embodiment of the disclosure, hole center coordinate information, hole direction vector information and hole cross section shape information in virtual cooling hole information are taken as references, a unified and general computing mode is adopted for various geometric position conditions, and the added source item grids are determined. In addition, at the cooling hole in-process of arranging, to some design principles that aeroengine combustion chamber flame tube cooling hole was arranged, compare the mode that current simple array was arranged and was added manual adjustment, realized automatic identification and the function of arranging, improve the degree of automation that engineering design used. Referring to fig. 2 and 3, the hole center coordinate information contains coordinates of the hole center. The hole direction vector information includes hole direction vectors of the cooling holes, and in the present embodiment, a vertical wall surface direction vector and an inclined hole deflection angle are used as the hole direction vectors, where the inclined hole deflection angle is a three-dimensional space included angle and is composed of three-dimensional space included angles. The hole cross section shape information and the hole center coordinate information can be used for determining the shape line point cloud coordinates of the cooling hole; if the cross section of the hole is a circular hole, the diameter can be directly given by the information of the cross section of the hole; if the hole cross-sectional shape is a hole of other shape, the hole cross-sectional shape information may be a characteristic size or a function curve that enables determination of the shape and size. Wherein the shape line point cloud coordinates are coordinates of a shape line point cloud. The contour point cloud refers to a set of points forming the cross section shape of the hole, and the plane of the contour point cloud passes through the hole center of the cooling hole and is vertical to the direction of the cooling hole; the point cloud-shaped line is a line formed by connecting adjacent points in the point cloud of the shape line.

According to hole center coordinate information, hole direction vector information, hole cross section shape information and flame tube wall thickness in the virtual cooling hole information, the shape line point cloud coordinate and the hole depth of the virtual cooling hole can be determined; according to the shape line point cloud coordinates and the hole depths of the cooling holes, an inlet face grid of the virtual cooling holes, an outlet face grid of the virtual cooling holes and a solid region heat sink area grid can be determined.

The fluid-solid-thermal coupling simulation can simulate a temperature field, a speed field and the like, and can be used for observing whether the wall surface of the flame tube is cooled by the air holes to cause temperature reduction, the effect of the speed field of air flow inflow and outflow at the wall surface and the like.

In one embodiment, generating virtual cooling hole information of the flame tube wall mesh model of the non-porous structure according to the hole arrangement requirement of the cooling holes comprises the following steps: acquiring hole arrangement requirements and hole basic information of the cooling holes;

referring to fig. 4 and 5, the hole arrangement of the cooling holes includes: the maximum number of holes per row and the stretching or rotating arrangement. The hole basic information of the cooling holes includes hole center position coordinate information, hole direction vector information, and hole cross-sectional shape information of the first hole in each row of holes. As shown in fig. 4, the arrangement of the rotation includes an included angle between adjacent holes, a rotation axis, and a rotation direction; as shown in fig. 5, the arrangement of the stretching includes a stretching direction and a distance between adjacent holes.

Determining hole center coordinate information, hole direction vector information and hole cross section shape information of all holes according to hole center position coordinate information, hole direction vector information, hole cross section shape information of a first hole in each row of holes and the stretching or rotating arrangement mode of cooling holes; when the system uses the method, only a small amount of information needs to be input, the added source item grids can be determined, and flow-solid-heat coupling simulation is completed through a source item method based on the added source item grids.

In one embodiment, referring to fig. 6 to 9, generating virtual cooling hole information of a flame tube wall mesh model with a non-porous structure according to the hole arrangement requirement of the cooling holes comprises:

y11, acquiring hole arrangement requirements and hole basic information of the cooling holes, wherein the hole basic information comprises center position coordinate information, hole direction vector information and hole cross section shape information of a first hole in each row of holes;

y12, if the hole arrangement requirement comprises the requirement of avoiding the main combustion holes and the mixing holes, generating hole center coordinate information and expansion shape line point cloud coordinate information of each cooling hole according to the hole basic information and the hole arrangement requirement of the cooling holes;

y13, judging whether the cooling holes meet a first preset condition and a second preset condition, and generating the virtual cooling hole information based on hole basic information of the cooling holes meeting the first preset condition and the second preset condition;

wherein, judging whether the cooling hole meets the first preset condition and the second preset condition comprises:

if the hole center point and the point cloud of the expansion line of the cooling hole are both in the solid domain grid of the flame tube wall, judging that the cooling hole meets a first preset condition, and otherwise, judging that the cooling hole does not meet the first preset condition;

when the side wall surface mesh of each main combustion hole or each blending hole is orthographically projected to the point cloud plane of the cooling hole, if a side wall surface mesh projection graph exists on the point cloud plane of the cooling hole, the side wall surface mesh projection graph is inside an expansion point cloud line of the cooling hole or is intersected with the expansion point cloud line of the cooling hole, and the projection length of the side wall surface mesh projection does not exceed half of the hole depth of the cooling hole, the cooling hole is judged not to meet a second preset condition, and otherwise, the cooling hole is judged to meet the second preset condition; and the side wall surface grid projection graph is a graph formed by orthographic projection of the side wall surface grid of the main combustion hole or the mixing hole to a point cloud plane of the cooling hole.

Referring to fig. 3, the expansion point cloud with lines is formed by expanding a point cloud with lines outward by a set expansion radius increment with the center of a hole as a dot, and the expansion point cloud with lines is formed by expanding a point cloud with lines outward by a set expansion radius increment; wherein the increment of the radius of expansion is obtained according to the proximity distance required for avoiding the main combustion hole or the mixing hole.

Wherein, referring to fig. 6 and 7, the flame tube wall solid domain mesh is a mesh of flame tube wall solid domains; the point cloud plane is a plane where the shape line point cloud is located; the side wall surface mesh projection graph is a graph formed when the side wall surface mesh of the main combustion hole or the mixing hole is orthographically projected to the point cloud plane of the cooling hole; the projection length of the side wall surface mesh projection is the vertical distance between the side wall surface mesh of the main burning hole or the mixing hole and the point cloud plane. It can be appreciated that if a cooling hole has a hole center point or some expanded shape line point cloud point not within the flame tube wall solid domain mesh, then the cooling hole does not satisfy the first predetermined condition, i.e., the hole is an internal hole.

According to the technical scheme, when the hole arrangement requirement includes the requirement of avoiding the main combustion holes and the mixing holes, the required cooling holes are automatically and accurately determined to be not the inner holes, the intersected holes and the adjacent holes according to whether the cooling holes meet the first preset condition and the second preset condition or not (namely, the cooling holes meeting the first preset condition and the second preset condition are not the inner holes, and the cooling holes meeting the second preset condition are not the intersected holes and the adjacent holes), the hole information of other cooling holes not belonging to the inner holes, the intersected holes and the adjacent holes generates the final hole information of the virtual cooling holes, and the avoiding function of the main combustion holes and the mixing holes is automatically realized. Referring to fig. 9, an effect of automatically identifying and deleting original cooling holes according to the embodiment of the present disclosure is schematically illustrated, wherein the "identified and deleted original cooling holes" in fig. 9 are internal holes, intersecting holes, or adjacent holes of the "main combustion holes or blending holes" in fig. 9.

According to the technical scheme of the embodiment of the disclosure, the arrangement positions of the cooling holes on the wall surface are judged by using a geometric algorithm, main combustion holes and mixing holes in any shapes are automatically avoided within a specified distance, and the efficiency is improved.

In one embodiment, referring to fig. 10 and 11, generating the virtual cooling hole information based on hole basic information of a cooling hole satisfying a first preset condition and satisfying a second preset condition includes:

y131: if the hole arrangement requirement comprises a requirement that the intersecting holes and the adjacent holes are replaced by the vertical holes, generating the vertical holes based on the cooling holes which meet the first preset condition and do not meet the second preset condition;

y132: generating a vertical hole center point and an expanded line point cloud according to the information of the vertical hole, judging whether the vertical hole meets a first preset condition and a second preset condition, if so, replacing the original cooling hole with the vertical hole, otherwise, deleting the cooling hole, not replacing the vertical hole, and generating hole information of all the cooling holes based on all the original cooling holes and the vertical holes which meet the first preset condition and the second preset condition.

The hole information of the vertical hole includes hole center position coordinate information, hole direction vector information, and hole cross-sectional shape information, which are vertical holes. The hole information of the cooling hole includes hole center position coordinate information, hole direction vector information, and hole cross-sectional shape information of the cooling hole.

According to the technical scheme of the embodiment of the disclosure, when the hole arrangement requirement comprises a requirement that a vertical hole replaces an intersected hole and an adjacent hole, the required cooling holes belong to the intersected hole and the adjacent hole (namely, the cooling holes meeting the first preset condition and not meeting the second preset condition) are automatically, quickly and accurately determined according to whether the cooling holes meet the first preset condition and the second preset condition, the hole information of the vertical hole replaces the hole information of the intersected hole and the hole information of the adjacent hole, the vertical hole is judged, and if the vertical hole meets the first preset condition and meets the second preset condition, the vertical hole information is reserved so as to achieve the purpose that the vertical hole replaces the intersected hole and the adjacent hole. Referring to fig. 11, a schematic view of the effect of replacing the original cooling holes with vertical holes is shown, wherein the "identified and deleted original cooling holes" are internal holes corresponding to the "primary combustion holes or blend holes" in the figure. The "vertical hole instead of the original cooling hole" in fig. 11 is an intersecting hole or an adjacent hole to the "main combustion hole or the blend hole" in fig. 11.

According to the embodiment of the disclosure, when the cooling inclined holes and the main combustion holes are geometrically interfered, the function of replacing the original inclined holes with the vertical holes perpendicular to the wall surface can be automatically realized, and the arrangement principle and the rule of the divergent cooling holes in engineering application are comprehensively considered.

In one embodiment, determining to add a source item mesh based on the virtual cooling hole information comprises:

orthographic projection is carried out on each inlet wall surface grid to the point cloud plane of the virtual cooling hole; if the projection graph of the inlet wall surface mesh is inside the point cloud line of the virtual cooling hole or is intersected with the point cloud line of the virtual cooling hole, and the projection length of the inlet wall surface mesh projection is not more than half of the depth of the virtual cooling hole, determining the inlet wall surface mesh as the inlet surface mesh of the virtual cooling hole;

and (3) orthographically projecting each outlet wall surface grid to a point cloud plane of the virtual cooling hole, and if the outlet wall surface grid projection graph is inside or intersected with a point cloud line of the virtual cooling hole and the projection length of the outlet wall surface grid projection is not more than half of the hole depth of the virtual cooling hole, determining the outlet wall surface grid as the outlet surface grid of the virtual cooling hole.

The inlet wall surface grid is a wall surface grid positioned on the inlet side of the cooling hole in the flame cylinder wall grid, and the outlet wall surface grid is a wall surface grid positioned on the outlet side of the cooling hole in the flame cylinder wall grid; the projection length of the entrance wall surface mesh projection is the maximum value of the distance between each vertex of the entrance wall surface mesh and the point cloud plane; the projection length of the exit wall surface mesh projection is the maximum value of the distance between each vertex of the exit wall surface mesh and the point cloud plane.

And judging each grid in the solid area of the wall of the flame tube, and if the central point of the grid is positioned in the virtual hole, determining the grid as the solid area heat sink area grid of the virtual cooling hole. The virtual hole interior can be a column region formed by projecting a point cloud line of the hole along a hole direction vector and the opposite direction of the hole direction vector, wherein the projection distance does not exceed half of the hole depth.

In the technical scheme of the embodiment, whether the inlet wall surface mesh is the inlet surface mesh of the virtual cooling hole is determined according to the condition that the projection length of the inlet wall surface mesh projection graph is not more than half of the hole depth of the virtual cooling hole in the inside of the point cloud-shaped line of the virtual cooling hole or is intersected with the point cloud-shaped line of the virtual cooling hole; determining whether the outlet wall surface mesh is the outlet surface mesh of the virtual cooling hole according to the condition that the projection graph of the outlet wall surface mesh is inside or intersected with the point cloud line of the virtual cooling hole, and the projection length of the outlet wall surface mesh projection is not more than half of the depth of the virtual cooling hole; and determining the solid region heat sink region grid according to the condition that the grid central point of the flame tube wall solid region grid is positioned in the virtual hole. Such that the computing device determines the add source item grid according to the method quickly and accurately.

In one embodiment, the method further comprises: calculating the effective flow area of the virtual cooling hole based on the added source item grid;

referring to fig. 12 and 13, calculating the effective flow area of the virtual cooling holes based on adding the grid of source terms includes: and projecting the added source item mesh onto a point cloud plane of a corresponding virtual cooling hole to form an added source item mesh projection graph, sequentially connecting mesh internal vertexes of the added source item mesh projection graph, setting intersection points and point cloud points of shape lines to form an overlapped polygon, and calculating an effective flow area according to the overlapped polygon, wherein the set intersection points are the intersection points of mesh edges of the added source item mesh projection graph and the point cloud shape lines.

Completing a fluid solid thermal coupling simulation based on the grid of added source items comprises: and completing the fluid-solid-thermal coupling simulation based on the addition of the source item grid and the effective flow area.

Determining the effective flow area of the cooling hole according to the area of an overlapped polygon formed by sequentially connecting the internal vertexes of the grids of the source item added grid projection graph and the set intersection points and the shape line point cloud points; compared with the prior art that the actual inlet and outlet areas of the cooling holes are replaced by the inlet surface grid area and the outlet surface grid area in the added source item grid, the technical scheme of the disclosure adopts the effective flow area, improves the calculation precision, and adopts the effective flow area calculation method in the embodiment, so that the calculation efficiency is higher, and the computer operation is convenient. Here, it can be known that, as shown in fig. 12, the mesh internal vertex is a vertex falling inside the point cloud shaped line in the mesh projection graph vertex. FIG. 13 is a schematic diagram comparing the effect of characterizing the cross-sectional area of holes of an embodiment of the present disclosure with the prior art, wherein the shaded portion on the left is the cross-sectional area of holes characterized by the area of the grid in the prior art, the shaded portion on the right is the cross-sectional area of holes characterized by the effective flow area in the present example, and the cross-sectional lines of cooling holes in the diagram characterize the actual cross-section of the holes; from the drawing, it can be intuitively obtained that, compared with the prior art in which the cross-sectional area of the hole is represented by the area of the grid, the cross-sectional area of the hole is represented by the effective flow area in the embodiment, and the precision of representing the cross-sectional area of the hole is higher.

In the embodiment, the effective flow area is used for replacing the grid area, so that the sum of the effective flow areas is always equal to the actual cross-sectional area of the hole no matter the size of the grid, and the error of calculating the flow of the virtual hole caused by different grid scales is avoided.

In one embodiment, referring to fig. 14, the method further comprises: determining an inlet data collection area and an outlet data collection area of a virtual cooling hole based on the virtual cooling hole information;

and judging the gas flow direction of the virtual cooling hole according to the flow-thermodynamic data, and calculating and obtaining the source item areal density of a control equation, wherein the fluid domain grid control equation comprises at least one of a continuity equation, a momentum equation and an energy equation. Calculating a source term of a grid energy equation of the solid domain heat sink area according to the flow, the temperature and the wall temperature (acquired by the solid domain heat sink area) of the virtual hole air hole;

completing the fluid-solid-thermal coupling simulation based on adding the source item grid and the effective flow area comprises the following steps: calculating a first source term of a fluid domain addition source term grid based on the virtual cooling hole gas flow direction, the source term areal density and the grid effective flow area, and adding the first source term to the fluid domain grid control equation. And calculating a second source item of the solid area heat sink area grid energy equation according to the flow, the temperature and the wall temperature of the virtual holes and the air holes, and adding the second source item into the solid area heat sink area grid energy equation.

The method also comprises the step of determining a solid area heat sink area of the virtual cooling hole according to the virtual cooling hole information, the wall thickness and the like, wherein the solid area heat sink area is used for simulating the effect that the airflow flowing through the air hole brings away the heat of the wall plate.

Optionally, the fluid domain grid control equation includes a continuity equation, a momentum equation, and an energy equation.

Optionally, calculating the first source term of the fluid domain addition source term grid based on the virtual cooling hole gas flow direction, the source term areal density and the grid effective flow area includes: and acquiring the size of a third source item of the previous calculation step, calculating the size of a fourth source item of the current calculation step of adding the source item grid according to the gas flow direction of the virtual cooling hole, the area density of the source item and the effective flow area of the grid, and calculating the first source item by a relaxation factor method according to the size of the third source item and the size of the fourth source item. Calculating a second source term of the solid region heat sink region grid energy equation according to the virtual orifice and air hole flow, the temperature and the wall temperature comprises: and calculating the size of a fifth source item of the current calculation step of the solid region heat sink region grid according to the flow, the temperature and the wall temperature of the virtual holes and the air holes, acquiring the size of a sixth source item of the previous calculation step, and calculating a second source item by a relaxation factor method according to the size of the fifth source item and the size of the sixth source item.

According to the technical scheme, when the cooling holes are arranged on the wall surface of the flame tube, the cooling holes can be automatically arranged according to an avoiding principle, and the cooling holes can automatically avoid the main combustion holes and the mixing holes within a certain approaching distance.

According to the technical scheme, when the cooling holes are arranged on the wall surface of the flame tube, the inclined holes and the main combustion holes can be automatically judged according to the 'replacement principle', and if the vertical holes perpendicular to the wall surface meet the condition of no interference or no proximity, the vertical holes can replace the interference holes.

The technical scheme disclosed by the invention can ensure that the flow calculated by the virtual holes is not influenced by the grid dimension when the mass flow source item is added on the grid with any size and dimension. However, the grid size needs to be ensured not to be too large, and a situation that one grid covers a plurality of holes of the entrances and exits occurs. Through tests, the technical scheme disclosed by the invention uses grid models with grid scales of 0.1mm, 0.2mm, 0.5mm and 1mm, and compared with the scheme with a porous structure, the errors are respectively 1.80%, 1.84%, 1.02% and 0.16%.

According to the technical scheme, under the condition that certain calculation precision is guaranteed, the porous structure is replaced by the non-porous flat plate structure, the grid model manufacturing process can be greatly simplified, the number of grids is reduced, and calculation time is greatly saved. In an example, the minimum grid number of the real holes is 388.87 ten thousand, and the grid surface is not closed at the position near the small air holes, so that the grid cannot be drawn. In the technical scheme of the disclosure, the grid quantities of the 0.2mm, 0.5mm and 1mm scale grids are 277.60 ten thousand, 75.66 ten thousand and 52.33 ten thousand respectively. The number of grids of the source item method scheme can be reduced to about 1/5 to 1/7 of the original number.

The technical scheme disclosed by the invention can be suitable for the polygonal convex mesh model with any number of sides and any size. However, the grid size needs to be ensured not to be too large, and a situation that one grid covers a plurality of holes of the entrances and exits occurs. The cooling hole cross-sectional shape may be any convex shape not limited to a circle. Under different grid sizes, the sum of the effective flow areas of the grids in each hole is always equal to the cross-sectional area of the hole.

According to the technical scheme, the combustion chamber fluid-solid-heat coupling simulation calculation with the divergent cooling structure can be realized.

Example 2:

referring to fig. 15, the fluid-solid-thermal coupling simulation apparatus based on the virtual pore grid self-identification technology includes:

the virtual cooling hole information confirmation module 1 is used for generating virtual cooling hole information of a flame cylinder wall mesh model with a non-porous structure according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement;

an adding source item mesh determining module 2, configured to determine an adding source item mesh based on the virtual cooling hole information, where the adding source item mesh includes an inlet face mesh of virtual cooling holes, an outlet face mesh of virtual cooling holes, and a solid-domain heat sink region mesh;

and the computing module 3 is used for completing fluid-solid-thermal coupling simulation computation based on the added source item grid and adding source items of a mark network.

In one embodiment, generating virtual cooling hole information of the flame tube wall mesh model of the non-porous structure according to the hole arrangement requirement of the cooling holes comprises the following steps: acquiring hole arrangement requirements and hole basic information of cooling holes;

wherein the hole arrangement requirement of the cooling hole comprises: the maximum number of holes per row and the stretching or rotating arrangement. The hole basic information of the cooling holes includes hole center position coordinate information, hole direction vector information, and hole cross-sectional shape information of the first hole in each row of holes. The rotary arrangement mode comprises an included angle between adjacent holes, a rotary shaft and a rotary direction; the stretching arrangement mode comprises the stretching direction and the space between adjacent holes.

Determining hole center coordinate information, hole direction vector information and hole cross section shape information of all holes according to hole center position coordinate information, hole direction vector information, hole cross section shape information of a first hole in each row of holes and the stretching or rotating arrangement mode of cooling holes; when the system uses the method, only a small amount of information needs to be acquired, the added source item grid can be determined, and further the flow-solid-heat coupling simulation is completed through a source item method based on the added source item grid.

In one embodiment, generating virtual cooling hole information of the flame tube wall mesh model of the non-porous structure according to the hole arrangement requirement of the cooling holes comprises the following steps:

acquiring hole arrangement requirements and hole basic information of cooling holes, wherein the hole basic information comprises hole center position coordinate information, hole direction vector information and hole cross section shape information of a first hole in each row of holes;

if the hole arrangement requirements comprise the requirements of avoiding the main combustion holes and the mixing holes, generating hole center coordinate information and expansion line point cloud coordinate information of each cooling hole according to the hole basic information and the hole arrangement requirements of the cooling holes;

if the hole center point and the point cloud of the expansion line of the cooling hole are both in the solid domain grid of the flame tube wall, judging that the cooling hole does not meet a first preset condition, and otherwise, judging that the cooling hole meets the first preset condition;

In one embodiment, generating the virtual cooling hole information based on hole basic information of a cooling hole satisfying a first preset condition and satisfying a second preset condition includes:

if the hole arrangement requirement comprises a requirement that the intersecting holes and the adjacent holes are replaced by vertical holes, generating the vertical holes based on the cooling holes which meet the first preset condition and do not meet the second preset condition;

and generating a vertical hole center point and an expanded line point cloud according to the information of the vertical hole, judging whether the vertical hole meets a first preset condition and a second preset condition, if so, replacing the original cooling hole with the vertical hole, otherwise, deleting the cooling hole and not replacing the vertical hole. And generating hole information of all the cooling holes based on all the original cooling holes and the vertical holes which accord with the first preset condition and the second preset condition.

In one embodiment, determining to add a grid of source items based on the virtual cooling hole information comprises:

orthographically projecting each entrance wall surface grid to a point cloud plane of a virtual cooling hole, and if an entrance wall surface grid projection graph is inside a point cloud line of the virtual cooling hole or is intersected with the point cloud line of the virtual cooling hole, and the projection length of the entrance wall surface grid projection is not more than half of the depth of the virtual cooling hole, determining the entrance wall surface grid as the entrance surface grid of the virtual cooling hole;

and (3) orthographically projecting each outlet wall surface grid to the point cloud plane of the virtual cooling hole, and if the projection graph of the outlet wall surface grid is inside or intersected with the point cloud line of the virtual cooling hole and the projection length of the outlet wall surface grid projection is not more than half of the hole depth of the virtual cooling hole, determining the outlet wall surface grid as the outlet surface grid of the virtual cooling hole.

And judging each grid in the solid area of the wall of the flame tube, and if the central point of the grid is positioned in the virtual hole, determining the grid as the solid area heat sink area grid of the virtual cooling hole. And the inside of the virtual hole indicates that the cloud-shaped line of the hole point is projected along the direction vector of the hole and the opposite direction of the hole, and the projection distance does not exceed half of the depth of the hole, so that the inside of a formed columnar area is formed.

In one embodiment, the apparatus further comprises an effective flow area calculation module for calculating an effective flow area of the virtual cooling hole based on adding the grid of source items;

calculating the effective flow area of the virtual cooling hole based on the grid of added source items includes: projecting the added source item mesh onto a point cloud plane of a corresponding virtual cooling hole to form an added source item mesh projection graph, sequentially connecting mesh internal vertexes of the added source item mesh projection graph, a set intersection point and a line point cloud point to form an overlapped polygon, and calculating an effective flow area according to the overlapped polygon, wherein the set intersection point is an intersection point of a mesh edge of the added source item mesh projection graph and a point cloud line;

completing a flow solid thermal coupling simulation based on the adding source item grid comprises: performing a fluid-solid-thermal coupling simulation based on adding a source term grid and the effective flow area.

In one embodiment, the apparatus further comprises an acquisition zone confirmation module,

the acquisition area confirmation module is used for:

determining an inlet data collection area and an outlet data collection area of a virtual cooling hole based on the virtual cooling hole information;

and judging the gas flow direction of the virtual cooling holes according to the flow-thermodynamic data, and calculating and obtaining the source term areal density of a fluid domain grid control equation, wherein the control equation comprises at least one of a continuity equation, a momentum equation and an energy equation. And calculating the source term of the solid region heat sink region grid energy equation according to the flow, the temperature and the wall temperature (collected by the solid region heat sink region) of the virtual holes and the air holes.

Completing the fluid-solid-thermal coupling simulation based on adding the source item grid and the effective flow area comprises the following steps: and calculating an adding source item of the fluid domain adding source item grid based on the gas flow direction of the virtual cooling hole, the source item areal density and the grid effective flow area, and adding the source item to the fluid domain grid control equation. And calculating a source item of an energy equation of the solid area heat sink area grid according to the flow, the temperature and the wall temperature of the virtual holes and the air holes (collected by the solid area heat sink area), and adding the source item into the energy equation of the solid area heat sink area grid.

In one embodiment, calculating the added source term for the fluid domain added source term grid based on the virtual cooling hole gas flow direction, the source term areal density, and the grid effective flow area comprises: and calculating the size of a first source item of the current calculation step of the added source item grid according to the gas flow direction of the virtual cooling hole, the area density of the source item and the effective flow area of the grid, acquiring the size of a second source item of the previous calculation step, and calculating the added source item through a relaxation factor method according to the size of the first source item and the size of the second source item. Calculating the source term of the solid domain heat sink area grid energy equation according to the virtual hole and air hole flow, the temperature and the wall temperature comprises the following steps: calculating the size of a first source item of a current calculation step of the solid region heat sink region grid according to the flow, the temperature and the wall temperature of the virtual holes and the air holes, obtaining the size of a second source item of a previous calculation step, and calculating an added source item through a relaxation factor method according to the size of the first source item and the size of the second source item.

The principle and effect of the apparatus of the present embodiment can be referred to the principle and effect of the method in embodiment 1, and the description of the present embodiment will not be repeated.

Example 3:

a computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any of the embodiment 1.

Example 4:

computer apparatus comprising a memory and a processor, the memory storing a computer program, the processor implementing the steps of the method of embodiment 1 when executing the computer program.

Example 5:

the embodiment relates to a flow-solid-thermal coupling simulation method based on a virtual hole grid self-identification technology, and with reference to fig. 16 to 18, the method comprises the following steps:

s1) automatically arranging cooling holes on the wall surface of the flame tube and transmitting and writing geometric position information. And according to the cooling hole arrangement requirement and the basic information of the cooling holes, calculating related calculation geometric algorithms such as hole geometric information generation, projection judgment position relation and the like, and outputting and writing a cooling hole information file containing all the cooling holes meeting the requirement. As shown in fig. 17, the specific implementation method includes the following steps:

s11) reading hole basic information of the cooling holes and setting hole arrangement requirements of the cooling holes.

The hole basic information of the cooling hole includes:

(1) The center point position coordinates and hole direction vector information of the first hole in each row of holes are shown in fig. 2. The hole direction vector information includes a vertical wall surface direction vector and an inclined hole deflection angle. The inclined hole deflection angle is a three-dimensional space included angle and is composed of three direction included angles.

(2) Pore cross-sectional shape (convex) information. The hole cross-sectional shape information is primarily used to generate the shaping line point cloud coordinates, which is a collection of a series of points that make up the shape of the hole cross-section, as shown in fig. 3. If the hole is a circular hole, the diameter is given; in the case of other shaped holes, it is sufficient to give a characteristic dimension or function curve that can be determined in shape and size.

The hole arrangement requirements of the cooling holes include:

(1) Maximum number of holes per row of holes.

(2) The hole arrangement mode is as follows: stretched or rotated. The stretching arrangement gives the stretching direction and the distance between adjacent holes in the same row of holes (the distance between adjacent holes is different and can be given in an array form). The rotation arrangement gives an equation of the rotation axis, the rotation direction and the included angle between adjacent holes in the same row of holes (the included angle between adjacent holes is different and can be given in an array form). The arrangement is shown in fig. 4 and 5.

(3) And if the hole arrangement needs to avoid the main combustion hole and the mixing hole, the adjacent distance between the main combustion hole and the mixing hole is avoided.

(4) Whether it is necessary to replace the inclined holes interfering with the main combustion hole and the dilution hole with vertical holes.

S12) judging whether hole avoiding is needed. If hole avoiding is needed, executing the step S13; and (5) avoiding the hole, and executing the step S14.

S13) array generation of hole center position and expansion line point cloud coordinate information. According to the hole arrangement requirements of the cooling holes and the basic information of the holes, center points of all the holes are arranged on the wall of the flame tube in an array mode in a stretching or rotating mode, and coordinate data of the center points are stored, as shown in fig. 4 and 5. And calculating the point cloud coordinates of the expansion shape line of each hole according to the set hole avoiding proximity distance, inclined hole deflection angle, hole center point coordinates and hole cross section shape information (circular hole radius or other shape functions), as shown in fig. 3. And the geometric position relation between the grid projection and the expansion holes is judged subsequently, so that the condition that the cooling holes are not interfered with the main combustion holes but are too close to the main combustion holes is eliminated, and the method is suitable for main combustion holes and mixing holes in any shapes.

S14) array generation of hole center position information. According to the hole arrangement requirements of the holes and the basic information of the holes, central points of all the holes are arranged on the wall of the flame tube in an array mode in a stretching or rotating mode, and coordinate data of the central points are stored, as shown in fig. 4 and 5.

S15) judging and marking the inner hole, the intersected hole and the adjacent hole. The marking process is as follows:

(1) And scanning the solid domain grids of the wall of each flame tube, and judging whether the central point of each hole and each point in the expanded point cloud are positioned in the solid grids. If there is a center point or some point cloud point that is not within the solid mesh, then the hole is marked.

(2) And scanning side wall surface grids of each main combustion hole or each blending hole, respectively carrying out orthographic projection on the point cloud plane of each unmarked cooling hole, marking the hole when the projection graph of the side wall surface grids is inside or intersected with the cloud line of the expansion point of a certain cooling hole and the projection length is not more than half of the hole depth, and unmarking when the projection length is more than half of the hole depth, as shown in fig. 7. The hole depth calculation is obtained by calculating the wall thickness and the hole direction vector. The process of determining the geometric position relationship by the surface mesh projection is shown in fig. 8.

S16) judging whether the inclined hole is interfered or close to the inclined hole or not, and replacing the inclined hole by a straight hole. If the replacement is needed, executing the step S17; no replacement is required and step S19 is executed.

S17) replacing the intersecting hole and the adjacent hole with a vertical hole. And (5) aiming at the cooling hole marked in the step (2) of S15, calculating expansion point cloud coordinate information of the vertical hole by using information such as a hole center point, a vertical wall surface direction vector, a hole avoiding approach distance, a hole cross section shape and the like. And simultaneously, deleting expanded point cloud coordinate information of the original inclined hole.

S18) judging intersecting holes and adjacent holes among the vertical holes replaced by the marks. Marking is performed according to the step S15.

And S19) transmitting and writing the coordinate information of the center position of the unmarked hole, the hole direction vector and the hole cross section line information, and storing the information in a file form to finish the transmission and writing of the full-hole information file.

And S2) grid identification and effective flow area calculation. Reading the cooling hole information file generated in the step S1), generating cooling hole geometric information, and performing calculation such as grid point projection, geometric position judgment and the like to realize the addition of the identifier of the source item grid. And constructing an overlapping region of the projection grid and the cross section of the cooling hole by using the geometric information of the cooling hole and the grid projection information, and calculating the area of the overlapping region as the effective flow area of the identification grid. At the same time, data acquisition regions for the virtual well entrance and exit are determined. The specific implementation method is shown in fig. 18, and includes the following steps:

s21) reading in the full-hole information file. Before adding the source item, reading all hole information files generated in S1), including hole center point coordinates, hole direction vectors (the direction is from the hole inlet to the hole outlet), and hole cross section geometric information (such as the diameter of a circular hole).

And S22) generating a contour point cloud coordinate of each hole according to the read file information and the step S13 (the approach distance is 0).

And S23) projection of the vertex of the wall surface mesh marks the entrance mesh. And scanning all the inlet wall surface grids, projecting the top point of each grid to the point cloud plane of each hole, judging whether the projection grids are in the point cloud line or are intersected, and if the projection grids are in the point cloud line or are intersected, and the projection distance is not more than half of the hole depth, marking the grids as the inlet surface grids of the holes. Similarly, the exit face grid of each well is labeled.

S24) calculating the effective circulation area of the mark grid. After the marked mesh is projected on the point cloud-shaped line plane, the internal vertex of the mesh, the intersection point of the projected mesh edge and the point cloud-shaped line, and the point cloud points of the shape line are sequentially identified and connected in sequence to form an overlapped polygon of the projected mesh and the point cloud-shaped line, and the area, namely the effective flow area, is calculated, as shown in fig. 12.

S25) searching for the identification acquisition area and the heat sink area. And calculating a hole inlet data acquisition area, a hole outlet data acquisition area and a solid area heat sink area according to the data such as the hole depth, the hole cross section size and the like. The solid zone heat sink region was used to simulate the effect of the air flow flowing through the air holes carrying away the heat of the wall plate. FIG. 14 shows the program identification effect of the data acquisition zone and the solid heat sink zone.

And S3) acquiring and calculating flow and thermodynamic data of the virtual hole. On the basis of flow field data calculated by a solver and before each calculation step starts, average values of pressure, temperature and the like of each hole data acquisition area are obtained, the flow direction of virtual holes is judged, the source item size of flow, heat exchange quantity and the like of each virtual hole is calculated, and the source item surface density is calculated. The specific implementation method is shown in fig. 19, and includes the following steps:

s31) collecting and processing flow data, thermodynamic data, etc. And collecting static pressure, speed, temperature and other data of the inlet and the outlet of the hole, and obtaining the data of the grid of the collection area by weighted average according to the volume.

S32) judging the flow direction according to the calculation model, and calculating the virtual hole parameters. According to a theoretical flow formula and a flow coefficient (a constant or an empirical formula), judging the flow direction of gas in the virtual hole according to the pressure difference between two sides of the wall plate by data collected by the collecting area, calculating the flow speed, the flow and the density of a flow surface (mass of fluid flowing on a unit cross-sectional area in unit time) and the like of the gas flowing through the virtual hole, calculating the convection heat transfer coefficient, and combining the flow data to obtain the heat transfer rate of the gas flowing through the virtual hole and the fixed wall for calculating the size of an energy equation source item of a solid area heat sink area grid. And calculating to obtain the source item surface density of the grid control equation of each fluid domain at the inlet and the outlet of each hole.

S33) displaying and monitoring the virtual hole to calculate related parameters. And synchronously displaying the average value, the maximum value and the minimum value of the parameters such as the flow coefficient, the pore flow, the flow velocity in the pores and the like of all the pores calculated in each iteration step, and displaying whether the pores have the counter-flow phenomenon or not.

And S4) adding source items to the identification grid by a relaxation factor method. And calculating the size of the added source item of each identification grid of the heat sink area of the fluid area and the solid area respectively according to the flow direction of the virtual holes, the surface density of the source item, the effective circulation area of the grid, the heat exchange rate of airflow flowing through the virtual holes and the solid wall and the like, calculating the size of the source item by combining the grid in the previous step, performing source item calculation of a relaxation factor method, obtaining the size of the added source item in the calculation step, and storing. And adding the calculated source term into a corresponding control equation to complete the calculation of the source term method. The specific implementation method is shown in fig. 20, and includes the following steps:

s41) calculating the size of the source term of the correlation equation. And calculating the source term size of related equations such as a grid continuity equation, a momentum equation, an energy equation and the like of each marking fluid domain according to the flow areal density, the flow velocity and other related acquired data calculated in the step S32). The airflow flow of each grid is the product of the flow area density of the virtual air hole and the effective flow area of the grid. And calculating the heat exchange rate of the airflow flowing through the virtual holes and the solid wall obtained in the step S32) on the solid area heat sink area grid by dividing the volume of the grid, and taking the heat exchange rate as the size of an energy equation source item of the solid area heat sink area grid.

S42) the grid source item size of the last step is called, and the relaxation factor method source item calculation is carried out. And (4) carrying out weighted summation on the size of the source item obtained by the calculation in the step S41) and the size of the source item obtained by the calculation in the previous step by using a relaxation factor, namely:

s41) Source item relaxation factor + previous Source item (1-relaxation factor)

As the source item size is calculated at this step, the source item is slowly applied, making the calculation more stable.

S43) saving the size of the calculated source item in the step, and returning the value to a corresponding control equation to finish the addition and calculation of the source item method.

In the description herein, reference to the description of the term "one embodiment/manner," "some embodiments/manners" or "specific examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment/manner or example is included in at least one embodiment/manner or example of the application. In this specification, the schematic representations of the terms used above are not necessarily intended to be the same embodiment/mode or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments/modes or examples. Furthermore, the various embodiments/aspects or examples and features of the various embodiments/aspects or examples described in this specification can be combined and combined by one skilled in the art without conflicting therewith.

Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present application, "plurality" means at least two, e.g., two, three, etc., unless explicitly specified otherwise.

It will be understood by those skilled in the art that the foregoing embodiments are merely for clarity of illustration of the disclosure and are not intended to limit the scope of the disclosure. Other variations or modifications may be made to those skilled in the art, based on the above disclosure, and still be within the scope of the present disclosure.

Claims

1. The fluid-solid-heat coupling simulation method based on the virtual hole grid self-identification technology is characterized by comprising the following steps of:

generating virtual cooling hole information of a non-porous flame tube wall mesh model according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement;

determining an addition source item mesh based on the virtual cooling hole information, the addition source item mesh comprising an inlet face mesh of virtual cooling holes, an outlet face mesh of virtual cooling holes, and a solid domain heat sink region mesh;

completing a fluid-solid-thermal coupling simulation based on the grid of added source items;

the method further comprises the following steps: calculating the effective circulation area of the virtual cooling hole based on the adding source item grid;

calculating the effective flow area of the virtual cooling hole based on the grid of added source items includes: projecting the added source item mesh onto a point cloud plane of a corresponding virtual cooling hole to form an added source item mesh projection graph, sequentially connecting mesh internal vertexes of the added source item mesh projection graph, a set intersection point and a shape line point cloud point to form an overlapped polygon, and calculating an effective flow area according to the overlapped polygon, wherein the set intersection point is an intersection point of a mesh edge of the added source item mesh projection graph and a point cloud shape line;

completing a fluid solid thermal coupling simulation based on the grid of added source items comprises: performing a fluid-solid-thermal coupling simulation based on adding a source item grid and the effective flow area.

2. The method of claim 1, wherein generating virtual cooling hole information for the non-porous structured flame tube wall mesh model based on the hole arrangement requirements of the cooling holes comprises:

acquiring hole arrangement requirements and hole basic information of cooling holes, wherein the hole basic information comprises: hole center position coordinate information, hole direction vector information and hole cross section shape information of a first hole in each row of holes;

if the hole arrangement requirements comprise the requirements for avoiding the main combustion holes and the mixing holes, generating hole center coordinate information and expansion line point cloud coordinate information of each cooling hole according to the hole basic information and the hole arrangement requirements of the cooling holes, judging whether the cooling holes meet a first preset condition and a second preset condition, and generating virtual cooling hole information based on the hole basic information of the cooling holes meeting the first preset condition and the second preset condition;

3. The method of claim 2, wherein generating the virtual cooling hole information based on hole basic information of cooling holes satisfying a first preset condition and satisfying a second preset condition comprises:

generating a vertical hole center point and an expanded line point cloud according to the information of the vertical hole, judging whether the vertical hole meets a first preset condition and a second preset condition, if so, replacing the original cooling hole with the vertical hole, otherwise, deleting the cooling hole, not replacing the vertical hole, and generating hole information of all the cooling holes based on all the original cooling holes and the vertical holes which meet the first preset condition and the second preset condition.

4. The method of claim 1, wherein determining to add a source item mesh based on the virtual cooling hole information comprises:

orthographically projecting each outlet wall surface grid to a point cloud plane of a virtual cooling hole, and if the outlet wall surface grid projection graph is inside or intersected with a point cloud line of the virtual cooling hole and the projection length of the outlet wall surface grid projection does not exceed half of the hole depth of the virtual cooling hole, determining the outlet wall surface grid as the outlet surface grid of the virtual cooling hole;

and judging each grid in the solid area of the wall of the flame tube, and if the central point of the grid is positioned in the virtual hole, determining the grid as the solid area heat sink area grid of the virtual cooling hole.

5. The method of claim 1, wherein the method further comprises: determining an inlet data collection area and an outlet data collection area of a virtual cooling hole based on the virtual cooling hole information;

judging the gas flow direction of the virtual cooling hole according to flow-thermodynamic data, and calculating and obtaining the source item surface density of a fluid domain grid control equation, wherein the control equation comprises at least one of a continuity equation, a momentum equation and an energy equation;

completing the fluid-solid-thermal coupling simulation based on adding the source item grid and the effective flow area comprises the following steps: calculating a first source term of a fluid domain added source term grid based on the gas flow direction of the virtual cooling holes, the area density of the source term and the effective flow area of the grid, adding the first source term to a fluid domain grid control equation, calculating a second source term of a solid domain heat sink area grid energy equation according to the flow rate, the temperature and the wall temperature of the virtual holes and air holes, and adding the second source term to the energy equation of the solid domain heat sink area grid.

6. The method of claim 5, wherein calculating the added source term for the fluid domain added source term mesh based on the virtual cooling hole gas flow direction, the source term areal density, and the mesh effective flow area comprises: calculating the size of a third source item of a current calculation step of adding a source item grid according to the gas flow direction of the virtual cooling hole, the area density of the source item and the effective circulation area of the grid, acquiring the size of a fourth source item of a previous calculation step, and calculating the first source item through a relaxation factor method according to the size of the third source item and the size of the fourth source item; calculating the source term of the solid region heat sink region grid energy equation according to the virtual hole and air hole flow, the temperature and the wall temperature comprises the following steps: and calculating the size of a fifth source item of the current calculation step of the solid region heat sink region grid according to the flow, the temperature and the wall temperature of the virtual holes and the air holes, acquiring the size of a sixth source item of the previous calculation step, and calculating the second source item by a relaxation factor method according to the size of the fifth source item and the size of the sixth source item.

7. A flow-solid-thermal coupling simulation device based on a virtual hole grid self-recognition technology is characterized by comprising the following components:

the virtual cooling hole information confirmation module is used for generating virtual cooling hole information of the flame cylinder wall grid model with the non-porous structure according to the hole arrangement requirement of the cooling holes; the virtual cooling hole information comprises hole center coordinate information, hole direction vector information and hole cross section shape information of virtual cooling holes meeting the hole arrangement requirement;

a computing module for completing a fluid-solid-thermal coupling simulation based on the grid of added source items;

the device also comprises an effective circulation area calculation module, a virtual cooling hole calculation module and a virtual cooling hole calculation module, wherein the effective circulation area calculation module is used for calculating the effective circulation area of the virtual cooling hole based on the adding source item grid;

calculating the effective flow area of the virtual cooling hole based on the grid of added source terms includes: projecting the added source item mesh onto a point cloud plane of a corresponding virtual cooling hole to form an added source item mesh projection graph, sequentially connecting mesh internal vertexes of the added source item mesh projection graph, a set intersection point and a shape line point cloud point to form an overlapped polygon, and calculating an effective flow area according to the overlapped polygon, wherein the set intersection point is an intersection point of a mesh edge of the added source item mesh projection graph and a point cloud shape line;

completing a fluid solid thermal coupling simulation based on the grid of added source items comprises: performing a fluid-solid-thermal coupling simulation based on adding a source term grid and the effective flow area.

8. Computer-readable storage medium, characterized in that it stores a computer program which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 6.

9. Computer arrangement comprising a memory and a processor, characterized in that the memory stores a computer program which, when executed by the processor, carries out the steps of the method according to any one of claims 1 to 6.