CN106066912A - A kind of generation method of multi partition structured grid - Google Patents
A kind of generation method of multi partition structured grid Download PDFInfo
- Publication number
- CN106066912A CN106066912A CN201610372959.7A CN201610372959A CN106066912A CN 106066912 A CN106066912 A CN 106066912A CN 201610372959 A CN201610372959 A CN 201610372959A CN 106066912 A CN106066912 A CN 106066912A
- Authority
- CN
- China
- Prior art keywords
- grid
- face
- subregion
- structured
- current bay
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Generation (AREA)
Abstract
The invention discloses a kind of generation method of multi partition structured grid, step is: set up the set of all grids subregion to be generated；Determine starting sector；Six faces to current bay carry out the generation of structured grid respectively, it is thus achieved that the structured surface grid of current bay, and generate the structured bodies grid of current bay；Face F [0] according to current bay generating structure grid, finds the adjacent sectors at the place, face that with this face four summits coincide, generating structure grid on F [0] face of its adjacent sectors；Adjacent sectors is considered as current bay, repeats previous step, until partition boundaries；Sequentially above step is repeated in remaining five face of starting sector；Travel through all of all subregions, if all generating structure grids on all, proceed to next step；Otherwise, above step is repeated；Surface grids according to each subregion is internally generated structured bodies grid at its subregion.The present invention improves the quality of structured grid in whole computational fields subregion.
Description
Technical field
The present invention relates to a kind of computer-aided engineering (CAE), Fluid Mechanics Computation (CFD) and numerical simulation software skill
Art, the structured grid being specifically related to multi partition flow field in Fluid Mechanics Computation generates method.
Background technology
Grid generation technique is one of key technology in hydrodynamics numerical simulation, directly affects the convergence of numerical computations
Property, it is resolved that the final precision of numerical result and the efficiency of the process of calculating.Due to the flowing run in engineering and heat transfer
Problem mostly occurs in complex region, thus in irregular area, the generation of grid is Fluid Mechanics Computation and numerical heat transfer
In an important field of research.
For on the whole, flowing and heat transfer problem numerical computations in use grid be broadly divided into structured grid and
Unstructured grid, owing to structured grid has, mess generation speed is fast, quality good, data structure is simple, and to space
The region that curved surface uses the method for parametrization or spline interpolation to simulate is smooth, and the advantage such as actual model is easier access to,
Structured grid is widely adopted.
In actual numerical value simulation process, the mess generation for complex region is essential, as irregular in zoning,
Border excessively bends, smooth degree is poor, and CFD worker has developed multiblock technique and zoning method for calculating, adopts in this way
Time, first according to monnolithic case feature, whole domain being divided into a few sub regions, each subregion all uses the knot of routine
Structure grid comes discrete, and the discrete equation in usual each region all each solves respectively.Its advantage is reduction of mess generation
Difficulty；Different mesh-densities can be chosen in different regions；It is easy to use parallel algorithm to solve the algebraically side in each piece
Journey.
Using multi partition structured grid generation method when complex region is generated grid, it has a problem in that: (l) is to institute
The grid having segmentation block to generate is difficult to form unified numbering so that Numerical Simulation Program must be designed to carry out polylith meter
The form calculated, this can increase more programming effort amount undoubtedly；(2) even if multiple segmentation blocks being formed unified numbering and entering
Row overall calculation, but owing to the grid on interface between block and block often becomes the most rough, thus computational accuracy is had relatively
Big impact.Method set forth in the present invention has carried out Unified number to the grid surface of all subregions, and can make subregion interface
On the ground transition of grid continuous and derivable, thus improve the total quality of grid and the precision of numerical computations.
Summary of the invention
It is an object of the invention to provide a kind of generation method of multi partition structured grid, to solve existing multi partition structure
Change the rough problem of grid on grid generation technique caused subregion interface.
For achieving the above object, the present invention is by the following technical solutions:
A kind of generation method of multi partition structured grid, comprises the following steps:
Step one: set up the set of all grids subregion to be generated；
Step 2: determine a starting sector；
Step 3: six face F [i] to the current bay chosen, wherein i=0,1,2,3,4,5, carry out structuring respectively
The generation of grid, it is thus achieved that the structured surface grid of current bay, and generate the structured bodies grid of current bay；
Step 4: according to the face F [0] of current bay generating structure grid, finds four summits with this face mutually overlapping
The adjacent sectors at the place, face closed, utilizes unlimited interpolation method, generating structure grid on F [0] face of this adjacent sectors；
Step 5: adjacent sectors is considered as current bay, repeats step 4, until partition boundaries；
Step 6: sequentially remaining five face F [i] to starting sector, wherein i=1,2,3,4,5, repeat step 4 and step
Rapid five；
Step 7: travel through all of all subregions, if all generating structure grids on all, proceeds to step 8；No
Then, step 2 is repeated to step 6；
Step 8: utilize unlimited interpolation method, produces structured bodies net according to the surface grids of each subregion inside its subregion
Lattice, grid produces end-of-job.
Step 2 determining, the situation of starting sector is divided into following five kinds:
(1) there is structured bodies grid in certain subregion, then specifying this subregion is starting sector, performs step 4 to step
Rapid eight；
(2) there is structured grid on the face of certain subregion, then specifying this subregion is starting sector, performs step 3 extremely
Step 8；
(3) there is grid node on the limit of certain subregion, then specifying this subregion is starting sector, performs step 3 to step
Rapid eight；
(4) specified mesh definition on the limit of certain subregion, then specifying this subregion is starting sector, performs step 3 to step
Rapid eight；
(5) there is not any gridding information, then specifying first subregion in set is starting sector, performs step 2 extremely
Step 8；
The priority orders of these five kinds of situations is (1)～(5).
In step 4, on F [0] face of adjacent sectors, the step of generating structure grid is:
(1) search for six adjacent sectors of current bay, find four summits on F [0] face with current bay mutually overlapping
The subregion at the place, face closed is adjacent sectors to be looked for, if this adjacent sectors does not exists, proceeds to step 6；
(2) search for six faces of this adjacent sectors, find identical with the normal direction in current bay F [0] face and divide away from current
F [0] face in the face in district F [0] face, i.e. adjacent sectors, is designated as OPP [0] face；
(3) if not having grid on OPP [0] face, then unlimited interpolation method generating structure grid on this face is utilized；If OPP
[0] there is structured grid on face, then proceed to step 5.
In step 5, partition boundaries includes:
(1) current bay is positioned at the border of partitioned set；
(2), i.e. there is, on OPP [0] face, the grid generated in F [0] face of the adjacent sectors of current bay.
In the present invention, one of them child partition in computational fields subregion starts to generate grid, expands with spreading over a whole area from one point formula
To whole computational fields multi partition.
In the present invention, nodes and node arrangement mode in partition boundaries are inputted by program, it is achieved repair grid
Change and control.
The invention has the beneficial effects as follows: avoided by the method first expanded to subregion about again by child partition gridding
Grid interstitial content and the unmatched situation of grid distribution situation on adjacent subarea interface, it is ensured that net between each child partition
The nodal information of lattice block transmits like clockwork so that the grid on subregion interface can transition glossily, be greatly improved
The quality of structured grid in whole computational fields subregion.
Accompanying drawing explanation
Fig. 1 is the flow chart of the inventive method；
Fig. 2 is the schematic diagram determining starting sector；
Fig. 3 is the schematic diagram of the adjacent sectors corresponding with current bay F [0] face；
Fig. 4 is the schematic diagram again establishing current bay；
Fig. 5 a and 5b is the schematic diagram of the partition boundaries in step 5.
Detailed description of the invention
Below in conjunction with the accompanying drawings and detailed description of the invention, the present invention is done further description, make inventive feature and
Advantage is clearer.
It is illustrated in figure 1 the generation method of a kind of multi partition structured grid of the present invention, comprises the steps:
Step one: set up the set of all grids subregion to be generated；
Step 2: determine a starting sector.As in figure 2 it is shown, on A block, there is structured bodies grid, so determine A block
For starting sector；
Step 3: six face F [i] to the current bay chosen, wherein i=0,1,2,3,4,5, carry out structuring respectively
The generation of grid, it is thus achieved that the structured surface grid of current bay, and utilize unlimited interpolation method to generate the structured bodies of current bay
Grid, body grid on A block in Fig. 2, this step can be omitted；
Step 4: according to the face F [0] of current bay generating structure grid, it may be judged whether there are four with this face
The adjacent sectors at the place, face that summit coincides；
If not existing, the most as shown in Figure 5 a, current bay has been in the border of whole computational fields.At this moment, step 6 is proceeded to；
If existing, utilizing unlimited interpolation method, generating structure grid on F [0] face of its adjacent sectors, concrete operations are such as
Under:
(1) search for six adjacent sectors of current bay, find four summits on F [0] face with current bay mutually overlapping
The subregion at the place, face closed is adjacent sectors to be looked for；
(2) search for six faces of this adjacent sectors, find identical with the normal direction in current bay F [0] face and divide away from current
F [0] face in the face in district F [0] face, i.e. adjacent sectors, is designated as OPP [0] face；
(3) judge the most to have existed on OPP [0] face structured grid:
If there is no grid on OPP [0] face, then utilize unlimited interpolation method generating structure grid on this face.Such as Fig. 3 institute
Showing, the adjacent sectors at the place, face coincided with the summit, four, F [0] face of A block is B block, and F [0] face of B block is designated as OPP [0]
Face, this face does not has grid, so utilizing unlimited interpolation method generating structure grid on this face；
If there is structured grid on OPP [0] face, but the non-grid that step is generated before, then proceed to step 5.
Step 5: adjacent sectors is considered as current bay, repeats step 4, until partition boundaries.As shown in Figure 4, B block
It is considered as current bay, finds the adjacent sectors at the place, face coincided with summit, four, block F [0] face of B, i.e. C block, find the F of C block
[0] face, i.e. OPP [0] face, utilize unlimited interpolation method to produce grid on this OPP [0] face.
If there is the structured grid generated on OPP [0] face, as shown in Figure 5 b, C block is current bay, its adjacent point
District is D block, and now D block is in partition boundaries, utilizes unlimited interpolation method, in F [0] face of D block, i.e. generates net on OPP [0] face
Lattice.When D block is considered as current bay, its adjacent sectors is A block, and the F of A block [0] face, i.e. exist on OPP [0] face and generated
Structured grid, i.e. its grid has produced in A block, and now the mess generation in F [0] direction is complete, proceeds to step 6；
Step 6: sequentially remaining five face F [i] to starting sector, wherein i=1,2,3,4,5, repeat step 4 and step
Rapid five；
Step 7: travel through all of all subregions, if all generating structure grids on all, proceeds to step 8；No
Then, step 2 is repeated to step 6；
Step 8: utilize unlimited interpolation method, produces structured bodies net according to the surface grids of each subregion inside its subregion
Lattice, grid produces end-of-job.
The invention provides a kind of method that multi partition structured grid entirety generates, implement the side of this technical scheme
Method and approach are a lot, and the above provide only embodiment generally.Nodes in partition boundaries and node
Arrangement mode can be specified by program, it is achieved control and the amendment to grid.Ingredient the clearest and the most definite in the present embodiment
Realized by prior art.
Claims (6)
1. the generation method of a multi partition structured grid, it is characterised in that: comprise the following steps:
Step one: set up the set of all grids subregion to be generated；
Step 2: determine a starting sector；
Step 3: six face F [i] to the current bay chosen, wherein i=0,1,2,3,4,5, carry out structured grid respectively
Generation, it is thus achieved that the structured surface grid of current bay, and generate the structured bodies grid of current bay；
Step 4: according to the face F [0] of current bay generating structure grid, finds with this face four summits to coincide
The adjacent sectors at place, face, utilizes unlimited interpolation method, generating structure grid on F [0] face of this adjacent sectors；
Step 5: adjacent sectors is considered as current bay, repeats step 4, until partition boundaries；
Step 6: sequentially remaining five face F [i] to starting sector, wherein i=1,2,3,4,5, repeat step 4 and step
Five；
Step 7: travel through all of all subregions, if all generating structure grids on all, proceeds to step 8；Otherwise,
Repeat step 2 to step 6；
Step 8: utilize unlimited interpolation method, produces structured bodies grid, net according to the surface grids of each subregion inside its subregion
Lattice produce end-of-job.
The generation method of multi partition structured grid the most according to claim 1, it is characterised in that: step 2 has determined
The situation of beginning subregion is divided into following five kinds:
(1) there is structured bodies grid in certain subregion, then specifying this subregion is starting sector, performs step 4 to step
Eight；
(2) there is structured grid on the face of certain subregion, then specifying this subregion is starting sector, performs step 3 to step
Eight；
(3) there is grid node on the limit of certain subregion, then specifying this subregion is starting sector, performs step 3 to step
Eight；
(4) specified mesh definition on the limit of certain subregion, then specifying this subregion is starting sector, performs step 3 to step
Eight；
(5) there is not any gridding information, then specifying first subregion in set is starting sector, performs step 2 to step
Eight；
The priority orders of these five kinds of situations is (1)～(5).
The generation method of multi partition structured grid the most according to claim 1, it is characterised in that: in step 4, in phase
On F [0] face of adjacent subregion, the step of generating structure grid is:
(1) search for six adjacent sectors of current bay, find four summits on the F with current bay [0] face to coincide
The subregion at place, face is adjacent sectors to be looked for, if this adjacent sectors does not exists, proceeds to step 6；
(2) search for six faces of this adjacent sectors, find identical with the normal direction in current bay F [0] face and away from current bay F
[0] F [0] face in the face in face, i.e. adjacent sectors, is designated as OPP [0] face；
(3) if not having grid on OPP [0] face, then unlimited interpolation method generating structure grid on this face is utilized；If OPP [0] face
On there is structured grid, then proceed to step 5.
The generation method of multi partition structured grid the most according to claim 1, it is characterised in that: subregion limit in step 5
Boundary includes:
(1) current bay is positioned at the border of partitioned set；
(2), i.e. there is, on OPP [0] face, the grid generated in F [0] face of the adjacent sectors of current bay.
The generation method of multi partition structured grid the most according to claim 1, it is characterised in that: in computational fields subregion
One of them child partition start to generate grid, be extended to whole computational fields multi partition with spreading over a whole area from one point formula.
The generation method of multi partition structured grid the most according to claim 1, it is characterised in that: the joint in partition boundaries
Count and node arrangement mode is inputted by program, it is achieved amendment and the control to grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610372959.7A CN106066912A (en) | 2016-05-31 | 2016-05-31 | A kind of generation method of multi partition structured grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610372959.7A CN106066912A (en) | 2016-05-31 | 2016-05-31 | A kind of generation method of multi partition structured grid |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106066912A true CN106066912A (en) | 2016-11-02 |
Family
ID=57421530
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610372959.7A Pending CN106066912A (en) | 2016-05-31 | 2016-05-31 | A kind of generation method of multi partition structured grid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106066912A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108470081A (en) * | 2018-01-17 | 2018-08-31 | 天津大学 | A kind of supersonic boundary layer multi-Block Grid positioning and quick flow field interpolation method |
| CN109766578A (en) * | 2018-12-10 | 2019-05-17 | 南京航空航天大学 | A kind of three-dimensional wing icing rear surface reconstructing method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103838852A (en) * | 2014-03-13 | 2014-06-04 | 北京大学工学院南京研究院 | Method for rapidly searching for multi-block structured grid docking relations |
| CN104408773A (en) * | 2014-12-02 | 2015-03-11 | 北京大学工学院南京研究院 | Method for interpolating structured grid non-matching interface |
-
2016
- 2016-05-31 CN CN201610372959.7A patent/CN106066912A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103838852A (en) * | 2014-03-13 | 2014-06-04 | 北京大学工学院南京研究院 | Method for rapidly searching for multi-block structured grid docking relations |
| CN104408773A (en) * | 2014-12-02 | 2015-03-11 | 北京大学工学院南京研究院 | Method for interpolating structured grid non-matching interface |
Non-Patent Citations (4)
| Title |
|---|
| Y. SU 等: "Automatic hexahedral mesh generation for multi-domain composite models using a hybrid projective grid-based method", 《COMPUTER-AIDED DESIGN》 * |
| 姚彦忠 等: "复杂区域上的一种结构网格生成方法", 《计算物理》 * |
| 张永华: "叶轮机CFD网格生成", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 徐伟祖: "叶轮机CFD网格预处理研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108470081A (en) * | 2018-01-17 | 2018-08-31 | 天津大学 | A kind of supersonic boundary layer multi-Block Grid positioning and quick flow field interpolation method |
| CN108470081B (en) * | 2018-01-17 | 2021-10-12 | 天津大学 | Method for positioning multiple grids of supersonic velocity boundary layer and quickly interpolating flow field |
| CN109766578A (en) * | 2018-12-10 | 2019-05-17 | 南京航空航天大学 | A kind of three-dimensional wing icing rear surface reconstructing method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102306396B (en) | Three-dimensional entity model surface finite element mesh automatic generation method | |
| CN104036095B (en) | Coupling high accuracy complex appearance flow field fast algorithm based on Region Decomposition | |
| CN106650147A (en) | Continuum structure non-probability topologicaloptimization method based on bounded uncertainty | |
| CN104881582A (en) | Air quality prediction method and device | |
| US9183328B2 (en) | Method and apparatus for modeling interactions of the fluid with system boundaries in fluid dynamic systems | |
| CN102129517B (en) | High-precision two-phase fluid interface capturing method | |
| CN103106256A (en) | Gray model (GM) (1,1) prediction method of orthogonal interpolation based on Markov chain | |
| CN111523270B (en) | Improved topological optimization post-processing method for continuum structure | |
| CN105760602A (en) | Total flow field numerical simulation method for finite volume weighted essentially non-oscillatory scheme | |
| CN101329700B (en) | Method for simulating fluid flow | |
| EP2828777A1 (en) | Boundaries in fluid dynamic systems | |
| CN103838852A (en) | Method for rapidly searching for multi-block structured grid docking relations | |
| Williamschen et al. | Parallel anisotropic block-based adaptive mesh refinement algorithm for three-dimensional flows | |
| CN103116698A (en) | GM (1, 1) model prediction method based on cubic spline | |
| CN110990926B (en) | Urban surface building hydrodynamic simulation method based on area correction rate | |
| CN106066912A (en) | A kind of generation method of multi partition structured grid | |
| CN105677995A (en) | Method for numerical solution of fuzzy steady state heat conduction problem based on full grid point collocation theory | |
| CN106844889B (en) | Quadrilateral structure grid division method based on non-uniform block principle | |
| US9087165B2 (en) | Automatic extremum detection on a surface mesh of a component | |
| CN103970715A (en) | Mesh generation computing method and device | |
| CN104598600B (en) | A kind of parallel analysis of digital terrain optimization method based on distributed memory | |
| CN104992046A (en) | Computing system and method of fluid mechanics | |
| CN104809297A (en) | Electromagnetic force density transferring method used among special-shaped grids in magnetic field-structure field coupling calculation | |
| CN114595644A (en) | High-precision processing method for virtual hierarchical boundary in tree grid lattice boltzmann method | |
| CN111274739B (en) | Search method and system for computational fluid dynamics analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161102 |
|
| RJ01 | Rejection of invention patent application after publication |