CN115952752B - Data processing method based on flow field simulation robustness calculation under enhanced inclined grid - Google Patents
Data processing method based on flow field simulation robustness calculation under enhanced inclined grid Download PDFInfo
- Publication number
- CN115952752B CN115952752B CN202310227478.7A CN202310227478A CN115952752B CN 115952752 B CN115952752 B CN 115952752B CN 202310227478 A CN202310227478 A CN 202310227478A CN 115952752 B CN115952752 B CN 115952752B
- Authority
- CN
- China
- Prior art keywords
- target
- pressure
- field
- value
- cell
- 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.)
- Active
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The application discloses a data processing method based on flow field simulation robustness calculation under an enhanced inclined grid, which is used for matching an interpolation method applied to an internal interface of a target body with an interpolation method of an outlet boundary, so that the value on the internal interface and the value on the outlet boundary are simultaneously considered for inclination correction of the grid, and the interpolated value is more fit. By adopting the method provided by the application, the robustness can be improved, and the data processing effect can be further improved.
Description
Technical Field
The application belongs to the field of simulation data processing research, and particularly relates to a data processing method based on flow field simulation robustness calculation under an enhanced inclined grid.
Background
The simulation data processing is based on a simulation theory, uses a computer system and physical effect equipment as tools, builds and runs a model according to a target body, and realizes and reforms a researched object, so that the simulation data processing is an information type selectable technology generated in the progress of an industrialized society to an informatization society.
In the related art, when a SIMPLE algorithm is adopted to solve low-speed flow in the process of processing simulation data, a grid with larger skewness is often encountered, and the grid can cause calculation to be unable to converge. Poor robustness results in a less than ideal data processing effect.
Disclosure of Invention
In order to solve the defects in the prior art, the application provides a data processing method based on flow field simulation robustness calculation under an enhanced inclined grid, which is used for matching an interpolation method applied to an internal interface of a target body with an interpolation method applied to an outlet boundary, so that the value on the internal interface and the value on the outlet boundary simultaneously consider the inclination correction of the grid, and the interpolated value is more fit. By adopting the method provided by the application, the robustness can be improved, and the data processing effect can be further improved.
The technical effects to be achieved by the application are realized by the following scheme:
in a first aspect, the present specification provides a data processing method based on flow field simulation robustness calculation under an enhanced inclined grid, the method comprising:
obtaining a model of a target body, wherein the model represents an internal interface and an outlet boundary surface of the target body; at least part of the flow field formed by the model is divided into a plurality of grid cells; one part of the grid cells is an inclined grid, and the other part is a non-inclined grid;
setting an initial velocity field and an initial pressure field;
correcting the inclined grid as a target, and performing momentum interpolation on the initial velocity field to obtain a calculation model of a target velocity field represented by a target pressure field;
Establishing a pressure correction equation by adopting a calculation model of the target speed field;
calculating the pressure correction equation to obtain a pressure correction intermediate value;
interpolation processing is carried out on the pressure correction intermediate value belonging to the outlet boundary surface, so as to obtain a pressure correction target value;
obtaining a target speed field output by a calculation model of the target speed field based on a target pressure field obtained by correcting the target value according to the pressure;
and carrying out data processing based on the target speed field and the target pressure field.
In an alternative embodiment of the present specification, the formula of the calculation model of the target speed field is expressed as:
wherein: f represents the value of the quantity to which the interface belongs on an interface, wherein the interface comprises the internal interface and/or the outlet boundary surface; the upper dash represents the average value; />Representing the value of the pressure gradient at the interface; />Representing the known coefficients; />Representing a normal direction of the interface; />Representing a target speed field obtained by the iterative computation; />Indicating the previous timeIteratively calculating the average value of the historical target speed field; p represents the target pressure field; m represents the number of iteration steps.
In an alternative embodiment of the present disclosure, for a target unit, the value of the pressure gradient at the interface is calculated by the following formula; wherein the target cell is any one of the grid cells:
/>
Wherein: c represents the target unit; n represents an adjacent cell, wherein the adjacent cell is a grid cell adjacent to the target cell; p (P) N Representing the pressure value at the center of the adjacent cell; pc represents a pressure value of the target cell center;representation->Is a projection point of (2); />Representation->Is a projection point of (2); />Representing the pressure gradient at the centers of the adjacent cells; />A pressure gradient representing the center of the target cell; />A vector representing a center point of the adjacent cell; />A vector representing a center point of the target cell; />A vector representing an interface center point between the target cell and the neighboring cell; />And the first distance is the projection distance of the distance from the central point of the target unit to the central point of the interface in the normal direction, and the second distance is the projection distance of the distance from the central point of the adjacent unit to the central point of the interface in the normal direction.
In an alternative embodiment of the present disclosure, for the target unit, the value of the pressure gradient at the interface is calculated by the following formula:
wherein: />A direction vector representing the center of the target cell to the center of the adjacent cell.
In an alternative embodiment of the present specification, the formula of the pressure correction equation is expressed as:
in (1) the->Is a diffusion term; p' is a pressure correction intermediate value; u (U) * Is the initial pressure field.
In an alternative embodiment of the present disclosure, the interpolation is performed on the pressure correction intermediate value that belongs to the outlet boundary surface, and the interpolation calculation formula adopted when obtaining the pressure correction target value is:
wherein d represents an interpolation distance from the center of the target unit to the center of the outlet boundary surface; />Representing the pressure correction target value; />A gradient field representing a pressure correction target value; />Is a pressure correction target value for the center of the known cell.
In an alternative embodiment of the present disclosure, obtaining the target velocity field output by the calculation model of the target velocity field based on the target pressure field obtained by correcting the target value according to the pressure includes:
solving a calculation model of the target speed field based on the target pressure field to obtain a speed correction target value;
and correcting the initial speed field by adopting the speed correction target value to obtain a target speed field.
In an alternative embodiment of the present specification, the target velocity field is calculated by the following formula:
Wherein: u (U) * Is the target velocity field; u (U) *-1 Is the historical target speed field obtained in the last iteration; u' is the speed correction target value.
In an alternative embodiment of the present disclosure, before obtaining the target velocity field output by the calculation model of the target velocity field based on the target pressure field obtained by correcting the target value according to the pressure, the method further includes:
and correcting the initial pressure field by adopting the pressure correction target value to obtain a target pressure field.
In an alternative embodiment of the present specification, the target pressure field is calculated by the following formula:
wherein: p (P) * Is the target pressure field; p (P) *-1 Is the historical target pressure field obtained by the last iteration; p' is the pressure correction target value.
The present specification provides an apparatus for data processing based on enhanced under-diagonal grid flow field simulation robustness calculation for implementing the method in the first aspect, the apparatus comprising:
the model acquisition module is configured to: obtaining a model of a target body, wherein the model represents an internal interface and an outlet boundary surface of the target body; at least part of the flow field formed by the model is divided into a plurality of grid cells; one part of the grid cells is an inclined grid, and the other part is a non-inclined grid;
A setting module configured to: setting an initial velocity field and an initial pressure field;
the calculation model generation module of the target speed field is configured to: correcting the inclined grid as a target, and performing momentum interpolation on the initial velocity field to obtain a calculation model of a target velocity field represented by a target pressure field;
the pressure correction equation establishment module is configured to: establishing a pressure correction equation by adopting a calculation model of the target speed field;
a pressure correction intermediate value determination module configured to: calculating the pressure correction equation to obtain a pressure correction intermediate value;
a pressure correction target value determination module configured to: interpolation processing is carried out on the pressure correction intermediate value belonging to the outlet boundary surface, so as to obtain a pressure correction target value;
a target speed field determination module configured to: obtaining a target speed field output by a calculation model of the target speed field based on a target pressure field obtained by correcting the target value according to the pressure;
a data processing module configured to: and carrying out data processing based on the target speed field and the target pressure field.
In a third aspect, the present specification provides an electronic device comprising:
A processor; and
a memory arranged to store computer executable instructions which, when executed, cause the processor to perform the method of the first aspect.
In a fourth aspect, the present description provides a computer-readable storage medium storing one or more programs, which when executed by an electronic device comprising a plurality of application programs, cause the electronic device to perform the method of the first aspect.
Drawings
In order to more clearly illustrate the embodiments of the application or the prior art solutions, the drawings which are used in the description of the embodiments or the prior art will be briefly described below, it being obvious that the drawings in the description below are only some of the embodiments described in the present application, and that other drawings can be obtained according to these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a flow chart of a data processing method based on flow field simulation robustness calculation under an enhanced inclined grid in an embodiment of the application;
FIG. 2 is a flow chart of a data processing method based on enhanced under-grid flow field simulation robustness calculation in an alternative embodiment of the application;
FIG. 3 is a schematic view of a center point projection of an unstructured grid cell in an embodiment of the present application;
FIG. 4 is a 2D simplified schematic of a sloped mesh at the exit boundary surface in an embodiment of the present application;
FIG. 5 is a schematic structural diagram of a data processing device based on flow field simulation robustness calculation under an enhanced inclined grid in an embodiment of the application;
fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The application will be described in further detail below with reference to the drawings by means of specific embodiments. Wherein like elements in different embodiments are numbered alike in association. In the following embodiments, numerous specific details are set forth in order to provide a better understanding of the present application. However, one skilled in the art will readily recognize that some of the features may be omitted, or replaced by other elements, materials, or methods in different situations. In some instances, related operations of the present application have not been shown or described in the specification in order to avoid obscuring the core portions of the present application, and may be unnecessary to persons skilled in the art from a detailed description of the related operations, which may be presented in the description and general knowledge of one skilled in the art.
Furthermore, the described features, operations, or characteristics of the description may be combined in any suitable manner in various embodiments. Also, various steps or acts in the method descriptions may be interchanged or modified in a manner apparent to those of ordinary skill in the art. Thus, the various orders in the description and drawings are for clarity of description of only certain embodiments, and are not meant to be required orders unless otherwise indicated.
The numbering of the components itself, e.g. "first", "second", etc., is used herein merely to distinguish between the described objects and does not have any sequential or technical meaning. The term "coupled" as used herein includes both direct and indirect coupling (coupling), unless otherwise indicated.
Various non-limiting embodiments of the present application are described in detail below with reference to the attached drawing figures. A data processing method based on flow field simulation robustness calculation under an enhanced inclined grid in the specification, as shown in fig. 1, comprises the following steps:
s100: a model of the target volume is acquired.
The model of the present specification shows the internal interface and the outlet boundary surface of the object. At least part of the flow field formed by the model is divided into a number of grid cells. One part of the grid cells is a tilted grid, and the other part is a non-tilted grid.
The model in this specification may be a CAD file, for example. The manner of acquisition in this step may be determined according to actual requirements. For example, the model may be obtained by file import, or may be obtained by modeling.
The object in this description is a component of which during operation there is an airflow through at least part of its structure. By way of example, the target may be an engine (e.g., an engine of an aircraft), a turbine, or the like. When the airflow passes through the target body, the space provided by the periphery and/or the interior of the target body for the airflow to flow is a flow field. The interface inside the flow field is an internal interface; the interface at which the gas flows out of the flow field is the outlet boundary surface.
The present specification does not limit the specific shape of the mesh unit, and the mesh unit may be tetrahedral, hexahedral, or the like, for example. In the same model, there may be grid cells of different shapes. It can be seen that a model divides a grid cell into a plurality of grid cells.
In an alternative embodiment of the present specification, the grid cells in the present specification are unstructured grids. Alternatively, a specific implementation of this step may be to import the grid cells that make up the model into a solver.
S102: an initial velocity field and an initial pressure field are set.
In an alternative embodiment of the present description, the initial velocity field and/or the initial pressure field may be set upon triggering by a person; in another alternative embodiment of the present description, the initial velocity field and/or the initial pressure field is automatically set by the executing body of the present application.
S104: and correcting the inclined grid as a target, and performing momentum interpolation on the initial velocity field to obtain a calculation model of the target velocity field expressed by the target pressure field.
The target pressure field in the specification is a pressure field actually adopted when data processing is performed, and the target pressure field has a certain difference compared with the initial pressure field. The target speed field in the present specification is a speed field actually used in data processing, and there is a certain difference between the target speed field and the initial speed field.
In this step, the target pressure field is not yet obtained, and the target pressure field may be characterized by a computational model of the target pressure field. The computational model in this specification can be characterized in terms of formulas. After the target pressure field is obtained in the subsequent step, the target velocity field can be further obtained based on the target pressure field.
In an alternative embodiment of the present disclosure, a SIMPLE algorithm may be used to calculate a momentum equation based on a given initial velocity field, initial pressure field, to obtain a computational model of the target velocity field.
In an alternative embodiment of the present description, the pre-processing of the grid cells is also performed before performing this step.
S106: and establishing a pressure correction equation by adopting the calculation model of the target speed field.
From this step, the target pressure field needs to be obtained by calculation. The pressure correction equation in this step is used to output a pressure correction intermediate value.
S108: and obtaining a pressure correction intermediate value through calculation of the pressure correction equation.
After the pressure correction intermediate value is obtained, it may be processed to obtain a pressure correction target value, and further a target pressure field is obtained based on the pressure correction target value.
In an alternative embodiment of the present disclosure, the pressure correction equation may be discretely solved according to the divergence theorem, so as to obtain the pressure correction intermediate value.
S110: interpolation processing is carried out on the pressure correction intermediate value belonging to the outlet boundary surface, so as to obtain a pressure correction target value;
Based on the steps, the interpolation method on the outlet boundary is matched with the interpolation method of the internal interface in the momentum correction process, so that the values on the internal interface and the outlet boundary simultaneously consider the inclination correction of the grid, and the calculation accuracy is ensured.
S112: and obtaining a target speed field output by a calculation model of the target speed field based on the target pressure field obtained by correcting the target value according to the pressure.
After the target pressure field is obtained, the calculation model of the target velocity field can be calculated based on the target pressure field, and the target velocity field can be obtained.
S114: and carrying out data processing based on the target speed field and the target pressure field.
The data relating to the flow field can be processed by the method in the specification under the condition of permitting. By way of example, the data processed by the methods in this specification may be pneumatic data. In this specification, the process of generating a target velocity field and/or generating a target pressure field may also be regarded as a kind of data processing.
SIMPLE type algorithms solve for low flow, often encounter grid cells with large skewness, which can result in computation failure to converge. The main reason for the failure to converge is that interpolation is performed on the interface in the calculation process, when the grid unit is inclined, the calculation error caused by the inclination of the grid unit is ignored by a general interpolation method, especially the interpolation is ignored at the position of the outlet boundary surface, so that the interpolation value on the interface has larger deviation from the true value, and therefore the flow field simulation cannot converge.
Since flow correction of the inlet and outlet is performed in the SIMPLE algorithm, the correction method mainly performs a specific process on the value of the outlet boundary surface, and thus the accuracy of the variable value on the outlet boundary surface is particularly important. In order to obtain a flow field with calculated convergence, it is necessary to perform tilt correction interpolation on grid cells with tilt, however, in general, tilt correction of interface interpolation is usually ignored, and even for values on the outlet boundary surface, values of the internal cells are duplicated without performing interface interpolation, so that abnormal values occur in numerical simulation.
At present, for the grid unit with inclination, interpolation at the outlet boundary surface is often ignored, so that the value of an internal interface and the value of the outlet boundary surface do not reach effective fit, the value after interpolation of the internal interface considers the inclination correction of the grid unit, but the value on the outlet boundary surface is not subjected to the inclination correction, the value deviation on the outlet boundary surface is larger in the flow correction process, the flow field calculation result is abnormal, the convergence cannot be well achieved, and the robustness is reduced.
According to the data processing method based on the flow field simulation robustness calculation under the enhanced inclined grid, the interpolation method applied to the internal interface of the target body is matched with the interpolation method of the outlet boundary, so that the value on the internal interface and the value on the outlet boundary are simultaneously considered for the inclination correction of the grid, and the interpolated value is more matched. By adopting the method provided by the application, the robustness can be improved, and the data processing effect can be further improved.
To further increase robustness in the data processing process and improve the data processing effect, a further alternative embodiment of the present application will now be described, with an exemplary, at least partially alternative embodiment having a flow as shown in fig. 2.
Before describing alternative embodiments, several concepts are described.
SIMPLE algorithm: semi-Implicit Method for Pressure Linked Equations, a Semi-implicit algorithm that solves the pressure coupling equation.
SIMPLE type algorithm: based on the SIMPLE algorithm, and with corresponding improvement to the defects of the SIMPLE algorithm, a series of algorithms are generated, and the algorithms are collectively called as the SIMPLE algorithms.
Momentum equation: the SIMPLE algorithm solves the equation for the velocity field, which is derived from the principle of conservation of momentum of the fluid.
Pressure correction equation: the equation for solving the pressure correction value in the SIMPLE algorithm reflects the conservation of mass of the fluid.
Grid cell: because the partial differential equation set needs to be solved in calculating the flow field, the current method for solving the equation adopts a numerical solution, namely the whole calculation domain is discrete, and the continuous flow field is decomposed into a plurality of small units, namely grids.
Grid skewness: because the shape of the flow field in reality is irregular, grid cells with larger offset angles often appear when grids are divided, and the grid cells are scattered, so that the solving of the equation can be adversely affected, and the convergence of the equation is not facilitated. Such grid cells are also referred to as tilting grids.
Boundary conditions: the change rule of the solved variable or the derivative thereof along time and position on the boundary of the solving area is the precondition that the control equation has a definite solution. Region boundaries such as exit boundary surfaces.
Internal interface: refers to the law of variation of the solved variables or their derivatives over time and position over the non-boundaries of the solved area.
In an alternative embodiment, the formulation of the calculation model of the target speed field in step S104 is shown in the following formula (one).
In the formula (one): f represents the value of the quantity to which the interface belongs on an interface, wherein the interface comprises the internal interface and/or the outlet boundary surface; the upper dash represents the average value; />Representing the value of the pressure gradient at the interface; />Representing the known coefficients; />Representing a normal direction of the interface; />Representing a target speed field obtained by the iterative computation; / >An average value of the historical target speed field obtained by the previous iterative calculation is represented; p represents the target pressure field; m represents the number of iteration steps.
In the calculation process of the pressure correction equation, velocity interpolation on the unit surface is used, and if the interpolation is incorrect, the solving of the pressure value or the velocity value can be caused to generate oscillation. The problem of velocity interpolation can be solved by performing momentum interpolation on the updated velocity field by the formula (I).
Any one of the grid cells is taken as a target cell.
In an alternative embodiment of the present disclosure, the pressure gradient is valued at the interface for the target unitThe method can be calculated by the following formula (II):
in the formula (II): />A direction vector representing the center of the target cell to the center of the adjacent cell; p (P) N Representing the pressure value at the center of the adjacent cell; p (P) C Representing the pressure value at the center of the target unit。
The pressure gradient solution realized by the formula (II) has the characteristics of simplicity and high efficiency. However, the correction capability of the inclined grid cells is poor, so that the robustness is still improved.
In a further alternative embodiment of the present disclosure, the pressure gradient is valued at the interface for the target unit The calculation can be performed by the following formulas (three) to (six): />Formula (III)>Formula (IV)>Formula (five)Formula (six)
Wherein: c represents the target unit; n represents an adjacent cell, wherein the adjacent cell is a grid cell adjacent to the target cell;representation->Is a projection point of (2); />Representation->Is a projection point of (2); />Representing the pressure gradient at the centers of the adjacent cells; />A pressure gradient representing the center of the target cell; />A vector representing a center point of the adjacent cell; />A vector representing a center point of the target cell; />A vector representing an interface center point between the target cell and the neighboring cell; />And the first distance is the projection distance of the distance from the central point of the target unit to the central point of the interface in the normal direction, and the second distance is the projection distance of the distance from the central point of the adjacent unit to the central point of the interface in the normal direction.
The calculation of the correction amount of the inclined grid cell is added through the pressure gradient solving realized from the formula (III) to the formula (six), and the advantage is that the calculation result obtained by interpolation of the opposite formula (II) can have a relaxation effect in the face of the grid cell with larger inclination, so that the calculation of the inclined correction amount is not excessively large, and the stability of the calculation process is ensured.
In an alternative embodiment, the formulation of the mesh pressure correction equation in step S106 is shown as the following equation (seven).
Formula (seven)
In the method, in the process of the invention,is a diffusion term; p' is a pressure correction intermediate value; u (U) * Is the initial pressure field.
As shown in fig. 3, considering that the improvement of interpolation by momentum alone is still insufficient, if the interpolation method is not still adopted for calculating the boundary surface value of the outlet, the overall robustness is difficult to be improved.
The pressure correction amount can be obtained by discrete solving of the pressure correction equation. The right end of the above equation is the velocity field after momentum interpolation, which is a known source term. The left end is a diffusion term which can be discretely solved according to the divergence theorem.
In different boundary conditions, since the SIMPLE algorithm is adopted, the flow correction of the outlet and the inlet exists in the calculation process, and the variable value on the outlet boundary surface is mainly adjusted in the correction process, so that the influence of the grid inclination at the outlet boundary surface is larger than the influence of the grid inclination at other positions, and the inclination grid of the outlet boundary surface is briefly shown in the following figure 4.
In the related art, the calculation of the pressure correction value P' is generally performed by directly assigning the value inside the cell to the boundary, without performing interpolation processing, as follows:
In (1) the->And->Respectively the center point position of the two outlet boundary surfaces.
As is apparent from the schematic diagram of figure 4,and->The positions of the two boundary surface center points are not identical to the distance of the cell center C, but it is obviously unreasonable to assign values to them without distinction.
Thus, the calculation of the value at the exit boundary is better done by interpolationAnd (5) calculating. By a known value of the cell center CAnd carrying out gradient interpolation according to the distances from the cell center C to the centers of different boundary surfaces. In an alternative embodiment, the interpolation calculation formula adopted for obtaining the pressure correction target value in step S110 is shown in the following formula (eight).
Formula (eight)
Wherein d represents an interpolation distance from the center of the target unit to the center of the outlet boundary surface;representing the pressure correction target value; />A gradient field representing a pressure correction target value; />Is a pressure correction target value for the center of the known cell.
Through the interpolation calculation formula shown in the formula (eight), different interpolation distances between the centers of the outlet boundary surfaces and the unit center C are provided, so that different pressure correction values can be obtained, the inclined grid is corrected to a certain extent, and the calculated value of the outlet boundary surfaces is more reasonable.
Based on the description, the interpolation method on the outlet boundary is creatively matched with the interpolation method of the internal interface in the momentum correction process, so that the values on the internal interface and the outlet boundary surface simultaneously consider the inclination correction of the grid, and the calculation accuracy is ensured.
According to the processing method of the present specification, even if there is a case where the mesh is inclined, interpolation accuracy on the interface can be satisfied. And when the grid unit is orthogonal and has no inclination, the calculation result of the inclination correction amount is 0, and the original calculation result is still not influenced. Therefore, the method can greatly improve abnormal values of calculation caused by grid inclination, and further improves the robustness of the algorithm.
After solving the obtained pressure correction value P ', the gradient field v P' of the pressure correction value can be solved. In an alternative embodiment, the process of obtaining the target speed field in step S112 may be: and solving a calculation model of the target speed field based on the target pressure field to obtain a speed correction target value. And correcting the initial speed field by adopting the speed correction target value to obtain a target speed field.
The speed correction target value can be obtained by the following formula (nine).
Formula (nine)
Where U' is a speed correction target value.
On the basis of obtaining the speed correction target value, the target speed field can be obtained by the following formula (ten).
Formula (ten)
Wherein: u (U) * Is the target velocity field; u (U) *-1 Is the historical target speed field obtained in the last iteration; u' is the speed correction target value.
For the pressure field, before the target speed field output by the calculation model of the target speed field is obtained, the initial pressure field is corrected by adopting the pressure correction target value to obtain the target pressure field. Then, the target pressure field is calculated by the following formula (eleven).
Formula (eleven)
Wherein: p (P) * Is the target pressure field; p (P) *-1 Is the historical target pressure field obtained by the last iteration; p' is the purpose of pressure correctionAnd (5) marking.
After the target speed field and the target pressure field are obtained, equation residuals are calculated, if preset residual requirements (in an alternative embodiment, manual rules can be met) are met, calculation is finished, and a flow field result is output; if the residual requirement is not satisfied, returning to the step S102 for iterative computation until the residual requirement is satisfied. Optionally, the converged flow field is then post-processed for display.
The method aims at solving the problem of non-convergence of calculation caused by large grid skewness in the existing SIMPLE algorithm flow field calculation process, creatively matches and uses two modes of internal interface interpolation and outlet boundary interpolation, and successfully solves the problem. Compared with the prior art, the method improves the robustness of the method calculation.
Based on the same thought, the embodiment of the specification also provides a data processing device which corresponds to part of the process shown in fig. 1 and is based on the flow field simulation robustness calculation under the enhanced inclined grid.
As shown in fig. 5, a data processing apparatus in the present specification based on enhanced under-diagonal grid flow field simulation robustness calculation may include one or more of the following modules:
the model acquisition module 500 is configured to: obtaining a model of a target body, wherein the model represents an internal interface and an outlet boundary surface of the target body; at least part of the flow field formed by the model is divided into a plurality of grid cells; one part of the grid cells is an inclined grid, and the other part is a non-inclined grid;
a setting module 502 configured to: setting an initial velocity field and an initial pressure field;
a calculation model generation module 504 of the target velocity field is configured to: correcting the inclined grid as a target, and performing momentum interpolation on the initial velocity field to obtain a calculation model of a target velocity field represented by a target pressure field;
The pressure correction equation establishment module 506 is configured to: establishing a pressure correction equation by adopting a calculation model of the target speed field;
the pressure correction intermediate value determination module 508 is configured to: calculating the pressure correction equation to obtain a pressure correction intermediate value;
the pressure correction target value determination module 510 is configured to: interpolation processing is carried out on the pressure correction intermediate value belonging to the outlet boundary surface, so as to obtain a pressure correction target value;
a target speed field determination module 512 configured to: obtaining a target speed field output by a calculation model of the target speed field based on a target pressure field obtained by correcting the target value according to the pressure;
a data processing module 514 configured to: and carrying out data processing based on the target speed field and the target pressure field.
In an alternative embodiment of the present description, the calculation model of the target speed field is formulated as:
wherein: f represents the value of the quantity to which the interface belongs on an interface, wherein the interface comprises the internal interface and/or the outlet boundary surface; the upper dash represents the average value;representing the value of the pressure gradient at the interface; />Representing the known coefficients; / >Representing a normal direction of the interface; />Representing a target speed field obtained by the iterative computation; />An average value of the historical target speed field obtained by the previous iterative calculation is represented; p represents the target pressure field.
In an alternative embodiment of the present disclosure, the calculation model generation module 504 of the target velocity field is specifically configured to: the value of the pressure gradient on the interface is calculated by the following formula; wherein the target cell is any one of the grid cells:
/> wherein: c represents the target unit; n represents an adjacent cell, wherein the adjacent cell is a grid cell adjacent to the target cell; />Representing the pressure value at the center of the adjacent cell; />A pressure value representing the center of the target cell; />Representation->Is a projection point of (2); />Representation->Is a projection point of (2); />Representing the pressure gradient at the centers of the adjacent cells; />A pressure gradient representing the center of the target cell; />A vector representing a center point of the adjacent cell; />A vector representing a center point of the target cell; />A vector representing an interface center point between the target cell and the neighboring cell; />And the first distance is the projection distance of the distance from the central point of the target unit to the central point of the interface in the normal direction, and the second distance is the projection distance of the distance from the central point of the adjacent unit to the central point of the interface in the normal direction.
In an alternative embodiment of the present disclosure, the calculation model generation module 504 of the target velocity field is specifically configured to: for a target unit, the value of the pressure gradient on the interface is calculated by the following formula:
wherein: />A direction vector representing the center of the target cell to the center of the adjacent cell. />
In an alternative embodiment of the present description, the formula of the pressure correction equation is expressed as:
in (1) the->Is a diffusion term; p' is the pressure correction intermediate value.
In an alternative embodiment of the present disclosure, the pressure correction target value determination module 510 is specifically configured to: performing interpolation processing on the pressure correction intermediate value belonging to the outlet boundary surface, and obtaining an interpolation calculation formula adopted when the pressure correction target value is obtained, wherein the interpolation calculation formula is as follows:
wherein d represents an interpolation distance from the center of the target unit to the center of the outlet boundary surface; />Representing the pressure correction target value; />A gradient field representing a pressure correction target value; />Is a pressure correction target value for the center of the known cell.
In an alternative embodiment of the present disclosure, the target speed field determination module 512 is specifically configured to: solving a calculation model of the target speed field based on the target pressure field to obtain a speed correction target value; and correcting the initial speed field by adopting the speed correction target value to obtain a target speed field.
In an alternative embodiment of the present disclosure, the target velocity field is calculated by the following formula:wherein: u (U) * Is the target velocity field; u (U) *-1 Is the historical target speed field obtained in the last iteration; u' is the speed correction target value.
In an alternative embodiment of the present specification, the apparatus further comprises a target pressure field determination module configured to: and correcting the initial pressure field by adopting the pressure correction target value to obtain a target pressure field.
In this specification an alternative practiceIn an embodiment, the target pressure field determination module is specifically configured to:wherein: p (P) * Is the target pressure field; p (P) *-1 Is the historical target pressure field obtained by the last iteration; p' is the pressure correction target value.
Fig. 6 is a schematic structural view of an electronic device according to an embodiment of the present application. Referring to fig. 6, at the hardware level, the electronic device includes a processor, and optionally an internal bus, a network interface, and a memory. The Memory may include a Memory, such as a Random-Access Memory (RAM), and may further include a non-volatile Memory (non-volatile Memory), such as at least 1 disk Memory. Of course, the electronic device may also include hardware required for other services.
The processor, network interface, and memory may be interconnected by an internal bus, which may be an ISA (Industry Standard Architecture ) bus, a PCI (Peripheral ComponentInterconnect, peripheral component interconnect standard) bus, or EISA (Extended Industry Standard Architecture ) bus, among others. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, only one bi-directional arrow is shown in FIG. 6, but not only one bus or type of bus.
And the memory is used for storing programs. In particular, the program may include program code including computer-operating instructions. The memory may include memory and non-volatile storage and provide instructions and data to the processor.
The processor reads the corresponding computer program from the nonvolatile memory to the memory and then runs the computer program to form a data processing method based on the flow field simulation robustness calculation under the enhanced inclined grid on a logic level. And the processor is used for executing the program stored in the memory and particularly used for executing any data processing method based on the flow field simulation robustness calculation under the enhanced inclined grid.
The data processing method based on the flow field simulation robustness calculation under the enhanced inclined grid disclosed in the embodiment shown in fig. 1 of the present application can be applied to a processor (namely, a deletion control module in the present specification) or implemented by the processor. The processor may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software. The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; but also digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (ApplicationSpecific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present application may be embodied directly in the execution of a hardware decoding processor, or in the execution of a combination of hardware and software modules in a decoding processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in a memory, and the processor reads the information in the memory and, in combination with its hardware, performs the steps of the above method.
The electronic device may also execute a data processing method based on flow field simulation robustness calculation under the enhanced inclined grid in fig. 1, and implement the functions of the embodiment shown in fig. 1, which is not described herein.
The embodiment of the application also provides a computer readable storage medium, which stores one or more programs, the one or more programs include instructions, which when executed by an electronic device including a plurality of application programs, enable the electronic device to execute a method executed by a data processing method based on flow field simulation robustness calculation under an enhanced diagonal grid in the embodiment shown in fig. 1, and is specifically used for executing any one of the foregoing data processing methods based on flow field simulation robustness calculation under an enhanced diagonal grid.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Disks (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises the element.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The foregoing is merely exemplary of the present application and is not intended to limit the present application. Various modifications and variations of the present application will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. which come within the spirit and principles of the application are to be included in the scope of the claims of the present application.
Claims (4)
1. The data processing method based on flow field simulation robustness calculation under the enhanced inclined grid is characterized by comprising the following steps:
obtaining a model of a target body, wherein the model represents an internal interface and an outlet boundary surface of the target body; at least part of the flow field formed by the model is divided into a plurality of grid cells; one part of the grid cells is an inclined grid, and the other part is a non-inclined grid;
setting an initial velocity field and an initial pressure field;
correcting the inclined grid as a target, and performing momentum interpolation on the initial velocity field to obtain a calculation model of a target velocity field represented by a target pressure field;
establishing a pressure correction equation by adopting a calculation model of the target speed field;
calculating the pressure correction equation to obtain a pressure correction intermediate value;
interpolation processing is carried out on the pressure correction intermediate value belonging to the outlet boundary surface, so as to obtain a pressure correction target value;
obtaining a target speed field output by a calculation model of the target speed field based on a target pressure field obtained by correcting the target value according to the pressure;
performing data processing based on the target speed field and the target pressure field;
Wherein, the formula of the calculation model of the target speed field is expressed as:
wherein: f represents the value of the quantity to which the interface belongs on an interface, wherein the interface comprises the internal interface and/or the outlet boundary surface; the upper dash represents the average value;representing the value of the pressure gradient at the interface; d (D) f Representing the known coefficients; />Representing a normal direction of the interface; />Representing a target speed field obtained by the iterative computation; />An average value of the historical target speed field obtained by the previous iterative calculation is represented; p represents the target pressure field; m represents the iteration step number;
the formula of the pressure correction equation is expressed as:
in the method, in the process of the invention,is a diffusion term; p' is a pressure correction intermediate value; u (U) * Is the initial pressure field.
2. The method of claim 1, wherein the value of the pressure gradient at the interface for a target cell is calculated by the following formula; wherein the target cell is any one of the grid cells:
wherein: c represents the target unit; n represents an adjacent cell, wherein the adjacent cell is a grid cell adjacent to the target cell; p (P) N Representing the pressure value at the center of the adjacent cell; p (P) C A pressure value representing the center of the target cell; r is (r) N′ R represents N Is a projection point of (2); r is (r) C′ R represents C Is a projection point of (2);representing the pressure gradient at the centers of the adjacent cells; />Representing the targetA pressure gradient at the center of the cell; r is (r) N A vector representing a center point of the adjacent cell; r is (r) C A vector representing a center point of the target cell; r is (r) f A vector representing an interface center point between the target cell and the neighboring cell; and a represents a first distance and a second distance, and then the first distance and the second distance are compared to take a minimum value, wherein the first distance is a projection distance of a distance from a central point of the target unit to a central point of an interface in a normal direction, and the second distance is a projection distance of a distance from a central point of an adjacent unit to the central point of the interface in the normal direction.
3. The method of claim 2, wherein the value of the pressure gradient at the interface for a target cell is calculated by the following formula:
wherein:a direction vector representing the center of the target cell to the center of the adjacent cell.
4. The method according to claim 2, wherein the interpolation process is performed on the pressure correction intermediate value belonging to the outlet boundary surface to obtain the pressure correction target value by using an interpolation calculation formula:
Wherein d represents an interpolation distance from the center of the target unit to the center of the outlet boundary surface; p (P) b ′ f Representing the pressure correction target value;a gradient field representing a pressure correction target value; p (P) C ' isThe pressure correction target value at the center of the cell is known.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310227478.7A CN115952752B (en) | 2023-03-10 | 2023-03-10 | Data processing method based on flow field simulation robustness calculation under enhanced inclined grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310227478.7A CN115952752B (en) | 2023-03-10 | 2023-03-10 | Data processing method based on flow field simulation robustness calculation under enhanced inclined grid |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115952752A CN115952752A (en) | 2023-04-11 |
| CN115952752B true CN115952752B (en) | 2023-09-12 |
Family
ID=87282821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310227478.7A Active CN115952752B (en) | 2023-03-10 | 2023-03-10 | Data processing method based on flow field simulation robustness calculation under enhanced inclined grid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115952752B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116245049B (en) * | 2023-05-11 | 2023-07-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Node type non-structural grid boundary correction method, device, equipment and medium |
| CN116384165B (en) * | 2023-06-05 | 2023-08-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Ultra-relaxation processing method and device for enhancing computing efficiency and robustness |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108319793A (en) * | 2018-02-12 | 2018-07-24 | 大连理工大学 | A kind of low speed based on CFD tilts into the water the analysis method of free movement |
| CN111814422A (en) * | 2020-09-15 | 2020-10-23 | 北京智芯仿真科技有限公司 | Interpolation matrix construction method and device of integrated circuit electromagnetic simulation multiple grid method |
| DE102021120727A1 (en) * | 2021-03-23 | 2022-09-29 | Huazhong University Of Science And Technology | Parallel calculation method and system for the numerical simulation of casting flow fields |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7107188B2 (en) * | 2003-01-08 | 2006-09-12 | Schlumberger Technology Corporation | Digital pressure derivative method and program storage device |
-
2023
- 2023-03-10 CN CN202310227478.7A patent/CN115952752B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108319793A (en) * | 2018-02-12 | 2018-07-24 | 大连理工大学 | A kind of low speed based on CFD tilts into the water the analysis method of free movement |
| CN111814422A (en) * | 2020-09-15 | 2020-10-23 | 北京智芯仿真科技有限公司 | Interpolation matrix construction method and device of integrated circuit electromagnetic simulation multiple grid method |
| DE102021120727A1 (en) * | 2021-03-23 | 2022-09-29 | Huazhong University Of Science And Technology | Parallel calculation method and system for the numerical simulation of casting flow fields |
Non-Patent Citations (1)
| Title |
|---|
| 非结构CFD并行网格变形算法及其应用;高翔;《中国博士学位论文全文数据库 (基础科学辑)》(第2期);A004-15 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115952752A (en) | 2023-04-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115952752B (en) | Data processing method based on flow field simulation robustness calculation under enhanced inclined grid | |
| RU2519331C2 (en) | Method of calculating physical value, numerical analysis method, programme for calculating physical value, numerical analysis programme, device for calculating physical value and numerical analysis device | |
| Braaten et al. | A study of recirculating flow computation using body-fitted coordinates: consistency aspects and mesh skewness | |
| CN110851929B (en) | Two-dimensional leaf-type optimization design method and device based on self-adaptive grid | |
| JP6049746B2 (en) | Method and apparatus for modeling fluid and system boundary interactions in fluid dynamic systems | |
| US20200104708A1 (en) | Training apparatus, inference apparatus and computer readable storage medium | |
| Patrikalakis et al. | Intersection problems | |
| CN109002632B (en) | Liquid cooling system simulation method and device | |
| CN110490302B (en) | Neural network compiling and optimizing method and device and related products | |
| Pringuey et al. | High order schemes on three-dimensional general polyhedral meshes—Application to the level set method | |
| US6256599B1 (en) | Method for the representation of physical phenomena extending in a bi- or tridimensional spatial domain through semistructured calculation grid | |
| CN116384165B (en) | Ultra-relaxation processing method and device for enhancing computing efficiency and robustness | |
| Smith | Enclosure methods for systems of polynomial equations and inequalities | |
| CN115563840B (en) | Simulation method and system for reducing cascade error of unit matrix and related equipment | |
| Luo | A finite volume method based on weno reconstruction for compressible flows on hybrid grids | |
| CN114611289A (en) | Method and device for determining reference line in automatic driving and electronic equipment | |
| CN114429534A (en) | Three-dimensional model processing method and device, electronic equipment and storage medium | |
| CN110717291A (en) | Welding structure deformation simulation method, device, equipment and storage medium | |
| Moore et al. | Efficient energy evaluations for active B-Spline/NURBS surfaces | |
| Zagitov et al. | Automatic Block-Structured Grid Generation in Turbo Machine Blade Passages by TurboR&D. Mesher Software | |
| CN117150951B (en) | Pump equipment three-dimensional flow field calculation acceleration method | |
| CN117669429B (en) | Fluid simulation method and device for rotary machine and computer storage medium | |
| Matsuno | Improvement and assessment of an arbitrary-high-order time-accurate algorithm | |
| Zhu et al. | NURBS Modeling and Curve Interpolation Optimization of 3D Graphics | |
| JP2009282748A (en) | Design support device |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |