A fluid variable gradient calculation method, device, equipment and storage medium

Abstract

The present application discloses a fluid variable gradient calculation method, device, equipment and storage medium, which relates to the field of fluid mechanics technology, including: when iteratively discretizing the fluid mechanics control equation of the fluid region, obtaining the body-centered independent variable value of the non-structured grid unit in the current iteration round region; determining the face-centered independent variable value of the boundary surface between the non-structured grid units based on the body-centered independent variable value, and determining the independent variable gradient of the unit using the volume and face-centered independent variable value of the unit; determining the new face-centered independent variable value of the boundary surface using the face-centered position vector of the boundary surface, the body-centered independent variable value of the non-structured grid unit located on the left and right of the boundary surface, the independent variable gradient restriction coefficient and the independent variable gradient; determining the target grid independent variable gradient based on the volume and new face-centered independent variable value of the non-structured grid unit, so as to determine the variable distribution of the fluid region using the target grid independent variable gradient. The compactness of the gradient calculation is efficiently maintained by implicitly expanding the unit gradient calculation template.

Description

Technical Field

[0001] The present invention relates to the field of fluid mechanics technology, and in particular to a fluid variable gradient calculation method, device, equipment and storage medium. Background Art

[0002] At present, the fluid region is spatially discretized using structured grids or unstructured grids (see Figure 1 ), that is, by iteratively and discretely calculating the values ​​of independent variables such as density and pressure at the center of the unstructured grid unit, the variable distribution of the entire fluid area can be analyzed.

[0003] In the process of spatial discretization, one of the key steps is to calculate the independent variable value of the boundary surface, and the determination of the independent variable value of the center of the surface requires the gradient of the independent variable of the unstructured grid unit. However, the existing gradient calculation scheme has poor adaptability and low convergence, and will further have an adverse effect on the analysis of the independent variable value of the center of the surface and the variable distribution of the entire fluid region.

[0004] Therefore, how to efficiently and robustly calculate the independent variable gradients of unstructured grid cells is a problem that needs to be solved urgently. Summary of the invention

[0005] In view of this, the purpose of the present invention is to provide a fluid variable gradient calculation method, device, equipment and storage medium, which can efficiently maintain the compactness of gradient calculation through implicit expansion unit gradient calculation template, thereby avoiding the introduction of more additional calculations and additional parallel value transmission, and improving the calculation convergence and the ability to adapt to complex shapes. The specific scheme is as follows:

[0006] In a first aspect, the present application provides a fluid variable gradient calculation method, comprising:

[0007] In the process of iterative discretization of the fluid mechanics governing equations of the fluid region, obtaining the body-centered independent variable value of the unstructured grid cell in the fluid region in the current iteration round;

[0008] Determine the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determine the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value;

[0009] Determine a new face center independent variable value corresponding to the boundary surface by using the face center position vector of the boundary surface, the body center independent variable value of the unstructured grid unit located on the left and right of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient;

[0010] The target grid independent variable gradient corresponding to each of the non-structured grid cells is determined based on the volume of the non-structured grid cell and the corresponding new face-centre independent variable value, so as to determine the variable distribution of the fluid region by using the target grid independent variable gradient.

[0011] Optionally, obtaining the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round includes:

[0012] The body-centered independent variable values ​​of the unstructured grid cells in the fluid region in the current iteration round are obtained; the body-centered independent variable values ​​include density values, velocity values, and pressure values.

[0013] Optionally, determining the face-centered independent variable value of the boundary surface between the unstructured grid units based on the body-centered independent variable value includes:

[0014] For any boundary surface between the unstructured grid units, the face-centred independent variable value of the current boundary surface is determined by taking an arithmetic average based on the body-centred independent variable values ​​of the unstructured grid units located on the left and right sides of the current boundary surface.

[0015] Optionally, determining the independent variable gradient of each of the unstructured grid cells by using the volume of the unstructured grid cells and the surface normal vector of the corresponding boundary surface and the surface center independent variable value includes:

[0016] For any of the non-structured grid cells, obtaining the body-centered independent variable value of the non-structured grid cell adjacent to the current non-structured grid cell to obtain an adjacent grid independent variable value;

[0017] The independent variable gradient of the current unstructured grid cell is determined based on the adjacent grid independent variable values, the volume of the current unstructured grid cell, the surface normal vector of the boundary surface adjacent to the current unstructured grid cell, and the surface center independent variable value.

[0018] Optionally, determining the target grid independent variable gradient corresponding to each of the unstructured grid cells based on the volume of the unstructured grid cell and the corresponding new face center independent variable value comprises:

[0019] For any of the unstructured grid cells, a target grid independent variable gradient of the current unstructured grid cell is determined based on the volume of the current unstructured grid cell, the adjacent grid independent variable values, the surface normal vector of the boundary surface adjacent to the current unstructured grid cell, and the new face center independent variable value.

[0020] Optionally, the determining of a new face center independent variable value corresponding to the boundary surface by using the face center position vector of the boundary surface, the body center independent variable value of the unstructured grid unit located on the left and right sides of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient comprises:

[0021] For any of the boundary surfaces, obtaining the body-centered independent variable values ​​and the independent variable gradients of the unstructured grid cells located on the left and right sides of the current boundary surface;

[0022] Obtaining the independent variable gradient restriction coefficient and the body center position vector of the unstructured grid unit located on the left and right sides of the current boundary surface;

[0023] A new face center independent variable value of the current boundary surface is determined based on the face center position vector of the current boundary surface, the body center independent variable values ​​of the unstructured grid units located on the left and right of the current boundary surface, the independent variable gradient, the independent variable gradient restriction coefficient and the body center position vector.

[0024] Optionally, after determining the target grid independent variable gradient corresponding to each of the unstructured grid cells based on the volume of the unstructured grid cell and the corresponding new face center independent variable value, the method further includes:

[0025] For any of the boundary surfaces, interpolation is performed based on the target grid independent variable gradient, the body center independent variable value, the body center position vector and the face center position vector of the current boundary surface of the unstructured grid cells located on the left and right of the current boundary surface to determine the target face center independent variable value of the current boundary surface.

[0026] In a second aspect, the present application provides a fluid variable gradient calculation device, comprising:

[0027] An independent variable value acquisition module is used to obtain the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round during the iterative discretization of the fluid mechanics control equation of the fluid region;

[0028] A gradient calculation module, used to determine the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determine the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value;

[0029] An independent variable value determination module, used to determine a new face-center independent variable value corresponding to the boundary surface by using the face-center position vector of the boundary surface, the body-center independent variable values ​​of the unstructured grid cells located on the left and right sides of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient;

[0030] The target gradient determination module is used to determine the target grid independent variable gradient corresponding to each of the non-structured grid cells based on the volume of the non-structured grid cells and the corresponding new face-center independent variable value, so as to determine the variable distribution of the fluid region using the target grid independent variable gradient.

[0031] In a third aspect, the present application provides an electronic device, including:

[0032] Memory, used to store computer programs;

[0033] The processor is used to execute the computer program to implement the steps of the aforementioned fluid variable gradient calculation method.

[0034] In a fourth aspect, the present application provides a computer-readable storage medium for storing a computer program, which implements the steps of the aforementioned fluid variable gradient calculation method when executed by a processor.

[0035] It can be seen that in the present application, in the process of iterative discretization of the fluid mechanics control equations of the fluid region, the body-centered independent variable values ​​of the unstructured grid cells in the fluid region in the current iteration round are obtained; the face-centered independent variable values ​​of the boundary surfaces between the unstructured grid cells are determined based on the body-centered independent variable values, and the independent variable gradients of each of the unstructured grid cells are determined using the volume of the unstructured grid cells and the corresponding surface normal vector of the boundary surface and the face-centered independent variable values; the new face-centered independent variable values ​​corresponding to the boundary surface are determined using the face-center position vector of the boundary surface, the body-centered independent variable values ​​of the unstructured grid cells located on the left and right of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient; the target grid independent variable gradients corresponding to each of the unstructured grid cells are determined based on the volume of the unstructured grid cells and the corresponding new face-center independent variable values, so as to determine the variable distribution of the fluid region using the target grid independent variable gradients. That is, the present application first uses the body-centered independent variable value and the face-centered independent variable value to determine the independent variable gradient of the unstructured grid unit, and then determines the new face-centered independent variable value based on the body-centered independent variable value, the independent variable gradient restriction coefficient and the independent variable value gradient on the left and right sides of the boundary surface, and uses the new face-centered independent variable value to calculate the body-centered independent variable gradient of the unstructured grid unit again to obtain the corresponding target grid independent variable gradient. In this way, the compactness of the gradient calculation can be efficiently maintained by implicitly expanding the unit gradient calculation template, thereby avoiding the introduction of more additional calculations and additional parallel value transmission, and improving the calculation convergence and the ability to adapt to complex shapes. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] 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 embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0037] Figure 1 A schematic diagram of using an unstructured grid to discretize the entire fluid region provided for this application;

[0038] Figure 2 A flow chart of a fluid variable gradient calculation method provided in this application;

[0039] Figure 3 A specific flow numerical calculation schematic diagram provided for this application;

[0040] Figure 4 A schematic diagram of a template for calculating the independent variable value of a boundary surface center provided in this application;

[0041] Figure 5 A schematic diagram of an independent variable gradient calculation template for a C0 unit provided in this application;

[0042] Figure 6 A schematic diagram of the structure of a fluid variable gradient calculation device provided in this application;

[0043] Figure 7 A structural diagram of an electronic device provided for this application. DETAILED DESCRIPTION

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

[0045] In the process of spatial discretization, one of the key steps is to calculate the independent variable value of the face center of the boundary surface, and the determination of the face center independent variable value requires the gradient of the independent variable of the unstructured grid unit. However, the existing gradient calculation scheme has poor adaptability and low convergence, and will further have an adverse effect on the analysis of the face center independent variable value and the variable distribution of the entire fluid region. To this end, the present application provides a fluid variable gradient calculation device solution, which can efficiently maintain the compactness of the gradient calculation through an implicit expansion unit gradient calculation template, thereby avoiding the introduction of more additional calculations and additional parallel value transmission, and improving the calculation convergence and the ability to adapt to complex shapes.

[0046] See also Figure 2 As shown, an embodiment of the present invention discloses a fluid variable gradient calculation device method, comprising:

[0047] Step S11, in the process of iterative discretization of the fluid mechanics governing equations of the fluid region, obtaining the body-centered independent variable values ​​of the unstructured grid cells in the fluid region in the current iteration round.

[0048] In this embodiment, the unstructured second-order finite volume method can be used for iterative discretization to discretize the control equation into a set of algebraic equations. The independent variable value of the center of the unstructured grid unit at a certain moment is obtained by an iteration step. Figure 3 As shown, obtain The value of the body-centered independent variable of the unstructured grid cell.

[0049] also, Figure 3 In represents the boundary surface between two mesh cells, for example, Represents a grid cell The boundary surface between Represents a grid cell The boundary surface between them, and so on. That is, the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round is obtained; the body-centered independent variable value includes but is not limited to the density value, the velocity value and the pressure value.

[0050] Step S12: determining the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determining the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value.

[0051] In this embodiment, regarding the gradient calculation, first, for the boundary surface between any of the unstructured grid cells, the arithmetic average is performed based on the body-centered independent variable values ​​of the unstructured grid cells located on the left and right of the current boundary surface to determine the face-centered independent variable value of the current boundary surface. Then, for any of the unstructured grid cells, the body-centered independent variable value of the unstructured grid cell adjacent to the current unstructured grid cell is obtained to obtain the adjacent grid independent variable value; the independent variable gradient of the current unstructured grid cell is determined based on the adjacent grid independent variable value, the volume of the current unstructured grid cell, and the face normal vector of the boundary surface of the current unstructured grid cell and the face-centered independent variable value. Combined with Figure 4 As shown below, The calculation process of the independent variable gradient is introduced by taking the pressure gradient calculation of the unit as an example. First, calculate Boundary surfaces around mesh cells The face center pressure value can be calculated by taking the arithmetic mean of the body center pressure values ​​of the units on the left and right sides of the boundary surface. The specific expression can be as follows:

[0052] ;

[0053] In the formula, for The pressure value of the grid cell; for The pressure value of the grid cell; for The pressure value of the grid cell; For the boundary surface The pressure value; For the boundary surface The pressure value; For the boundary surface pressure value.

[0054] Then, combine Figure 5 As shown, based on the adoption A mesh cell and its adjacent mesh cells , , The body center pressure is calculated The pressure gradient value of the unit ,in , , is the adjacent boundary surface of C0 grid cell The surface normal vector of the face center (determined by multiplying the area of ​​the boundary surface and the face center normal vector), for The volume of the grid cell, its independent variable gradient The calculation formula can be shown as follows:

[0055] .

[0056] Step S13, using the center position vector of the boundary surface, the center independent variable values ​​of the unstructured grid units located on the left and right of the boundary surface, the independent variable gradient restriction coefficient and the independent variable gradient to determine a new center independent variable value corresponding to the boundary surface.

[0057] It should be understood that after determining the independent variable gradient, the body-centered independent variable value of the units on the left and right sides of the non-structured unit boundary surface and the gradient of the body-centered independent variable value are used to obtain the new independent variable value of the boundary surface center through arithmetic averaging. That is, for any of the boundary surfaces, the body-centered independent variable value and the independent variable gradient of the non-structured grid units located on the left and right sides of the current boundary surface are obtained; the independent variable gradient restriction coefficient and the body-centered position vector of the non-structured grid units located on the left and right sides of the current boundary surface are obtained; based on the center position vector of the current boundary surface, the body-centered independent variable value of the non-structured grid units located on the left and right sides of the current boundary surface, the independent variable gradient, the independent variable gradient restriction coefficient and the body-centered position vector, the new center independent variable value of the current boundary surface is determined. It can be understood that the independent variable gradient restriction coefficient includes but is not limited to the density gradient restriction coefficient, the velocity gradient restriction coefficient and the pressure gradient restriction coefficient.

[0058] Based on the calculation in step S12 The pressure gradient of the grid cell is calculated again Adjacent boundary surfaces of mesh cells The center pressure value of:

[0059] ;

[0060] In the formula, , , Boundary surface The face center position vector, for The center position vector of the grid cell. The grid cells are , , The grid cells are surrounded by a state where there are no sources or sinks. The center pressure value of the adjacent boundary surface of the grid unit cannot exceed , , The center pressure value of the grid unit is obtained based on this principle Pressure gradient limit coefficient of single grid element . Use limiter to limit Boundary surfaces of mesh cells The center pressure value of the surface, and considering that the boundary surface can be respectively from the left ( ) and on the right ( ) interpolation, we can get the following expression:

[0061] ;

[0062] In the formula, , They are , The pressure gradient of the grid cell, , They are , Pressure gradient limit factor for grid cells.

[0063] Afterwards, in order to improve the convergence and efficiency of the calculation, the calculation template is expanded to correct The independent variable gradient of the grid cell. Specifically, the arithmetic mean of the center pressure value and the pressure gradient value of the cells on the left and right sides of the boundary surface is used to recalculate the boundary surface. The face center pressure value is used to obtain the new face center pressure value , , .

[0064] .

[0065] Step S14: determining target grid independent variable gradients corresponding to each of the non-structured grid cells based on the volumes of the non-structured grid cells and the corresponding new face-centre independent variable values, so as to determine the variable distribution of the fluid region using the target grid independent variable gradients.

[0066] Specifically, in this embodiment, after determining the new face center independent variable value of the boundary surface, the independent variable gradient of the unstructured grid unit is recalculated. That is, for any of the unstructured grid units, the target grid independent variable gradient of the current unstructured grid unit is determined based on the volume of the current unstructured grid unit, the adjacent grid independent variable value, the surface normal vector of the boundary surface adjacent to the current unstructured grid unit, and the new face center independent variable value. Taking the grid unit as an example, the new center pressure value obtained in step S13 is recalculated. Pressure gradient across a grid cell:

[0067] ;

[0068] Combination Figure 3 As shown, it can be seen that this embodiment calculates Unstructured grid cells are used when calculating the pressure gradient of a grid cell. , , , .calculate Unstructured grid cells are used when calculating the pressure gradient of a grid cell. , , , .calculate Unstructured grid cells are used when the pressure gradient of the cell is , , , ,calculate Unstructured grid cells are used when the pressure gradient of the cell is , , , That is, The grid cell pressure gradient calculation template is composed of unstructured grid cells , , , Extension to unstructured grid cells , , , , , , , , , ,Template expansion can effectively enhance its adaptability to complex shapes and improve the convergence of ,computation.

[0069] Furthermore, in this embodiment, after determining the target grid independent variable gradient corresponding to each of the unstructured grid cells based on the volume of the unstructured grid cells and the corresponding new face-center independent variable value, for any of the boundary surfaces, interpolation is performed based on the target grid independent variable gradient, the body-center independent variable value, the body-center position vector and the face-center position vector of the current boundary surface of the unstructured grid cells located on the left and right of the current boundary surface to determine the target face-center independent variable value of the current boundary surface.

[0070] In summary, this embodiment uses the body-centered independent variable values ​​of the left and right units of the non-structured unit boundary surface and the gradient of the body-centered independent variable values ​​on the basis of calculating the variable gradient of the non-structured grid unit by arithmetic averaging, obtains the new independent variable value of the boundary surface center by arithmetic averaging, and uses the new independent variable value of the boundary surface center to recalculate the body-centered independent variable gradient of the non-structured grid unit, so as to introduce more adjacent calculation unit variable influences in the gradient calculation, implicitly expand the unit gradient calculation template in this process, maintain the compactness of the gradient calculation, and thus avoid introducing more additional calculations and additional parallel value transmission. The gradient calculated based on the new method has better calculation convergence and the ability to adapt to complex shapes.

[0071] It can be seen that in the embodiment of the present application, in the process of iterative discretization of the fluid mechanics control equations of the fluid region, the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round is obtained; the face-centered independent variable value of the boundary surface between the unstructured grid units is determined based on the body-centered independent variable value, and the independent variable gradient of each unstructured grid unit is determined using the volume of the unstructured grid unit and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value; the new face-centered independent variable value corresponding to the boundary surface is determined using the face-center position vector of the boundary surface, the body-centered independent variable value of the unstructured grid unit located on the left and right of the boundary surface, the independent variable gradient restriction coefficient and the independent variable gradient; based on the volume of the unstructured grid unit and the corresponding new face-center independent variable value, the target grid independent variable gradient corresponding to each of the unstructured grid units is determined, so as to determine the variable distribution of the fluid region using the target grid independent variable gradient. That is, the present application first uses the body-centered independent variable value and the face-centered independent variable value to determine the independent variable gradient of the unstructured grid unit, and then determines the new face-centered independent variable value based on the body-centered independent variable value, the independent variable gradient restriction coefficient and the independent variable value gradient on the left and right sides of the boundary surface, and uses the new face-centered independent variable value to calculate the body-centered independent variable gradient of the unstructured grid unit again to obtain the corresponding target grid independent variable gradient. In this way, the compactness of the gradient calculation can be efficiently maintained by implicitly expanding the unit gradient calculation template, thereby avoiding the introduction of more additional calculations and additional parallel value transmission, and improving the calculation convergence and the ability to adapt to complex shapes.

[0072] See also Figure 6 As shown, the embodiment of the present application also discloses a fluid variable gradient calculation device, including:

[0073] An independent variable value acquisition module 11 is used to obtain the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round during the iterative discretization of the fluid mechanics control equation of the fluid region;

[0074] A gradient calculation module 12, for determining the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determining the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding face normal vector of the boundary surface and the face-centered independent variable value;

[0075] An independent variable value determination module 13 is used to determine a new face center independent variable value corresponding to the boundary surface by using the face center position vector of the boundary surface, the body center independent variable value of the unstructured grid unit located on the left and right of the boundary surface, the independent variable gradient restriction coefficient and the independent variable gradient;

[0076] The target gradient determination module 14 is used to determine the target grid independent variable gradient corresponding to each of the non-structured grid cells based on the volume of the non-structured grid cells and the corresponding new face center independent variable value, so as to determine the variable distribution of the fluid region using the target grid independent variable gradient.

[0077] Among them, for more specific working processes of the above-mentioned modules, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be repeated here.

[0078] It can be seen that the present application first uses the body-centered independent variable value and the face-centered independent variable value to determine the independent variable gradient of the unstructured grid unit, and then determines the new face-centered independent variable value based on the body-centered independent variable value, the independent variable gradient restriction coefficient and the independent variable value gradient on the left and right sides of the boundary surface, and uses the new face-centered independent variable value to calculate the body-centered independent variable gradient of the unstructured grid unit again to obtain the corresponding target grid independent variable gradient. In this way, the compactness of the gradient calculation can be efficiently maintained by implicitly expanding the unit gradient calculation template, thereby avoiding the introduction of more additional calculations and additional parallel value transmission, and improving the calculation convergence and the ability to adapt to complex shapes.

[0079] In some specific embodiments, the independent variable value acquisition module 11 may specifically include:

[0080] The body-centered independent variable acquisition unit is used to obtain the body-centered independent variable values ​​of the unstructured grid units in the fluid area in the current iteration round; the body-centered independent variable values ​​include density values, velocity values ​​and pressure values.

[0081] In some specific embodiments, the gradient calculation module 12 may specifically include:

[0082] The face-centre independent variable value acquisition unit is used to determine the face-centre independent variable value of the current boundary surface based on the arithmetic average of the body-centre independent variable values ​​of the non-structured grid cells located on the left and right of the current boundary surface for any boundary surface between the non-structured grid cells.

[0083] In some specific embodiments, the gradient calculation module 12 may specifically include:

[0084] An adjacent grid independent variable value obtaining unit is used to obtain the body-centered independent variable value of the non-structured grid cell adjacent to the current non-structured grid cell for any of the non-structured grid cells, so as to obtain an adjacent grid independent variable value;

[0085] The gradient calculation unit is used to determine the independent variable gradient of the current unstructured grid unit based on the adjacent grid independent variable value, the volume of the current unstructured grid unit, and the surface normal vector of the boundary surface adjacent to the current unstructured grid unit and the face center independent variable value.

[0086] In some specific embodiments, the target gradient determination module 14 may specifically include:

[0087] The target gradient determination unit is used to determine the target grid independent variable gradient of the current unstructured grid cell for any of the unstructured grid cells based on the volume of the current unstructured grid cell, the adjacent grid independent variable values, the surface normal vector of the boundary surface adjacent to the current unstructured grid cell, and the new face center independent variable value.

[0088] In some specific embodiments, the independent variable value determination module 13 can be specifically used to obtain the body-centered independent variable value and the independent variable gradient of the non-structured grid cells located on the left and right of the current boundary surface for any of the boundary surfaces; obtain the independent variable gradient restriction coefficient and the body-centered position vector of the non-structured grid cells located on the left and right of the current boundary surface; and determine the new face-centered independent variable value of the current boundary surface based on the face-centered position vector of the current boundary surface, the body-centered independent variable value of the non-structured grid cells located on the left and right of the current boundary surface, the independent variable gradient, the independent variable gradient restriction coefficient and the body-centered position vector.

[0089] In some specific embodiments, the fluid variable gradient calculation device can also be used to interpolate, for any of the boundary surfaces, based on the target grid independent variable gradient, the body center independent variable value, the body center position vector and the face center position vector of the current boundary surface of the unstructured grid units located on the left and right of the current boundary surface to determine the target face center independent variable value of the current boundary surface.

[0090] Furthermore, the present application also discloses an electronic device. Figure 7 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram cannot be regarded as any limitation on the scope of use of the present application.

[0091] Figure 7 A schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the fluid variable gradient calculation method disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0092] In this embodiment, the power supply 23 is used to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device, and the communication protocol it follows is any communication protocol that can be applied to the technical solution of the present application, and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0093] In addition, the memory 22, as a carrier for storing resources, can be a read-only memory, a random access memory, a disk or an optical disk, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0094] The operating system 221 is used to manage and control the hardware devices and computer program 222 on the electronic device 20, and can be Windows Server, Netware, Unix, Linux, etc. In addition to the computer program that can be used to complete the fluid variable gradient calculation method performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include computer programs that can be used to complete other specific tasks.

[0095] Furthermore, the present application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the fluid variable gradient calculation method disclosed above is implemented. The specific steps of the method can refer to the corresponding contents disclosed in the above embodiments, and will not be repeated here.

[0096] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0097] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0098] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0099] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0100] The technical solution provided by the present application is introduced in detail above. Specific examples are used in this article to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for general technicians in this field, according to the idea of ​​the present application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A fluid variable gradient calculation method, characterized in that: include: In the process of iterative discretization of the fluid mechanics governing equations of the fluid region, obtaining the body-centered independent variable value of the unstructured grid cell in the fluid region in the current iteration round; Determine the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determine the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value; Determine a new face center independent variable value corresponding to the boundary surface by using the face center position vector of the boundary surface, the body center independent variable value of the unstructured grid unit located on the left and right of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient; The target grid independent variable gradient corresponding to each of the non-structured grid cells is determined based on the volume of the non-structured grid cell and the corresponding new face-centre independent variable value, so as to determine the variable distribution of the fluid region by using the target grid independent variable gradient.

2. The fluid variable gradient calculation method according to claim 1, characterized in that: The step of obtaining the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round includes: The body-centered independent variable values ​​of the unstructured grid cells in the fluid region in the current iteration round are obtained; the body-centered independent variable values ​​include density values, velocity values, and pressure values.

3. The fluid variable gradient calculation method according to claim 1, characterized in that: The step of determining the face-centered independent variable value of the boundary surface between the unstructured grid units based on the body-centered independent variable value comprises: For any boundary surface between the unstructured grid units, the face-centred independent variable value of the current boundary surface is determined by taking an arithmetic average based on the body-centred independent variable values ​​of the unstructured grid units located on the left and right sides of the current boundary surface.

4. The fluid variable gradient calculation method according to claim 1, characterized in that: The method of determining the independent variable gradient of each of the unstructured grid cells by using the volume of the unstructured grid cells and the surface normal vector of the corresponding boundary surface and the surface center independent variable value comprises: For any of the non-structured grid cells, obtaining the body-centered independent variable value of the non-structured grid cell adjacent to the current non-structured grid cell to obtain an adjacent grid independent variable value; The independent variable gradient of the current unstructured grid cell is determined based on the adjacent grid independent variable values, the volume of the current unstructured grid cell, the surface normal vector of the boundary surface of the current unstructured grid cell, and the surface center independent variable value.

5. The fluid variable gradient calculation method according to claim 4, characterized in that: The step of determining the target grid independent variable gradient corresponding to each of the non-structured grid cells based on the volume of the non-structured grid cell and the corresponding new face center independent variable value comprises: For any of the unstructured grid cells, a target grid independent variable gradient of the current unstructured grid cell is determined based on the volume of the current unstructured grid cell, the adjacent grid independent variable values, the surface normal vector of the boundary surface adjacent to the current unstructured grid cell, and the new face center independent variable value.

6. The fluid variable gradient calculation method according to any one of claims 1 to 5, characterized in that: The method of determining a new face center independent variable value corresponding to the boundary surface by using the face center position vector of the boundary surface, the body center independent variable value of the unstructured grid unit located on the left and right sides of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient comprises: For any of the boundary surfaces, obtaining the body-centered independent variable values ​​and the independent variable gradients of the unstructured grid cells located on the left and right sides of the current boundary surface; Obtaining the independent variable gradient restriction coefficient and the body center position vector of the unstructured grid unit located on the left and right sides of the current boundary surface; A new face center independent variable value of the current boundary surface is determined based on the face center position vector of the current boundary surface, the body center independent variable values ​​of the unstructured grid units located on the left and right of the current boundary surface, the independent variable gradient, the independent variable gradient restriction coefficient and the body center position vector.

7. The fluid variable gradient calculation method according to claim 1, characterized in that: After determining the target grid independent variable gradients corresponding to each of the non-structured grid cells based on the volume of the non-structured grid cells and the corresponding new face center independent variable values, the method further includes: For any of the boundary surfaces, interpolation is performed based on the target grid independent variable gradient, the body center independent variable value, the body center position vector and the face center position vector of the current boundary surface of the unstructured grid cells located on the left and right of the current boundary surface to determine the target face center independent variable value of the current boundary surface.

8. A fluid variable gradient calculation device, characterized in that: include: An independent variable value acquisition module is used to obtain the body-centered independent variable value of the unstructured grid unit in the fluid region in the current iteration round during the iterative discretization of the fluid mechanics control equation of the fluid region; A gradient calculation module, used to determine the face-centered independent variable value of the boundary surface between the unstructured grid cells based on the body-centered independent variable value, and determine the independent variable gradient of each of the unstructured grid cells using the volume of the unstructured grid cell and the corresponding surface normal vector of the boundary surface and the face-centered independent variable value; An independent variable value determination module, used to determine a new face-center independent variable value corresponding to the boundary surface by using the face-center position vector of the boundary surface, the body-center independent variable values ​​of the unstructured grid cells located on the left and right sides of the boundary surface, the independent variable gradient restriction coefficient, and the independent variable gradient; The target gradient determination module is used to determine the target grid independent variable gradient corresponding to each of the non-structured grid cells based on the volume of the non-structured grid cells and the corresponding new face-center independent variable value, so as to determine the variable distribution of the fluid region using the target grid independent variable gradient.

9. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the fluid variable gradient calculation method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: Used to store a computer program, which, when executed by a processor, implements the fluid variable gradient calculation method as described in any one of claims 1 to 7.

Images