CN115906709A - Fluid pressure field prediction method and device - Google Patents

Fluid pressure field prediction method and device Download PDF

Info

Publication number
CN115906709A
CN115906709A CN202211670441.3A CN202211670441A CN115906709A CN 115906709 A CN115906709 A CN 115906709A CN 202211670441 A CN202211670441 A CN 202211670441A CN 115906709 A CN115906709 A CN 115906709A
Authority
CN
China
Prior art keywords
grid
target
boundary
target area
determining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211670441.3A
Other languages
Chinese (zh)
Inventor
段志伟
许伟程
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Yundao Zhizao Technology Co ltd
Original Assignee
Beijing Yundao Zhizao Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Yundao Zhizao Technology Co ltd filed Critical Beijing Yundao Zhizao Technology Co ltd
Priority to CN202211670441.3A priority Critical patent/CN115906709A/en
Publication of CN115906709A publication Critical patent/CN115906709A/en
Pending legal-status Critical Current

Links

Images

Abstract

The application provides a method and a device for predicting a fluid pressure field, wherein the method comprises the following steps: acquiring a target area file containing the boundary type and the label of each grid in the target area; determining a plurality of connected domains in the target area according to the boundary type between adjacent grids; the boundary type is used for describing whether the adjacent grids are connected or not; determining the reference number and the reference pressure of one reference point of each connected domain; and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target region. According to the method and the device, the boundary type of each grid in the target area file is identified, the multiple connected domains of the target area are determined according to the boundary types of the adjacent grids, and the predicted pressure of the target area is determined according to the reference point and the reference pressure of each connected domain, so that the technical problem that the predicted pressure is inaccurate or cannot be predicted due to the fact that the target area containing the multiple connected domains can only depend on one reference point and reference pressure in the prior art is solved.

Description

Fluid pressure field prediction method and device
Technical Field
The present disclosure relates to the field of fluid mechanics, and more particularly, to a method and an apparatus for predicting a fluid pressure field.
Background
In actual engineering, predicting the pressure of a fluid typically solves partial differential equations by finite volume methods. The principle of the finite volume method is to discretize the prediction region into a finite number of "grid cells", each of which has a "finite small" volume and is interconnected. In CFD (predictive Fluid Dynamics) software, it is common to randomly select a reference point inside a prediction region where "grid cells" are connected to each other, and to specify a reference pressure of the reference point. If a limited number of grid cells in the prediction area are not all connected with each other, but a plurality of grid cells exist in the prediction area and respectively form different 'connected domains', the fluid pressure of the target area cannot be predicted through one reference point and the reference pressure.
Disclosure of Invention
In view of the above, an object of the present application is to provide at least a prediction of a fluid pressure field, which determines a plurality of connected domains of a target area according to boundary types of adjacent grids by identifying boundary types of each grid in a target area file, so as to determine a predicted pressure of the target area according to a reference point and a reference pressure of each connected domain, thereby solving a technical problem in the prior art that a predicted pressure determined for a target area including a plurality of connected domains only according to one reference point and reference pressure is inaccurate, and achieving a technical effect of improving accuracy of determining the predicted pressure.
The application mainly comprises the following aspects:
in a first aspect, an embodiment of the present application provides a method for predicting a fluid pressure field, where the method includes: acquiring a target area file containing the boundary type and the label of each grid in the target area; determining a plurality of connected domains in the target area according to the boundary type between adjacent grids; the boundary type is used for describing whether adjacent grids are connected or not; determining the label and reference pressure of a reference point of each connected domain; and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area.
Optionally, determining a plurality of connected domains in the target region according to the boundary type between the adjacent grids includes: randomly selecting one grid from all grids of the target area file as a first grid; determining a target connected domain containing the first grid and marking the target connected domain according to the boundary type of the first grid; determining whether the number of grids of the marked target connected domain is the same as the number of grids of the target area; if the number of the grids of the marked target connected domain is different from the number of the grids of the target area, randomly selecting a second grid which does not belong to the target connected domain in the target area; taking the second grid as a new first grid, jumping to a boundary type according to the first grid, determining a target connected domain containing the first grid, marking the target connected domain and continuing to execute the process until the grid number of all marked target connected domains is the same as that of the target area; all marked target connected domains are taken as a plurality of connected domains in the target area.
Optionally, the boundary type includes connected and disconnected, and determining a target connected domain including the first mesh and marking the target connected domain according to the boundary type of the first mesh includes: taking the first grid as a target grid; determining whether a boundary of which the boundary type is connected exists in the target grid; and if the boundary type of the target grid is not a connected boundary, taking the target grid as a target connected domain containing the first grid and marking the target connected domain.
Optionally, after determining whether the target mesh has a boundary with a connected boundary type, the method further includes: if the boundary type of the target grid is a connected boundary, determining the boundary with the connected boundary type in the target grid as a first boundary; determining a third grid corresponding to the first boundary except the target grid; determining whether the third grid is a present target grid; if the third grid is not the appeared target grid, the third grid is taken as a new target grid, and the third grid is jumped to determine whether boundaries with communicated boundary types exist in all boundaries of the target grid to continue execution until the third grid is a repeated target grid; and combining all the target grids into a target connected domain, taking the target connected domain as a target connected domain containing the first grid, and marking the target connected domain.
Optionally, determining coordinates and reference pressure of one reference point of each connected domain comprises: randomly determining a grid in each connected domain as a reference grid; and taking the midpoint of the reference grid as a reference point, and taking the preset reference pressure as the reference pressure of the reference point.
Optionally, the boundary type includes non-connectivity, and the method further includes: determining whether a target area file has a grid containing a non-connected boundary type; if the target area file contains grids with non-connected boundary types, determining a plurality of connected domains in the target area according to the boundary types between adjacent grids; and if the target area file does not contain grids with the boundary types of non-connected grids, determining all grids in the target area file as connected domains of the target area.
Optionally, the method further comprises: randomly determining one grid as a reference grid in all grids of the target area; and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking the preset reference pressure as the reference pressure of the reference point.
In a second aspect, an embodiment of the present application further provides a prediction apparatus for a fluid pressure field, where the prediction apparatus includes: the acquisition module is used for acquiring a target area file containing the boundary type and the label of each grid in the target area; the first determining module is used for determining a plurality of connected domains in the target area according to the boundary type between the adjacent grids; the boundary type is used for describing whether the adjacent grids are connected or not; a second determining module for determining the reference number and the reference pressure of one reference point of each connected domain; and the third determination module is used for substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area.
In a third aspect, an embodiment of the present application further provides an electronic device, including: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory communicating via the bus when the electronic device is running, the machine-readable instructions being executed by the processor to perform the steps of the method for predicting a fluid pressure field in the first aspect or any one of the possible embodiments of the first aspect.
In a fourth aspect, the present application further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of predicting the fluid pressure field in the first aspect or any one of the possible implementation manners of the first aspect.
The embodiment of the application provides a method and a device for predicting a fluid pressure field, wherein the method comprises the following steps: acquiring a target area file containing the boundary type and the label of each grid in the target area; determining a plurality of connected domains in the target area according to the boundary type between adjacent grids; the boundary type is used for describing whether the adjacent grids are connected or not; determining the reference number and the reference pressure of one reference point of each connected domain; and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area. According to the method and the device, the boundary type of each grid in the target area file is identified, the multiple connected domains of the target area are determined according to the boundary types of the adjacent grids, and therefore the predicted pressure of the target area is determined according to the reference point and the reference pressure of each connected domain, the technical problem that the predicted pressure determined by the target area comprising the multiple connected domains only according to one reference point and the reference pressure is inaccurate in the prior art is solved, and the technical effect of improving the accuracy of the determined predicted pressure is achieved.
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
Fig. 1 shows a flowchart of a method for predicting a fluid pressure field according to an embodiment of the present disclosure.
Fig. 2 is a flow chart illustrating another fluid pressure field prediction method provided in the embodiment of the present application.
Fig. 3 shows one of the schematic diagrams of the target area provided by the embodiment of the present application.
Fig. 4 shows a second schematic diagram of the target area provided in the embodiment of the present application.
Fig. 5 shows a third schematic diagram of the target area provided in the embodiment of the present application.
Fig. 6 shows a fourth schematic diagram of the target area provided by the embodiment of the present application.
Fig. 7 is a functional block diagram of a device for predicting a fluid pressure field according to an embodiment of the present disclosure.
Fig. 8 shows a schematic structural diagram of an electronic device provided in an embodiment of the present application.
Detailed Description
To make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it should be understood that the drawings in the present application are for illustrative and descriptive purposes only and are not used to limit the scope of protection of the present application. Additionally, it should be understood that the schematic drawings are not necessarily drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of the present application. It should be understood that the operations of the flow diagrams may be performed out of order, and that steps without logical context may be performed in reverse order or concurrently. One skilled in the art, under the guidance of this application, may add one or more other operations to, or remove one or more operations from, the flowchart.
In addition, the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, as presented in the figures, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be obtained by a person skilled in the art without making any inventive step based on the embodiments of the present application, fall within the scope of protection of the present application.
In the prior art, a prediction area comprising a plurality of connected domains is still predicted by using a reference point and a reference pressure, so that the predicted fluid pressure cannot be obtained.
Based on this, embodiments of the present application provide a method and an apparatus for predicting a fluid pressure field, where a boundary type of each mesh in a target area file is identified, and multiple connected domains of a target area are determined according to the boundary type of an adjacent mesh, so as to determine a predicted pressure of the target area according to a reference point and a reference pressure of each connected domain, thereby solving a technical problem in the prior art that, for a target area including multiple connected domains, the predicted pressure is inaccurate or the pressure cannot be predicted only according to one reference point and one reference pressure, and achieving a technical effect of improving accuracy of determining the predicted pressure. The method comprises the following specific steps:
referring to fig. 1, fig. 1 is a flowchart illustrating a method for predicting a fluid pressure field according to an embodiment of the present disclosure. As shown in fig. 1, a method for predicting a fluid pressure field provided in an embodiment of the present application includes the following steps:
s101: and acquiring a target area file containing the boundary type and the label of each grid in the target area.
That is, a target area file is obtained, which contains the boundary type of each mesh in the target area, as well as the total number of boundaries of the target area and the label of each mesh. The target area here is the fluid pressure field.
The boundary type is used for describing whether the adjacent grids are connected or not, and the boundary type comprises connected and non-connected. The boundaries of the meshes refer to the contact areas between adjacent meshes within the target area. It is understood that for each mesh in the target area, a mesh is considered to be a boundary mesh of the target area if there is at least one boundary of the mesh that is not connected to other meshes outside the mesh.
For example, please refer to fig. 2, fig. 2 is a schematic diagram of a target area according to an embodiment of the present disclosure. As shown in fig. 2, the target area includes 24 grids, the grids numbered 1, 2, 3, 4, 5, 9, 13, 17, 8, 12, 16, 20, 21, 22, 23, and 24 are boundary grids of the target area, the grid numbered 1 is abbreviated as grid 1, and the subsequent grids are all abbreviated as such. Taking grid 1 as an example, grid 1 is adjacent to grids 2 and 5, and then grid 1 is considered to have only two boundaries, and the target area file includes: grid 1 has two boundaries, each of which is connected in boundary type. Taking grid 6 as an example, grid 6 is adjacent to grid 2, grid 5, grid 7 and grid 10, such that grid 6 has 4 boundaries, and the target region file contains: the grid 6 has 4 boundaries, one boundary being of non-connected boundary type and three boundary being of connected boundary type.
S102: determining a plurality of connected domains in the target region according to the boundary type between the adjacent grids.
Determining a plurality of connected domains in the target region according to the boundary type between the adjacent grids, comprising: randomly selecting one grid from all grids of the target area file as a first grid; determining a target connected domain containing the first grid and marking the target connected domain according to the boundary type of the first grid; determining whether the number of grids of the marked target connected domain is the same as the number of grids of the target area; if the number of the grids of the marked target connected domain is different from the number of the grids of the target area, randomly selecting a second grid which does not belong to the target connected domain in the target area; taking the second grid as a new first grid, jumping to a boundary type according to the first grid, determining a target connected domain containing the first grid, marking the target connected domain and continuing to execute the process until the grid number of all marked target connected domains is the same as that of the target area; all marked target connected domains are taken as a plurality of connected domains in the target area.
For example, as shown in fig. 2, if the first grid is grid 1, it is determined that the connected domain including grid 1 is composed of all grids in the target region except grid 10 and grid 11; marking the connected domain containing grid 1 as a first connected domain, wherein the first connected domain has 22 grids, and if the grid number of the first connected domain is determined to be less than that of the target area, randomly selecting a grid which does not belong to the first connected domain in the target area as grid 10; taking the grid 10 as a new first grid, and determining that a connected domain containing the grid 10 consists of the grid 10 and a grid 11 in the target area; grid 10 and grid 11 are combined and labeled as a second connected domain; and determining that the sum of the grid number of the first connected domain and the grid number of the second connected domain is equal to the grid number of the target area, and taking the first connected domain and the second connected domain as a plurality of connected domains of the target area.
Determining a target connected domain containing the first grid and marking the target connected domain according to the boundary type of the first grid, wherein the method comprises the following steps: taking the first grid as a target grid; determining whether a boundary of which the boundary type is connected exists in the target grid; and if the boundary type of the target grid is not a connected boundary, taking the target grid as a target connected domain containing the first grid and marking the target connected domain.
That is, the first mesh is used as the target mesh, and if the boundary types of all the boundaries of the target mesh are all non-connected, the target mesh is used as a connected domain of the target area, and the connected domain is marked.
After determining whether the target grid has a boundary with a connected boundary type, the method further includes: if the boundary type of the target grid is a connected boundary, determining the boundary with the connected boundary type in the target grid as a first boundary; determining a third grid corresponding to the first boundary except the target grid; determining whether the third grid is a present target grid; if the third grid is not the appeared target grid, the third grid is taken as a new target grid, and the third grid is jumped to determine whether a boundary with a connected boundary type exists in all boundaries of the target grid to continue execution until the third grid is a repeated target grid; and combining all the target grids into a target connected domain, taking the target connected domain as a target connected domain containing the first grid, and marking the target connected domain.
For example, as shown in fig. 2, if grid 1 is the first grid selected randomly, grid 1 is used as the target grid, and it is determined that there are two boundaries of which the boundary types are connected in the target grid, namely, a boundary adjacent to grid 5 and a boundary adjacent to grid 2. And then determining grid 5 and grid 2 as a third grid, determining that grid 5 and grid 2 are not repeated with grid 1, further taking grid 5 and grid 2 as new target grids, determining third grids corresponding to grid 5 and grid 2 respectively, wherein the third grids corresponding to grid 5 are grid 6 and grid 9, the third grids corresponding to grid 2 are grid 6 and grid 3, grid 6, grid 9 and grid 3 are not repeated with grid 5, grid 2 and grid 1, taking grid 6, grid 9 and grid 3 as new target grids, stopping querying until the third grids are the existing target grids, and combining all target grids into a first connected domain containing grid 1.
Illustratively, the number of grids of the target connected domain including grid 1 is different from the number of grids of the target region, and then grid 10 is determined as a new first grid, grid 10 is taken as a target grid, and grid 10 is determined to include 1 boundary type as a connected boundary; determining a third grid corresponding to the grid 10 as a grid 11, not repeating the grid 11 and the grid 10, further taking the grid 11 as a new target grid, and determining that the grid 11 comprises 1 boundary with a connected boundary type; and determining that the third grid corresponding to the grid 11 is the grid 10, and the grid 10 is the appeared target grid, stopping querying, and combining all the target grids, namely the grid 10 and the grid 11, into a second connected domain containing the grid 10.
S103: and determining the corresponding index and reference pressure of one reference point of each connected domain.
Determining coordinates and reference pressure of a reference point for each connected domain, comprising: randomly determining a grid in each connected domain as a reference grid; and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking the preset reference pressure as the reference pressure of the reference point.
The preset reference pressure may be set by a user, or may be a preset pressure value designated in the CFD software.
That is, one grid may be randomly determined as a reference grid in each connected domain by the user; taking the midpoint of the reference grid as a reference point, and taking the mark of the reference grid as a mark corresponding to the reference point; and taking the preset reference pressure set by the user as the reference pressure of the reference point. Or, a grid is randomly determined in each connected domain by CFD software to be used as a reference grid; taking the midpoint of the reference grid as a reference point, and taking the mark number of the reference grid as the mark number corresponding to the reference point; and taking a preset pressure value in the CFD software as a reference pressure of a reference point.
S104: and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target region.
Figure M_221223104156870_870464001
(1)
Figure M_221223104156932_932949001
(2)
Figure M_221223104157011_011996001
(3)
The formulas (1) to (3) are a Navier-Stokes equation, the formula (1) is a continuity equation, the formula (2) is a momentum conservation equation, and the formula (3) is an energy conservation equation.
In the formulas (1) to (3),
Figure M_221223104157074_074534001
is the fluid density of the target area>
Figure M_221223104157090_090657002
Is time, is>
Figure M_221223104157106_106291003
Is the fluid speed of the target area>
Figure M_221223104157137_137533004
Is the fluid pressure of the target area>
Figure M_221223104157153_153157005
Is the first fluid viscosity index of the target area>
Figure M_221223104157168_168780006
Is the fluid temperature of the target area>
Figure M_221223104157184_184413007
Is a second fluid viscosity index of the target area and +>
Figure M_221223104157217_217108008
,/>
Figure M_221223104157248_248896009
Is acceleration of gravity, </or>
Figure M_221223104157264_264496010
Is the total enthalpy of the target region and>
Figure M_221223104157295_295761011
,/>
Figure M_221223104157326_326990012
is the fluid enthalpy of the target region>
Figure M_221223104157342_342611013
Is kinetic energy and->
Figure M_221223104157373_373919014
,/>
Figure M_221223104157389_389481015
Is the heat diffusion coefficient of the fluid of the target area, < > is >>
Figure M_221223104157422_422171016
Is the fluid stress tensor of the target region.
For the case where the target region is incompressible, the fluid density of the target region is constant, and for most incompressible cases, the energy conservation equation is generally not considered, so the navier-stokes equation can be reduced to:
Figure M_221223104157438_438283001
(4)
Figure M_221223104157469_469554001
(5)
in the formula (4) and the formula (5),
Figure M_221223104157500_500804001
is a moving viscosity index of the target area and->
Figure M_221223104157532_532064002
In the formula (4) and the formula (5), the unknown quantity is
Figure M_221223104157547_547697001
And &>
Figure M_221223104157578_578949002
All in the momentum equation. In contrast to the compressible continuity equation, the incompressible continuity equation, equation (4), becomes the constraint of velocity and cannot be directly solved for fluid density. Therefore, in order to solve, the form written as laplace's equation, i.e., pressure poisson's equation, is required.
The pressure poisson equation is:
Figure M_221223104157594_594571001
(6)
in the formula (6), the first and second groups,
Figure M_221223104157622_622855001
is the first coefficient of the pressure poisson equation>
Figure M_221223104157654_654626002
Is the second coefficient of the pressure Poisson equation, n is the number of grids in the target area, and->
Figure M_221223104157670_670244003
Refers to the pressure of the grid labeled n in the target area.
That is, after the label and the reference pressure corresponding to the reference point are determined, the predicted pressure of each grid in the target area can be obtained according to the pressure poisson equation, and the predicted pressure of each grid is used as the predicted pressure of the target area.
The method further comprises the following steps: determining whether a target area file has a grid containing a non-connected boundary type; if the target area file contains grids with non-connected boundary types, determining a plurality of connected domains in the target area according to the boundary types between adjacent grids; and if the target area file does not contain grids with the boundary types of non-connected grids, determining all grids in the target area file as connected domains of the target area.
That is, if no mesh in the target area file contains a boundary whose boundary type is non-connected, it is determined that there are no multiple connected domains in the target area, and all meshes of the target area are combined to be one connected domain.
The method further comprises the following steps: randomly determining a grid in all grids of the target area as a reference grid; and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking the preset reference pressure as the reference pressure of the reference point.
That is, if the target area only contains one connected domain, acquiring that a grid randomly determined by the user in all grids of the target area is used as a reference grid; and taking the midpoint of the reference grid as a reference point, taking the label of the reference grid as the label of the reference point, and taking the preset reference pressure set by the user as the reference pressure of the reference point. Or, the CFD software randomly determines one mesh among all meshes of the target area as a reference mesh; and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking the preset pressure value of the CFD software as the reference pressure of the reference point.
When calculating the pressure at the boundary of the target area, different boundary conditions are defined, and are generally classified into a first type of boundary condition and a second type of boundary condition. The first type of boundary condition is the pressure value at the specified boundary, or the boundary is the specified pressure boundary, i.e. the specified pressure boundary is specified
Figure M_221223104157701_701500001
,/>
Figure M_221223104157717_717128002
Is the pressure value of the grid labeled i located at the boundary of the target region. A second type of boundary condition is the assignment of a pressure derivative value at the boundary, for example pressure Poisson's equation, i.e. assignment of->
Figure M_221223104157748_748385003
Wherein->
Figure M_221223104157795_795246004
Refers to the normal gradient of the pressure of the grid i in the preset x-axis direction, <' >>
Figure M_221223104157827_827957005
Refers to the normal gradient of the pressure of the grid i in the direction of the preset y-axis, is>
Figure M_221223104157859_859227006
Refers to the normal gradient of the pressure of the grid i in the direction of the predetermined z-axis, is->
Figure M_221223104157890_890454007
Is directed to a predetermined x-axis directionIs predetermined normal gradient of (4), is greater than or equal to>
Figure M_221223104157906_906099008
Refers to a predetermined normal gradient in the direction of a predetermined y-axis>
Figure M_221223104157937_937357009
Refers to a predetermined normal gradient to a predetermined z-axis direction. The boundary conditions of all the boundaries of the target area may be all the first type of boundary conditions, may be all the second type of boundary conditions, or may include both the first and second types of boundary conditions.
In the process of solving the pressure poisson equation, if all the boundary conditions of the pressure poisson equation are the second type of boundary conditions, the equation cannot be uniquely solved. A common situation in CFD is wall conditions. It is then necessary to pick a reference point inside the calculation area and to specify a reference pressure for this reference point.
There are several typical situations in the target area, for example, please refer to fig. 3, and fig. 3 is one of the schematic diagrams of the target area. As shown in FIG. 3, the target regions are all connected together, the boundary types of all the grids in the target region file are connected, and the target region has more than one boundary condition as the boundary of the first type of boundary condition, i.e. p 1 Corresponding boundary, p 1 =p ambient ,p ambient As an ambient pressure, or may be a user-specified pressure. At this time, it is not necessary to specify the reference point and the reference pressure.
For example, please refer to fig. 4, fig. 4 is a second schematic diagram of the target area. As shown in FIG. 4, the target areas are all connected together, and the boundary types of all the grids in the target area file are connected. The boundary conditions of all the boundaries of the target region are the second type of boundary conditions, and when the target region is considered to be a connected domain, a reference point and a reference pressure of the reference point are required to be specified to calculate a pressure Poisson equation, such as p 2 Pointed black dots as reference points, p 2 Is a reference pressure.
Referring to FIG. 5, FIG. 5 shows a target areaIn the third schematic diagram, as shown in FIG. 5, the boundary condition of the target area having more than one boundary is the specified pressure boundary condition, i.e. p 3 Corresponding boundary, p 3 =p ambient ,p ambient As ambient pressure or may be a user-specified pressure. But at least one connected domain surrounded by the second type of boundary condition exists in the target area, and the connected domain is not communicated with the fluid outside the connected domain. It is now necessary to specify a reference point and a reference pressure, for example p, in each communication field 31 Pointed black dots as reference points, p 31 Is a reference pressure.
Illustratively, referring to fig. 6, fig. 6 is a fourth schematic diagram of the target area, as shown in fig. 6, the boundary conditions of all the boundaries of the target area are wall surfaces or second type boundary conditions, and at least one connected domain surrounded by the second type boundary conditions exists inside the target area, and the connected domain is not communicated with the fluid outside the connected domain. It is then necessary to specify a reference point and a reference pressure, for example p, in each communication field 41 Pointed black dots as reference points, p 41 Is a reference pressure. The portion outside the connected component that belongs to the target area is also a connected component and requires the specification of a reference point and a reference pressure, such as p 42 Pointed black dots as reference points, p 42 Is a reference pressure.
The boundary conditions of all the boundaries of the target area are defaulted to be wall surfaces or are the second type of boundary conditions, namely the boundary conditions of the boundaries of the target area are not considered. The boundary of the target area corresponds to a boundary condition which is divided into a first type of boundary condition and a second type of boundary condition; the boundaries of all grids in the target area correspond to boundary types which are divided into connected and non-connected.
Based on the same application concept, a fluid pressure field prediction device corresponding to the fluid pressure field prediction method provided by the above embodiment is also provided in the embodiment of the present application, and as the principle of solving the problem of the device in the embodiment of the present application is similar to the fluid pressure field prediction method provided by the above embodiment of the present application, the implementation of the device may refer to the implementation of the method, and repeated details are omitted.
Fig. 7 is a functional block diagram of a fluid pressure field predicting apparatus according to an embodiment of the present disclosure, as shown in fig. 7. The fluid pressure field prediction device 10 includes: an acquisition module 101, a first determination module 102, a second determination module 103, and a third determination module 104.
An obtaining module 101, configured to obtain a target area file including a boundary type and a label of each mesh in a target area; a first determining module 102, configured to determine a plurality of connected domains in the target region according to a boundary type between adjacent grids; the boundary type is used for describing whether adjacent grids are connected or not; a second determination module 103 for determining the reference number and reference pressure of one reference point of each connected domain; and a third determining module 104, configured to substitute the reference number and the reference pressure corresponding to each connected component into a pressure poisson equation, and determine a predicted pressure of the target region.
Based on the same application concept, referring to fig. 8, a schematic structural diagram of an electronic device provided in the embodiment of the present application is shown, where the electronic device 20 includes: a processor 201, a memory 202 and a bus 203, wherein the memory 202 stores machine-readable instructions executable by the processor 201, and when the electronic device 20 is running, the processor 201 communicates with the memory 202 via the bus 203, and the machine-readable instructions are executed by the processor 201 to perform the steps of the method for predicting the fluid pressure field according to any of the above embodiments.
In particular, the machine readable instructions, when executed by the processor 201, may perform the following: acquiring a target area file containing the boundary type and the label of each grid in the target area; determining a plurality of connected domains in the target area according to the boundary type between adjacent grids; the boundary type is used for describing whether adjacent grids are connected or not; determining the reference number and the reference pressure of one reference point of each connected domain; and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area.
Based on the same application concept, the present application further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program performs the steps of the method for predicting a fluid pressure field provided by the foregoing embodiment.
Specifically, the storage medium can be a general storage medium, such as a mobile disk, a hard disk, or the like, when a computer program on the storage medium is executed, the method for predicting the fluid pressure field can be executed, the boundary type of each grid in the target area file is identified, and a plurality of connected domains of the target area are determined according to the boundary type of the adjacent grid, so that the predicted pressure of the target area is determined according to the reference point and the reference pressure of each connected domain, the technical problem that the predicted pressure determined for the target area including a plurality of connected domains only according to one reference point and reference pressure is inaccurate in the prior art is solved, and the technical effect of improving the accuracy of the determined predicted pressure is achieved.
It can be clearly understood by those skilled in the art that, for convenience and simplicity of description, the specific working process of the system and the apparatus described above may refer to the corresponding process in the foregoing method embodiment, and details are not described herein again. In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described apparatus embodiments are merely illustrative, and for example, the division of the units into only one type of logical function may be implemented in other ways, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not implemented. In addition, the shown or discussed coupling or direct coupling or communication connection between each other may be through some communication interfaces, indirect coupling or communication connection between devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of grid units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a grid device) to execute all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (10)

1. A method of predicting a fluid pressure field, the method comprising:
acquiring a target area file containing the boundary type and the label of each grid in the target area;
determining a plurality of connected domains in the target area according to the boundary type between adjacent grids; the boundary type is used for describing whether adjacent grids are connected or not;
determining the reference number and the reference pressure of one reference point of each connected domain;
and substituting the reference number and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area.
2. The method of claim 1, wherein determining the plurality of connected components in the target region according to the boundary type between the adjacent grids comprises:
randomly selecting one grid from all grids of the target area file as a first grid;
determining a target connected domain containing the first grid and marking the target connected domain according to the boundary type of the first grid;
determining whether the number of grids of the marked target connected domain is the same as the number of grids of the target area;
if the number of the grids of the marked target connected domain is different from the number of the grids of the target area, randomly selecting a second grid which does not belong to the target connected domain in the target area;
taking the second grid as a new first grid, jumping to a boundary type according to the first grid, determining a target connected domain containing the first grid, marking the target connected domain and continuing to execute until the grid number of all marked target connected domains is the same as that of the target area;
and taking all marked target connected domains as a plurality of connected domains in the target area.
3. The method of claim 2, wherein the boundary types include connected and non-connected, and wherein determining a target connected domain containing a first mesh and marking the target connected domain according to the boundary type of the first mesh comprises:
taking the first grid as a target grid;
determining whether a boundary of which the boundary type is connected exists in the target grid;
and if the target grid does not have a boundary with a connected boundary type, taking the target grid as a target connected domain containing the first grid and marking the target connected domain.
4. The method of claim 3, wherein after determining whether the target mesh has a boundary with a boundary type of connected, the method further comprises:
if the boundary type of the target grid is a connected boundary, determining the boundary with the connected boundary type in the target grid as a first boundary;
determining a third grid corresponding to the first boundary except the target grid;
determining whether the third mesh is the present target mesh;
if the third grid is not the appeared target grid, taking the third grid as a new target grid, and jumping to determine whether a boundary with a connected boundary type exists in all boundaries of the target grid to continue execution until the third grid is a repeated target grid;
and combining all the target grids into a target connected domain, serving as the target connected domain containing the first grid, and marking the target connected domain.
5. The method of claim 1, wherein determining coordinates and reference pressure for one reference point for each connected domain comprises:
randomly determining a grid in each connected domain as a reference grid;
and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking preset reference pressure as the reference pressure of the reference point.
6. The method of claim 1, wherein the boundary type comprises non-connectivity, the method further comprising:
determining whether a grid containing a boundary type which is non-connected exists in the target area file;
if the target area file contains grids with non-connected boundary types, determining a plurality of connected domains in the target area according to the boundary types between adjacent grids;
if the target area file does not contain grids with non-connected boundary types, all grids in the target area file are determined as connected domains of the target area.
7. The method of claim 6, wherein after determining all grids in the target area file as connected domains of the target area, the method further comprises:
randomly determining a grid in all grids of the target area as a reference grid;
and taking the midpoint of the reference grid as a reference point, taking the mark of the reference grid as the mark of the reference point, and taking preset reference pressure as the reference pressure of the reference point.
8. A prediction apparatus of a fluid pressure field, the prediction apparatus comprising:
the acquisition module is used for acquiring a target area file containing the boundary type and the label of each grid in the target area;
a first determining module, configured to determine a plurality of connected domains in the target region according to a boundary type between adjacent grids; the boundary type is used for describing whether adjacent grids are connected or not;
the second determining module is used for determining the label and the reference pressure of one reference point of each connected domain;
and the third determination module is used for substituting the reference numbers and the reference pressure corresponding to each connected domain into a pressure Poisson equation to determine the predicted pressure of the target area.
9. An electronic device, comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory communicating over the bus when the electronic device is operating, the machine-readable instructions when executed by the processor performing the steps of the method of predicting a fluid pressure field according to any one of claims 1 to 7.
10. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method of predicting a fluid pressure field according to any one of claims 1 to 7.
CN202211670441.3A 2022-12-26 2022-12-26 Fluid pressure field prediction method and device Pending CN115906709A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211670441.3A CN115906709A (en) 2022-12-26 2022-12-26 Fluid pressure field prediction method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211670441.3A CN115906709A (en) 2022-12-26 2022-12-26 Fluid pressure field prediction method and device

Publications (1)

Publication Number Publication Date
CN115906709A true CN115906709A (en) 2023-04-04

Family

ID=86474562

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211670441.3A Pending CN115906709A (en) 2022-12-26 2022-12-26 Fluid pressure field prediction method and device

Country Status (1)

Country Link
CN (1) CN115906709A (en)

Similar Documents

Publication Publication Date Title
US8682850B2 (en) Method of enhancing de-duplication impact by preferential selection of master copy to be retained
CN102681940B (en) Method for carrying out performance test on memory management subsystem of Linux operation system
WO2017020747A1 (en) Method and device for detecting slow disk
WO2017092480A1 (en) Data migration method and device
JP2009110346A (en) Management method of performance history, and management system of performance history
US10101188B2 (en) Fluid blob tracking for evaluation of multiphase flow simulations
JP5891680B2 (en) Power control apparatus, power control method, and power control program
CN105117489B (en) Database management method and device and electronic equipment
WO2008016493A1 (en) Detection of memory leaks
CN109634960B (en) Key value data storage method, device, equipment and storage medium
JP4208079B2 (en) Database server, program, recording medium, and control method
CN115906709A (en) Fluid pressure field prediction method and device
CN109656479A (en) A kind of method and device constructing memory command sequence
CN105117168A (en) Information processing method and electronic equipment
JP2014026440A (en) Simulation program, simulation device, and simulation method
CN112486988A (en) Data processing method, device, equipment and storage medium
JP2014174844A (en) Performance evaluation support program, performance evaluation support device, and performance evaluation support method
CN115630592B (en) Pressure initialization method and device for flow field
US9690545B2 (en) Floating-point calculation apparatus, program, and calculation apparatus
KR102181981B1 (en) Apparatus, method and computer program for sph based fluid analysis simulation
CN117216454B (en) Reliability assessment method and device based on fuzzy non-probability, terminal and storage medium
Dedering et al. On new challenges for CFD simulation in filtration
CN111487685B (en) Method and device for calculating resistivity of rock reservoir and storage medium
CN109388583B (en) Parameter detection method and related equipment
CN114564689A (en) Data processing method, data processing device, computer equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination