CN112364362B - Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction - Google Patents
Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction Download PDFInfo
- Publication number
- CN112364362B CN112364362B CN202011279002.0A CN202011279002A CN112364362B CN 112364362 B CN112364362 B CN 112364362B CN 202011279002 A CN202011279002 A CN 202011279002A CN 112364362 B CN112364362 B CN 112364362B
- Authority
- CN
- China
- Prior art keywords
- layer
- block structure
- grid
- layer block
- field variable
- 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 57
- 239000012530 fluid Substances 0.000 title claims abstract description 13
- 238000004088 simulation Methods 0.000 title claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 230000008859 change Effects 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 19
- 230000004069 differentiation Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
- G06F21/60—Protecting data
- G06F21/602—Providing cryptographic facilities or services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Bioethics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Software Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a parallel multi-layer self-adaptive local encryption method oriented to a fluid simulation direction, which comprises the following steps: step S1, constructing an initial multi-layer block structure layout according to the actual physical scale of the space domain and initial conditions, wherein the initial multi-layer block structure layout comprises: base layer block structures and other layer block structures; step S2, carrying out grid division layer by layer on the initial multi-layer block structure layout in the step S1; step S3, applying physical boundary conditions to the space domain after the step S2 layer-by-layer grid division; s4, reconstructing a multi-layer block structure layout based on the field variable gradient at intervals of a designated time step, so that a block structure of a finer layer can timely cover a region with severe field variable change; and S5, outputting a picture result of the moment at intervals of a designated time. The beneficial effects of the invention are as follows: the problem that the calculation accuracy is not reduced while the calculation grid quantity scale is reduced is solved by adopting a multi-layer block structure local encryption technology.
Description
Technical Field
The invention relates to the field of CAE encryption, in particular to a parallel multi-layer self-adaptive local encryption method oriented to a fluid simulation direction.
Background
The industrial field CAE generally uses a uniform grid or a variable-size grid with partial encryption to discrete the computational domain and performs the analog computation of the flow field variables based on the discrete grid. The uniform grid is simple to process in the calculation process, but huge calculation amount can be caused for the calculation problems of complex calculation domain and huge calculation under the condition of keeping the same calculation precision. The local encryption variable-size grid performs grid local encryption on a computation domain with complex and variable topological structure, and coarse grid approximation processing is adopted on the computation domain with simple structure; the local encryption variable-size grid can reduce the size of the calculation grid and improve the calculation efficiency under the condition of keeping certain calculation precision, but the complexity can be improved when the grid is split, and meanwhile, a parallel calculation architecture cannot be adopted efficiently in the flow calculation process, so that the improvement of the calculation efficiency is limited greatly.
Disclosure of Invention
The invention overcomes the defects in the prior art and provides a parallel multi-layer self-adaptive local encryption method facing the fluid simulation direction.
The aim of the invention is achieved by the following technical scheme.
A parallel multi-layer self-adaptive local encryption method facing a fluid simulation direction comprises the following steps:
step S1, constructing an initial multi-layer block structure layout according to the actual physical scale of the space domain and initial conditions, wherein the initial multi-layer block structure layout comprises the following steps: base layer block structures and other layer block structures;
step S2, carrying out grid division layer by layer on the initial multi-layer block structure layout in the step S1;
step S3, physical boundary conditions are applied to the space domain after the step S2 layer-by-layer grid division;
s4, reconstructing a multi-layer block structure layout based on the field variable gradient at intervals of a designated time step, so that a block structure of a finer layer can timely cover a region with severe field variable change;
and S5, outputting a picture result of the moment at intervals of a designated time step.
Further, in the step S1, the base layer block structure and the region with severe field variable variation are marked according to different field variable variations, respectively, to obtain a multi-layer block structure layout.
Further, in the step S1, the base layer block structure covers a complete calculation domain when applied to flow field calculation.
Further, in the step S2, the layer-by-layer meshing includes: and adopting a layer-by-layer division mode aiming at each layer of block structure, and carrying out grid division on the current layer of block structure by using a fixed grid size aiming at each layer of block structure from a base layer to the layer-by-layer sequence of other layer of block structures.
Further, in the step S3, the applying physical boundary conditions includes:
the first is the physical boundary condition of the described physical problem;
the second type is boundary conditions corresponding to fine meshes at the junction of the block structure of the thinner layer and the coarse meshes intersected below.
Further, applying a first type of boundary condition comprises the steps of:
firstly, obtaining boundary condition types on different grids of a basic layer block structure according to the application field;
the respective type of boundary condition is then applied to the corresponding grid of the base layer block structure.
Further, applying a second type of boundary condition, comprising the steps of:
for the boundary condition corresponding to the boundary grid of the block structure of the thinner layer, interpolation acquisition of the field variable values of the coarse grid intersected below the block structure of the thinner layer is adopted.
Further, the specific calculation equation for solving the field variable in the step S3 is as follows:
wherein ρ represents density, c represents specific heat capacity, T represents temperature, λ thermal conductivity is time T, x, y, z are space coordinates, and Q is an internal heat source;is a partial differential operator.
The above equation is shown in the space discrete process as follows, taking one dimension as an example:
in the method, in the process of the invention,is the temperature at position i at time n,/-, for example>Is the temperature at the i-1 position at time n, +.>Is the temperature of the i+1 position at time n, delta x is the space step size, o [ (delta x) 2 ]Is the second order of the spatial step size infinitesimal.
The time differentiation adopts Euler method, namely forward explicit differential format, as follows:
in the above-mentioned method, the step of,for the temperature value of the i node position at the time n+1, Δt is the time step, omicron (Δt) is Δt first order, less, and so on.
The discretized temperature field variable equation is as follows:
the variables in the above formula have the meanings described above.
Further, the specific steps for solving the field variable in the step S3 are as follows:
step S31, solving the discrete field variable of the ith node at the (n+1) th moment for the basic layer block structure
Step S32, after the solving of the basic layer block structure is finished, assigning the grid intersected by the first layer and the basic layer block structure as a grid field variable value corresponding to the basic layer block structure;
step S33, solving the field variables from the first layer to the topmost layer by analogy, and meanwhile, directly giving adjacent coarse layers after interpolating the field variables obtained by calculating the fine layers.
The beneficial effects of the invention are as follows:
different from the traditional local encryption method, the method adopts parallel multi-layer, block structure and self-adaptive local encryption method to reconstruct the calculation process of the flow field variable in the industrial CAE field, and improves the calculation efficiency by nearly two orders of magnitude while keeping the same calculation precision;
the method aims to solve the problem that the calculation accuracy is not reduced while the calculation grid quantity scale is reduced by adopting a multi-layer block structure local encryption technology, and simultaneously adopts an MPI parallel mechanism in the calculation process, so that the calculation efficiency is improved by two orders of magnitude compared with the traditional method.
Drawings
FIG. 1 is a schematic diagram of the encryption method of the present invention;
FIG. 2 is a schematic diagram of an initial multi-layer block structure layout;
FIG. 3 is a schematic diagram of an internal boundary interpolation process;
FIG. 4 is a diagram of a field variable control equation solving process under a multi-layer block structure architecture;
FIG. 5 is a schematic diagram of the dependency of a multi-layer block structure on MPI calculation threads;
FIG. 6 is a schematic diagram of a dynamic reconfiguration of a multi-layer block structure layout.
Detailed Description
The technical scheme of the invention is further described by specific examples.
The parallel multi-layer self-adaptive local encryption method facing to the fluid simulation direction as shown in fig. 1 comprises the following steps:
step S1, constructing an initial multi-layer block structure layout according to the actual physical scale of the space domain and initial conditions, wherein the initial multi-layer block structure layout comprises the following steps: base layer block structures and other layer block structures;
in the step S1, marking the basic layer block structure and the severe field variable change area according to different field variable changes to obtain a multi-layer block structure layout;
in the step S1, when applied to flow field calculation, the base layer block structure covers a complete calculation domain;
step S2, carrying out grid division layer by layer on the initial multi-layer block structure layout in the step S1; in the step S2, the layer-by-layer meshing includes: adopting a layer-by-layer division mode aiming at each layer of block structure, and carrying out grid division on the current layer of block structure according to a fixed grid size aiming at each layer of block structure in a layer-by-layer sequence from a base layer to other layers of block structures;
step S3, physical boundary conditions are applied to the space domain after the step S2 layer-by-layer grid division; in the step S3, the applying physical boundary conditions includes:
the first is the physical boundary condition of the described physical problem;
the second type is boundary conditions corresponding to fine meshes at the junction of the block structure of the thinner layer and the coarse meshes intersected below.
Further, applying a first type of boundary condition comprises the steps of:
firstly, obtaining boundary condition types on different grids of a basic layer block structure according to the application field;
the respective type of boundary condition is then applied to the corresponding grid of the base layer block structure.
Further, applying a second type of boundary condition, comprising the steps of:
for the boundary condition corresponding to the boundary grid of the block structure of the thinner layer, interpolation acquisition of the field variable values of the coarse grid intersected below the block structure of the thinner layer is adopted.
Further, the specific calculation equation for solving the field variable in the step S3 is as follows:
wherein ρ represents density, c represents specific heat capacity, T represents temperature, λ thermal conductivity is time T, x, y, z are space coordinates, and Q is an internal heat source;is a partial differential operator.
The above equation is shown in the space discrete process as follows, taking one dimension as an example:
in the method, in the process of the invention,is the temperature at position i at time n,/-, for example>Is the temperature at the i-1 position at time n, +.>Is the temperature of the i+1 position at time n, delta x is the space step size, o [ (delta x) 2 ]Is the second order of the spatial step size infinitesimal.
The time differentiation adopts Euler method, namely forward explicit differential format, as follows:
in the above-mentioned method, the step of,for the temperature value of the i node position at the time n+1, Δt is the time step, omicron (Δt) is Δt first order, less, and so on.
The discretized temperature field variable equation is as follows:
the variables in the above formula have the meanings described above.
S4, reconstructing a multi-layer block structure layout based on the field variable gradient at intervals of a designated time step, so that a block structure of a finer layer can timely cover a region with severe field variable change;
the specific steps for solving the field variable in the step S3 are as follows:
step S31, solving the discrete field variable of the ith node at the (n+1) th moment for the basic layer block structure
Step S32, after the solving of the basic layer block structure is finished, assigning the grid intersected by the first layer and the basic layer block structure as a grid field variable value corresponding to the basic layer block structure;
step S33, solving the field variables from the first layer to the topmost layer by analogy, and meanwhile, directly giving adjacent coarse layers after interpolating the field variables obtained by calculating the fine layers.
And S5, outputting a picture result of the moment at intervals of a designated time step.
Example 1
The flow field is taken as an example to illustrate the implementation steps of the method.
Step S1, constructing an initial multi-layer block structure layout
The initial multi-layer and block structure layout is determined according to the initial conditions, and as shown in the round tube round flow calculation example shown in fig. 2, the flow state change around the round tube is generally severe, so that marking is required before calculation according to the outer contour of the round tube, and the multi-layer block structure layout is constructed. The specific construction results are shown in fig. 1, wherein the basic layer block structure covers the complete calculation domain, the first layer block structure covers a certain range around the circular tube, and the second layer block structure tightly surrounds the outer contour of the circular tube.
Step S2, carrying out grid division layer by layer on the initial multi-layer block structure layout in the step S1;
starting from the base layer, each layer performs meshing on the current intra-layer block structure with a fixed mesh size.
Step S3, physical boundary conditions are applied to the space domain after the step S2 layer-by-layer grid division;
the multi-layer block structure layout has two types of boundary conditions, one is the physical boundary condition of the described physical problem, and the other is the boundary condition corresponding to the fine grid at the junction of the finer layer and the lower teaching coarse layer grid.
For the first type of boundary conditions, the boundary conditions need to be added to the basic layer boundary grid, the inlet is a speed boundary condition, the outlet is pressure, the speed boundary condition is applied to the inlet of the corresponding basic layer boundary grid, and the pressure boundary condition is applied to the outlet.
For the second type of boundary conditions, namely boundary conditions corresponding to boundary grids corresponding to a thinner layer of block structures, interpolation of field variable values of coarse grids intersected below the boundary conditions is adopted, and a specific interpolation process is shown in fig. 3.
The specific calculation equation for solving the field variables is as follows:
wherein ρ represents density, c represents specific heat capacity, T represents temperature, λ thermal conductivity is time T, x, y, z is space coordinates, and Q is internal heatA source;is a partial differential operator.
The above equation is shown in the space discrete process as follows, taking one dimension as an example:
in the method, in the process of the invention,is the temperature at position i at time n,/-, for example>Is the temperature at the i-1 position at time n, +.>Is the temperature of the i+1 position at time n, delta x is the space step size, o [ (delta x) 2 ]Is the second order of the spatial step size infinitesimal.
The time differentiation adopts Euler method, namely forward explicit differential format, as follows:
in the above-mentioned method, the step of,for the temperature value of the i node position at the time n+1, Δt is the time step, omicron (Δt) is Δt first order, less, and so on.
The discretized temperature field variable equation is as follows:
the variables in the above formula have the meanings described above.
The specific steps of solving are as follows, and the specific flow is shown in fig. 4:
step S31, the foundation layer solves discrete field variables of the ith node at the n+1th moment;
step S32, after the solution of the base layer is finished, assigning the grid intersected by the first layer and the base layer as a grid field variable value corresponding to the base layer;
step S33, solving the field variables from the first layer to the topmost layer by analogy, and meanwhile, calculating the fine layer to obtain the field variable interpolation and then directly giving the field variable interpolation to the adjacent coarse layer;
the whole solving process is based on a parallel method of a message passing Mechanism (MPI), each layer of each block structure is distributed to a specific calculation thread according to a certain algorithm, and data communication between each calculation block is completed by the MPI.
The dependency of the multi-layer block structure on the MPI calculation thread is shown in FIG. 5, where the cuboid representations of the same color are owned by the same calculation thread, and vice versa.
S4, reconstructing a multi-layer block structure layout based on the field variable gradient at intervals of a designated time step, so that a block structure of a finer layer can timely cover a region with severe field variable change;
as shown in fig. 6, the field variable is velocity in the figure, i.e. the multi-layer block structure layout performs marker reconstruction from the velocity gradient. The time interval for layout reconstruction in the figure is 4000 time steps.
In the process of solving the field variable control equation, in order to encrypt the local grid of the field variable region with severe change in time, the multi-layer block structure layout needs to be reconfigured every a plurality of time periods, so that the block structure with finer layers can timely cover the field variable region with severe change. The whole block structure layout reconstruction process is self-adaptive, namely marking is carried out according to the gradient value of the field variable, marking is carried out when the gradient value is changed greatly, and the marked area is covered by a finer layer block structure during layout reconstruction.
And S5, outputting a picture result of the moment at intervals of a designated time step.
The foregoing describes one embodiment of the present invention in detail, but the description is only a preferred embodiment of the present invention and should not be construed as limiting the scope of the invention. All equivalent changes and modifications within the scope of the present invention are intended to be covered by the present invention.
Claims (6)
1. A parallel multi-layer self-adaptive local encryption method facing to a fluid simulation direction is characterized in that: the method comprises the following steps:
step S1, constructing an initial multi-layer block structure layout according to the actual physical scale of the space domain and initial conditions, wherein the initial multi-layer block structure layout comprises the following steps: base layer block structures and other layer block structures;
step S2, carrying out grid division layer by layer on the initial multi-layer block structure layout in the step S1;
step S3, physical boundary conditions are applied to the space domain after the step S2 layer-by-layer grid division;
in the step S3, the applying physical boundary conditions includes:
the first is the physical boundary condition of the described physical problem;
the second type is a boundary condition corresponding to a fine grid at the junction of a block structure of a thinner layer and a coarse grid intersected below;
applying a second type of boundary condition comprising the steps of:
for the boundary condition corresponding to the boundary grid of the thinner one-layer block structure, adopting interpolation acquisition of field variable values of the intersected coarse grid below the thinner one-layer block structure;
s4, reconstructing a multi-layer block structure layout by taking the field variable gradient as a criterion at intervals of a designated time step, so that a block structure of a finer layer can timely cover a region with severe field variable change;
step S5, outputting a picture result at a moment corresponding to the current time step at intervals of a designated time step;
the specific calculation equation for solving the field variable in the step S3 is as follows:
wherein ρ represents density, c represents specific heat capacity, T represents temperature, λ thermal conductivity is time T, x, y, z are space coordinates, and Q is an internal heat source;is a partial differential operator;
the above equation is shown in the space discrete process as follows, taking one dimension as an example:
in the method, in the process of the invention,is the temperature at position i at time n,/-, for example>Is the temperature at the i-1 position at time n, +.>The temperature of the i+1 position at time n, Δx is the space step size, o [ (Δx) 2 ]Is the second order infinitesimal of the space step;
the time differentiation adopts Euler method, namely forward explicit differential format, as follows:
in the above-mentioned method, the step of,for the temperature value of the i node position at the time n+1, deltat is the time step, o (deltat) is deltat first-order anecdotal, and the other analogy is that;
the discretized temperature field variable equation is as follows:
wherein the variables are as defined above.
2. The parallel multi-layer adaptive partial encryption method oriented to the fluid simulation direction according to claim 1, wherein the method comprises the following steps: in the step S1, the base layer block structure and the field variable severe region are marked according to different field variable changes, so as to obtain a multi-layer block structure layout.
3. The parallel multi-layer adaptive partial encryption method facing the fluid simulation direction according to claim 1 or 2, wherein: in the step S1, the base layer block structure covers the complete calculation domain when applied to flow field calculation.
4. The parallel multi-layer adaptive partial encryption method oriented to the fluid simulation direction according to claim 1, wherein the method comprises the following steps: in the step S2, the layer-by-layer meshing includes: and adopting a layer-by-layer division mode aiming at each layer of block structure, and carrying out grid division on the current layer of block structure by using a fixed grid size aiming at each layer of block structure from a base layer to the layer-by-layer sequence of other layer of block structures.
5. The parallel multi-layer adaptive partial encryption method oriented to the fluid simulation direction according to claim 1, wherein the method comprises the following steps: applying a first type of boundary condition comprising the steps of:
firstly, obtaining boundary condition types on different grids of a basic layer block structure according to the application field;
the respective type of boundary condition is then applied to the corresponding grid of the base layer block structure.
6. The parallel multi-layer adaptive partial encryption method oriented to the fluid simulation direction according to claim 1, wherein the method comprises the following steps: the specific steps for solving the field variable in the step S3 are as follows:
step S31, solving the discrete field variable of the ith node at the (n+1) th moment for the basic layer block structure
Step S32, after the solving of the basic layer block structure is finished, assigning the grid intersected by the first layer and the basic layer block structure as a grid field variable value corresponding to the basic layer block structure;
step S33, solving the field variables from the first layer to the topmost layer by analogy, and meanwhile, directly giving adjacent coarse layers after interpolating the field variables obtained by calculating the fine layers.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011279002.0A CN112364362B (en) | 2020-11-16 | 2020-11-16 | Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011279002.0A CN112364362B (en) | 2020-11-16 | 2020-11-16 | Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112364362A CN112364362A (en) | 2021-02-12 |
| CN112364362B true CN112364362B (en) | 2023-12-29 |
Family
ID=74514948
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011279002.0A Active CN112364362B (en) | 2020-11-16 | 2020-11-16 | Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112364362B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114741983B (en) * | 2022-03-01 | 2023-10-13 | 华南理工大学 | Multi-layer fluid sloshing simulation method based on OpenFOAM |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101833605A (en) * | 2010-04-29 | 2010-09-15 | 浙江工业大学 | Control method for precision processing of micro-flow path abrasive flow based on flow volume mould |
| EP2287757A1 (en) * | 2009-08-20 | 2011-02-23 | Fujitsu Limited | Multilevel-Multigrid simulation techniques |
| CN103970989A (en) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | Immersing boundary flow field calculation method based on fluid/solid interface consistency |
| CN106227951A (en) * | 2016-07-26 | 2016-12-14 | 南京航空航天大学 | It is applicable to the two-dimentional wind pulse computational methods of the super-huge cooling tower of supergage |
| CN109063248A (en) * | 2018-06-26 | 2018-12-21 | 广东工业大学 | It is a kind of for laser impact intensified gas-liquid-solid coupling calculation |
| CN110309543A (en) * | 2019-05-31 | 2019-10-08 | 中国航天空气动力技术研究院 | A kind of simulation flow design method of multi-component fluid interface movement |
| CN110489777A (en) * | 2019-07-02 | 2019-11-22 | 哈尔滨理工大学 | A kind of array hole electric spark ultrasonic Compound Machining flow field simulation method |
| CN111324947A (en) * | 2020-01-21 | 2020-06-23 | 中国石油化工股份有限公司 | Oil reservoir numerical simulation non-uniform angular point grid optimization method based on attribute distribution |
| CN111444614A (en) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | Flow field reconstruction method based on graph convolution |
| CN111581858A (en) * | 2020-03-24 | 2020-08-25 | 国网甘肃省电力公司经济技术研究院 | Design method of heat preservation and insulation system for power transmission and transformation foundation |
| CN111859766A (en) * | 2020-07-28 | 2020-10-30 | 深圳拳石科技发展有限公司 | Lagrange integral point finite element numerical simulation system and method of variable calculation domain |
| CN111859645A (en) * | 2020-07-09 | 2020-10-30 | 南京理工大学 | Improved MUSL format material dot method for shock wave solving |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1617309B1 (en) * | 2004-07-15 | 2011-01-12 | Fujitsu Limited | Simulation technique with local grid refinement |
| CN106295120B (en) * | 2016-07-26 | 2019-10-18 | 上海无线电设备研究所 | A kind of electromagnetic characteristic of scattering quick calculation method |
| EP3740886A1 (en) * | 2018-03-16 | 2020-11-25 | Siemens Aktiengesellschaft | Topology optimization with design-dependent loads and boundary conditions for multi-physics applications |
-
2020
- 2020-11-16 CN CN202011279002.0A patent/CN112364362B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2287757A1 (en) * | 2009-08-20 | 2011-02-23 | Fujitsu Limited | Multilevel-Multigrid simulation techniques |
| CN101833605A (en) * | 2010-04-29 | 2010-09-15 | 浙江工业大学 | Control method for precision processing of micro-flow path abrasive flow based on flow volume mould |
| CN103970989A (en) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | Immersing boundary flow field calculation method based on fluid/solid interface consistency |
| CN106227951A (en) * | 2016-07-26 | 2016-12-14 | 南京航空航天大学 | It is applicable to the two-dimentional wind pulse computational methods of the super-huge cooling tower of supergage |
| CN109063248A (en) * | 2018-06-26 | 2018-12-21 | 广东工业大学 | It is a kind of for laser impact intensified gas-liquid-solid coupling calculation |
| CN110309543A (en) * | 2019-05-31 | 2019-10-08 | 中国航天空气动力技术研究院 | A kind of simulation flow design method of multi-component fluid interface movement |
| CN110489777A (en) * | 2019-07-02 | 2019-11-22 | 哈尔滨理工大学 | A kind of array hole electric spark ultrasonic Compound Machining flow field simulation method |
| CN111324947A (en) * | 2020-01-21 | 2020-06-23 | 中国石油化工股份有限公司 | Oil reservoir numerical simulation non-uniform angular point grid optimization method based on attribute distribution |
| CN111581858A (en) * | 2020-03-24 | 2020-08-25 | 国网甘肃省电力公司经济技术研究院 | Design method of heat preservation and insulation system for power transmission and transformation foundation |
| CN111444614A (en) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | Flow field reconstruction method based on graph convolution |
| CN111859645A (en) * | 2020-07-09 | 2020-10-30 | 南京理工大学 | Improved MUSL format material dot method for shock wave solving |
| CN111859766A (en) * | 2020-07-28 | 2020-10-30 | 深圳拳石科技发展有限公司 | Lagrange integral point finite element numerical simulation system and method of variable calculation domain |
Non-Patent Citations (4)
| Title |
|---|
| 光滑粒子流体动力学算法研究及其在流体力学中的数值模拟应用;王永鹏;《光滑粒子流体动力学算法研究及其在流体力学中的数值模拟应用》;全文 * |
| 唐军.《大重型数控静压转台设计及承载能力研究》.原子能出版社,2019,第62-66页. * |
| 立体旋流筛板空间气相流场的CFD模拟;王娜;张少峰;王德武;唐猛;张伟;刘燕;;化学工程(第02期);全文 * |
| 自适应结构网格上扩散方程隐式时间积分算法及其应用;徐小文;莫则尧;刘青凯;安恒斌;;计算物理(第05期);第684-691页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112364362A (en) | 2021-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111709171B (en) | Isogeometric solving and heat dissipation topology generation method for heat flow strong coupling problem | |
| CN111737835B (en) | Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method | |
| Zhang et al. | A hybrid boundary node method | |
| CN102819647B (en) | A kind of heterogeneous material random microscopic structure finite element modeling method | |
| CN111832204B (en) | Non-iterative topology optimization method of heat flow coupling structure | |
| WO2022000132A1 (en) | Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface | |
| US9183328B2 (en) | Method and apparatus for modeling interactions of the fluid with system boundaries in fluid dynamic systems | |
| CN112364362B (en) | Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction | |
| He et al. | An unsupervised learning approach to solving heat equations on chip based on auto encoder and image gradient | |
| CN115048901B (en) | Quantum layout optimization method and device and computer readable storage medium | |
| CN115758841A (en) | Finite element reduced modeling method for digital twinning application-oriented transformer temperature field | |
| CN114969860A (en) | Automatic hexahedron non-structural grid generation method | |
| CN111339688A (en) | Method for solving rocket simulation model time domain equation based on big data parallel algorithm | |
| Sitharam et al. | Corner-sharing tetrahedra for modeling micro-structure | |
| CN114077801B (en) | Particle data feature visualization method based on NURBS | |
| CN106066912A (en) | A kind of generation method of multi partition structured grid | |
| CN111239210A (en) | Capacitance tomography complex flow type data set establishing method | |
| CN114707291B (en) | Statics calculation method based on novel grid-free stable distribution points | |
| CN116579151B (en) | Non-uniform lattice structure optimization design method based on MMC framework | |
| CN115719045A (en) | Structured multi-layer multi-block nested background Cartesian grid rapid automatic generation method | |
| CN112632825B (en) | Electrostatic field smooth finite element numerical algorithm based on finite element super-convergence | |
| CN114004115A (en) | Novel rapid three-dimensional pipe network half-gas thermal coupling evaluation method for complex cooling blades | |
| Carré et al. | Multigrid strategies for CFD problems on non-structured meshes | |
| Liu et al. | Feature-aligned poly-square mapping of large-scale 2D geometries for semi-structured quad mesh generation | |
| Escobar et al. | The meccano method for isogeometric solid modeling |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Zhonglin Inventor after: Zhang Wei Inventor after: Wang Zhongwei Inventor before: Li Zhonglin Inventor before: Zhang Wei Inventor before: Guo Zhipeng Inventor before: Wang Zhongwei |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |