CN102156806A - Method for quickly acquiring information of unstructured gaseous phase grid with solid particles - Google Patents
Method for quickly acquiring information of unstructured gaseous phase grid with solid particles Download PDFInfo
- Publication number
- CN102156806A CN102156806A CN2011100688677A CN201110068867A CN102156806A CN 102156806 A CN102156806 A CN 102156806A CN 2011100688677 A CN2011100688677 A CN 2011100688677A CN 201110068867 A CN201110068867 A CN 201110068867A CN 102156806 A CN102156806 A CN 102156806A
- Authority
- CN
- China
- Prior art keywords
- grid
- gas phase
- small
- scale
- scale structure
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for quickly acquiring information of an unstructured gaseous phase grid with solid particles. The method comprises the following steps: (1) establishing a large structured grid and a small structured grid to cover the unstructured gaseous phase grid; (2) scanning the gaseous phase grid, and carrying out statistic registration on the gaseous phase grid covered by each large structured grid; (3) establishing an array for the relationship of the gaseous phase grids when a minimum value is acquired in calculating the distance between the small structured grid and the gaseous phase grid thereof in the large structured grid, namely the gaseous phase grid of the small structure grid; and (4) in subsequent calculation, determining the small structured grid with the particles according to the particle positions, and acquiring the information of the particles in the unstructured gaseous phase grid by utilizing the array established in the step (3). By utilizing the method, while the numeric value simulation accuracy of gaseous phase and solid phase is ensured, and huge calculation amount caused by searching the gaseous phase grid with the particles in the traditional method is eliminated, so that the numeric value calculation efficiency is greatly improved.
Description
Technical field
The invention belongs to Gas-solid Two-phase Flow Computer Numerical Simulation field, be specifically related to a kind of method of obtaining solid phase particles non-structure gas phase gridding information of living in the gas-solid two-phase numerical simulation fast.
Background technology
Gas-solid Two-phase Flow extensively is present in fields such as the energy, chemical industry, environmental protection, and along with the significantly raising of computer speed, the computer numerical simulation of Gas-solid Two-phase Flow is also fast-developing thereupon, and has given play to more and more important effect in production reality.
In the Gas-solid Two-phase Flow Computer Numerical Simulation, a kind of effective measure is with gas phase field and solid phase particles field separate processes, adopt Euler's method and Lagrangian method to carry out numerical simulation respectively, can fully improve the accuracy of counting yield and increase numerical simulation like this.In the counteractive solution procedure, must obtain the information of particle local gas phase grid of living in to gas phase flow field calculating suffered gas drag force of discrete particle and particle, as local gas flow rate etc.For the non-structured grid of complex appearance solid, judge the residing current grid of particle how accurately, fast, be a difficult problem.Solve if employing is carried out the average method of data to a plurality of non-structure gas phase grids, then this process will be sacrificed a large amount of gas phase field information inevitably, thereby reduce computational accuracy; Yet take in each time step, all to find the solution particle distance between all non-structure gas phase grid element center points in the zoning, obtain the method for the residing current non-structure gas phase grid of particle, can bring extremely huge calculated amount again, thereby greatly reduce computing velocity.If can find a kind of method of determining solid phase particles place non-structure gas phase grid fast, the efficient of gas-solid two-phase numerical simulation is improved greatly.
Summary of the invention
Goal of the invention: at the deficiencies in the prior art, the purpose of this invention is to provide a kind of method of obtaining solid phase particles non-structure gas phase gridding information of living in the gas-solid two-phase numerical simulation fast, this method has been removed the huge calculated amount of bringing because of search particle non-structure gas phase grid of living in the classic method, thereby has improved the speed and the precision of gas-solid two-phase numerical evaluation greatly.
Technical scheme: in order to realize the foregoing invention purpose, the present invention adopts the technical scheme of the method for obtaining solid phase particles non-structure gas phase gridding information of living in fast as follows:
Step 1) is calculated in the non-structured grid of being divided in gas phase, sets up large-scale, small-sized two nested structure grids respectively whole zoning is covered; Wherein, the large scale structure mesh scale is greater than non-structure gas phase mesh scale, and the small-scale structure mesh scale is less than non-structure gas phase mesh scale;
Step 2) scanning all non-structure gas phase grids in the zoning, the sequence number of the non-structure gas phase grid that each large scale structure grid is covered is added up registration, sets up array in order to follow-up calculating use;
Step 3) is determined the large scale structure grid at each small-scale structure grid place, according to step 2) the middle array of setting up, find the solution the distance of interior all the non-structure gas phase grid element center points of central point large scale structure grid residing of small-scale structure grid with it, wherein minor increment person is the residing non-structure gas phase grid of this small-scale structure grid, afterwards the non-structure gas phase grid sequence number under each small-scale structure grid is added up registration, set up array and use in order to follow-up calculating;
Step 4) is determined its residing small-scale structure grid according to particle position in the calculating of each time step, thereby obtains the residing non-structure gas phase gridding information of this particle by the array of setting up in the step 3).
In step 1), the large scale structure mesh scale is greater than non-structure gas phase mesh scale, its yardstick is big more, then the large scale structure meshes number of Jian Liing is few more, the contained non-structure gas phase meshes number of single large scale structure grid is many more, thereby number of times of all non-structure gas phase grid element center point distances is just many more in the central point of finding the solution the small-scale structure grid in the step 3) large scale structure grid residing with it, in numerical evaluation, the large scale structure mesh scale should be soft according to the above principle and computing machine, hardware condition is set.
In step 1), the small-scale structure mesh scale is less than non-structure gas phase mesh scale, its yardstick is more little, then the small-scale structure meshes number of Jian Liing is many more, the number of the small-scale structure grid that single non-structure gas phase grid is contained is many more, thereby it is accurate more to obtain the residing non-structure gas phase gridding information of particle in the step 4), the number of times of finding the solution all non-structure gas phase grid element center point distances in the central point large scale structure grid residing with it of small-scale structure grid in the step 3) is just many more, in numerical evaluation, the small-scale structure mesh scale is soft according to the above principle and computing machine, hardware condition is set.
In step 3), in the time of near the small-scale structure grid is positioned at large scale structure grid edge, calculate for two dimension, the expanded range of finding the solution small-scale structure grid and non-structure gas phase grid element center point distance is to 8 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid.
In step 3), in the time of near the small-scale structure grid is positioned at large scale structure grid edge, for three-dimensional computations, the expanded range of finding the solution small-scale structure grid and non-structure gas phase grid element center point distance is to 26 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid.
Beneficial effect:The method of determining the non-structured grid of solid phase particles place gas phase fast of the present invention, non-structured grid be need not to average processing, but at first by cutting apart zoning, neighborhood search and calculating small-scale structure grid-non-structure gas phase grid element center point minimum distance method and determine the residing non-structure gas phase grid of small-scale structure grid, in iterative computation, directly obtain the residing non-structure gas phase gridding information of this particle afterwards, thereby strengthened the accuracy of numerical simulation by solid phase particles small-scale structure grid of living in; In addition, the present invention need not to find the solution particle distance between the non-structure gas phase grid element center point in the zoning, has therefore removed the calculated amount of bringing because of search particle non-structure gas phase grid of living in, from and improved the speed of calculating greatly.
Description of drawings
Fig. 1 is that gas phase is calculated the non-structured grid divided and and the large scale structure grid synoptic diagram that covered.
Fig. 2 is that gas phase is calculated the non-structured grid divided and and the small-scale structure grid synoptic diagram that covered.
Fig. 3 is the residing non-structure gas phase grid synoptic diagram of each small-scale structure grid of search before the iterative computation.
Fig. 4 determines this particle non-structure gas phase grid synoptic diagram of living in by the residing small-scale structure grid of discrete particle in the numerical evaluation.
Embodiment
The method of obtaining solid phase particles non-structure gas phase gridding information of living in the gas-solid two-phase numerical simulation fast of the present invention may further comprise the steps:
(1) calculates in the non-structured grid of being divided in gas phase, set up large-scale, small-sized two nested structure grids respectively whole zoning is covered.Wherein, the large scale structure mesh scale is greater than non-structure gas phase mesh scale, its yardstick is big more, then the large scale structure meshes number of Jian Liing is few more, the contained non-structure gas phase meshes number of single large scale structure grid is many more, thereby it is just many more hereinafter to find the solution the number of times of all non-structure gas phase grid element center point distances in the central point large scale structure grid residing with it of small-scale structure grid in the step (3).In addition, the small-scale structure mesh scale is less than non-structure gas phase mesh scale, its yardstick is more little, then the small-scale structure meshes number of Jian Liing is many more, the number of the small-scale structure grid that single non-structure gas phase grid is contained is many more, thereby hereinafter step (4) to obtain the residing non-structure gas phase gridding information of particle accurate more, but the number of times of hereinafter finding the solution all non-structure gas phase grid element center point distances in the central point large scale structure grid residing with it of small-scale structure grid in the step (3) is just many more;
(2) scanning all non-structure gas phase grids in the zoning, the sequence number of the non-structure gas phase grid that each large scale structure grid is covered is added up registration, sets up array in order to follow-up calculating use;
(3) determine the large scale structure grid at each small-scale structure grid place, according to the array of setting up in the step (2), find the solution the distance of interior all the non-structure gas phase grid element center points of central point large scale structure grid residing of small-scale structure grid with it, wherein minor increment person is the residing non-structure gas phase grid of this small-scale structure grid, non-structure gas phase grid sequence number under each small-scale structure grid is added up registration, set up array and use in order to follow-up calculating.
In the time of near the small-scale structure grid is positioned at large scale structure grid edge, calculate for two dimension, the expanded range of finding the solution small-scale structure grid and non-structure gas phase grid element center point distance is to 8 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid.And, then should expand 26 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid to for three-dimensional computations.
(4) in the calculating of each time step, determine its residing small-scale structure grid according to particle position, by step (3) thus in the array set up obtain the residing non-structure gas phase gridding information of this particle.
The present invention will be further explained below in conjunction with specific embodiment.
The numerical simulation of embodiment 1 gas phase flow field adopts two dimension or three-dimensional to carry out usually, now narrates with two-dimensional case, and three-dimensional case can the rest may be inferred:
(1) calculates the non-structured grid of being divided for gas phase, set up large-scale, small-sized two nested structure grids at first respectively whole zoning is covered.Wherein, the large scale structure mesh scale is greater than non-structure gas phase mesh scale, and the small-scale structure mesh scale is less than non-structure gas phase mesh scale, as shown in Figure 1, 2.Elliptic region is that the non-structured grid that gas phase calculating is divided is found the solution the zone, and large and small two kinds of squares are the two nested structure grid overlay areas of being set up; (2) scanning all non-structured grids in the zoning, the sequence number of the non-structure gas phase grid that each large scale structure grid is covered is added up registration, sets up array in order to follow-up calculating use; (3) determine the large scale structure grid at each small-scale structure grid place, according to the array of setting up in the step (2), find the solution the distance of interior all the non-structured grid central points of central point large scale structure grid residing of small-scale structure grid with it, wherein minor increment person is the residing non-structured grid of this small-scale structure grid, non-structure gas phase grid sequence number under each small-scale structure grid is added up registration, set up array and use in order to follow-up calculating; (4) consider when the small-scale structure grid is positioned near the large scale structure grid edge, its residing non-structure gas phase grid might produce in the adjacent mesh of large scale structure grid, as shown in Figure 3, the central point of small-scale structure grid a non-structure gas phase grid I of living in is arranged in the adjacent large scale structure grid B of this small-scale structure grid large scale structure grid of living in A.Therefore, need the required expanded range of particle and non-structure gas phase grid element center point distance of finding the solution in above-mentioned (3) to adjacent 8 the large scale structure grids (, expanding 26 adjacent large scale structure net regions to) of this structured grid for three-dimensional computations; (5) in the calculating of each time step, determine its residing small-scale structure grid according to particle position, by step (3) thus in the array set up obtain the residing non-structure gas phase gridding information of this particle.As shown in Figure 4, particle 1 is arranged in small-scale structure grid a, and small-scale structure grid a is arranged in non-structure gas phase grid I, thereby can determine that particle 1 is arranged in non-structure gas phase grid I.
The large scale structure grid that method utilization of the present invention is set up is determined small-scale structure grid location in space fast, and makes the hunting zone to non-structure gas phase grid dwindle greatly.With 200,000 non-structure gas phase grid computing space is an example, if the large scale structure grid number that adopts is 2,000, on average, each large scale structure grid covers 100 non-structure gas phase grids, and then calculating each small-scale structure grid can be from 200 to the number of times of non-structured grid central point distance, reduce to for 000 time 900 times (reducing to 2,700 times under the three-dimensional situation).What is more important when judging the residing non-structure gas phase grid of particle, can directly be determined according to the residing small-scale structure grid of this particle.That is to say, when adopting classic method, only be in single time step, 1 particle number of times that need search for non-structure gas phase grid of living in just reached 200,000 times.And adopt method of the present invention, this extremely huge calculated amount to save, thereby improved counting yield greatly.
Claims (5)
1. method of obtaining solid phase particles non-structure gas phase gridding information of living in fast, it is characterized in that this method may further comprise the steps: step 1) is calculated in the non-structured grid of being divided in gas phase, sets up large-scale, small-sized two nested structure grids respectively whole zoning is covered; Wherein, the large scale structure mesh scale is greater than non-structure gas phase mesh scale, and the small-scale structure mesh scale is less than non-structure gas phase mesh scale; Step 2) scanning all non-structure gas phase grids in the zoning, the sequence number of the non-structure gas phase grid that each large scale structure grid is covered is added up registration, sets up array in order to follow-up calculating use; Step 3) is determined the large scale structure grid at each small-scale structure grid place, according to step 2) the middle array of setting up, find the solution the distance of interior all the non-structure gas phase grid element center points of central point large scale structure grid residing of small-scale structure grid with it, wherein minor increment person is the residing non-structure gas phase grid of this small-scale structure grid, afterwards the non-structure gas phase grid sequence number under each small-scale structure grid is added up registration, set up array and use in order to follow-up calculating; Step 4) is determined its residing small-scale structure grid according to particle position in the calculating of each time step, thereby obtains the residing non-structure gas phase gridding information of this particle by the array of setting up in the step 3).
2. the method for obtaining solid phase particles non-structure gas phase gridding information of living in fast according to claim 1, it is characterized in that in the step 1), the large scale structure mesh scale is greater than non-structure gas phase mesh scale, its yardstick is big more, then the large scale structure meshes number of Jian Liing is few more, the contained non-structure gas phase meshes number of single large scale structure grid is many more, thereby number of times of all non-structure gas phase grid element center point distances is just many more in the central point of finding the solution the small-scale structure grid in the step 3) large scale structure grid residing with it, in numerical evaluation, the large scale structure mesh scale should be soft according to the above principle and computing machine, hardware condition is set.
3. the method for obtaining solid phase particles non-structure gas phase gridding information of living in fast according to claim 1, it is characterized in that: in the step 1), the small-scale structure mesh scale is less than non-structure gas phase mesh scale, its yardstick is more little, then the small-scale structure meshes number of Jian Liing is many more, the number of the small-scale structure grid that single non-structure gas phase grid is contained is many more, thereby it is accurate more to obtain the residing non-structure gas phase gridding information of particle in the step 4), the number of times of finding the solution all non-structure gas phase grid element center point distances in the central point large scale structure grid residing with it of small-scale structure grid in the step 3) is just many more, in numerical evaluation, the small-scale structure mesh scale is soft according to the above principle and computing machine, hardware condition is set.
4. the method for obtaining solid phase particles non-structure gas phase gridding information of living in fast according to claim 1, it is characterized in that: in the step 3), in the time of near the small-scale structure grid is positioned at large scale structure grid edge, calculate for two dimension, the expanded range of finding the solution small-scale structure grid and non-structure gas phase grid element center point distance is to 8 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid.
5. the method for obtaining solid phase particles non-structure gas phase gridding information of living in fast according to claim 1, it is characterized in that: in the step 3), in the time of near the small-scale structure grid is positioned at large scale structure grid edge, for three-dimensional computations, the expanded range of finding the solution small-scale structure grid and non-structure gas phase grid element center point distance is to 26 adjacent large scale structure grids of the residing large scale structure grid of this small-scale structure grid.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100688677A CN102156806A (en) | 2011-03-22 | 2011-03-22 | Method for quickly acquiring information of unstructured gaseous phase grid with solid particles |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100688677A CN102156806A (en) | 2011-03-22 | 2011-03-22 | Method for quickly acquiring information of unstructured gaseous phase grid with solid particles |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102156806A true CN102156806A (en) | 2011-08-17 |
Family
ID=44438300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100688677A Pending CN102156806A (en) | 2011-03-22 | 2011-03-22 | Method for quickly acquiring information of unstructured gaseous phase grid with solid particles |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102156806A (en) |
-
2011
- 2011-03-22 CN CN2011100688677A patent/CN102156806A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104853434A (en) | Indoor positioning method based on SVM and K mean value clustering algorithm | |
CN103247041B (en) | A kind of dividing method of the cloud data of the many geometric properties based on local sampling | |
CN107392875A (en) | A kind of cloud data denoising method based on the division of k neighbours domain | |
CN101620270B (en) | Wireless location method based on cluster-fusion | |
CN105547301A (en) | Indoor map construction method and device based on geomagnetism | |
CN103136393B (en) | A kind of areal coverage computing method based on stress and strain model | |
CN105549005B (en) | A kind of dynamic object direction of arrival tracking based on mesh generation | |
CN104869541A (en) | Indoor positioning tracking method | |
CN107703480A (en) | Mixed kernel function indoor orientation method based on machine learning | |
CN103220777A (en) | Mobile device positioning system | |
CN104038901A (en) | Indoor positioning method for reducing fingerprint data acquisition workload | |
CN104933291B (en) | Altimetry Data mean recovery time production method based on net―function | |
CN106646453A (en) | Doppler radar target tracking method based on predicted value measurement conversion | |
CN105101090B (en) | A kind of node positioning method of environmental monitoring wireless sense network | |
CN102325370A (en) | High-precision three-dimensional positioner for wireless sensor network node | |
CN104036095A (en) | Regional-decomposition based high-precision coupling fast-calculation method for complex-shape flow field | |
CN111551895A (en) | Motion radiation source TDOA and FDOA positioning method based on weighted multidimensional scale and Lagrange multiplier technology | |
CN104182571A (en) | Kriging interpolation method based on Delaunay and GPU | |
CN109661030A (en) | Unknown object location algorithm in wireless sensor network based on dynamic grid | |
CN105334497A (en) | Three-dimensional positioning method and device | |
CN104535067A (en) | Method for quickly calculating arrival time of pulse signal based on sector search | |
CN103260238A (en) | Mobile node positioning method based on speed sampling zone in wireless sensor network | |
CN101964024B (en) | Method for fast determining gas phase unstructured grid of solid particle | |
Wang et al. | A bluetooth location method based on kNN algorithm | |
Guo et al. | A localization algorithm for underwater wireless sensor networks based on ranging correction and inertial coordination |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110817 |