CN112084724B - A CFD Local Fast Calculation Method - Google Patents
A CFD Local Fast Calculation Method Download PDFInfo
- Publication number
- CN112084724B CN112084724B CN202010936081.1A CN202010936081A CN112084724B CN 112084724 B CN112084724 B CN 112084724B CN 202010936081 A CN202010936081 A CN 202010936081A CN 112084724 B CN112084724 B CN 112084724B
- Authority
- CN
- China
- Prior art keywords
- calculation
- cfd
- block
- local
- calculation example
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000013507 mapping Methods 0.000 claims abstract description 12
- 239000000463 material Substances 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims 1
- 239000013589 supplement Substances 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 10
- 239000012530 fluid Substances 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 description 12
- 230000000903 blocking effect Effects 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种CFD局部快速计算方法，涉及计算流体动力学中的局部快速计算方法。解决现有技术存在浪费计算资源，限制了优化迭代的速度等问题，提供一种CFD局部快速计算方法，包括映射：将完成计算的原始算例的CFD计算结果映射到未进行计算的新算例中；分块：对新算例进行分块处理，将新算例与原始算例几何结构存在差异的部分进行分离，获得分块算例；算例处理：将获得的分块算例进行处理，获得用于独立计算的分块算例；计算：调用CFD软件计算获得的分块算例，完成CFD局部快速计算。本发明将庞大复杂的完整算例拆分，只对需要的局部进行计算，大大降低计算资源，加快单次迭代计算速度，使得局部优化可以得到及时的反馈，大幅减少迭代周期。
A CFD local fast calculation method relates to a local fast calculation method in computational fluid dynamics. To solve the problems of wasting computing resources and limiting the speed of optimization iterations in the prior art, a local fast calculation method for CFD is provided, including mapping: mapping the CFD calculation result of the original calculation example that has completed the calculation to the new calculation example that has not been calculated. Medium; block: perform block processing on the new example, separate the parts with different geometric structures from the original example, and obtain block examples; example processing: process the obtained block examples , obtain the block calculation example for independent calculation; Calculation: call the block calculation example obtained by CFD software calculation to complete the local fast calculation of CFD. The invention splits the huge and complex complete calculation example, and only calculates the required part, which greatly reduces the calculation resources, speeds up the calculation speed of a single iteration, enables the local optimization to obtain timely feedback, and greatly reduces the iteration cycle.
Description
Technical Field
The invention relates to a local fast calculation method in computational fluid dynamics.
Background
Computational Fluid Dynamics (CFD for short) is based on hydrodynamics and computer science, and an approximate solution of a Fluid control equation is obtained by starting from a calculation method and utilizing the quick calculation capability of a computer, so that the Computational Fluid Dynamics is an important technical means for researching related theories of hydrodynamics and engineering application.
Because problems in actual engineering usually cannot be analyzed and the cost for making a model is high, a great amount of CFD is usually relied on to carry out numerical simulation in the product design process, and the accuracy and the speed of CFD calculation results directly determine the length of a design iteration period. Therefore, how to improve the calculation speed and accuracy of CFD has been a key area of research.
The product structure in actual engineering is often extremely complex, in order to ensure accuracy, a sufficiently fine physical model and a sufficiently dense grid are required to be used for numerical simulation, so that the calculation speed is greatly reduced, and if some calculation problems are met, the design iteration cycle is also prolonged; conversely, if the model is simplified, the computation speed may be increased, but the accuracy will be affected. Especially, when local optimization iteration is performed, even if the model is changed a little, the same long numerical simulation time is needed, and the time is wasted.
In addition, because the required computing resources are often very huge, in order to save resources, the optimization of different parts of the product needs to unify the iteration cycle and calculate uniformly, so that the local optimization works in an elbow manner, and the working efficiency is influenced.
Therefore, when the product is optimized in local detail, the traditional calculation method wastes calculation resources seriously, limits the speed of optimization iteration, and needs a new flexible and quick calculation method aiming at local parts.
Disclosure of Invention
The invention aims to solve the problems of wasting computing resources, limiting the speed of optimization iteration and the like in the prior art, and provides a CFD local rapid computing method, which mainly aims at the complex original model and the algorithm which finish CFD computing, and realizes the CFD local rapid computing of a new model and the new algorithm after the local geometric structure of the complex original model and the complex original algorithm is modified.
A CFD local rapid calculation method is realized by the following steps:
step one, mapping;
mapping the CFD calculation result of the original example which is calculated to a new example which is not calculated;
step two, partitioning;
carrying out blocking processing on the new example obtained in the step one, and separating the part of the new example, which has a difference with the geometric structure of the original example, to obtain a blocking example;
step three, example processing;
processing the blocking calculation example obtained in the step two to obtain a blocking calculation example for independent calculation;
step four, calculating;
and calling the block calculation example obtained in the CFD software calculation step three to finish the CFD local quick calculation.
The invention has the beneficial effects that: the invention splits a large and complex complete example, and only calculates the required local part, thereby greatly reducing the calculation resource, accelerating the single iterative calculation speed, leading the local optimization to obtain the timely feedback and greatly reducing the iterative period;
the invention can make the local optimization work more flexibly and dispersedly, and avoid the situation that a plurality of local optimization works are mutually locked when synchronously carried out;
the method of the invention applies a mapping interpolation method to provide a reasonable initial field for CFD calculation, can accelerate calculation convergence and reduce iteration times;
the invention can realize fast self-defined calculation region division, and is convenient for users to adjust the block calculation region so as to give consideration to calculation speed and accuracy.
Drawings
Fig. 1 is a schematic diagram of a CFD local fast computation method according to the present invention.
Detailed Description
In a first embodiment, the present embodiment is described with reference to fig. 1, and a method for local fast CFD calculation is implemented by the following steps:
1. mapping;
and (3) performing interpolation mapping on the CFD calculation result of the original calculation example which is subjected to calculation into a new calculation example which is not subjected to calculation by adopting an interpolation mapping tool, and performing reasonable assignment on the grid area which is not subjected to mapping by adopting a custom tool, thereby providing an accurate boundary condition and a reasonable initial field for the calculation of the subsequent block calculation example.
Since the new and old algorithms only have slight geometrical structure difference, the calculation results are almost the same in the areas far away from the geometrical difference, so the new algorithm can directly use the calculation results of the original algorithm in most areas, and the areas do not need further calculation. The new example thus has areas that need to be calculated and areas that do not need to be calculated.
2. Block processing;
the new algorithm is subjected to parallel block processing by using a parallel block tool, a local geometric structure of the new algorithm, which is different from the geometric structure of the original algorithm, can be accurately positioned by using a self-defined block method in the process, and the local geometric structure is separated from the complete structure, so that the corresponding local grid and the corresponding physical field are separated into the block algorithm which cannot be calculated yet.
In this embodiment, the grid and the physical field of the operator are divided into a plurality of parallel sub-operators before parallel computation is performed, and each computer thread only needs to compute a single parallel sub-operator, so that the parallel computation efficiency is higher than that of single-core computation. But the parallel sub-algorithm cannot independently calculate because the parallel sub-algorithm needs to transmit data through a newly-built boundary.
Based on the above principle, the new example is subjected to parallel blocking processing, and is divided into an area requiring calculation and an area not requiring calculation, and the former is called a blocking example. At the moment, the block calculation example is not different from the parallel sub-calculation example, the boundary condition is connected with other parallel sub-calculation examples, and the files such as physical properties are lacked, so that independent calculation cannot be carried out.
3. Example processing;
the grid and the physical field of the block calculation example which can not be calculated are transplanted into an independent calculation example file, the boundary conditions are modified, and necessary files such as the gravity acceleration, the material physical property, the turbulence model and the like are supplemented to enable the block calculation example to become the block calculation example which can be independently calculated.
4. Calculating;
and calling CFD software to calculate a block calculation example to obtain a calculation result and finish the local rapid calculation of the CFD.
In a second specific embodiment, the present embodiment is an application example of the CFD local fast calculation method according to the first specific embodiment:
the present embodiment takes CFD simulation calculation of a turbine guide vane as an example. After the original example A is calculated, the local geometry of the original example A is modified (one air film hole at the leaf back is deleted) to generate a new example B, and the CFD local calculation steps of the new example B are as follows:
1. the calculation result of the example A is subjected to interpolation mapping to the example B, the process calls a tool mapFields of OpenFOAM software, and a user-defined tool is adopted to reasonably assign a physical field corresponding to a grid area which cannot be mapped by the example B, so that an accurate boundary condition and a reasonable initial field are provided for the calculation of the subsequent block-based example;
2. the method comprises the following steps of performing parallel block processing on an operator B, calling a tool decomposePar of OpenFOAM software in the process, and realizing automatic self-defined block division based on the tool, so that leaf back regions modified by the operator A and the operator B can be accurately positioned, and the local geometric structure is separated from a complete structure, so that a corresponding local grid and a corresponding physical field are separated into a block operator C which cannot be calculated;
3. transplanting a blocking example C, modifying the boundary conditions of the blocking example C, supplementing necessary files such as material thermophysical properties, turbulence model parameters and the like, wherein the example C can carry out independent calculation and the process is mainly realized by an automatic script;
4. and (4) calculating an example C, completing the local rapid calculation of the CFD, and calling OpenFOAM software in the process.
In the actual operation of the calculation method described in this embodiment, only the example positions of the example a and the example B and the parameters required for the custom partitioning need to be filled in the file, and only the written script needs to be run to realize the above-described whole flow, so that the CFD local fast calculation is realized.
In this embodiment, the number of grids in the example B is about 1315 ten thousand, if the calculation occupies about 10G of memory according to the conventional calculation method, the single-step iteration time is about 62.5 s; however, by adopting the method of the embodiment, only the calculation example C is needed, the number of grids is about 145 ten thousand, the memory occupied by calculation is about 1.2G, the single-step iteration time is about 4.6s, and the method of the embodiment has obvious advantages compared with the traditional method.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (5)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010936081.1A CN112084724B (en) | 2020-09-08 | 2020-09-08 | A CFD Local Fast Calculation Method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010936081.1A CN112084724B (en) | 2020-09-08 | 2020-09-08 | A CFD Local Fast Calculation Method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112084724A CN112084724A (en) | 2020-12-15 |
| CN112084724B true CN112084724B (en) | 2021-06-22 |
Family
ID=73731572
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010936081.1A Active CN112084724B (en) | 2020-09-08 | 2020-09-08 | A CFD Local Fast Calculation Method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112084724B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102750414A (en) * | 2012-06-25 | 2012-10-24 | 江苏大学 | Method for testing relationship between mesh quality of centrifugal pump and computational accuracy |
| CN104573152A (en) * | 2013-10-22 | 2015-04-29 | 北京临近空间飞行器系统工程研究所 | Method for shortening CFD (computation fluid dynamics) numerical-simulation calculating time on basis of pre-operation |
| CN109446660A (en) * | 2018-10-31 | 2019-03-08 | 潍柴动力股份有限公司 | A kind of engine inlets automatic optimization method and device |
| US10684078B1 (en) * | 2019-05-24 | 2020-06-16 | Macrae Technologies, Inc. | Method for stabilizing thermal conduction of block coolers with cast-in coolant pipes |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107992664B (en) * | 2017-11-27 | 2021-04-13 | 重庆大学 | A Calculation Method of Commutation Transformer Electricity-Thermal Coupled Field Based on Improved Meshless Method |
| CN110851987B (en) * | 2019-11-14 | 2022-09-09 | 上汽通用五菱汽车股份有限公司 | Method, apparatus and storage medium for predicting calculated duration based on acceleration ratio |
-
2020
- 2020-09-08 CN CN202010936081.1A patent/CN112084724B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102750414A (en) * | 2012-06-25 | 2012-10-24 | 江苏大学 | Method for testing relationship between mesh quality of centrifugal pump and computational accuracy |
| CN104573152A (en) * | 2013-10-22 | 2015-04-29 | 北京临近空间飞行器系统工程研究所 | Method for shortening CFD (computation fluid dynamics) numerical-simulation calculating time on basis of pre-operation |
| CN109446660A (en) * | 2018-10-31 | 2019-03-08 | 潍柴动力股份有限公司 | A kind of engine inlets automatic optimization method and device |
| US10684078B1 (en) * | 2019-05-24 | 2020-06-16 | Macrae Technologies, Inc. | Method for stabilizing thermal conduction of block coolers with cast-in coolant pipes |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112084724A (en) | 2020-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105760624B (en) | It is a kind of support extensive three dimensional integrated circuits heat emulation and Thermal design | |
| CN103869755B (en) | A kind of machining code feed speed that adjusts is to ensure the method that lathe output is smooth | |
| CN103279627B (en) | A kind of heat-machinery based on finite element-wearing and tearing coupling analysis method for numerical simulation | |
| CN104200045B (en) | A Parallel Calculation Method for Distributed Hydrodynamic Model of Large-Scale Water System in Watershed | |
| CN110826188B (en) | GPU acceleration-based natural gas pipe network hydraulic parameter simulation method | |
| CN108009349B (en) | A Computational Grid Optimization Drawing Method for Two-dimensional Water Quality Numerical Simulation Model of Rivers | |
| CN116258042B (en) | Large-scale heat transfer heterogeneous parallel simulation method based on DDM | |
| CN107045571A (en) | A kind of die casting simulation control method based on cloud computing | |
| US11602908B1 (en) | Method of mesh generation for resin transfer molding process | |
| CN113392472B (en) | OpenMP parallel disturbance domain updating method for aircraft aerodynamic characteristic simulation | |
| CN114004115A (en) | Novel rapid three-dimensional pipe network half-gas thermal coupling evaluation method for complex cooling blades | |
| CN112084724B (en) | A CFD Local Fast Calculation Method | |
| CN102063539A (en) | Inertial platform residual stress release simulating method based on finite element | |
| CN112861374A (en) | Multi-physical coupling simulation processing method, device and equipment based on pre-controller | |
| CN109033733A (en) | A kind of finite element Region Decomposition improves the parallel method in SSORPCG solution temperature field | |
| CN114398728B (en) | Mold temperature simulation analysis method considering cooling water temperature change | |
| CN110442935B (en) | Production beat design and optimization method for automatic equipment | |
| CN108733869B (en) | A large-scale three-dimensional integrated circuit partition method and device | |
| US20120316842A1 (en) | Automatic extremum detection on a surface mesh of a component | |
| CN102930589A (en) | Method for correcting intersecting line of unstructured Cartesian grid | |
| CN111666644A (en) | Performance analysis method for parameterized CFD aerostatic bearing | |
| CN111090917A (en) | Automatic analysis method for whole vehicle air intake and exhaust system based on computational fluid dynamics | |
| CN114169100B (en) | Efficient design optimization method and system for super-large variable impeller machinery and application | |
| CN116644561A (en) | Multi-scale numerical simulation method for whole-stage cascade of gas turbine | |
| CN112182771B (en) | Data processing method based on numerical simulation, storage medium and electronic device |
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 |