[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111008489B - Grid unit particle number optimization method in rarefied airflow numerical simulation - Google Patents

Grid unit particle number optimization method in rarefied airflow numerical simulation Download PDF

Info

Publication number
CN111008489B
CN111008489B CN202010159460.4A CN202010159460A CN111008489B CN 111008489 B CN111008489 B CN 111008489B CN 202010159460 A CN202010159460 A CN 202010159460A CN 111008489 B CN111008489 B CN 111008489B
Authority
CN
China
Prior art keywords
particles
particle
set range
grid
grid unit
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
Application number
CN202010159460.4A
Other languages
Chinese (zh)
Other versions
CN111008489A (en
Inventor
赵杰
陈灏
蒋光南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai suochen Information Technology Co., Ltd
Original Assignee
Shanghai Suochen Information Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Suochen Information Technology Co ltd filed Critical Shanghai Suochen Information Technology Co ltd
Priority to CN202010159460.4A priority Critical patent/CN111008489B/en
Publication of CN111008489A publication Critical patent/CN111008489A/en
Application granted granted Critical
Publication of CN111008489B publication Critical patent/CN111008489B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for optimizing the particle number of grid units in rarefied airflow numerical simulation, which comprises the following steps: in the process of carrying out numerical simulation on the rarefied airflow by adopting the DSMC, judging whether the number of simulation particles in any grid unit is within a set range, and if the number of the simulation particles is within the set range, optimizing the number of particles in the grid unit is not needed; if the number of the simulation particles exceeds the upper limit value of the set range, repeatedly executing particle fusion operation according to particle fusion criteria so that the number of the simulation particles in each grid unit is in the set range; if the number of the simulated particles is lower than the lower limit value of the set range, the particle separation operation is repeatedly executed according to the particle separation criterion, so that the number of the simulated particles in each grid cell is in the set range. And controlling the particle number of each grid cell in the analog domain, so that the simulated particle number in each grid cell does not change greatly along with the cell density and the cell size, and balancing the distribution of the simulated particles on different grid cells.

Description

Grid unit particle number optimization method in rarefied airflow numerical simulation
Technical Field
The invention relates to the technical field of rarefied airflow numerical simulation, in particular to a grid unit particle number optimization method in rarefied airflow numerical simulation.
Background
In recent years, the method of numerically simulating the lean airflow using the monte carlo Direct Simulation Method (DSMC) has attracted more and more attention due to the failure of the continuous medium assumption of the N-S equation in the lean airflow region.
DSMC is a discrete particle simulation method of probability theory. The DSMC method uses a small number of simulated particles to represent a large number of real gas molecules in a computer, the space coordinates, the speed, the internal energy and the like of the simulated particles are stored in the computer and change along with time due to the movement, the collision and the interaction with a boundary of the particles, and a molecular speed distribution function is restored by sampling the simulated particles.
Generally, for a lean gas stream supersonic flow, the fluid density increases significantly on the shock wave; furthermore, if the wall temperature is low, the fluid density will increase further towards the boundary layer region. Assuming uniform grid cells are used in the simulation, the number of particles per cell in the simulation will increase from free-flow cells to boundary layer cells due to the significant increase in flow density from free-flow to boundary layer regions. Thus, if the free-flow cell has about 10-20 particles, the cells near the boundary layer may have more than several hundred particles. As previously mentioned, this would make the simulation time long and inefficient. In contrast, if this number is too small, the simulation results may be inaccurate due to the inability to solve for collisions in each analog cell, according to physical problems.
Disclosure of Invention
Aiming at the problems and the defects in the prior art, the invention provides a grid unit particle number optimization method in the rarefied airflow numerical simulation, which carries out self-adaptive adjustment and optimization on the number of simulated particles in a grid unit, thereby obviously improving the precision and the efficiency of the rarefied airflow numerical calculation.
The invention solves the technical problems through the following technical scheme:
the invention provides a method for optimizing the particle number of grid cells in rarefied airflow numerical simulation, which is characterized by comprising the following steps of:
in the process of carrying out numerical simulation on the rarefied airflow by adopting the DSMC, judging the relation between the number of simulation particles in any grid unit and a set range to judge whether to execute particle fusion operation or particle separation operation, if the number of the simulation particles is in the set range, optimizing the number of the particles of the grid unit is not needed, namely, the particle fusion operation or the particle separation operation is not executed; if the number of the simulation particles exceeds the upper limit value of the set range, repeatedly executing particle fusion operation according to particle fusion criteria so that the number of the simulation particles in each grid unit is in the set range; if the number of the simulated particles is lower than the lower limit value of the set range, the particle separation operation is repeatedly executed according to the particle separation criterion, so that the number of the simulated particles in each grid cell is in the set range.
Preferably, in the numerical simulation of the lean airflow by adopting the DSMC, the particles with the particle mass larger than the average particle mass in any grid unit are defined as large particles, the particles with the particle mass smaller than or equal to the average particle mass are defined as small particles, the particles in the grid unit are numbered, and the number NB and the number NS of the large particles and the small particles in the grid unit are counted;
the particle fusion criterion is used for judging the relationship between the number of simulation particles containing large particles and small particles in the grid unit and a set range so as to judge whether to execute particle fusion operation, and comprises the following steps:
if the number NB + NS of the simulated particles exceeds the upper limit value of the set range, then:
a) if NB is larger than or equal to the upper limit value of the set range, adjusting the grid size to ensure that the number NB of large particles in the adjusted grid unit is smaller than the upper limit value of the set range;
a1) if the number of the simulated particles in the adjusted grid unit is within the set range, the particle fusion operation or the particle separation operation is not executed;
a2) if the number of the simulated particles in the adjusted grid unit still exceeds the upper limit value of the set range, judging whether the distance between any two small particles is smaller than the set distance and whether the particle fusion judgment condition is met or not in the adjusted grid unit
Figure 650387DEST_PATH_IMAGE001
If so, fusing two small particles into one large particle, otherwise, adjusting the grid size, renumbering the particles in the adjusted grid unit, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle fusion or separation operation again;
b) if NB is less than the upper limit value of the set range, in the grid unit, the distance between any two small particles is judged to be less than the set distance and whether the particle fusion judgment condition is satisfied
Figure 490167DEST_PATH_IMAGE002
If so, fusing the two small particles into one large particle, otherwise, adjusting the size of the grid, renumbering the particles in the grid unit after adjustment, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle fusion or separation operation again;
in steps a2) and b), judging whether the number of the simulated particles in the grid cell after the two small particle fusion operations are executed is within a set range every time the two small particle fusion operations are executed, and if not, continuing the two small particle fusion operations;
wherein,
Figure 908510DEST_PATH_IMAGE003
is the average mass of large particles around two small particles,
Figure 200951DEST_PATH_IMAGE004
is the average mass of the two small particles,
Figure 995731DEST_PATH_IMAGE005
is the number of large particles surrounding two small particles,
Figure 690018DEST_PATH_IMAGE006
is the number of small particles that are present,
Figure 279262DEST_PATH_IMAGE007
is an artificial threshold;
particle fusion follows conservation of mass and momentum:
Figure 324579DEST_PATH_IMAGE008
wherein m is1、m2Mass of two small particles, v1、v2The velocities of the two small particles are respectively, m and v are respectively the mass and velocity of the fused large particle, and the position of the fused large particle is the midpoint of the positions of the two small particles before fusion.
Preferably, the mean free path of small particles is set at a distance of 0.5 times.
Preferably, the determining the particle separation criterion determines the relationship between the number of simulated particles containing large particles and small particles in the grid cell and the set range to determine whether to perform the particle separation operation includes:
if the number NB + NS of the simulated particles is lower than the lower limit value of the set range:
c) if NB =0, adjusting the grid size of the grid unit so that the number NB of macro-particles in the adjusted grid unit is more than 0;
c1) if the number of the simulated particles in the adjusted grid unit is within the set range, the particle fusion operation or the particle separation operation is not executed;
c2) if the adjusted grid cellIf the number of middle simulation particles is still lower than the lower limit value of the set range, judging whether any single large particle meets the particle separation judgment condition or not in the adjusted grid unit
Figure 719788DEST_PATH_IMAGE009
If so, separating a single large particle into two small particles, otherwise, adjusting the grid size, and judging whether to execute particle fusion or separation operation;
d) if NB is greater than 0, in the grid unit, judging whether any single large particle meets the particle separation judgment condition
Figure 206264DEST_PATH_IMAGE009
If so, separating a single large particle into two small particles, otherwise, adjusting the grid size, and judging whether to execute particle fusion or separation operation;
in steps c2) and d), each time the macro-particle separating operation is executed, judging whether the number of the simulated particles in the grid unit after the macro-particle separating operation is executed is within a set range, and if not, continuing the macro-particle separating operation;
wherein,
Figure 294306DEST_PATH_IMAGE010
is the mass of the individual large particles,
Figure 499022DEST_PATH_IMAGE011
is the average mass of small particles surrounding a single large particle,
Figure 697922DEST_PATH_IMAGE012
is the number of large particles and is,
Figure 38905DEST_PATH_IMAGE013
is the number of small particles surrounding a single large particle,
Figure 297848DEST_PATH_IMAGE014
is an artificial threshold;
particle separation follows conservation of mass and momentum:
Figure 986931DEST_PATH_IMAGE015
wherein,
Figure 723943DEST_PATH_IMAGE016
the mass of the two separated small particles is respectively,
Figure 247328DEST_PATH_IMAGE017
respectively the speed of the two separated small particles,
Figure 614855DEST_PATH_IMAGE018
the mass and speed of the large particles before separation are respectively, and the positions of the two small particles after separation are as follows:
Figure 590901DEST_PATH_IMAGE019
wherein p is1、p2The positions of two small particles after separation are respectively shown, p is the position of a single large particle before separation, and h is a set threshold value.
Preferably, h is taken as the radius r of the separated small particles.
Preferably, the setting range is 10-20.
On the basis of the common knowledge in the field, the above preferred conditions can be combined randomly to obtain the preferred embodiments of the invention.
The positive progress effects of the invention are as follows:
1. and controlling the particle number of each grid cell in the analog domain, so that the simulated particle number in each grid cell does not change greatly along with the cell density and the cell size, and balancing the distribution of the simulated particles on different grid cells.
2. The calculation amount of numerical simulation is controlled, the simulation efficiency is improved, and the simulation precision is improved.
Drawings
FIG. 1 is a schematic view of particle fusion according to a preferred embodiment of the present invention.
FIG. 2 is a schematic view of particle separation according to a preferred embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
The embodiment provides a method for optimizing the particle number of grid cells in a rarefied airflow numerical simulation, which comprises the following steps of:
in the numerical simulation of the thin airflow by adopting the DSMC, judging the relationship between the number of simulation particles in any grid unit and a set range (such as 10-20) to judge whether to execute particle fusion operation or particle separation operation, if the number of the simulation particles is within the set range (10-20), optimizing the number of particles of the grid unit is not needed, namely, the particle fusion operation or the particle separation operation is not executed; if the number of the simulated particles exceeds the upper limit value (20) of the set range, repeatedly executing particle fusion operation according to particle fusion criteria so that the number of the simulated particles in each grid cell is within the set range (10-20); if the number of the simulated particles is lower than the lower limit value (10) of the set range, the particle separation operation is repeatedly executed according to the particle separation criterion so that the number of the simulated particles in each grid cell is within the set range (10-20).
In the numerical simulation of the thin gas flow by adopting the DSMC, the particles with the particle mass larger than the average particle mass in any grid unit are defined as large particles, the particles with the particle mass smaller than or equal to the average particle mass are defined as small particles, the particles in the grid unit are numbered, and the number NB and the number NS of the large particles and the small particles in the grid unit are counted.
(1) The particle fusion criterion is used for judging the relationship between the number of simulation particles containing large particles and small particles in the grid unit and a set range so as to judge whether to execute particle fusion operation, and comprises the following steps:
fusing small particles in a region with more large particles to obtain larger particles. When large particles are mostly present around two small particles and the two small particles are very close to each other, the particles are fused, the new fused particle mass is the sum of the two small particle masses, and the position is the midpoint of the positions of the two particles before fusion.
If the number NB + NS of the simulated particles exceeds the upper limit value of the set range, then:
a) if NB is larger than or equal to the upper limit value of the set range, adjusting the grid size to enable the number NB of large particles in the adjusted grid unit to be smaller than the upper limit value of the set range.
For example: if NB + NS in the grid cell is 24 (20 set upper limits), where NB is 22, the number NB of large particles in the grid cell exceeds 20 set upper limits, and at this time, the number of simulated particles in the grid cell cannot be within the set range even if the particle fusion operation is performed. At this time, the grid size of the grid unit needs to be adjusted, so that the number NB of large particles in the adjusted grid unit is less than the upper limit value 20 of the set range, and at this time, the number of simulated particles in the grid unit may be in the set range through the particle fusion operation.
a1) And if the number of the simulated particles in the adjusted grid unit is within the set range, not executing the particle fusion operation or the particle separation operation. For example: if the number NB of large particles in the adjusted grid cell is 18 and the number of small particles in the adjusted grid cell is 2, the number of simulated particles in the adjusted grid cell is 20, and is within the set range, and then it is not necessary to perform the particle fusion operation or the particle separation operation.
a2) If the number of the simulated particles in the adjusted grid cell still exceeds the upper limit value of the set range, judging whether the distance between any two small particles is smaller than the set distance (the small particle mean free path which is 0.5 times) and whether the particle fusion judgment condition is met in the adjusted grid cell
Figure 69287DEST_PATH_IMAGE020
If so, two granulesAnd (3) merging the particles into a large particle (see figure 1), otherwise, adjusting the grid size, renumbering the particles in the grid unit after adjustment, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle merging or separation operation again.
For example: if the number NB of large particles in the adjusted grid unit is 19 and the number of small particles is 2, the number 21 of simulated particles in the adjusted grid unit still exceeds the upper limit value 20 of the set range, at this time, it is determined whether the distance between any two small particles is smaller than the set distance and the particle fusion determination condition is satisfied, if so, two small particles are fused into one large particle, so that the number of simulated particles in the adjusted grid unit is 20 and is within the set range, and then, there is no need to perform the particle fusion operation or the particle separation operation.
b) If NB is less than the upper limit value of the set range, in the grid unit, the distance between any two small particles is judged to be less than the set distance and whether the particle fusion judgment condition is satisfied
Figure 650441DEST_PATH_IMAGE020
If so, fusing two small particles into one large particle (see fig. 1), otherwise, adjusting the grid size, renumbering the particles in the adjusted grid unit, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle fusion or separation operation again.
In steps a2) and b), each time two small particle fusion operations are executed, whether the number of simulated particles in the grid cell after the two small particle fusion operations are executed is within a set range is judged, and if not, the two small particle fusion operations are continued.
Wherein,
Figure 985608DEST_PATH_IMAGE003
is the average mass of large particles around two small particles,
Figure 714529DEST_PATH_IMAGE004
is the average mass of the two small particles,
Figure 731027DEST_PATH_IMAGE005
is the number of large particles surrounding two small particles,
Figure 229004DEST_PATH_IMAGE006
is the number of small particles that are present,
Figure 672755DEST_PATH_IMAGE007
is an artificial threshold (γ = 1).
Particle fusion follows conservation of mass and momentum:
Figure 888973DEST_PATH_IMAGE008
wherein m is1、m2Mass of two small particles, v1、v2The velocities of the two small particles are respectively, m and v are respectively the mass and velocity of the fused large particle, and the position of the fused large particle is the midpoint of the positions of the two small particles before fusion.
(2) The particle separation criterion is used for judging the relationship between the number of simulation particles containing large particles and small particles in the grid unit and a set range so as to judge whether to execute particle separation operation, and comprises the following steps:
when large particles are mixed in a small particle region, a separation mechanism opposite to fusion is added. When there are many cases of small particles around a single large particle, the large particle is divided into smaller particles. The velocity and position of the two small particles separated are determined from conservation of energy and conservation of momentum.
If the number NB + NS of the simulated particles is lower than the lower limit (10) of the set range:
c) if NB =0, adjusting the grid size of the grid unit so that the number of macro-particles NB in the adjusted grid unit is greater than 0. If NB =0, the particle separating operation cannot be performed, and the simulated particle data in the mesh cell cannot be made to be within the set range, so the particle separating operation is performed on the premise that NB > 0.
c1) And if the number of the simulated particles in the adjusted grid unit is within the set range, not executing the particle fusion operation or the particle separation operation.
c2) If the number of the simulated particles in the adjusted grid unit is still lower than the lower limit value of the set range, judging whether any single large particle meets the particle separation judgment condition or not in the adjusted grid unit
Figure 709161DEST_PATH_IMAGE009
If so, separating a single large particle into two small particles, otherwise, adjusting the grid size, and judging whether to execute particle fusion or separation operation;
d) if NB is greater than 0, in the grid unit, judging whether any single large particle meets the particle separation judgment condition
Figure 61645DEST_PATH_IMAGE009
If so, separating a single large particle into two small particles, otherwise, adjusting the grid size, and judging whether to execute particle fusion or separation operation;
in steps c2) and d), each time the macro-particle separating operation is executed, judging whether the number of the simulated particles in the grid unit after the macro-particle separating operation is executed is within a set range, and if not, continuing the macro-particle separating operation;
wherein,
Figure 4193DEST_PATH_IMAGE010
is the mass of the individual large particles,
Figure 379811DEST_PATH_IMAGE011
is the average mass of small particles surrounding a single large particle,
Figure 800428DEST_PATH_IMAGE012
is the number of large particles and is,
Figure 213611DEST_PATH_IMAGE013
is the number of small particles surrounding a single large particle,
Figure 61481DEST_PATH_IMAGE014
is an artificial threshold;
particle separation follows conservation of mass and momentum:
Figure 924395DEST_PATH_IMAGE015
wherein,
Figure 148703DEST_PATH_IMAGE016
the mass of the two separated small particles is respectively,
Figure 475779DEST_PATH_IMAGE017
respectively the speed of the two separated small particles,
Figure 697813DEST_PATH_IMAGE018
the mass and speed of the large particles before separation are respectively, and the positions of the two small particles after separation are as follows:
Figure 110339DEST_PATH_IMAGE021
wherein p is1、p2The positions of two small particles after separation, p the position of a single large particle before separation, h a set threshold, the size of h generally affects the distance between two small particles after separation, and is not recommended to be too large, and herein, the radius r of a small particle is taken.
The invention provides a grid unit particle number optimization method in rarefied airflow numerical simulation, which is used for adaptively adjusting and optimizing the number of simulation particles in each grid unit in a simulation domain, so that the number of the simulation particles in each grid unit is in a set range, the number of the simulation particles in each grid unit is not greatly changed along with the unit density and the unit size, the distribution of the simulation particles on different grid units is balanced, and the accuracy and the efficiency of rarefied airflow numerical calculation can be obviously improved.
While specific embodiments of the invention have been described above, it will be appreciated by those skilled in the art that these are by way of example only, and that the scope of the invention is defined by the appended claims. Various changes and modifications to these embodiments may be made by those skilled in the art without departing from the spirit and scope of the invention, and these changes and modifications are within the scope of the invention.

Claims (5)

1. A method for optimizing the particle number of grid cells in a rarefied airflow numerical simulation is characterized by comprising the following steps of:
in the process of carrying out numerical simulation on the rarefied airflow by adopting the DSMC, judging the relation between the number of simulation particles in any grid unit and a set range to judge whether to execute particle fusion operation or particle separation operation, if the number of the simulation particles is in the set range, optimizing the number of the particles of the grid unit is not needed, namely, the particle fusion operation or the particle separation operation is not executed; if the number of the simulation particles exceeds the upper limit value of the set range, repeatedly executing particle fusion operation according to particle fusion criteria so that the number of the simulation particles in each grid unit is in the set range; if the number of the simulated particles is lower than the lower limit value of the set range, repeatedly executing particle separation operation according to the particle separation criterion so as to enable the number of the simulated particles in each grid unit to be in the set range;
in the numerical simulation of the thin airflow by adopting the DSMC, defining the particles with the particle mass larger than the average particle mass in any grid unit as large particles and the particles with the particle mass smaller than or equal to the average particle mass as small particles, numbering the particles in the grid unit, and counting the number NB and NS of the large particles and the small particles in the grid unit;
the particle fusion criterion is used for judging the relationship between the number of simulation particles containing large particles and small particles in the grid unit and a set range so as to judge whether to execute particle fusion operation, and comprises the following steps:
if the number NB + NS of the simulated particles exceeds the upper limit value of the set range, then:
a) if NB is larger than or equal to the upper limit value of the set range, adjusting the grid size to ensure that the number NB of large particles in the adjusted grid unit is smaller than the upper limit value of the set range;
a1) if the number of the simulated particles in the adjusted grid unit is within the set range, the particle fusion operation or the particle separation operation is not executed;
a2) if the number of the simulated particles in the adjusted grid unit still exceeds the upper limit value of the set range, judging whether the distance between any two small particles is smaller than the set distance and whether the particle fusion judgment condition is met or not in the adjusted grid unit
Figure FDA0002460992140000011
If so, fusing two small particles into one large particle, otherwise, adjusting the size of the grid, renumbering the particles in the grid unit after adjustment, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle fusion or separation operation again;
b) if NB is less than the upper limit value of the set range, in the grid unit, the distance between any two small particles is judged to be less than the set distance and whether the particle fusion judgment condition is satisfied
Figure FDA0002460992140000021
If so, fusing the two small particles into one large particle, otherwise, adjusting the size of the grid, renumbering the particles in the grid unit after adjustment, counting the number NB of the large particles and the number NS of the small particles, and judging whether to execute particle fusion or separation operation again;
in steps a2) and b), judging whether the number of the simulated particles in the grid cell after the two small particle fusion operations are executed is within a set range every time the two small particle fusion operations are executed, and if not, continuing the two small particle fusion operations;
wherein m isbigIs the average mass of large particles around two small particles, msmaIs the average mass of two small particles, NbigIs the number of large particles surrounding two small particles, NsmaIs the number of small particles, γ is the artificial threshold;
particle fusion follows conservation of mass and momentum: m is m1+m2,m1=m2,v=v1=v2
Wherein m is1、m2The mass of each of the two small particles,v1、v2the velocities of the two small particles are respectively, m and v are respectively the mass and velocity of the fused large particle, and the position of the fused large particle is the midpoint of the positions of the two small particles before fusion.
2. The method for optimizing the particle count of grid cells in a lean flow numerical simulation of claim 1, wherein a mean free path of small particles is set at a distance of 0.5 times.
3. The method of claim 1, wherein the particle separation criterion for determining the relationship between the number of simulated particles containing large particles and small particles in the grid cell and a set range to determine whether to perform a particle separation operation comprises:
if the number NB + NS of the simulated particles is lower than the lower limit value of the set range:
c) if NB is 0, adjusting the grid size of the grid unit so that the number NB of large particles in the adjusted grid unit is more than 0;
c1) if the number of the simulated particles in the adjusted grid unit is within the set range, the particle fusion operation or the particle separation operation is not executed;
c2) if the number of the simulated particles in the adjusted grid unit is still lower than the lower limit value of the set range, judging whether any single large particle meets the particle separation judgment condition or not in the adjusted grid unit
Figure FDA0002460992140000031
If so, separating a single large particle into two small particles, otherwise, adjusting the grid size, and judging whether to execute particle fusion or separation operation;
d) if NB>0, judging whether any single large particle meets the particle separation judgment condition or not in the grid unit
Figure FDA0002460992140000032
If so, the single large particle is separated into two small particles, otherwise, the net is adjustedGrid size, judging whether to execute particle fusion or separation operation;
in steps c2) and d), each time the macro-particle separating operation is executed, judging whether the number of the simulated particles in the grid unit after the macro-particle separating operation is executed is within a set range, and if not, continuing the macro-particle separating operation;
wherein m'bigIs the mass of individual macroparticles, m'smaIs the average mass of small particles surrounding a single large particle, N'bigIs the number of macroparticles, N'smaIs the number of small particles surrounding a single large particle, γ' is the artificial threshold;
particle separation follows conservation of mass and momentum: m '═ m'1+m′2,m′1=m′2
Figure FDA0002460992140000033
Wherein m'1、m′2Respectively the mass v 'of the two separated small particles'1、v′2The speeds of the two small particles after separation are respectively, m 'and v' are respectively the mass and the speed of the large particle before separation, and the positions of the two small particles after separation are as follows:
Figure FDA0002460992140000034
Figure FDA0002460992140000035
Figure FDA0002460992140000036
wherein p is1、p2The positions of two small particles after separation are respectively shown, p is the position of a single large particle before separation, and h is a set threshold value.
4. The method of optimizing the particle count of grid cells in a lean flow numerical simulation of claim 3, wherein h is taken as the radius r of the separated small particles.
5. The method for optimizing the particle count of grid cells in a lean airflow numerical simulation of claim 1, wherein the setting range is 10 to 20.
CN202010159460.4A 2020-03-10 2020-03-10 Grid unit particle number optimization method in rarefied airflow numerical simulation Active CN111008489B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010159460.4A CN111008489B (en) 2020-03-10 2020-03-10 Grid unit particle number optimization method in rarefied airflow numerical simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010159460.4A CN111008489B (en) 2020-03-10 2020-03-10 Grid unit particle number optimization method in rarefied airflow numerical simulation

Publications (2)

Publication Number Publication Date
CN111008489A CN111008489A (en) 2020-04-14
CN111008489B true CN111008489B (en) 2020-07-14

Family

ID=70121018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010159460.4A Active CN111008489B (en) 2020-03-10 2020-03-10 Grid unit particle number optimization method in rarefied airflow numerical simulation

Country Status (1)

Country Link
CN (1) CN111008489B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107942309A (en) * 2017-10-19 2018-04-20 上海无线电设备研究所 Ultrahigh speed Electromagnetic Scattering of Target quick calculation method in a kind of rarefied atmosphere

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930102B (en) * 2012-10-31 2015-04-22 西安空间无线电技术研究所 Particle combination method in secondary electron multiplication simulation of microwave part
CN105574335A (en) * 2015-12-16 2016-05-11 武汉科技大学 Parallel stochastic simulation method for predicting concentration of atmospheric particulate matters
CN106650064B (en) * 2016-12-09 2019-07-26 华东师范大学 A kind of condensation phenomenon emulation mode based on particle model
CN110848065B (en) * 2019-11-13 2021-07-02 北京工业大学 Method for automatically recognizing diesel oil spray crushing process and automatically realizing spray continuous calculation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107942309A (en) * 2017-10-19 2018-04-20 上海无线电设备研究所 Ultrahigh speed Electromagnetic Scattering of Target quick calculation method in a kind of rarefied atmosphere

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
动态网格的DSMC方法在GPU上的并行;文敏华等;《计算机科学与探索》;20130219;摘要、第472页至479页 *
稀薄流到连续流的气体运动论统一算法研究;李志辉等;《空气动力学学霸》;20030928;第21卷(第3期);第255页志266页 *

Also Published As

Publication number Publication date
CN111008489A (en) 2020-04-14

Similar Documents

Publication Publication Date Title
Sørensen An analytic model of wind-blown sand transport
Krumholz et al. Embedding Lagrangian sink particles in Eulerian grids
Rice et al. An experimental study of multiple grain‐size ejecta produced by collisions of saltating grains with a flat bed
Wise et al. ENZO+ MORAY: radiation hydrodynamics adaptive mesh refinement simulations with adaptive ray tracing
Edwards et al. The surface statistics of a granular aggregate
Oh et al. Computations of high-speed, high Knudsen number microchannel flows
JP5280863B2 (en) Piecewise meshing and neighborhood search for object interaction simulation
Theofanous et al. Shock dispersal of dilute particle clouds
CN115793470B (en) Parameter control method of ore separator and ore separator
CN111008489B (en) Grid unit particle number optimization method in rarefied airflow numerical simulation
Pieri et al. Anisotropic galactic outflows and enrichment of the intergalactic Medium. I. Monte Carlo simulations
CN113406957B (en) Mobile robot autonomous navigation method based on immune deep reinforcement learning
Flomenbom Clustering in anomalous files of independent particles
CN109829975A (en) A kind of variable orifice porosity porous media construction method
CN105117537A (en) Weight comparison based resampling method for particle filter system
Chen et al. Research on an improved PSO algorithm with dual self-adaptation and dual variation
SCOTT et al. Comparison of limiters in flux-split algorithms for Euler equations
Joshi et al. Particle-resolved study of the onset of turbulence
CN115496133A (en) Density data stream clustering method based on self-adaptive online learning
CN107341849B (en) Rapid real-time smoke simulation algorithm
CN113570120A (en) Wind power prediction method and device based on improved particle swarm optimization
Li et al. The potential of the height-to-width ratio as an indicator to characterize the ventilation performance of a canyon in the urban context
Zeng et al. Adaptive particle swarm optimization guided by acceleration information
Elling Carbuncles as self-similar entropy solutions
Schruff et al. Numerical simulation of pore fluid flow and fine sediment infiltration into the riverbed

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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 201204 building 27, Lane 676, Wuxing Road, Pudong New Area, Shanghai

Patentee after: Shanghai suochen Information Technology Co., Ltd

Address before: No.27, Lane 676, Wuxing Road, Pudong New Area, Shanghai, 201204

Patentee before: SHANGHAI SUOCHEN INFORMATION TECHNOLOGY Co.,Ltd.