Interpolation method, device and equipment for fluid-solid coupling interface in engine and medium

Abstract

The invention relates to a fluid-solid coupling interface interpolation method, device and equipment in an engine and a medium. Comprising the following steps: acquiring an interpolation surface node and a to-be-interpolated surface node; grouping the interpolation surface nodes to obtain a plurality of interpolation surface node subsets; distributing each interpolation surface node subset and the to-be-interpolated surface node to different processes by adopting a parallel computing strategy; in each process, adaptively setting a support radius, and then solving to obtain a radial basis interpolation function; broadcasting a radial basis interpolation function, and performing physical field interpolation on surface nodes to be interpolated; in an interpolation process, judging whether a surface node to be interpolated is a single-grouping influence area or a double-grouping influence area, and if the surface node to be interpolated is the single-grouping influence area, calculating an interpolation result by adopting a single-grouping interpolation function; and if the influence area is a double-grouping influence area, weighting by adopting two adjacent radial basis interpolation functions to obtain an interpolation result. According to the method, both precision and efficiency are considered, and accurate data interaction of the large-scale fluid-solid coupling interface can be ensured.

Description

Technical Field

[0001] The invention relates to the technical field of fluid-solid coupling simulation, and in particular to a method, device, equipment and medium for fluid-solid coupling interface interpolation in an engine. Background Art

[0002] Today, air-breathing aircraft technology has become a priority technology, and its core is engine technology. Due to the needs of low-cost space transportation and other aspects, wide-range Mach number, long-endurance reusable engine technology has become a hot issue. The current engine design mainly relies on a large number of tests, but there are still problems such as the lack of high-temperature long-range test equipment, high test costs, and lack of full-scale global data. Long-endurance, reusable engines face severe thermal safety issues. Long-term ultra-high temperatures, limited wall cooling flow, and especially repeated use require the thermal interface to remain intact and the structure to be free of fatigue damage. These characteristics make the study of the fluid-solid thermal coupling characteristics of the engine a necessary part of the design.

[0003] Numerical calculation is currently a feasible research method, which strongly supports the research, development and design of scramjet engines. In the field of fluid-solid thermal coupling research, a multi-physics coupling simulation framework such as flow, structural heat transfer and structural deformation has been established, and a variety of coupling strategies and coupling methods have been proposed. A key factor in fluid-solid thermal multi-field coupling is the efficiency and accuracy of fluid-solid interface data transmission. Currently, local interpolation methods such as mapping point interpolation method and weighted residual method, as well as global interpolation methods represented by radial basis function interpolation method and Shepard method have been proposed. When these traditional methods are applied to fluid-solid interface data interaction of large-scale complex configurations, there are still problems of low efficiency and lack of accuracy. Summary of the invention

[0004] Based on this, it is necessary to provide a method, device, equipment and medium for interpolation of fluid-solid coupling interface in an engine with high efficiency and high precision, which is suitable for the internal flow of the engine, in order to solve the above technical problems.

[0005] A method for interpolation of fluid-solid coupling interface in an engine, the method comprising: Generate engine fluid domain mesh and solid domain mesh, obtain the coordinates of interpolation surface nodes of the fluid domain and the coordinates of surface nodes to be interpolated of the solid domain; The interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes to obtain a plurality of interpolation surface node subsets; a parallel computing strategy is adopted to allocate each of the interpolation surface node subsets to different processes, and the surface nodes to be interpolated are evenly distributed to each process; In each process, the support radius of the interpolation surface node subset is adaptively set, and then a linear algebraic equation system is constructed according to the radial basis interpolation function, and the equation system is solved to obtain the radial basis interpolation function based on the interpolation surface node subset; The radial basis interpolation function in each process is broadcast to all processes, and then the radial basis interpolation function of the interpolation surface node subset is used to interpolate the physical field of the surface node to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, two adjacent radial basis interpolation functions are weighted to obtain the interpolation result.

[0006] A fluid-solid coupling interface interpolation device in an engine, the device comprising: An initialization module is used to generate a fluid domain mesh and a solid domain mesh of the engine, and obtain the coordinates of the interpolation surface nodes of the fluid domain and the coordinates of the surface nodes to be interpolated of the solid domain; A node grouping and parallel allocation module is used to group the interpolation surface nodes according to the longest direction of the fluid domain surface nodes to obtain a plurality of interpolation surface node subsets; a parallel computing strategy is adopted to allocate each of the interpolation surface node subsets to different processes, and the surface nodes to be interpolated are evenly distributed to each process; A radial basis interpolation function calculation module is used to adaptively set the support radius of the interpolation surface node subset in each process, then construct a linear algebraic equation system according to the radial basis interpolation function, and solve the equation system to obtain a radial basis interpolation function based on the interpolation surface node subset; The interpolation module is used to broadcast the radial basis interpolation function in each process to all processes, and then use the radial basis interpolation function of the interpolation surface node subset to interpolate the physical field of the surface node to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, two adjacent radial basis interpolation functions are weighted to obtain the interpolation result.

[0007] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the fluid-solid coupling interface interpolation method in the engine when executing the computer program.

[0008] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the fluid-solid coupling interface interpolation method in the engine.

[0009] The above-mentioned method, device, equipment and medium for interpolation of fluid-solid coupling interface in the engine generate the engine fluid domain grid and solid domain grid to obtain the coordinates of the interpolation surface nodes of the fluid domain and the coordinates of the surface nodes to be interpolated of the solid domain; the interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes to obtain a number of interpolation surface node subsets; a parallel computing strategy is adopted to allocate each interpolation surface node subset to different processes, and the surface nodes to be interpolated are evenly distributed to each process; in each process, the support radius of the interpolation surface node subset is adaptively set, and then according to the radial basis function, the support radius of the interpolation surface node subset is calculated. The interpolation function constructs a linear algebraic equation group and solves the equation group to obtain a radial basis interpolation function based on a subset of interpolation surface nodes; the radial basis interpolation function in each process is broadcast to all processes, and then the radial basis interpolation function of the subset of interpolation surface nodes is used to interpolate the physical field of the surface nodes to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, the interpolation result is obtained by weighting two adjacent radial basis interpolation functions.

[0010] The present invention groups the interpolation surface nodes and assigns each group to different processes, so that the complexity of solving the subsequent linear algebraic equations is approximately reduced from Reduce to , significantly improves the computational efficiency, among which, Represents the number of boundary nodes. By adaptively setting the support radius of the interpolation surface node subset, it is ensured that the interpolation function has high accuracy in the local area, and at the same time, the positive definiteness problem of the linear system caused by the excessive support radius is avoided. When interpolating, by judging whether the surface node to be interpolated is a single-group influence area or a double-group influence area, different interpolation calculation methods are adopted to reduce unnecessary calculation overhead, while ensuring the interpolation accuracy, and further improve the interpolation efficiency; in addition, in the double-group influence area, two adjacent interpolation functions are used for weighted interpolation, which ensures a smooth transition of the interpolation result at the junction of the groups. The present invention takes into account both accuracy and efficiency, is particularly suitable for large-scale complex configurations, and can ensure accurate interaction of large-scale fluid-solid coupling interface data. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the structures shown in these drawings without paying creative work.

[0012] Figure 1It is a schematic flow chart of a fluid-solid coupling interface interpolation method in an engine in one embodiment; Figure 2 A schematic diagram of a framework of a fluid-solid coupling interface interpolation method in an engine in one embodiment; Figure 3 A schematic diagram of interpolation surface node grouping in one embodiment, wherein: Figure 3 (a) is a 3D stereoscopic view. Figure 3 (b) is a 2D top view; Figure 4 is a graph of simple fluid-solid interface interpolation results in one embodiment, where: Figure 4 (a) is the original interpolation surface and physical quantity p The distribution map of Figure 4 (b) is the distribution diagram of physical quantities after interpolation. Figure 4 (c) is the relative error distribution diagram between the exact value and the interpolation result. Figure 4 (d) is the relative error scatter plot; Figure 5 This is a diagram of the interpolation results of the fluid-solid interface of a complex 1 / 4 configuration axisymmetric engine in one embodiment, where: Figure 5 (a) is the heat flux distribution on the engine surface obtained by computational fluid dynamics (CFD). Figure 5 (b) is the heat flux distribution on the solid surface obtained by interpolation using the method proposed in the present invention; Figure 6 is a structural block diagram of a fluid-solid coupling interface interpolation device in an engine in one embodiment; Figure 7 FIG. 4 is a diagram showing the internal structure of a computer device in one embodiment.

[0013] The purpose, features and advantages of the present invention will be further described with reference to the accompanying drawings in conjunction with the embodiments. DETAILED DESCRIPTION

[0014] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0015] It can be understood that the technical solutions between the various embodiments of the present invention can be combined with each other, but it must be based on the fact that ordinary technicians in the field can implement it. When the combination of technical solutions is contradictory or cannot be implemented, it should be deemed that such combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0016] The following will describe the implementation of the present invention in detail with reference to the accompanying drawings in the embodiment diagram of the present invention.

[0017] Example 1 This embodiment discloses a method for interpolating a fluid-solid coupling interface in an engine. By grouping interpolation surface nodes and assigning each group to a different process, the complexity of solving the subsequent linear algebraic equations is approximately reduced from Reduce to , significantly improves the computational efficiency, among which, Represents the number of boundary nodes. By adaptively setting the support radius of the interpolation surface node subset, it is ensured that the interpolation function has high accuracy in the local area, and at the same time, the positive definiteness problem of the linear system caused by the excessive support radius is avoided. When interpolating, by judging whether the surface node to be interpolated is a single-group influence area or a double-group influence area, different interpolation calculation methods are adopted to reduce unnecessary calculation overhead, while ensuring the interpolation accuracy, and further improve the interpolation efficiency; in addition, in the double-group influence area, two adjacent interpolation functions are used for weighted interpolation, which ensures a smooth transition of the interpolation result at the junction of the groups. The present invention takes into account both accuracy and efficiency, is particularly suitable for large-scale complex configurations, and can ensure accurate interaction of large-scale fluid-solid coupling interface data.

[0018] like Figure 1 and Figure 2 As shown, the fluid-solid coupling interface interpolation method in the engine provided by this embodiment includes the following steps: Step 201 : Generate a fluid domain mesh and a solid domain mesh of the engine, and obtain the coordinates of the interpolation surface nodes of the fluid domain and the coordinates of the surface nodes to be interpolated of the solid domain.

[0019] Step 202, grouping the interpolation surface nodes according to the longest direction of the fluid domain surface nodes to obtain a number of interpolation surface node subsets; using a parallel computing strategy, allocating each interpolation surface node subset to a different process, and evenly distributing the surface nodes to be interpolated to each process.

[0020] Step 203, in each process, adaptively set the support radius of the interpolation surface node subset, then construct a linear algebraic equation system according to the radial basis interpolation function, and solve the equation system to obtain the radial basis interpolation function based on the interpolation surface node subset.

[0021] Step 204, broadcast the radial basis interpolation function in each process to all processes, and then use the radial basis interpolation function of the interpolation surface node subset to interpolate the physical field of the surface node to be interpolated; during the interpolation process, determine whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, use a single group interpolation function to calculate the interpolation result; if it is a double-group influence area, use two adjacent radial basis interpolation functions to weight the interpolation result.

[0022] In the specific implementation process of step 201, the engine fluid domain mesh is generated by CFD, and the coordinates and heat flux density of the fluid domain surface nodes are output; at the same time, computational solid mechanics (CSD) generates the engine solid domain mesh and outputs the solid wall node coordinates. Then, the wall node coordinates and heat flux density are read through the fluid-solid coupling interpolation program.

[0023] Usually, the radial basis function used for fluid-structure interaction interpolation is Wendland's C 2 Function, Wendland's C 2 The function is a kind of compact function, and its specific form is: (1) In the formula, , represents the position vector of the point to be interpolated, Indicates The position vectors of the control points, Indicates the support radius.

[0024] In the specific implementation process of step 202, MPI and OpenMP are used to implement a parallel architecture for parallel computing. The number of interpolation surface node groups, the number of parallel cores, and the overlap coefficient between groups are set, the interpolation surface nodes are grouped according to the longest direction and assigned to different MPI processes, and the surface nodes to be interpolated are evenly grouped into different processes. For the surface nodes to be interpolated, they are evenly distributed to each process according to the node index. It can be understood that the parallel architecture generally includes sub-processes and main processes. In the process of this embodiment, both the main process and the sub-process participate in the calculation.

[0025] Specifically, the interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes to obtain several interpolation surface node subsets, including: Set the number of groups and overlap coefficient between groups ;generally Set to an integer multiple of the number of parallel cores to balance the load; coefficient f Set to a number between 0 and 1. Assume that the longest direction of the nodes on the surface of the fluid domain is direction, based on the overlap coefficient between groups , the surface of the fluid domain is The direction is divided into equal lengths copy, get A subset of interpolated surface nodes.

[0026] Among them, based on the overlap coefficient between groups , the surface of the fluid domain is The direction is divided into equal lengths The grouping rules are: Calculate the group determination integer Determining decimals with grouping : (2) (3) In the formula, floor is the rounding function; Represents the coordinates of the node along the maximum scale direction; For all interpolation surface nodes x The minimum value of the coordinates; is the value of all interpolation surface nodes. x The maximum value of the coordinates.

[0027] The interpolation surface nodes are grouped, when , , the interpolation surface node is located in the first interpolation surface node subset.

[0028] when , When the interpolation surface nodes are only located at the In a subset of interpolation surface nodes; when and When the interpolation surface nodes are only located at the In a subset of interpolation surface nodes; when and When the interpolation surface nodes are simultaneously located at and In a subset of interpolation surface nodes; when and When the interpolation surface nodes are simultaneously located at and A subset of interpolation surface nodes.

[0029] in, Characterizes the size of the overlapping area between grouped nodes, indicates completely non-overlapping groups, Indicates that all groups overlap each other, and the width of the group area is ,in, According to experience, the preferred A value of 0.4 can get better results. After grouping, the nodes are evenly distributed to different processes. The purpose of setting the overlapping area is to eliminate the discontinuity of the interpolation results near the boundary caused by grouping.

[0030] In the specific implementation process of step 203, the support radius of the interpolation surface node subset is adaptively set, and the expression is: (4) In the formula, For all interpolation surface nodes The minimum value of the coordinates; is the value of all interpolation surface nodes. The maximum value of the coordinates; Indicates the adjustable multiple; Indicates that the surface of the fluid domain is The number of equal lengths in the direction.

[0031] It is worth mentioning that Generally, a number greater than 1 is taken to ensure that all nodes in the group contribute; but it should not be too large, otherwise it will lead to the positive definiteness problem of the linear system. Preferably, Take 3.

[0032] In each process, a radial basis interpolation function is constructed according to the subset of interpolation surface nodes assigned to it, and the expression is as follows: (5) In order to solve the coefficients of the interpolation function , so that the interpolation function is exactly equal to its physical quantity at the nodes of the interpolation surface, and a system of linear algebraic equations is constructed. The linear algebraic equations in the dimensional space are expressed as: (6) In the formula, express The physical field vector at the position; Indicates The weight coefficient vector of control points; represents the position vector; Indicates The position vector of the control point; Represents control points Control Points The impact of Represents control points The position vector of Represents control points The position vector of Represents control points Corresponding to The unknown coefficients of dimension, Represents control points The first Dimension value. Control point With control point Maximum , is the number of nodes in the interpolation surface group.

[0033] LDL using the open source Eigen library T Cholevsky decomposition method, solving the coefficients of the interpolation function based on a system of linear algebraic equations , and obtain the radial basis interpolation function based on the subset of interpolation surface nodes.

[0034] In the specific implementation process of step 204, the radial basis interpolation function in each process is broadcast to all processes, mainly broadcasting the weight coefficient vector of the control point and the position vector of the control point , and then substitute the coordinates of the surface nodes to be interpolated into the formula (5) Interpolate the physical field. During the interpolation process, OpenMP is used to open multiple threads to further accelerate the process.

[0035] If the coordinates of the nodes on the surface to be interpolated are If , then it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. The specific method is: Get the coordinates of the surface nodes to be interpolated , based on a subset of interpolated surface nodes Loop, during the loop, the index is , Cycle from 1 to .

[0036] If satisfied , it means that the surface nodes to be interpolated are affected by and The common influence of the groups is the double-group influence area, which directly jumps out of the loop.

[0037] If not satisfied , it means that the surface node to be interpolated is a single-group influence area. When it is a single-group influence area, it is also necessary to determine whether it is affected by the first or second The influence of a subset of interpolation surface nodes, if not satisfied, In 1 to The judgment is repeated in a cycle.

[0038] The single-group impact area is divided into three cases: when , it is only affected by the first subset of interpolation surface nodes.

[0039] when , it is only affected by A subset of interpolated surface nodes is affected.

[0040] when , which is only indexed by The interpolation surface node subset affects In 1 to between.

[0041] If it is a single-group influence area, the radial basis interpolation function of a single interpolation surface node subset is used for interpolation. If it is a double-group influence area, the interpolation result is obtained by weighting two adjacent radial basis interpolation functions. The calculation expression is: (7) (8) In the formula, Represents the final interpolation result obtained by combining two subsets of interpolation surface nodes; and Respectively indicate that the index is and The physical field is calculated using the radial basis interpolation function of a subset of interpolation surface nodes; represents the weighting coefficient; Indicates that the index is Interpolation surface nodes in a subset of interpolation surface nodes The minimum value of the coordinates; Indicates that the index is Interpolation surface nodes in a subset of interpolation surface nodes The maximum value of the coordinate. The interpolation result is obtained by weighting two adjacent radial basis interpolation functions, which can ensure a smooth transition between the interpolation results in the overlapping area (the area affected by two groups at the same time) and the surrounding area.

[0042] Finally, the interpolation results are collected and sent to the main process, which then outputs them. This completes the interface data transmission and completes the fluid-solid coupling interface interpolation in the engine.

[0043] The fluid-solid coupling interface interpolation method in the engine provided by the present invention is to group the interpolation surface nodes into The RBF interpolation function is constructed independently based on each group using MPI and OpenMP parallel technology, which significantly improves the efficiency; and an overlapping area is set between adjacent groups to ensure a smooth transition of the interpolation result. The method proposed in the present invention takes into account both accuracy and efficiency, ensuring the accuracy of fluid-solid coupling data interaction in large-scale complex configurations.

[0044] It is worth noting that although Figure 2 The framework diagram of the CPC has been developed parallel processes, but in fact the number of processes opened can be equal to or less than the number of interpolation surface node subsets , that is, each process processes multiple interpolation surface node subsets. When the number of processes is equal to the number of interpolation surface node subsets When , each process is assigned a subset of interpolation surface nodes for calculation; when the number of processes is less than the number of interpolation surface node subsets When , each process is assigned more than two surface node subsets. Preferably, the number of surface node subsets can be Set it to an integer multiple of the number of processes to facilitate even distribution and thus ensure parallel load balancing.

[0045] In one embodiment, Figure 4 As shown, a simple fluid-solid interface interpolation result diagram is provided to verify the effectiveness of the method proposed in the present invention.

[0046] The original interpolation surface point cloud and physical quantities are generated by functions. The coordinates are evenly generated in the range of 0 to 1. The coordinates are also uniformly generated in the range of 0 to 1, with 10 values. The coordinates are generated using the following function: (9) Surface physical quantity p Function generation using coordinates: (10) Original interpolation surface nodes and physical quantities p The distribution of Figure 4 (a) shown.

[0047] The surface point cloud to be interpolated is generated in a similar way, except that and 30 points are used for the generation of coordinates. The surface physical quantity is obtained by interpolation using the method provided by the present invention. p The distribution of Figure 4 (b) It can be seen that the distribution of the physical quantity after interpolation is the same as that of the original surface.

[0048] In order to quantitatively compare the reliability of the interpolation results, the physical quantities of the nodes on the interpolated surface are calculated using formula (10) and regarded as the exact values. , and compare it with the interpolation result to calculate the relative error: (11) Figure 4 (c) and Figure 4 (d) shows the relative error distribution diagram and relative error scatter diagram respectively, and it can be seen that the maximum error is less than 1.6%, and the errors of most points are less than 0.2%, which proves the effectiveness of the interpolation method provided by the present invention.

[0049] In one embodiment, Figure 5 As shown, a fluid-solid interface interpolation result diagram of a complex 1 / 4 configuration axisymmetric engine is provided, verifying the effectiveness of the method provided by the present invention on large-scale configurations.

[0050] The static temperature of the inlet flow at the engine inlet is 945.8K, the static pressure of the inlet flow is 77100Pa, the Mach number is 2.5, the wall temperature is 300K, and ethylene fuel is used. The heat flux of the engine wall is calculated by CFD as follows: Figure 5 (a) is shown. The method provided by the present invention is used to transfer data, and the heat flux on the wall of the solid domain is obtained as follows Figure 5 (b) It can be seen that a consistent wall heat flux distribution is obtained.

[0051] Although this embodiment Figure 1 The steps in the process are shown in sequence as indicated by the arrows, but the steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of the steps, and the steps can be executed in other orders. Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0052] Example 2 Based on the method for interpolating the fluid-solid coupling interface in the engine in Example 1, this embodiment discloses an interpolation device for the fluid-solid coupling interface in the engine, such as Figure 6 As shown, the interpolation device for the fluid-solid coupling interface in the engine includes: an initialization module 401, a node grouping and parallel allocation module 402, a radial basis interpolation function calculation module 403 and an interpolation module 404, wherein: The initialization module 401 is used to generate the engine fluid domain mesh and the solid domain mesh, and obtain the coordinates of the interpolation surface nodes of the fluid domain and the coordinates of the surface nodes to be interpolated of the solid domain.

[0053] The node grouping and parallel allocation module 402 is used to group the interpolation surface nodes according to the longest direction of the fluid domain surface nodes to obtain a plurality of interpolation surface node subsets; adopt a parallel computing strategy to allocate each of the interpolation surface node subsets to different processes, and evenly distribute the surface nodes to be interpolated to each process.

[0054] The radial basis interpolation function calculation module 403 is used to adaptively set the support radius of the interpolation surface node subset in each process, then construct a linear algebraic equation group according to the radial basis interpolation function, and solve the equation group to obtain the radial basis interpolation function based on the interpolation surface node subset.

[0055] The interpolation module 404 is used to broadcast the radial basis interpolation function in each process to all processes, and then use the radial basis interpolation function of the interpolation surface node subset to interpolate the physical field of the surface node to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, two adjacent radial basis interpolation functions are weighted to obtain the interpolation result.

[0056] In this embodiment, the specific working process and working principle of the initialization module 401, the node grouping and parallel allocation module 402, the radial basis interpolation function calculation module 403 and the interpolation module 404 are the same as those in the method of Embodiment 1, so they are not described in detail in this embodiment. Each unit module can be implemented in whole or in part by software, hardware and a combination thereof, and each unit module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to each of the above unit modules.

[0057] Example 3 like Figure 7 The terminal device disclosed in this embodiment includes a transmitter, a receiver, a memory, and a processor. The transmitter is used to send instructions and data, the receiver is used to receive instructions and data, the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions stored in the memory to implement the method in the above-mentioned embodiment 1.

[0058] It should be noted that the above memory can be independent or integrated with the processor. When the memory is independently provided, the terminal device further includes a bus for connecting the memory and the processor.

[0059] Example 4 This embodiment discloses a computer-readable storage medium, in which computer-executable instructions are stored. When a processor executes the computer-executable instructions, the method in the above-mentioned embodiment 1 is implemented.

[0060] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0061] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0062] The above-mentioned embodiments only express several implementation modes of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the invention. It should be pointed out that, for a person of ordinary skill in the art, several modifications and improvements can be made without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the attached claims.

Claims

1. A method for interpolation of fluid-solid coupling interface in an engine, characterized in that: The method comprises: Generate engine fluid domain mesh and solid domain mesh, obtain the coordinates of interpolation surface nodes of the fluid domain and the coordinates of surface nodes to be interpolated of the solid domain; The interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes to obtain a plurality of interpolation surface node subsets; a parallel computing strategy is adopted to allocate each of the interpolation surface node subsets to different processes, and the surface nodes to be interpolated are evenly distributed to each process; In each process, the support radius of the interpolation surface node subset is adaptively set, and then a linear algebraic equation system is constructed according to the radial basis interpolation function, and the equation system is solved to obtain the radial basis interpolation function based on the interpolation surface node subset; The radial basis interpolation function in each process is broadcast to all processes, and then the radial basis interpolation function of the interpolation surface node subset is used to interpolate the physical field of the surface node to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, two adjacent radial basis interpolation functions are weighted to obtain the interpolation result.

2. The method for interpolation of fluid-solid coupling interface in an engine according to claim 1, characterized in that: The interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes to obtain several interpolation surface node subsets, including: Set the number of groups and overlap coefficient between groups ; Assume that the longest direction of the fluid domain surface node is direction, based on the overlap coefficient between groups , the surface of the fluid domain is The lengths of the directions are equally divided copy, get A subset of interpolated surface nodes.

3. The method for interpolation of fluid-solid coupling interface in an engine according to claim 2, characterized in that: Based on the overlap coefficient between groups , the surface of the fluid domain is The lengths of the directions are equally divided The grouping rules are: Calculate the group determination integer Determining decimals with grouping : ； ； In the formula, floor is the rounding function; Represents the coordinates of the node along the maximum scale direction; For all interpolation surface nodes x The minimum value of the coordinates; is the value of all interpolation surface nodes. x The maximum value of the coordinates; The interpolation surface nodes are grouped, when , When , the interpolation surface nodes are only located in the first interpolation surface node subset; when , When the interpolation surface nodes are only located at the In a subset of interpolation surface nodes; when and When the interpolation surface nodes are only located at the In a subset of interpolation surface nodes; when and When the interpolation surface nodes are simultaneously located at and In a subset of interpolation surface nodes; when and When the interpolation surface nodes are simultaneously located at and A subset of interpolation surface nodes.

4. The method for interpolation of fluid-structure coupling interface in an engine according to any one of claims 1 to 3, characterized in that: The support radius of the interpolation surface node subset is adaptively set, and the expression is: ； In the formula, For all interpolation surface nodes The minimum value of the coordinates; is the value of all interpolation surface nodes. The maximum value of the coordinates; Indicates the adjustable multiple; Indicates that the surface of the fluid domain is The number of equal lengths in the direction; among them, Represents the coordinate of the node along the maximum direction of the scale.

5. The method for interpolation of fluid-structure coupling interface in an engine according to claim 4, characterized in that: A linear algebraic equation system is constructed according to the radial basis interpolation function, and the linear algebraic equation system is solved to obtain a radial basis interpolation function based on the interpolation surface node subset, including: Each process constructs a radial basis interpolation function based on the subset of interpolation surface nodes assigned to it, and the expression is as follows: ； In order to solve the coefficients of the interpolation function , so that the interpolation function is exactly equal to its physical quantity at the nodes of the interpolation surface, and a system of linear algebraic equations is constructed. The linear algebraic equations in the dimensional space are expressed as: ； In the formula, express The physical field vector at the position; Indicates The weight coefficient vector of control points; represents the position vector; Indicates The position vector of the control point; Represents control points Control Points The impact of Represents control points Corresponding to The unknown coefficients of dimension, Represents control points The first Dimension value; Solve for coefficients of interpolation functions based on a system of linear algebraic equations , and obtain a radial basis interpolation function based on the interpolation surface node subset.

6. The method for interpolation of fluid-structure coupling interface in an engine according to claim 5, characterized in that: During the interpolation process, determining whether the surface node to be interpolated is a single-group influence area or a double-group influence area includes: Get the coordinates of the surface nodes to be interpolated , based on a subset of interpolated surface nodes Loop, during the loop, the index is , Cycle from 1 to ; If satisfied , it means that the surface node to be interpolated is affected by the and The common influence of the groups is the double-group influence area; If not satisfied , it means that the surface node to be interpolated is a single-group influence area.

7. The method for interpolation of fluid-structure coupling interface in an engine according to claim 6, characterized in that: If it is a double-grouped influence area, the interpolation result is obtained by weighting two adjacent radial basis interpolation functions. The calculation expression is: ； ； In the formula, Represents the final interpolation result obtained by combining two subsets of interpolation surface nodes; and Respectively indicate that the index is and The physical field is calculated using the radial basis interpolation function of a subset of interpolation surface nodes; represents the weighting coefficient; Indicates that the index is Interpolation surface nodes in a subset of interpolation surface nodes The minimum value of the coordinates; Indicates that the index is Interpolation surface nodes in a subset of interpolation surface nodes The maximum value of the coordinates.

8. A fluid-solid coupling interface interpolation device in an engine, characterized in that: The device comprises: An initialization module is used to generate a fluid domain mesh and a solid domain mesh of the engine, and obtain the coordinates of the interpolation surface nodes of the fluid domain and the coordinates of the surface nodes to be interpolated of the solid domain; A node grouping and parallel allocation module is used to group the interpolation surface nodes according to the longest direction of the fluid domain surface nodes to obtain a plurality of interpolation surface node subsets; a parallel computing strategy is adopted to allocate each of the interpolation surface node subsets to different processes, and the surface nodes to be interpolated are evenly distributed to each process; A radial basis interpolation function calculation module is used to adaptively set the support radius of the interpolation surface node subset in each process, then construct a linear algebraic equation system according to the radial basis interpolation function, and solve the equation system to obtain a radial basis interpolation function based on the interpolation surface node subset; The interpolation module is used to broadcast the radial basis interpolation function in each process to all processes, and then use the radial basis interpolation function of the interpolation surface node subset to interpolate the physical field of the surface node to be interpolated; during the interpolation process, it is determined whether the surface node to be interpolated is a single-group influence area or a double-group influence area. If it is a single-group influence area, a single group interpolation function is used to calculate the interpolation result; if it is a double-group influence area, two adjacent radial basis interpolation functions are weighted to obtain the interpolation result.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for interpolation of fluid-structure coupling interface in an engine as described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for interpolation of fluid-structure coupling interface in an engine as claimed in any one of claims 1 to 7 are implemented.

Images