CN107943747B - Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation - Google Patents
Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation Download PDFInfo
- Publication number
- CN107943747B CN107943747B CN201711143763.1A CN201711143763A CN107943747B CN 107943747 B CN107943747 B CN 107943747B CN 201711143763 A CN201711143763 A CN 201711143763A CN 107943747 B CN107943747 B CN 107943747B
- Authority
- CN
- China
- Prior art keywords
- temperature
- node
- heat conduction
- differential equation
- nodes
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004891 communication Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 2
- 238000004590 computer program Methods 0.000 abstract 1
- 238000000638 solvent extraction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000004088 simulation 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
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Feedback Control In General (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a method for automatically decomposing a multi-connected region based on a two-dimensional heat conduction differential equation. The method obtains the temperature field by utilizing the computer program, thereby automatically decomposing the multi-connected region, liberating a large amount of labor force, having applicability to any complex multi-connected region and having high efficiency.
Description
Technical Field
The invention relates to a method for automatically decomposing a multi-communication area based on a two-dimensional heat conduction differential equation, and belongs to the technical field of computational fluid mechanics numerical values.
Background
The first step of numerical simulation is grid division, a large number of multi-connected regions can appear in the practical application process, when the multi-connected regions are subjected to grid division, the regions need to be decomposed, and then each decomposed block is subjected to grid division. The quality of the partition directly affects the quality of the grid. The existing partitioning method is mainly manual partitioning and is decomposed according to the characteristic structure of the region. However, for different areas, the partitioning must be performed again, and this method will consume a lot of manpower and time in practical application, and the efficiency and processing capability cannot meet the requirement of real-time performance.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention aims to provide a method for automatically decomposing multiple connected regions based on a two-dimensional heat conduction differential equation, which can automatically divide any multiple connected region into regions and has the advantages of easiness in automatically processing complex connected regions, high efficiency and the like.
In order to achieve the above object, the present invention adopts the following technical solutions:
a method for automatically decomposing a multi-connected region based on a two-dimensional heat conduction differential equation is characterized by comprising the following steps:
step 1) carrying out discretization treatment on a multi-connected area needing to be partitioned;
step 2) setting a certain characteristic structure in the area as a constant heat source, wherein the temperature is constant T;
setting the initial temperature of other characteristic structures to be 0 ℃, and keeping the boundary temperature value to be 0;
step 3) based on two-dimensional steady-state heat conduction differential equation
Where t is temperature, τ is time, α is thermal diffusion coefficient, qvThe heat generated by an internal heat source is rho, the density and the specific heat are c; the temperature of each internal node can be obtained by a series expansion method or a thermal balance method, and the numerical solution of the two-dimensional steady-state heat conduction of the multi-communication area under the first class boundary condition is calculated to obtain the temperature distribution of the whole multi-communication area when the characteristic structure is used as a constant heat source;
step 4) repeating steps 2) and 3) until all the characteristic structures are used as constant heat sources to perform steps 2) and 3);
step 5) each node has different temperature values under different heat sources, and the connected area to which the node belongs is judged according to the maximum value of the node temperature;
the method for automatically decomposing the multi-connected area based on the two-dimensional heat conduction differential equation is characterized in that the specific content in the step 1) is as follows: a limited set of nodes is determined to replace the original continuous space, and the continuous space area is divided into a plurality of non-overlapping sub-areas.
The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation is characterized in that the sub-region division in the step 1) adopts a structured grid to solve: the solution area is divided into a plurality of sub-areas by a series of grid lines parallel to the coordinate axes, the intersection points of the grid lines are used as nodes, the positions of the nodes are represented by reference numerals X, Y of the nodes in two directions, the distance between two adjacent nodes is used as a step length and is recorded as delta X and delta Y, the delta X is the step length on the X axis, and the delta Y is the step length on the Y axis.
The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation is characterized in that the step 2) of setting a certain characteristic structure as a constant heat source means that the characteristic structure is mapped into a calculation domain, the characteristic structure is represented by a node, the position of the node is represented by reference numeral X, Y of the node in two directions, and the temperature of the nodes is set as a constant value T.
The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation is characterized in that the specific content in the step 3) is as follows: firstly, setting an error value epsilon according to an empirical value, carrying out iterative operation based on a two-dimensional heat conduction differential equation, wherein the temperature of each internal node is equal to the average value of the temperatures of four adjacent nodes at the upper, lower, left and right sides in the last iterative operation, the temperature value of a boundary node is constantly 0 ℃, comparing whether the error value of the temperature of each node and the temperature obtained in the last iterative operation is smaller than epsilon, if so, finishing the iterative operation, obtaining the temperature distribution of the whole multi-connected region when the characteristic structure is used as a constant heat source, and otherwise, continuing the iterative operation.
The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation is characterized in that the specific content in the step 5) is as follows: each node has different temperature values under different heat sources, and if the temperature value of the node under the Nth constant heat source is the maximum, the node belongs to the Nth communication area.
The invention achieves the following beneficial effects: the invention has applicability to any complex multi-connected area, is easy to process the complex connected area, has high efficiency, can meet the requirement of engineering, can automatically decompose the multi-connected area during processing, and liberates a large amount of labor force.
Drawings
FIG. 1 is a flow chart of an automatic decomposition of multiple connected regions into blocks method of the present invention;
FIG. 2 is a schematic view of the multiple connectivity regions of the present invention;
FIG. 3 is a schematic diagram of the computational domain discretization of the present invention;
fig. 4 (a), (b), and (c) are schematic diagrams of temperature distributions of regions in which three features are used as constant heat sources, respectively;
FIG. 5 is an exploded view of the multiple connectivity regions of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
The invention aims to solve the problems that when a multi-connected region is subjected to grid division, the multi-connected region needs to be partitioned, manual partitioning is long in time consumption, and workload is repeated and burdensome. With reference to the attached figure 1, the content of the method is as follows:
step 1) discretizing the multi-connected region shown in fig. 2, dividing the boundary line of the calculation domain into a plurality of sections by using a plurality of points, and connecting the points according to the shape of the calculation domain so that the calculation domain is divided into a plurality of sub-regions which are not overlapped with each other. The method adopts structured grid solving, namely a series of grid lines parallel to coordinate axes are used for dividing a solving area into a plurality of sub-areas, and intersection points of the grid lines are called nodes. The distance between two adjacent nodes is called step size and is marked as Δ X and Δ Y. In this example, as shown in fig. 3, the selected calculation field is a regular rectangular area, X is 100, and Y is 500.
And 2) taking the intersection points of the grid lines as space positions, namely nodes, of which the temperature values need to be determined. The position of the node is indicated by reference X, Y for the point in both directions. A feature structure is mapped into the calculation domain, the feature structure is represented by a plurality of nodes, the feature structure is also a regular rectangular area, coordinates of four vertexes of the feature structure are respectively (Xs, Ys) (Xe, Ye) (Xs, Ye), and Xe is larger than or equal to Xs, and Ye is larger than or equal to Yes. Setting the temperature value of the node representing the characteristic structure to be constant T, namely when X is larger than or equal to Xs, X is smaller than or equal to Xe, Y is larger than or equal to Ys, and Y is smaller than or equal to Ye, the values of the internal nodes are all constant T. The initial temperature values of the remaining nodes are all set to 0 degrees.
Step 3) first, an error value epsilon is set to 0.01. Based on two-dimensional steady-state heat conduction differential equation
For the internal nodes, the temperature of each internal node can be obtained by a series expansion method or a thermal equilibrium method
Wherein u [ x ]][y]Representing internal nodes [ x ] in the last iteration][y]Temperature value of, node x][y]Is [ x-1 ] as an adjacent node][y]、[x+1][y]、[x][y-1]、[x][y+1]；w[x][y]For the node [ x ] in the iteration operation][y]And for nodes on the boundary, its temperature value w x][y]Is always 0. Compare all nodes for | w [ x ]][y]-u[x][y]And if the value of | is less than epsilon, finishing the iterative operation, otherwise, continuing the iterative operation. The temperature distribution of the entire multiply connected region is finally obtained, as shown in fig. 4.
Step 4) repeating steps 2) and 3) until all the characteristic structures are used as constant heat sources to perform steps 2) and 3);
and 5) each node has different temperature values under different heat sources, and the connected area to which the node belongs is judged according to the maximum value of the temperature of the node, namely if the temperature value of the node under the Nth constant heat source is the maximum, the node belongs to the Nth connected area.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (6)
1. A method for automatically decomposing a multi-connected region based on a two-dimensional heat conduction differential equation is characterized by comprising the following steps:
step 1) carrying out discretization treatment on a multi-connected area needing to be partitioned;
step 2) setting a certain characteristic structure in the area as a constant heat source, wherein the temperature is constant T;
setting the initial temperature of other characteristic structures to be 0 ℃, and keeping the boundary temperature value to be 0;
step 3) based on two-dimensional steady-state heat conduction differential equation
Where t is temperature, τ is time, α is thermal diffusion coefficient, qvThe heat generated by an internal heat source is rho, the density is rho, the specific heat is c, and x and y are coordinates; the temperature of each internal node can be obtained by a series expansion method or a thermal balance method, and the numerical solution of the two-dimensional steady-state heat conduction of the multi-communication area under the first class boundary condition is calculated to obtain the temperature distribution of the whole multi-communication area when the characteristic structure is used as a constant heat source;
step 4) repeating steps 2) and 3) until all the characteristic structures are used as constant heat sources to perform steps 2) and 3);
and 5) each node has different temperature values under different heat sources, and the connected area to which the node belongs is judged according to the maximum value of the temperature of the node.
2. The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation as claimed in claim 1, wherein the specific contents in the step 1) are as follows: a limited set of nodes is determined to replace the original continuous space, and the continuous space area is divided into a plurality of non-overlapping sub-areas.
3. The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation as claimed in claim 2, wherein the sub-region division in the step 1) is solved by using a structured grid: the solution area is divided into a plurality of sub-areas by a series of grid lines parallel to the coordinate axes, the intersection points of the grid lines are used as nodes, the positions of the nodes are represented by reference numerals X, Y of the nodes in two directions, the distance between two adjacent nodes is used as a step length and is recorded as delta X and delta Y, the delta X is the step length on the X axis, and the delta Y is the step length on the Y axis.
4. The method according to claim 1, wherein the step 2) of setting a certain feature as a constant heat source means that the feature is mapped into a calculation domain, the feature is represented by a node, the position of the node is represented by reference numeral X, Y of the node in two directions, and the temperature of the nodes is set as a constant value T.
5. The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation as claimed in claim 1, wherein the specific contents in the step 3) are as follows: firstly, setting an error value epsilon according to an empirical value, carrying out iterative operation based on a two-dimensional heat conduction differential equation, wherein the temperature of each internal node is equal to the average value of the temperatures of four adjacent nodes at the upper, lower, left and right sides in the last iterative operation, the temperature value of a boundary node is constantly 0 ℃, comparing whether the error value of the temperature of each node and the temperature obtained in the last iterative operation is smaller than epsilon, if so, finishing the iterative operation, obtaining the temperature distribution of the whole multi-connected region when the characteristic structure is used as a constant heat source, and otherwise, continuing the iterative operation.
6. The method for automatically decomposing the multi-connected region based on the two-dimensional heat conduction differential equation as claimed in claim 1, wherein the specific contents in the step 5) are as follows: each node has different temperature values under different heat sources, and if the temperature value of the node under the Nth constant heat source is the maximum, the node belongs to the Nth communication area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711143763.1A CN107943747B (en) | 2017-11-17 | 2017-11-17 | Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711143763.1A CN107943747B (en) | 2017-11-17 | 2017-11-17 | Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107943747A CN107943747A (en) | 2018-04-20 |
| CN107943747B true CN107943747B (en) | 2021-01-08 |
Family
ID=61932715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711143763.1A Active CN107943747B (en) | 2017-11-17 | 2017-11-17 | Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107943747B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246764A (en) * | 2013-04-19 | 2013-08-14 | 哈尔滨工业大学 | Submersible electric machine with oil temperature field simulation method based on ansys workbench |
| CN106781187A (en) * | 2016-12-02 | 2017-05-31 | 武汉科技大学 | Scene of fire burning things which may cause a fire disaster area positioning method based on multidimentional system FM models |
| CN106886564A (en) * | 2017-01-03 | 2017-06-23 | 北京国能日新系统控制技术有限公司 | A kind of method and device that NWP wind energy collection of illustrative plates is corrected based on space clustering |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2639730C (en) * | 2007-09-21 | 2014-05-20 | Analytica Of Branford, Inc. | Improvements of finite differences methods |
| CN102129715B (en) * | 2011-03-24 | 2015-05-20 | 山东大学 | Generation method of quadrilateral grid of geometric model with any internal feature constraints |
| CN103324836B (en) * | 2013-05-31 | 2016-04-27 | 清华大学 | Based on the chip heat analysis method that 3D region decomposes |
| CN106874591B (en) * | 2017-02-10 | 2019-12-03 | 中冶华天南京工程技术有限公司 | A kind of calculation method of square billet heating process temperature distribution |
| CN107066737B (en) * | 2017-04-14 | 2019-05-17 | 北京科技大学 | A kind of two-dimentional staggered difference method for predicting hot rolling process plate belt temperature field |
-
2017
- 2017-11-17 CN CN201711143763.1A patent/CN107943747B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246764A (en) * | 2013-04-19 | 2013-08-14 | 哈尔滨工业大学 | Submersible electric machine with oil temperature field simulation method based on ansys workbench |
| CN106781187A (en) * | 2016-12-02 | 2017-05-31 | 武汉科技大学 | Scene of fire burning things which may cause a fire disaster area positioning method based on multidimentional system FM models |
| CN106886564A (en) * | 2017-01-03 | 2017-06-23 | 北京国能日新系统控制技术有限公司 | A kind of method and device that NWP wind energy collection of illustrative plates is corrected based on space clustering |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107943747A (en) | 2018-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Löhner | Matching semi-structured and unstructured grids for Navier-Stokes calculations | |
| CN111709171A (en) | Isogeometric solving and heat dissipation topology generation method for heat flow strong coupling problem | |
| CN111768502A (en) | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology | |
| CN108399290B (en) | Fracture flow MODFLOW simulation method based on fracture network connectivity | |
| CN112699623A (en) | High-precision heat flow calculation method based on unstructured grid regularized reconstruction technology | |
| CN111489447B (en) | Right-angle grid adaptive modeling method suitable for lattice Boltzmann method | |
| Krasnopolsky | The reordered BiCGStab method for distributed memory computer systems | |
| CN104182209A (en) | PETSc-based GCRO-DR algorithm parallel processing method | |
| CN111159939A (en) | Fin configuration topology optimization system and method | |
| CN102496168A (en) | Complicated river channel gridding method used for simulation of river channel hydrological numerical value | |
| Bahar et al. | Numerical solution of advection-diffusion equation using operator splitting method | |
| CN103729506B (en) | Complicated model complete hexahedron modeling and geometry remodeling and encryption method | |
| CN102902742A (en) | Spatial data partitioning method in cloud environment | |
| CN107153724B (en) | Chip temperature analysis method based on iterative algorithm | |
| CN107943747B (en) | Method for automatically decomposing multiple connected regions based on two-dimensional heat conduction differential equation | |
| CN111914447B (en) | Novel finite-volume multi-scale finite element method for simulating underground water solute transport | |
| CN102073796B (en) | Lattice traveling method for simulating solute three-dimensional transport process | |
| CN111339688B (en) | Method for solving rocket simulation model time domain equation based on big data parallel algorithm | |
| CN105224726A (en) | Structured grid Dynamic mesh is used for the method for unstrctured grid flow field calculation device | |
| CN105045767B (en) | A kind of method of immediate access and reading power system sparse matrix data | |
| Wang et al. | Pattern formation in a cross-diffusive Holling-Tanner model | |
| CN109299502B (en) | Two-dimensional numerical simulation method and system for continuous-discontinuous medium heat conduction | |
| CN111210522B (en) | Method for tracking streamline distribution in three-dimensional unstructured grid flow field by using FEM (finite element modeling) | |
| CN112364362A (en) | Parallel multilayer self-adaptive local encryption method facing fluid simulation direction | |
| CN111859646A (en) | Shock wave variable step length solving method based on B spline mapping function material point method |
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 |