Interpolation Method, Device, Equipment and Medium for Fluid-Structure Interaction Interface in Engine

Abstract

The present invention relates to an interpolation method, device, equipment and medium for the internal fluid-structure coupling interface of an engine. The method includes: obtaining interpolation surface nodes and surface nodes to be interpolated; grouping the interpolation surface nodes to obtain several subsets of interpolation surface nodes; adopting a parallel computing strategy to allocate each subset of interpolation surface nodes and the surface nodes to be interpolated to different processes; in each process, adaptively setting a support radius, and then obtaining a radial basis interpolation function by solving; broadcasting the radial basis interpolation function to perform interpolation of the physical field for the surface nodes to be interpolated; during the interpolation process, determining whether the surface nodes to be interpolated are in a single-group influence area or a double-group influence area. If it is a single-group influence area, the interpolation result is calculated using a single-group interpolation function; if it is a double-group influence area, the interpolation result is obtained by weighting two adjacent radial basis interpolation functions. The present invention takes into account both accuracy and efficiency, and can ensure accurate data interaction at the large-scale fluid-structure coupling interface.

Description

Technical Field

[0001] The present invention relates to the technical field of fluid-structure interaction simulation, and particularly to an interpolation method, device, equipment and medium for the fluid-structure interaction interface in an engine. Background Technique

[0002] Currently, the technology of air-breathing aircraft has become a technology with priority development, and its core is engine technology. Due to the needs of low-cost space transportation and other aspects, the engine technology with a wide range of Mach numbers and long endurance and reusable has become a hot issue. The current design of engines mainly relies on a large number of tests, but there are still problems such as the lack of high-temperature long-range test devices, the too high test cost, and the lack of full-scale global data. The long-endurance and reusable engines face severe thermal safety problems. The characteristics of long-time ultra-high temperature, limited wall cooling flow rate, especially the requirement of maintaining the integrity of the thermal interface and no fatigue damage inside the structure for repeated use, make the research on the fluid-structure-thermal coupling characteristics of engines an essential link in the design.

[0003] Numerical calculation is a currently practical research method, which strongly supports the research and development design of scramjet engines. In the field of fluid-structure-thermal coupling research, a multi-physical field coupling simulation framework for fluid 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 the multi-field coupling of fluid-structure-thermal is the efficiency and accuracy of the data transfer at the fluid-structure interface. Currently, local interpolation methods represented by the mapped point interpolation method and the weighted residual method, and global interpolation methods represented by the radial basis function interpolation method and the Shepard method have been proposed. These traditional methods still have problems of too low efficiency and inability to guarantee accuracy when applied to the data interaction at the fluid-structure interface of large-scale complex configurations. Summary of the Invention

[0004] Based on this, it is necessary to provide an interpolation method, device, equipment and medium for the fluid-structure interaction interface in an engine that is applicable to the internal flow of the engine and has high efficiency and high accuracy in view of the above technical problems.

[0005] An interpolation method for the fluid-structure interaction interface in an engine, the method includes:

[0006] Generate the grid of the engine fluid domain and the grid of the solid domain, and obtain the coordinates of the interpolation surface nodes in the fluid domain and the coordinates of the surface nodes to be interpolated in the solid domain;

[0007] Group the interpolation surface nodes according to the longest direction of the surface nodes in the fluid domain to obtain several subsets of interpolation surface nodes; adopt a parallel computing strategy, allocate each subset of interpolation surface nodes to different processes, and evenly distribute the surface nodes to be interpolated to each process;

[0008] In each process, adaptively set the support radius of the interpolation surface node subset, then construct a system of linear algebraic equations according to the radial basis interpolation function, and solve the system of equations to obtain the radial basis interpolation function based on the interpolation surface node subset;

[0009] Broadcast the radial basis interpolation functions 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 the weighted sum of two adjacent radial basis interpolation functions to obtain the interpolation result.

[0010] An interpolation device for the internal fluid-structure coupling interface of an engine, the device includes:

[0011] An initialization module, configured to generate an engine fluid domain grid and a solid domain grid, and obtain the coordinates of the interpolation surface nodes in the fluid domain and the coordinates of the surface nodes to be interpolated in the solid domain;

[0012] A node grouping and parallel allocation module, configured to group the interpolation surface nodes according to the longest direction of the fluid domain surface nodes to obtain several interpolation surface node subsets; adopt a parallel computing strategy to allocate each interpolation surface node subset to different processes, and evenly distribute the surface nodes to be interpolated to each process;

[0013] A radial basis interpolation function calculation module, configured to adaptively set the support radius of the interpolation surface node subset in each process, then construct a system of linear algebraic equations according to the radial basis interpolation function, and solve the system of equations to obtain the radial basis interpolation function based on the interpolation surface node subset;

[0014] An interpolation module, configured to broadcast the radial basis interpolation functions 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 the weighted sum of two adjacent radial basis interpolation functions to obtain the interpolation result.

[0015] A computer device, including a memory and a processor, the memory stores a computer program, and when the processor executes the computer program, the steps of the interpolation method for the internal fluid-structure coupling interface of the engine are implemented.

[0016] A computer-readable storage medium stores a computer program thereon, and when the computer program is executed by a processor, the steps of the interpolation method for the fluid-structure coupling interface in the engine are implemented.

[0017] For the above-mentioned interpolation method, device, equipment and medium for the fluid-structure coupling interface in the engine, by generating a fluid domain mesh and a solid domain mesh of the engine, the coordinates of the interpolation surface nodes in the fluid domain and the coordinates of the surface nodes to be interpolated in the solid domain are obtained; the interpolation surface nodes are grouped according to the longest direction of the surface nodes in the fluid domain to obtain several subsets of interpolation surface nodes; a parallel computing strategy is adopted to allocate each subset of interpolation surface nodes to different processes, and the surface nodes to be interpolated are evenly distributed to each process; in each process, the support radius of the subset of interpolation surface nodes 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 subset of interpolation surface nodes; the radial basis interpolation functions in each process are broadcast to all processes, and then the physical field of the surface nodes to be interpolated is interpolated using the radial basis interpolation functions of the subset of interpolation surface nodes; during the interpolation process, it is judged 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, the interpolation result is calculated using a single-group interpolation function; if it is a double-group influence area, the interpolation result is obtained by weighting two adjacent radial basis interpolation functions.

[0018] In the present invention, by grouping the interpolation surface nodes and allocating each group to different processes, the complexity of solving the subsequent linear algebraic equation system is approximately reduced from to , significantly improving the calculation efficiency, where represents the number of boundary nodes. By adaptively setting the support radius of the subset of interpolation surface nodes, it is ensured that the interpolation function has high precision in the local area, and at the same time, the problem of positive definiteness of the linear system caused by too large a support radius is avoided. When performing interpolation, 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, reducing unnecessary calculation overhead, further improving the interpolation efficiency while ensuring the interpolation accuracy; in addition, in the double-group influence area, two adjacent interpolation functions are weighted for interpolation, ensuring a smooth transition of the interpolation result at the junction of the groups. The present invention takes both accuracy and efficiency into account, is particularly suitable for large-scale complex configurations, and can ensure accurate data interaction at the large-scale fluid-structure coupling interface. Description of the Drawings

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on the structures shown in these drawings.

[0020] Figure 1 It is a schematic flowchart of the interpolation method for the internal fluid-structure coupling interface of the engine in an embodiment;

[0021] Figure 2 It is a schematic framework diagram of the interpolation method for the internal fluid-structure coupling interface of the engine in an embodiment;

[0022] Figure 3 It is a schematic diagram of the grouping of interpolation surface nodes in an embodiment, where Figure 3 (a) is a 3D stereoscopic view, Figure 3 (b) is a 2D top view;

[0023] Figure 4 It is a diagram of the interpolation result of the simple fluid-structure interface in an embodiment, where Figure 4 (a) is the original interpolation surface and the distribution diagram of physical quantities p of Figure 4 (b) is the distribution diagram of physical quantities after interpolation, Figure 4 (c) is the distribution diagram of the relative error between the exact value and the interpolation result, Figure 4 (d) is the scatter diagram of the relative error;

[0024] Figure 5 It is a diagram of the interpolation result of the fluid-structure interface of the complex 1 / 4 configuration axisymmetric engine in an embodiment, where Figure 5 (a) is the heat flux distribution on the engine surface obtained by computational fluid dynamics (CFD) calculation, Figure 5 (b) is the heat flux distribution on the solid surface interpolated by the method proposed in the present invention;

[0025] Figure 6 It is a structural block diagram of the interpolation device for the internal fluid-structure coupling interface of the engine in an embodiment;

[0026] Figure 7 It is an internal structure diagram of a computer device in an embodiment.

[0027] The realization of the object, functional features, and advantages of the present invention will be further described with reference to the embodiments and the accompanying drawings. Specific Embodiments

[0028] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the scope of protection of the present invention.

[0029] 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 ability of those of ordinary skill in the art to implement. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0030] Next, the embodiments of the present invention will be described in detail in conjunction with the accompanying drawings in the embodiments of the present invention.

[0031] Embodiment 1

[0032] This embodiment discloses a method for interpolating the fluid-structure interaction interface in an engine. By grouping the interpolation surface nodes and distributing each group to different processes, the complexity of solving the subsequent system of linear algebraic equations is approximately reduced from to , significantly improving the calculation efficiency, where represents the number of boundary nodes. By adaptively setting the support radius of the subset of interpolation surface nodes, it is ensured that the interpolation function has high accuracy in the local area, while avoiding the positive definiteness problem of the linear system caused by too large a support radius. When performing interpolation, by determining whether the surface node to be interpolated is in a single-group influence area or a double-group influence area, different interpolation calculation methods are adopted, reducing unnecessary calculation overhead, further improving the interpolation efficiency while ensuring the interpolation accuracy; in addition, in the double-group influence area, adjacent two interpolation functions are weighted for interpolation, ensuring a smooth transition of the interpolation result at the group boundary. The present invention takes into account both accuracy and efficiency, is particularly suitable for large-scale complex configurations, and can ensure accurate data interaction at the large-scale fluid-structure interaction interface.

[0033] As shown in Figure 1 and Figure 2 , the method for interpolating the fluid-structure interaction interface in an engine provided in this embodiment includes the following steps:

[0034] Step 201, generate the fluid domain grid and the solid domain grid of the engine, and obtain the coordinates of the interpolation surface nodes in the fluid domain and the coordinates of the surface nodes to be interpolated in the solid domain.

[0035] Step 202: Group the interpolated surface nodes according to the longest direction of the fluid domain surface nodes to obtain several subsets of interpolated surface nodes; adopt a parallel computing strategy, allocate each subset of interpolated surface nodes to different processes, and evenly distribute the surface nodes to be interpolated among the processes.

[0036] Step 203: In each process, adaptively set the support radius of the subset of interpolated surface nodes, then construct a system of linear algebraic equations according to the radial basis interpolation function, and solve this system of equations to obtain the radial basis interpolation function based on the subset of interpolated surface nodes.

[0037] Step 204: Broadcast the radial basis interpolation functions in each process to all processes, and then use the radial basis interpolation functions of the subsets of interpolated surface nodes to interpolate the physical field of the surface nodes 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, calculate the interpolation result using a single-group interpolation function; if it is a double-group influence area, obtain the interpolation result by weighting two adjacent radial basis interpolation functions.

[0038] In the specific implementation process of Step 201, generate the mesh of the engine fluid domain through CFD, and output the coordinates and heat flux density of the fluid domain surface nodes; at the same time, calculate the engine solid domain mesh by computational solid dynamics (CSD) and output the coordinates of the solid wall nodes. Then, read the wall node coordinates and heat flux density through the fluid-structure interaction interpolation program.

[0039] Generally, the radial basis function used for fluid-structure interaction interpolation is Wendland’s C 2 function, and Wendland’s C 2 function is a kind of compact function, and its specific form is:

[0040] (1)

[0041] In the formula, , represents the position vector of the point to be interpolated, represents the position vector of the th control point, represents the support radius.

[0042] 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 in 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 a parallel architecture generally includes child processes and a main process. In the processes of this embodiment, both the main process and the child processes participate in the calculation.

[0043] Specifically, the interpolation surface nodes are grouped according to the longest direction of the fluid domain surface nodes, obtaining several interpolation surface node subsets, including:

[0044] Setting the number of groups and the overlap coefficient between groups ; Usually set to an integer multiple of the number of parallel cores to achieve load balancing; the overlap coefficient between groups is set to a number between 0 and 1. Assuming that the longest direction of the fluid domain surface nodes is direction, then based on the overlap coefficient between groups , the fluid domain surface is evenly divided into equal-length parts in the direction, obtaining the interpolation surface node subset , and the interpolation surface node subset is expressed as .

[0045] Among them, based on the overlap coefficient between groups , the fluid domain surface is evenly divided into equal-length parts in the direction, and the grouping rule is:

[0046] Calculating the grouping determination integer and the grouping determination decimal :

[0047] (2)

[0048] (3)

[0049] In the formula, floor is the floor function; represents the coordinate of the node along the direction of the largest scale; is the minimum value of the coordinates of all interpolation surface nodes x ; is the x coordinate maximum value of all interpolation surface nodes.

[0050] Grouping the interpolation surface nodes, when , When, the interpolation surface nodes are located in the first subset of interpolation surface nodes.

[0051] When , When, the interpolation surface nodes are only located in the th subset of interpolation surface nodes;

[0052] When and When, the interpolation surface nodes are only located in the th subset of interpolation surface nodes;

[0053] When and When, the interpolation surface nodes are simultaneously located in the th and the th subsets of interpolation surface nodes;

[0054] When and When, the interpolation surface nodes are simultaneously located in the th and the th subsets of interpolation surface nodes.

[0055] Among them, , so as to determine which subset of interpolation surface nodes the interpolation surface nodes are specifically assigned to. The overlap coefficient between groups can be used to characterize the size of the overlapping area between grouped nodes. represents groups that do not overlap at all. represents that all groups overlap with each other, and the width of the grouped area is . Among them, . According to experience, it is preferable that takes the value of 0.4, and better results can be obtained. 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.

[0056] In the specific implementation process of step 203, the support radius of the subset of the interpolation surface nodes is adaptively set, and the expression is:

[0057] (4)

[0058] In the formula, is the minimum value of the coordinates of all interpolation surface nodes; is the coordinate maximum value of all interpolation surface nodes; represents an adjustable multiple; represents dividing the surface of the fluid domain according to The number of equal parts in terms of direction and length.

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

[0060] In each process, a radial basis interpolation function is constructed according to the subset of interpolation surface nodes assigned to itself, and the expression is as follows:

[0061] (5)

[0062] In order to solve the coefficients of the interpolation function , so that the interpolation function is exactly equal to its physical quantity at the interpolation surface nodes, a system of linear algebraic equations is constructed. The expression of the system of linear algebraic equations in the

[0063] (6)

[0064] In the formula, represents the physical field vector at the represents the weight coefficient vector of the th control point; represents the position vector; represents the th position vector of the control point; represents the influence of the control point on the control point , where represents the position vector of the control point , represents the position vector of the control point ; represents the unknown coefficient corresponding to the th control point in the th dimension, represents the th dimension value of the physical field at the control point . The control point and the control point can take a maximum of , being the number of nodes in the interpolation surface grouping.

[0065] Adopt the LDL T Cholesky decomposition method of the open-source Eigen library to solve the coefficients of the interpolation function, and obtain the radial basis interpolation function based on the subset of interpolation surface nodes.

[0066] In the specific implementation process of step 204, the radial basis interpolation functions in each process are broadcast to all processes, mainly broadcasting the weight coefficient vector of the control points and the position vector of the control points , and then substituting the coordinates of the surface nodes to be interpolated into in formula (5) for interpolation of the physical field. During the interpolation process, OpenMP is used to create multiple threads to further accelerate the process.

[0067] If the coordinates of the surface nodes to be interpolated are represented by , it is determined whether the surface nodes to be interpolated are in a single-group influence area or a double-group influence area. The specific method is as follows:

[0068] Obtain the coordinates of the surface nodes to be interpolated , and loop based on the subset of interpolation surface nodes . During the loop, the index is , looping from 1 in sequence to .

[0069] If is satisfied, it means that the surface nodes to be interpolated are affected by both and groups, and it is a double-group influence area, and the loop is directly exited.

[0070] If is not satisfied, it means that the surface nodes to be interpolated are in a single-group influence area. When it is in a single-group influence area, it is also necessary to determine whether it is affected by the first or the th subset of interpolation surface nodes. If not satisfied, is looped again between 1 and for judgment.

[0071] The single-group influence area is subdivided into three cases:

[0072] When , it is only affected by the first subset of interpolation surface nodes.

[0073] When , it is only affected by the th subset of interpolation surface nodes.

[0074] When , it is only affected by the subset of interpolation surface nodes with index , and at this time is between 1 and .

[0075] 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:

[0076] (7)

[0077] (8)

[0078] 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.

[0079] 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.

[0080] 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.

[0081] It is worth noting that although Figure 2 The framework diagram of the CPC has been developed parallel processes, but in practice the number of processes opened can be equal to or less than the subset of interpolation surface nodes , that is, each process processes multiple subsets of interpolation surface nodes. When the number of processes is equal to the subset of interpolation surface nodes When , each process is assigned a subset of interpolation surface nodes for calculation; when the number of processes is less than the subset of interpolation surface nodes When, two or more subsets of surface nodes are allocated to each process. Preferably, the subsets of surface nodes can be set as an integer multiple of the number of processes for easy average distribution, thus ensuring parallel load balancing.

[0082] In one embodiment, as Figure 4 shown, a simple fluid-solid interface interpolation result graph is provided to verify the effectiveness of the method proposed by the present invention.

[0083] Both the original interpolated surface point cloud and physical quantities are generated by functions. The coordinates of the interpolated surface nodes are uniformly generated 10 in the range from 0 to 1, and the coordinates are also uniformly generated 10 in the range from 0 to 1. The

[0084] (9)

[0085] The surface physical quantity p is generated by a function of the

[0086] (10)

[0087] The distribution of the original interpolated surface nodes and physical quantities p is as shown in Figure 4 (a).

[0088] The surface point cloud to be interpolated is generated in a similar manner, with only 30 points each used for the generation of the and coordinates. After interpolation using the method provided by the present invention, the distribution of the surface physical quantity p is as shown in 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.

[0089] To quantitatively compare the reliability of the interpolation results, the physical quantities of the surface nodes to be interpolated are calculated using formula (10) and regarded as the exact values , and the relative error is calculated by comparing with the interpolation results:

[0090] (11)

[0091] Figure 4 (c) and Figure 4 (d) respectively give the relative error distribution graph and the relative error scatter plot. It can be seen that the maximum error is less than 1.6%, and the errors of most points are less than 0.2%. This proves the effectiveness of the interpolation method provided by the present invention.

[0092] In one embodiment, as Figure 5As shown, a graph of the interpolation results at the fluid-structure interface of a complex 1 / 4 configuration axisymmetric engine is provided to verify the effectiveness of the method provided by the present invention on a large-scale configuration.

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

[0094] Although the steps in this embodiment Figure 1 are shown in sequence according to the indication of the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise clearly stated in this article, the execution of these steps has no strict order limit, and these steps can be executed in other orders. Moreover, Figure 1 at least a part of the steps in this embodiment 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 either, but can be executed alternately or in turn with at least a part of other steps or sub-steps or stages of other steps.

[0095] Embodiment 2

[0096] Based on the interpolation method at the fluid-structure coupling interface in the engine of Embodiment 1, this embodiment discloses an interpolation device at the fluid-structure coupling interface in the engine, as Figure 6 shown. The interpolation device at the fluid-structure 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, where:

[0097] The initialization module 401 is used to generate the fluid domain grid and the solid domain grid of the engine, and obtain the coordinates of the interpolation surface nodes in the fluid domain and the coordinates of the surface nodes to be interpolated in the solid domain.

[0098] The node grouping and parallel allocation module 402 is used to group the interpolation surface nodes according to the longest direction of the surface nodes in the fluid domain, obtaining several subsets of interpolation surface nodes; adopting a parallel computing strategy, distributing each subset of interpolation surface nodes to different processes, and evenly distributing the surface nodes to be interpolated to each process.

[0099] 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 system of linear algebraic equations according to the radial basis interpolation function, and solve the system of equations to obtain the radial basis interpolation function based on the interpolation surface node subset.

[0100] 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 to-be-interpolated surface nodes; during the interpolation process, it is judged whether the to-be-interpolated surface node 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.

[0101] In this embodiment, the specific working processes and working principles 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 of the method in Embodiment 1, so they will not be elaborated herein. Each of these unit modules can be implemented in whole or in part by software, hardware, and their combination. Each unit module can be embedded in the processor of the computer device in hardware form or be independent of it, or be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to each of these unit modules.

[0102] Embodiment 3

[0103] As Figure 7 shown, a terminal device disclosed in this embodiment includes a transmitter, a receiver, a memory, and a processor. Among them, 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 execution instructions, and the processor is used to execute the computer execution instructions stored in the memory to implement the method in Embodiment 1 above.

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

[0105] Embodiment 4

[0106] This embodiment discloses a computer-readable storage medium, in which computer execution instructions are stored. When the processor executes the computer execution instructions, the method in Embodiment 1 above is implemented.

[0107] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. 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 methods. Among them, any reference to a memory, storage, database, or other medium used in the embodiments provided in the present application can include non-volatile and / or volatile memories. Non-volatile memories can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memories can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0108] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, 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, it should be considered as the scope recorded in this specification.

[0109] The above-described embodiments merely represent several implementation manners of the present invention. The description thereof is relatively specific and detailed, but it cannot be construed as a limitation on the scope of the invention. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present invention, several modifications and improvements can still be made, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the appended 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; Broadcasting the radial basis interpolation function in each process to all processes, and then using 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, judging 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, using a single group interpolation function to calculate the interpolation result; if it is a double-group influence area, using two adjacent radial basis interpolation functions to weight the interpolation result; 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, where represents the overlap coefficient between groups; the interpolation surface node subset is expressed as .

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 Partition, get the interpolation surface node subset .

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 1, 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.

7. A fluid-solid coupling interface interpolation device in an engine, characterized in that: The method for interpolation of fluid-solid coupling interface in an engine according to any one of claims 1 to 6 is adopted, wherein 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.

8. 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 6 are implemented.

9. 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 6 are implemented.

Images