CN115455853A - Three-order compact reconstruction method based on finite volume method on unstructured grid - Google Patents
Three-order compact reconstruction method based on finite volume method on unstructured grid Download PDFInfo
- Publication number
- CN115455853A CN115455853A CN202211076411.XA CN202211076411A CN115455853A CN 115455853 A CN115455853 A CN 115455853A CN 202211076411 A CN202211076411 A CN 202211076411A CN 115455853 A CN115455853 A CN 115455853A
- Authority
- CN
- China
- Prior art keywords
- grid
- reconstruction
- cell
- grid cell
- derivative
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 106
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000010354 integration Effects 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 2
- 210000004027 cell Anatomy 0.000 description 160
- 238000012360 testing method Methods 0.000 description 9
- 239000012530 fluid Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 210000003888 boundary cell Anatomy 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a three-order compact reconstruction method based on a finite volume method on an unstructured grid, which comprises the following steps: constructing an unstructured grid; reading coordinates of grid nodes in the unstructured grid and connection relations among grid units; calculating an integral average value of a reconstruction basis function in a grid unit according to the grid node coordinates; calculating an average flow field according to a finite volume method; calculating a reconstruction relation weight coefficient according to density distribution in a flow field; constructing a first derivative relation matrix and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and the average flow field; and constructing a reconstruction polynomial second-order derivative term solving equation set, solving to obtain a reconstruction polynomial second-order derivative term, bringing the reconstruction polynomial second-order derivative term into the relation matrix of the first-order derivative and the second-order derivative, solving to obtain a reconstruction polynomial first-order derivative term, and finishing reconstruction.
Description
Technical Field
The disclosure belongs to the technical field of computational fluid mechanics, and particularly relates to a three-order compact reconstruction method based on a finite volume method on an unstructured grid.
Background
Computational fluid mechanics is a method for obtaining an approximate numerical solution of a flow field by dispersing a fluid mechanics control equation set and determining solution conditions through a numerical method on a computer. The high-precision numerical value dispersion method has an important effect on improving the precision of numerical solution and further expanding the application range of computational fluid mechanics. Due to the advantages of processing complex geometric shapes, realizing flexible automatic generation and the like, the computational fluid mechanics method based on the unstructured grid has become a hot spot of current research. The finite volume method is a numerical method which is most widely applied in computational fluid mechanics due to the advantages of automatically meeting the flow conservation characteristic, the convection term windward characteristic and the like. However, the finite volume method widely applied to the unstructured grid at present only has second-order precision, and cannot meet the requirements of high-precision and high-resolution scenes such as pneumatic acoustics, vortex-dominated flow, direct numerical simulation of turbulence, large vortex simulation and the like, and the development of a high-order numerical method is an effective method for solving related problems. The following important challenges are faced in constructing a high-order finite volume method calculation method based on an unstructured grid:
unstructured grids do not have regular topological structures in three spatial dimensions, and a high-order method cannot be independently constructed in each dimension like structured grids. For the finite volume method, the core of the method for constructing high-precision numerical values on an unstructured grid is to construct a high-order reconstruction polynomial function of a reconstruction variable in a grid unit by using an integral average value of the reconstruction variable on a reconstruction template, wherein pioneering work is a k-order accurate reconstruction method, namely the integral average value of the high-order reconstruction function in the unit in the reconstruction template unit is equal to an average value obtained by calculation of the finite volume method, each coefficient of a k-order distribution function is solved by using a least square method, and in addition, a student popularizes a weighting basic non-oscillation reconstruction method on the structured grid to the unstructured grid. Because each grid unit can only provide a reconstruction condition for the conservation of the integral average value in the grid unit, and the number of undetermined coefficients of the reconstruction polynomial is rapidly increased along with the increase of dimensionality and order, the two reconstruction methods, namely a k-order accurate reconstruction method and a weighting basic flutter-free reconstruction method, have the problem of overlarge reconstruction template.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a three-order compact reconstruction method based on a finite volume method on an unstructured grid, which does not need to solve an implicit algebraic equation set in the reconstruction process, and has small calculated amount and high solving efficiency; in addition, boundary conditions are considered in the reconstruction process, and third-order reconstruction accuracy with the same boundary grid cells and internal grid cells is achieved.
In order to achieve the above object, the present disclosure provides the following technical solutions:
a three-order compact reconstruction method based on a finite volume method on an unstructured grid comprises the following steps:
s100: constructing an unstructured grid;
s200: reading coordinates of grid nodes in the unstructured grid and connection relations among grid units;
s300: calculating an integral average value of a reconstruction basis function in a grid unit according to the grid node coordinates;
s400: calculating an average flow field according to a finite volume method;
s500: calculating a reconstruction relation weight coefficient according to density distribution in a flow field;
s600: constructing a first derivative relation matrix and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and the average flow field;
s700: and constructing a reconstruction polynomial second-order derivative term solving equation set, solving to obtain a reconstruction polynomial second-order derivative term, bringing the reconstruction polynomial second-order derivative term into the relation matrix of the first-order derivative and the second-order derivative, solving to obtain a reconstruction polynomial first-order derivative term, and finishing reconstruction.
Preferably, in step S300, the integrated average of the reconstruction basis functions in the grid cell is calculated by gaussian integration:
wherein omega i Is the volume of grid cell i, N is the number of Gaussian integration points, ω j Is the gaussian integral weight of the jth gaussian integral point,for the ith reconstructed basis function of grid cell i,to reconstruct basis functionsIntegral mean, cor, within grid cell i i j Is the coordinate of the jth gaussian integration point in grid cell i.
Preferably, in step S500, the reconstruction relation weight coefficient is calculated by the following formula:
wherein,is the weight of the adjacent cell nb of the grid cell i, nc is the number of cells adjacent to the i-face of the grid cell, ISS i j The smoothness ISS of each grid cell is the smoothness of the jth adjacent grid cell of grid cell i i Calculated by the following formula:
wherein ISS i To be the degree of smoothness of the grid cell i,is the average density in the jth cell adjacent to grid cell i,is the average density within grid cell i.
Preferably, in step S600, the first derivative and the second derivative relation matrix are expressed as follows:
wherein,
wherein,the weight coefficients of the reconstructed relation in j cells adjacent to the grid cell i,the mean value of the reconstructed variable in cell i calculated for the finite volume method,average value, x, of reconstructed variables in j cells adjacent to cell i calculated by finite volume method i,j,c Is the abscissa, y, of the center of the jth grid cell adjacent to grid cell i i,j,c J is 1-Nc as the ordinate of the center of the j-th grid cell adjacent to grid cell i, nc is the number of grid cells adjacent to grid cell i, x i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of the grid cell i,is the integral average of the function f within grid cell i.
Preferably, in step S700, the equation system for solving the second derivative term of the reconstruction polynomial is expressed as follows:
wherein,
wherein,is the kth row vector of the relationship matrix AuxMat of the jth cell adjacent to grid cell i,is the kth row vector of the relationship matrix AuxMat of grid cell i,the k-th component of the relationship matrix AuxVec for the j-th cell adjacent to grid cell i,the k-th component, x, of the relationship matrix AuxVec for grid cell i i,j,c Is the abscissa, y, of the center of the jth cell adjacent to grid cell i i,j,c Is the ordinate, x, of the center of the jth cell adjacent to grid cell i i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of grid cell i. j ranges from 1 to Nc, and k ranges from 1 to 2.
The present disclosure further provides a third-order compact reconstruction apparatus based on the finite volume method on unstructured grid, including:
the unstructured grid generating module is used for generating unstructured grids;
the preprocessing module is used for reading node coordinates in the unstructured grid and the connection relation between grid units;
the integral calculation module is used for calculating an integral average value of the reconstruction basis function in the grid unit according to the grid node coordinates;
the finite volume method calculating module is used for calculating an average flow field according to a finite volume method;
the weight coefficient calculation module is used for calculating a weight coefficient of a reconstruction relation according to density distribution in a flow field;
the reconstruction relation matrix building module is used for building a first derivative and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and an average flow field obtained by a finite volume method;
and the third-order reconstruction polynomial solving module is used for solving an equation set according to the constructed reconstruction polynomial second-order derivative term, solving to obtain a reconstruction polynomial second-order derivative term, and substituting the reconstruction polynomial second-order derivative term into a first-order derivative and second-order derivative relation matrix to solve to obtain a reconstruction polynomial first-order derivative term.
Preferably, the third-order reconstruction polynomial solving module includes:
the second order derivative item solving submodule is used for solving an algebraic equation system of the reconstruction polynomial second order derivative item to obtain a reconstruction polynomial second order derivative item;
and the first derivative item solving submodule is used for substituting the reconstructed polynomial second derivative item into the first derivative and second derivative item relation matrix to solve and obtain the reconstructed polynomial first derivative item.
Compared with the prior art, the beneficial effect that this disclosure brought does:
1. the problem of large templates of a high-order precision finite volume method on an unstructured grid is solved;
2. the explicit solution of the three-order reconstruction polynomial is realized, and compared with an implicit solution method, the calculation amount is small and the solution efficiency is high;
3. boundary conditions are considered in the reconstruction process, and full-field consistent third-order precision is achieved.
Drawings
FIG. 1 is a flow chart of a method provided by an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of a two-dimensional unstructured grid;
FIG. 3 is a diagram showing the correspondence between the coordinates of Gaussian integration points in the actual unit and the coordinates of a reference unit;
FIG. 4 is a graph of the numerical accuracy results of the reconstruction method obtained by test function testing;
fig. 5 is a streamline distribution diagram of the square cavity drive flow Re =400 calculated by the method.
Detailed Description
Specific embodiments of the present disclosure will be described in detail below with reference to fig. 1 to 5. While specific embodiments of the disclosure are shown in the drawings, it should be understood that the disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
It should be noted that certain terms are used throughout the description and claims to refer to particular components. As one skilled in the art will appreciate, various names may be used to refer to a component. This specification and claims do not intend to distinguish between components that differ in name but not function. In the following description and in the claims, the terms "include" and "comprise" are used in an open-ended fashion, and thus should be interpreted to mean "include, but not limited to. The description which follows is a preferred embodiment of the disclosure, but is made for the purpose of illustrating the general principles of the disclosure and not for the purpose of limiting the scope of the disclosure. The scope of the present disclosure is to be determined by the terms of the appended claims.
To facilitate an understanding of the embodiments of the present disclosure, the following detailed description is to be considered in conjunction with the accompanying drawings, and the drawings are not to be construed as limiting the embodiments of the present disclosure.
In one embodiment, as shown in fig. 1, a third-order compact reconstruction method based on finite volume on unstructured grid method includes the following steps:
s100: constructing an unstructured grid;
s200: reading coordinates of grid nodes in the unstructured grid and connection relations among grid units;
s300: calculating an integral average value of a reconstruction basis function in a grid unit according to the grid node coordinates;
s400: calculating an average flow field according to a finite volume method;
s500: calculating a weight coefficient of a reconstruction relation according to density distribution in a flow field;
s600: constructing a first derivative relation matrix and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and the average flow field;
s700: and constructing a reconstruction polynomial second-order derivative term solving equation set, solving to obtain a reconstruction polynomial second-order derivative term, bringing the reconstruction polynomial second-order derivative term into the relation matrix of the first-order derivative and the second-order derivative, solving to obtain a reconstruction polynomial first-order derivative term, and finishing reconstruction.
The above embodiments constitute a complete technical solution of the present disclosure. Compared with the prior art, the technical scheme of the embodiment has the following advantages: 1. the problem of large templates of a high-order precision finite volume method on an unstructured grid is solved; 2. the explicit solution of the three-order reconstruction polynomial is realized, and compared with an implicit solution method, the calculation amount is small and the solution efficiency is high; 3. boundary conditions are considered in the reconstruction process, and full-field consistent third-order precision is achieved.
In another embodiment, in step S100, an unstructured grid, such as the unstructured grid shown in fig. 2, is generated by an unstructured grid generating program, and the grid includes two grid cells, five grid edges, one internal grid edge and four boundary grid edges, and four grid nodes.
In another embodiment, in step S200, grid node coordinates of the unstructured grid are read, and in the unstructured grid shown in fig. 2, grid node n is 1 -n 4 Are numbered 1-4 in sequence, and grid cell E 1 、E 2 Are numbered 1 and 2, respectively, and the grid boundary e 1 -e 5 The numbers of the two are 1 to 5 in sequence.
The coordinates of the grid nodes are:
n 1 :(0,0)
n 2 :(1,0)
n 3 :(0,1)
n 4 :(1,1)
the connection relationship of the grid cells includes:
internal grid edge:
2 3 1 2
boundary grid edge:
4 3 2 0
3 1 1 0
1 2 1 0
2 4 2 0
the unit connection relation rule is as follows: the first two numbers represent the numbers of two end points of the grid boundary, wherein the direction of the first point pointing to the second point is the direction of the grid boundary; the last two numbers represent the numbers of two grid cells connected by a grid boundary formed by the first two grid points, wherein the normal direction inside the grid boundary points to the grid cell represented by the third number, and the normal direction outside the grid boundary points to the grid cell represented by the fourth number. Note: the units outside the computation domain are denoted by 0.
In another embodiment, in step S300, according to the two-dimensional Taylor expansion:
wherein u is i (x, y) is a reconstruction polynomial function within grid cell i, (x) i,c ,y i,c ) Is the coordinate of the center of the grid cell i,to reconstruct a polynomial u i The order derivatives at the center of the grid cell.
Taking integral average value of the reconstruction polynomial (1) in a grid unit i to obtain:
wherein omega i Is the volume of the grid cell i,is the integral average of the reconstructed polynomial in grid cell i.
Substituting formula (1) into formula (2) yields:
Substituting formula (3) into formula (1) for u in formula (1) i (x i,c ,y i,c ) The reconstruction polynomial is obtained as follows:
wherein,is the integrated average of the reconstructed variable in grid cell i,for the l-th undetermined coefficient of the reconstruction polynomial in grid cell i,for the ith reconstruction basis function of the reconstruction polynomial in grid cell i,is the integrated average of the basis functions in grid cell i,is a zero mean basis function, x i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of the grid cell i,is the integral average of the function f within grid cell i.
The reconstruction basis function in the grid unit has the specific form:
the integral average of the reconstructed basis functions in each grid cell is calculated by gaussian integration:
wherein omega i Is the volume of grid cell i, N is the number of Gaussian integration points, ω j Is the gaussian integral weight of the jth gaussian integral point,for the ith reconstructed basis function of grid cell i,is its integral mean, cor, within grid cell i i j Is the j-th Gaussian integration point coordinate in the grid unit i.
For the reference triangle composed of vertices 1' (0,0), 2' (1,0), 3' (0,1) given in fig. 3, the mapping relationship between the actual mesh cells and the reference cells is as follows:
wherein (x) 1 ,y 1 ),(x 2 ,y 2 ),(x 3 ,y 3 ) The physical coordinates of the three vertices 1,2,3 in the real cell, (x ', y') are the coordinates of the gaussian integration points a ', b', c ', d' in the reference cell, and (x, y) are the coordinates of the corresponding gaussian integration points a, b, c, d in the real cell.
The third-order gaussian integration point coordinates, i.e. the integration weights, in the reference unit are shown in table 1:
TABLE 1
In another embodiment, in step S400, the method of calculating the average flow field by using the finite volume method is a common method, and will not be described in detail, and the integral average value in the grid cell can be calculated by using the finite volume method
In another embodiment, in step S500, the reconstruction relation weight coefficient is calculated by the following formula:
wherein,is the weight of the neighboring cell nb of the grid cell i, nc is the number of cells adjacent to the face of the grid cell i,the smoothness ISS of each cell is the smoothness of the jth adjacent cell of the grid cell i i The calculation formula is as follows:
wherein ISS i To be the degree of smoothness of the unit i,is the integrated average of the density in the jth cell adjacent to cell i,is the integrated average of the density within cell i.
In another embodiment, in step S600, the integral average value of the cells adjacent to the i-plane of the grid cell is conserved, that is:
wherein nb is the number of the cell adjacent to the i-plane of the grid cell, Ω nb The volume of nb cells adjacent to the i-plane of the grid cell,the integral average of the reconstructed variable in nb cells adjacent to the i-plane of the grid cell.
Based on the weight coefficients and the integrated average of the reconstructed basis functions, equation (10) may be further expressed as:
wherein,is the integrated average of the reconstructed variable in grid cell i,for the integrated average of the reconstructed variable in grid cell nb adjacent to grid cell i,for the l-th undetermined coefficient of the reconstruction polynomial in grid cell i, (x) i,c ,y i,c ) Is the center coordinate of grid cell i, (x) nb,c ,y nb,c ) Is the center coordinate of the grid cell nb,as the integrated average of the function f within grid cell i,is the integral average of the function f within the grid cell nb.
Scanning all grid boundaries according to the connection relation among the grid units read in the step S200 to obtain an algebraic equation set constructed by a reconstructed variable integral mean conservation relation, wherein the algebraic equation set comprises:
wherein,
wherein,the weight coefficients of the reconstructed relation in j cells adjacent to the grid cell i,the mean value of the reconstructed variable in cell i calculated for the finite volume method,average value, x, of reconstructed variables in j cells adjacent to cell i calculated by finite volume method i,j,c Is the abscissa, y, of the center of the jth grid cell adjacent to grid cell i i,j,c J is 1-Nc as the ordinate of the center of the j-th grid cell adjacent to grid cell i, nc is the number of grid cells adjacent to grid cell i, x i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of the grid cell i,is the integral average of the function f within grid cell i.
For the boundary cell, the boundary condition also needs to be considered in the reconstruction process, and the general boundary condition is in the form as follows:
wherein u is i,b Is the value, c, of the reconstructed variable u at the boundary point b on the boundary cell i 1 ,c 2 ,c 3 The control parameter for the general boundary condition is determined according to the type of the boundary condition, for example, if a reconstruction variable U = U is specified on a boundary, c is determined in the reconstruction process of the cell i adjacent to the boundary 1 =0,c 2 =1,c 3 =U。
Since the reconstruction relation is provided through the boundary condition in the reconstruction process, the order of the reconstruction polynomial in the boundary unit does not need to be reduced, and therefore the third-order reconstruction precision with the same boundary unit and internal unit can be realized.
Based on the algebraic equation system constructed as above, the first and second derivatives relationship is expressed as follows:
wherein,
wherein,the weight coefficients for the reconstruction relationship in the j cells adjacent to grid cell i,the mean value of the reconstructed variable in cell i calculated for the finite volume method,average value, x, of reconstructed variables in j cells adjacent to cell i calculated by finite volume method i,j,c Is the abscissa, y, of the center of the jth grid cell adjacent to grid cell i i,j,c Is a grid celli the vertical coordinate of the center of the j-th grid cell adjacent to the grid cell i, j is 1-Nc, nc is the number of grid cells adjacent to the grid cell i, and x i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of the grid cell i,is the integral average of the function f within grid cell i.
the first and second derivative relationships may also be expressed as:
for the traditional high-order finite volume method, reconstruction is carried out only by utilizing the integral average value conservation relation (10) of the reconstruction variable in the reconstruction unit, for example, in the example, the two-dimensional third-order reconstruction polynomial (4) has five undetermined coefficientsAt least five grid cells need to be applied with the conservation relation (10) of integral mean value of the reconstruction variable, but for a two-dimensional triangular grid, each grid cell only has three adjacent grid cells, so that the information of non-surface adjacent grid cells is needed in the reconstruction process, namely the reconstruction process of the traditional high-order finite volume method is non-compact. The present invention provides not only the conservation condition of the integral average value of equation (10) but also the conservation condition of the first derivative of the reconstruction polynomial in the face-adjacent cells, so that the third order reconstruction can be completed only by using the information of the face-adjacent cells, and the specific steps are as described in the following embodiments.
In another embodiment, in step S700, an algebraic equation system for solving the second-order derivative term is first constructed from the first-order derivative conservation relation of the reconstructed polynomial.
The first derivative of the reconstruction polynomial (4) for x is in the form:
the conservation relation (10) of the integral average of the analog reconstruction variables can obtain the conservation relation of the reconstruction polynomial to the first derivative of x:
namely:
wherein nb is a cell number adjacent to the mesh cell i,for the jth undetermined coefficient of the reconstruction polynomial in unit i, (x) i,c ,y i,c ) Is the center coordinate of grid cell i, (x) nb,c ,y nb,c ) Is the center coordinate of the nb cell adjacent to grid cell i.
In conjunction with the first derivative and second derivative relationship matrix (15) obtained in step S600, equation (18) can be expressed as:
wherein,reconstruct the jth row vector of the matrix AuxMat for cell i,the jth component of the matrix AuxVec is reconstructed for cell i.
Due to the fact thatAndall the unknown, so the equation system formed by the formula (19) usually needs to solve all the units in the whole field simultaneously, namely implicit solution, the solution method has large calculation amount and low parallel solution efficiency, and in order to solve the problem, the reconstruction solution process is explicitly realized by the following method.
The second derivative term in equation (1) is approximated by a Taylor expansion, for example:
wherein,for the second cross derivative of the reconstructed polynomial at the center of the grid cell i,reconstructing the second cross derivative of the polynomial for the cell nb centers adjacent to grid cell i, (x) i,c ,y i,c ) Is the center coordinate of grid cell i, (x) nb,c ,y nb,c ) As the center coordinate of the nb cell adjacent to grid cell i, (x) nb,c -x i,c ) And (y) nb,c -y i,c ) Is the grid size h.
Namely:
similarly, for the other two second derivative terms:
substituting formula (21) and formula (22) into formula (1) yields:
in grid cell i (x-x) i,c )∝h,(y-y i,c ) Oc h, therefore, expanding equation (23) can be written as:
in order to remember that,
is a primitive polynomial u i Approximation of (x, y), equation (24) can be written as:
i.e. byApproximation u i (x, y), then still with third order reconstruction accuracy. Therefore, it can be concluded that the coefficients of the highest order terms of the reconstruction polynomials in adjacent grids are made equal without affecting the reconstruction accuracy.
The highest order terms, i.e. second derivative terms, of neighboring reconstruction polynomials are made equal, i.e.,
equation (19) can then be:
similarly, a first derivative conservation relation of the reconstruction polynomial to y is obtained:
scanning all internal grid edges to obtain an algebraic equation system satisfied by a second derivative term:
wherein,
wherein,is the kth row vector of the relationship matrix AuxMat of the jth cell adjacent to grid cell i,is the kth row vector of the relationship matrix AuxMat of grid cell i,the kth component of the relationship matrix AuxVec for the jth cell adjacent to grid cell i,the k-th component, x, of the relationship matrix AuxVec for grid cell i i,j,c Is the abscissa, y, of the center of the jth cell adjacent to grid cell i i,j,c Is the ordinate, x, of the center of the jth cell adjacent to grid cell i i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of grid cell i. j ranges from 1 to Nc, and k ranges from 1 to 2.
Solving the algebraic equation set to obtain second derivative termSubstituting into the relation matrix of the first derivative and the second derivative obtained in step S600 to obtain the first derivative term of the reconstruction polynomialAnd completing the reconstruction.
Next, the flow field calculation results based on the method of the present disclosure are given.
Dividing a 10 × 10, 14 × 14, 20 × 20, 28 × 28, 40 × 40 unstructured grid in a computational domain (0,1) × (0,1) and using an exponential function f 1 =x 2 e y As test function 1, with complex functionsAs a test function 2, the reconstruction method is used for reconstructing the test function, the difference between a reconstructed variable value and an accurate value reconstructed at the edge center of a grid is used as a reconstruction error, and the method is defined according to the error:
where Error is the numerical Error, h is the grid size, and n is the precision order.
As can be seen from the above equation, the slope of the logError-high curve represents the accuracy of the numerical method. Fig. 4 is a result diagram of numerical accuracy of a reconstruction method obtained through a test function test, in fig. 4, an average slope of a logError-high curve reconstructed by the method of the present disclosure is 3.01, that is, the method of the present disclosure does have third-order accuracy, compared with a classical second-order center method, in fig. 4, an average slope of a logError-high curve reconstructed by a second-order center format is 2.00, that is, the second-order center format has second-order reconstruction accuracy, and for a test function 1, to achieve a reconstruction error of the reconstruction method of the present disclosure on a 10 × 10 grid, the second-order center method requires a grid number of 71 × 71, and for a test function 2, to achieve a reconstruction error of the reconstruction method of the present disclosure on a 10 × 10 grid, the second-order center method requires a grid number of 26 × 26, so the method of the present disclosure has higher accuracy compared with the classical second-order center method, and the grid number required for achieving the same calculation residual error is less.
The reconstruction method is utilized to carry out numerical calculation on the incompressible square cavity driving flow calculation example under the condition that the Reynolds number (Re) is 400, 80 multiplied by 80 unstructured grids are divided into a calculation domain (0,1) × (0,1), four edges of a square cavity are non-slip boundaries, the upper boundary is a driving edge, and the speed distribution in the x direction is set as: u (x, 1) =16x 2 (1-x) 2 The other boundary speeds are all 0. Fig. 5 is a distribution diagram of a square cavity driving flow Re =400 streamline calculated by the method, and it can be seen that a large vortex in the middle of a flow field and small vortices at two base angles can be clearly captured by applying the method of the present invention. Table 2 shows the comparison of the vortex center coordinates calculated by the method with the literature results, as follows:
TABLE 2
The calculation result of the spectrum type high-precision algorithm for the calculation example shows that the vortex center position of the middle large vortex is in the range of (0.565-0.578,0.616-0.629), the vortex center position of the left bottom angle small vortex is in the range of (0.038-0.048,0.038-0.043), and the vortex center position of the right bottom angle small vortex is in the range of (0.897-0.902,0.103-0.114). In summary, it can be seen that the present disclosure can better identify the shape and location of vortices in a separation flow.
In conclusion, the three examples with different complexities verify the accuracy and effectiveness of the method of the third-order compact reconstruction method based on the finite volume method on the unstructured grid, and have feasibility and necessity of popularization and application.
The technical solutions provided by the present disclosure are described in detail with reference to specific embodiments, and the description of the embodiments is only for assisting understanding of the core ideas of the present disclosure, and a person skilled in the art may change the specific implementation and application scope according to the ideas of the present disclosure. Accordingly, the description should not be construed as limiting the disclosure.
Claims (7)
1. A three-order compact reconstruction method based on a finite volume method on an unstructured grid comprises the following steps:
s100: constructing an unstructured grid;
s200: reading coordinates of grid nodes in the unstructured grid and connection relations among grid units;
s300: calculating an integral average value of a reconstruction basis function in a grid unit according to the grid node coordinates;
s400: calculating an average flow field according to a finite volume method;
s500: calculating a weight coefficient of a reconstruction relation according to density distribution in a flow field;
s600: constructing a first derivative relation matrix and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and the average flow field;
s700: and constructing a reconstruction polynomial second-order derivative term solving equation set, solving to obtain a reconstruction polynomial second-order derivative term, bringing the reconstruction polynomial second-order derivative term into the relation matrix of the first-order derivative and the second-order derivative, solving to obtain a reconstruction polynomial first-order derivative term, and finishing reconstruction.
2. The method according to claim 1, wherein preferably, in step S300, the integrated average of the reconstruction basis functions in the grid cell is calculated by gaussian integration:
wherein omega i Is the volume of grid cell i, N is the number of Gaussian integration points, ω j Is the gaussian integral weight of the jth gaussian integral point,for the ith reconstructed basis function of grid cell i,to reconstruct basis functionsIntegral mean, cor, within grid cell i i j Is the coordinate of the jth gaussian integration point in grid cell i.
3. The method according to claim 1, wherein in step S500, the reconstruction relation weight coefficient is calculated by:
wherein,is the weight of the neighboring cell nb of the grid cell i, nc is the number of cells adjacent to the face of the grid cell i,the smoothness ISS of each grid cell is the smoothness of the jth adjacent grid cell of grid cell i i Calculated by the following formula:
4. The method according to claim 1, wherein in step S600, the first derivative and second derivative relation matrices are represented as follows:
wherein,
wherein,the weight coefficients of the reconstructed relation in j cells adjacent to the grid cell i,the mean value of the reconstructed variable in cell i calculated for the finite volume method,average value, x, of reconstruction variables in j cells adjacent to cell i calculated for finite volume method i,j,c Is the abscissa, y, of the center of the jth grid cell adjacent to grid cell i i,j,c For the ordinate of the center of the jth grid cell adjacent to grid cell i, j takes 1-Nc, nc being the number of grid cells adjacent to grid cell i, x i,c Is the abscissa, y, of the center of the grid cell i i,c Is the ordinate of the center of the grid cell i,is the integral average of the function f in the grid cell i;
5. The method of claim 1, wherein the solving of the equation set by the reconstruction polynomial second derivative term in step S700 is represented as follows:
wherein,
wherein,is the kth row vector of the relationship matrix AuxMat of the jth cell adjacent to grid cell i,for the kth row vector of the relation matrix AuxMat of grid cell i,the k-th component of the relationship matrix AuxVec for the j-th cell adjacent to grid cell i,the k-th component, x, of the relationship matrix AuxVec for grid cell i i,j,c Is the abscissa, y, of the center of the jth cell adjacent to grid cell i i,j,c Is the ordinate, x, of the center of the jth cell adjacent to grid cell i i,c Is the abscissa, y, of the center of the grid cell i i,c The value range of j is 1-Nc, and the value range of k is 1-2.
6. A three-order compact reconstruction device based on finite volume method on unstructured grid comprises:
the unstructured grid generating module is used for generating unstructured grids;
the preprocessing module is used for reading node coordinates in the unstructured grid and the connection relation between grid units;
the integral calculation module is used for calculating an integral average value of the reconstruction basis function in the grid unit according to the grid node coordinates;
the finite volume method calculating module is used for calculating an average flow field according to a finite volume method;
the weight coefficient calculation module is used for calculating a reconstruction relation weight coefficient according to density distribution in the flow field;
the reconstruction relation matrix building module is used for building a first derivative and a second derivative relation matrix based on the weight coefficient, the integral average value of the reconstruction basis function, the connection relation among the grid units and an average flow field obtained by a finite volume method;
and the third-order reconstruction polynomial solving module is used for solving an equation set according to the constructed reconstruction polynomial second-order derivative term, solving to obtain a reconstruction polynomial second-order derivative term, and substituting the reconstruction polynomial second-order derivative term into a first-order derivative and second-order derivative relation matrix to solve to obtain a reconstruction polynomial first-order derivative term.
7. The apparatus of claim 6, wherein the third order reconstruction polynomial solving module comprises:
the second order derivative item solving submodule is used for solving an algebraic equation system of the reconstruction polynomial second order derivative item to obtain a reconstruction polynomial second order derivative item;
and the first derivative item solving submodule is used for substituting the reconstructed polynomial second derivative item into the first derivative and second derivative item relation matrix to solve and obtain the reconstructed polynomial first derivative item.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211076411.XA CN115455853A (en) | 2022-09-02 | 2022-09-02 | Three-order compact reconstruction method based on finite volume method on unstructured grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211076411.XA CN115455853A (en) | 2022-09-02 | 2022-09-02 | Three-order compact reconstruction method based on finite volume method on unstructured grid |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115455853A true CN115455853A (en) | 2022-12-09 |
Family
ID=84300677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211076411.XA Pending CN115455853A (en) | 2022-09-02 | 2022-09-02 | Three-order compact reconstruction method based on finite volume method on unstructured grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115455853A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116822312A (en) * | 2023-08-30 | 2023-09-29 | 北京智芯微电子科技有限公司 | Post-processing display method, system and medium for high-order unstructured grid unit |
-
2022
- 2022-09-02 CN CN202211076411.XA patent/CN115455853A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116822312A (en) * | 2023-08-30 | 2023-09-29 | 北京智芯微电子科技有限公司 | Post-processing display method, system and medium for high-order unstructured grid unit |
CN116822312B (en) * | 2023-08-30 | 2024-01-26 | 北京智芯微电子科技有限公司 | Post-processing display method, system and medium for high-order unstructured grid unit |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110516388B (en) | Harmonic mapping-based curved surface discrete point cloud model circular cutter path generation method | |
Kukelova et al. | A clever elimination strategy for efficient minimal solvers | |
CN115455853A (en) | Three-order compact reconstruction method based on finite volume method on unstructured grid | |
Feng et al. | An adaptive sampling method for accurate measurement of aeroengine blades | |
Bai et al. | On the comparison of bilinear, cubic spline, and fuzzy interpolation techniques for robotic position measurements | |
Hanna et al. | Intersection of parametric surfaces by means of look-up tables | |
CN112862956A (en) | Human body and garment model collision detection and processing method based on HRBFs | |
CN111783343A (en) | Method for selecting direction template of unstructured grid and flow field solving method | |
CN116205153A (en) | Grid density optimization method based on flow field physical information | |
CN112001109A (en) | Method for realizing structural impact dynamics simulation by using regenerative nuclear particle algorithm | |
CN110134915B (en) | Method and device for processing magnetorheological polishing residence time | |
CN114347029B (en) | Model order reduction method for rapid simulation of pneumatic soft robot | |
CN112446163A (en) | Energy finite element topological optimization method based on parameterized level set | |
CN114996658A (en) | Projection-based hypersonic aircraft aerodynamic heat prediction method | |
CN114861315A (en) | Two-dimensional impeller profile optimization method based on machine learning | |
CN112214855A (en) | Planar region parameterization construction method based on frame field | |
CN111079353A (en) | Method and device for rapid and uniform grid division applied to complex fluid analysis | |
Wang et al. | Kirchhoff–Love shell representation and analysis using triangle configuration b-splines | |
CN115659619A (en) | Geometric topological optimization and additive manufacturing based integrated method | |
CN114399611A (en) | Polyhedral optimization method based on grid quality measurement characteristics | |
CN110555189A (en) | Spatial interpolation method based on reverse computing thinking | |
CN114332295A (en) | Method for mesh generation by using BP neural network | |
CN112632825B (en) | Electrostatic field smooth finite element numerical algorithm based on finite element super-convergence | |
CN102800100B (en) | Image segmentation method based on distance potential field and self-adaptive balloon force | |
CN113821885A (en) | Aero-engine blade flow field computing platform based on grid automatic generation technology |
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 |