CN110110446A - A kind of parallel embedded grids methodology and device of achievable overlay region rapid Optimum - Google Patents

A kind of parallel embedded grids methodology and device of achievable overlay region rapid Optimum Download PDF

Info

Publication number
CN110110446A
CN110110446A CN201910384986.XA CN201910384986A CN110110446A CN 110110446 A CN110110446 A CN 110110446A CN 201910384986 A CN201910384986 A CN 201910384986A CN 110110446 A CN110110446 A CN 110110446A
Authority
CN
China
Prior art keywords
component
interpolation
grid
point
boundary
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.)
Granted
Application number
CN201910384986.XA
Other languages
Chinese (zh)
Other versions
CN110110446B (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.)
Xian Aeronautics Computing Technique Research Institute of AVIC
Original Assignee
Xian Aeronautics Computing Technique Research Institute of AVIC
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 Xian Aeronautics Computing Technique Research Institute of AVIC filed Critical Xian Aeronautics Computing Technique Research Institute of AVIC
Priority to CN201910384986.XA priority Critical patent/CN110110446B/en
Publication of CN110110446A publication Critical patent/CN110110446A/en
Application granted granted Critical
Publication of CN110110446B publication Critical patent/CN110110446B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Image Analysis (AREA)

Abstract

The invention discloses the parallel embedded grids methodologies and device of a kind of achievable overlay region rapid Optimum.This method comprises: the initial interpolation boundary of the nested grid of one group parts of label, a group parts include at least one moving component and at least one stationary parts;The information of flow of the nested grid characterization all parts;The edge that initial interpolation boundary is pushed away to all parts is obtained positioned at the intermediate region interpolation boundary of the adjacent component, so that the overlapping region of the nested grid reaches minimum;Determine that interpolation relationship of the interpolation boundary in the grid of the adjacent component, the interpolation relationship indicate relative position of the interpolation boundary point of the A component in the adjacent component in the adjacent component in the grid of B component;According to interpolation relationship, the information of flow in the grid of the B component is passed to the interpolation boundary point of the A component.This method solve model assemblies there are problems that the Numerical Investigation of Unsteady Flow Field of relative motion.

Description

A kind of parallel embedded grids methodology and device of achievable overlay region rapid Optimum
Technical field
The invention belongs to the numerical simulation field Fluid Mechanics Computation (Computational Fluid Dynamics, CFD), It is related to the parallel embedded grids methodology and device of a kind of achievable overlay region rapid Optimum.
Background technique
The aircraft of relative motion is in flight course between component, can generate that turn to twist separation, Shock wave interaction etc. complicated non- Steady Flow, these Complex Flows gently then influence flight quality, heavy then influence flight safety.For such as capsule ejection, guided missile The Unsteady Flows problem such as dispensing, propeller rotational, embedded grids methodology are a kind of effective solution approach.
The thinking of nested grid is: model sub-unit independently generates grid, and grid can do rigid motion, different portions with component Part grid is overlapped, the couple solution in the form of interpolation.Nested grid independence is strong, and quality is high, and formation speed is fast, but from Moving, efficiently establishing nested grid interpolation relationship is the key that embedded grids methodology.
In domestic and international existing embedded grids methodology, it is primarily present two deficiencies.1. some embedded grids methodologies do not make weight Folded area's optimization, interpolation boundary are located near object plane, and it is not high that biggish size of mesh opening difference will lead to interpolation precision;Even if other are embedding Optimize in nested grid method comprising overlay region, but needs to find donor to all mesh points in overlay region in optimization process Unit, and complicated donor element finding algorithm causes overlay region optimization efficiency lower.2. some embedded grids methodologies are monokaryons Serial operation, when grid is larger, the CFD numerical simulation used time is difficult to meet demand;Even if other embedded grids methodologies are more Core is run parallel, but two process groups used with CFD solver, insufficient for the inheritance of solver data with existing, right When the solver transformation of non-nesting function, process grouping can be related to and more complex data are transmitted, change amount is larger.
Summary of the invention
Goal of the invention: the invention discloses the parallel embedded grids methodologies and dress of a kind of achievable overlay region rapid Optimum It sets, compensates for the main deficiency in domestic and international existing embedded grids methodology, solving model assembly, there are the non-of relative motion Steady flow field Numerical Simulation Problems.
In a first aspect, providing a kind of parallel embedded grids methodology of achievable overlay region rapid Optimum, comprising:
The initial interpolation boundary of the nested grid of one group parts of label, a group parts include at least one moving component With at least one stationary parts;The information of flow of the nested grid characterization all parts;
The edge that initial interpolation boundary is pushed away to all parts is obtained positioned at the intermediate region interpolation side of the adjacent component Boundary, so that the overlapping region of the nested grid reaches minimum;
Determine that interpolation relationship of the interpolation boundary in the grid of the adjacent component, the interpolation relationship are indicated in institute State relative position of the interpolation boundary point of A component in adjacent component in the adjacent component in the grid of B component;
According to interpolation relationship, the information of flow in the grid of the B component is passed to the interpolation boundary of the A component Point.
Further, the initial interpolation boundary of the nested grid of one group parts of label, comprising:
The auxiliary cartesian grid of ambient surface grids is created for each component, marks the cartesian grid of inclusion surface grids Unit is edge cells;
Judge whether the grid of A component intersects with the auxiliary cartesian grid of other component;The A component is a group parts Any of component;
If so, judging whether the grid of the A component intersects with the edge cells of other component;
If so, whether the object plane grid cell for judging that the grid lines of A component and the edge cells of other component include intersects;
If so, the mesh point on the grid lines of label A component is initial interpolation boundary point.
Further, the edge that initial interpolation boundary is pushed away to all parts is obtained positioned at the adjacent component Intermediate region interpolation boundary, comprising:
Wall surface distance is calculated for the mesh point of each component;
Chained list is promoted for the initial interpolation boundary creation of A component, and sets the hole face reference point of each chained list point as itself point, Hole identity distance is from being 0;The hole face is the non-streaming site that initial interpolation boundary surrounds;
Promote A component propulsion chained list, it is more each promote chained list point wall surface distance and hole identity distance from;
If hole identity distance is from being less than with a distance from wall surface, marking chained list point is hole point, and flow field point adjacent thereto is marked to be labeled as New interpolation boundary point;
The propulsion chained list is deleted, new propulsion chained list is established;Above-mentioned forward operation is repeated, until A component does not generate newly Until interpolation point.
Further, interpolation relationship of the determination interpolation boundary in the grid of the adjacent component, comprising:
The outer boundary mesh point of marker motion component is interpolation boundary point.
Judge interpolation boundary point in A component whether the mesh overlay with other component;
If so, finding out the donor element of interpolation boundary point by intersection algorithm, and calculate corresponding interpolation system Number, the donor element are the grid cell of the B component where interpolation boundary point.
Further, the method is realized by multiple task parallelisms.
Second aspect provides a kind of parallel nested grid device of achievable overlay region rapid Optimum, comprising:
Mark module, the initial interpolation boundary of the nested grid for one group parts of label, a group parts include extremely A few moving component and at least one stationary parts;The information of flow of the nested grid characterization all parts;
Module is pushed away, for initial interpolation boundary to be pushed away to the edge of all parts, is obtained positioned at the adjacent component Intermediate region interpolation boundary, so that the overlapping region of the nested grid reaches minimum;
Determining module, it is described to insert for determining interpolation relationship of the interpolation boundary in the grid of the adjacent component Value relationship indicates phase of the interpolation boundary point of the A component in the adjacent component in the adjacent component in the grid of B component To position;
Transfer module, for according to interpolation relationship, the information of flow in the grid of the B component to be passed to the portion A The interpolation boundary point of part.
Further, the mark module includes:
Creating unit marks inclusion veil for creating the auxiliary cartesian grid of ambient surface grids for each component The cartesian grid unit of lattice is edge cells;
Judging unit, for judging whether the grid of A component intersects with the auxiliary cartesian grid of other component;The A Component is any of group parts component;If so, judge the A component grid whether the edge cells with other component Intersection;If so, whether the object plane grid cell for judging that the grid lines of A component and the edge cells of other component include intersects;If It is that then to mark mesh point on the grid lines of A component be initial interpolation boundary point.
Further, the module that pushes away includes:
Computing unit, for calculating wall surface distance for the mesh point of each component;
Table unit is created, promotes chained list for the initial interpolation boundary creation for A component, and set the hole face of each chained list point Reference point is itself point, and hole identity distance is from being 0;The hole face is the non-streaming site that initial interpolation boundary surrounds;
Comparing unit is promoted, for promoting the propulsion chained list of A component, more each wall surface distance for promoting chained list point and hole Identity distance from;
Hole dot element is marked, if marking chained list point is hole point, label and its phase for hole identity distance from being less than with a distance from wall surface Adjacent flow field point is labeled as new interpolation boundary point;
Unit is deleted, for deleting the propulsion chained list, establishes new propulsion chained list;Above-mentioned forward operation is repeated, until the portion A Until part does not generate new interpolation point.
Further, the determining module includes:
Demarcation of boundary dot element, the outer boundary mesh point for marker motion component are interpolation boundary point;
Judge overlapped elements, for judge the interpolation boundary point in A component whether the mesh overlay with other component;
Query unit for if so, find out the donor element of interpolation boundary point by intersection algorithm, and calculates Corresponding interpolation coefficient, the donor element are the grid cell of the B component where interpolation boundary point.
The third aspect provides a kind of computer-readable storage medium, instruction, the finger is stored on the storage medium The step of any one of first aspect the method is realized when order is executed by processor.
The utility model has the advantages that 1. are capable of handling the nested grid of any subregion, it include single portion in either each process sub-grid Part or multiple components are stronger to the universality of grid division;
2. only parallelization executes nested net in associated process by judging the nest relation of different components in each sub-grid Links in lattice method, parallel communications amount are less;
3. it is separate that interpolation boundary can be quickly propelled by comparing the wall surface distance and hole face distance parameter of interpolation boundary point Object plane boundary, flow field interpolation precision are higher;
4. under conditions of not changing main flow, need to only add interface letter for the CFD solver of non-nesting grid function Number can have nested function, preferably with solver compatibility.
Detailed description of the invention
Fig. 1 is the broad flow diagram of the method for the present invention;
Fig. 2 (a) is using overlay region parallel optimization result schematic diagram of the invention;
Fig. 2 (b) is using Flow Field Calculation result schematic diagram of the invention;
Fig. 3 is using the process schematic that drops a bomb of the invention.
Specific embodiment
Realization process of the invention is discussed in detail with reference to the accompanying drawing.
The present invention relates to aircraft components there are problems that the Numerical Investigation of Unsteady Flow Field of relative motion, as capsule ejection, Guided missile dispensing, propeller rotational etc. provide a kind of parallel embedded grids methodology of achievable overlay region rapid Optimum.This method packet It includes: obtaining subregion nested grid and the parallel communications domain of CFD solver, and the initial interpolation boundary of line flag nested grid;Parallel Optimize nested grid overlay region;Parallel creation nested grid interpolation relationship;It is parallel during CFD solver flow field calculation to update stream The flow field variable of field interpolation boundary point.
The invention discloses a kind of parallel embedded grids methodology of achievable overlay region rapid Optimum, methods as shown in Figure 1 Steps are as follows:
1] subregion nested grid and the parallel communications domain of CFD solver are obtained.
It should be noted that being capable of handling tetrahedron, pyramid, trigone the present invention is based on the exploitation of non-structural nested grid The grid cells type such as column, hexahedron and polyhedron.
2] and the initial interpolation boundary of line flag nested grid.
2.1] judge in each process sub-grid whether inclusion surface grids.If so, then creating the auxiliary of ambient surface grids Cartesian grid, the cartesian grid unit of label inclusion surface grids are edge cells;If nothing does not create.
2.2] some component is selected as target component, collects the cartesian grid spatial domain of other component in all processes With edge cells information, make Intersection judgement.If target component sub-grid intersects with cartesian grid spatial domain in current process, And intersect with edge cells, then carry out next step operation;Otherwise 2.4 are jumped to].
2.3] it sends the target component grid lines intersected with edge cells in the corresponding process of edge cells, judges net Whether the object plane grid cell that ruling and edge cells include intersects.If intersection, the mesh point on indicia grid line is initial Interpolation boundary point;Otherwise, it does not mark.After all grid lines have differentiated, by corresponding mesh point information feedback into current process.
2.4] circulation other component executes step 2.2 as target component]~2.4], until components all in each process All complete the label on initial interpolation boundary.
Since the subregion sub-grid of nested grid is from CFD solver, the partitioning strategies of each solver are different, lead It causes to include one or more component grid in final each process sub-grid, not only between sub-grid there may be nest relation, And nest relation is also likely to be present inside sub-grid.When multi-process is parallel, it is complicated that there are nested logical relations, parallel communications amount Big problem.Therefore, the method for the present invention by circulation target component, and 2.2] and 2.3] in addition decision criteria, realize Simplify nesting logical relation, minimizes the purpose of parallel communications amount.
3] parallel optimization nested grid overlay region.
3.1] the wall surface distance of all mesh points is calculated.
3.1.1 some component] is selected as target component, is judged whether in each process sub-grid containing target component Physical boundary point and partition boundaries point.If so, then creating propulsion chained list, next step operation is carried out;Otherwise 3.1.4 is jumped to].
3.1.2] determine that current process target component promotes the wall surface reference point and wall surface distance of each point in chained list.If chain Table point is object plane boundary point, then its wall surface reference point is itself, and wall surface distance is 0;Otherwise, target component in all processes is traversed Object plane boundary point, calculate the minimal wall identity distance of chained list point from and determining its wall surface reference point.
3.1.3] to promote chained list as front is promoted, the wall surface of the consecutive points of each chained list point is determined according to adjacent principle Reference point and wall surface distance, promote layer by layer, finally determine the wall surface distance of each all mesh points of process goal component.
3.1.4] circulation other component executes step 3.1.1 as target component]~3.1.4], until all component netlists Lattice point all obtains wall surface distance.
3.2] based on hole identity distance from wall surface with a distance from optimization interpolation boundary.
3.2.1 it] selects some component as target component, judges the partition boundaries point of target component in each process sub-grid Whether nested grid overlay region is located at.If so, traversing the initial interpolation boundary point of target component in all processes, calculate each The most duck eye identity distance of partition boundaries point from, and determine the face Qi Dong reference point.
3.2.2] judge whether in each process sub-grid containing the initial interpolation boundary point of target component.If so, then creating Chained list is promoted, and determines that the hole face reference point of each chained list point is itself, hole identity distance is from being 0;Otherwise 3.2.4 is jumped to].
3.2.3] to promote chained list as front is promoted, the wall surface distance and hole identity distance of more each chained list point are from if hole For identity distance from being less than with a distance from wall surface, then marking chained list point is hole point, marks flow field point adjacent thereto labeled as new interpolation boundary point. Former chained list is deleted, new propulsion chained list is established, repeats 3.2.3], until target component does not generate new interpolation point.
3.2.4] circulation other component executes step 3.2.1 as target component]~3.2.4], until all components are inserted The optimization of value boundary finishes.
Existing embedded grids methodology, which is also disclosed, carries out overlay region optimization based on wall surface distance, but which needs Wall surface distance of all the points in other component inside overlay region, by comparing mesh point in own component and in other component Wall surface distance, the overlay region Lai Youhua, optimization efficiency are not satisfactory.And the present invention innovatively uses hole identity distance from substituting grid Wall surface distance of the point in other component, algorithm is simple, pertains only to promote the mesh point on front, under paralleling tactic, optimization Significant effect.
Fig. 2 (a) is using overlay region parallel optimization result schematic diagram of the invention.Using two-dimentional three Duan Yiwei models, net Lattice are million magnitudes, optimize nested grid overlay region using 8 task parallelisms, optimize the time less than 1 second, final nested grid is inserted It is worth boundary far from object plane boundary, the different component grid adjoining dimensions of interpolation area are conducive to improve flow field interpolation precision.
4] nested grid interpolation relationship is created parallel.
4.1] in each process, the outer boundary mesh point of marker motion component is interpolation boundary point.
4.2] select some component as target component, judge interpolation boundary point in current process target component whether with Other component mesh overlay in other processes.If overlapping, sends these interpolation points in corresponding process, passes through line face Intersection algorithm finds out the donor element of interpolation point, and calculates corresponding interpolation coefficient.
4.3] circulation other component executes step 4.2 as target component], until components all in each process all interpolation Boundary point determines donor element and interpolation coefficient.
Fig. 2 (b) is using Flow Field Calculation result schematic diagram of the invention.Mach number isopleth in figure is in interpolation area light It slips over and crosses, even if Mach number also matches preferably in the interpolation area of leading edge slat seam and trailing edge flap tail.
5] according to the interpolation relationship of nested grid, the information of flow of flow field interpolation boundary point is updated parallel for CFD solver.
Fig. 3 is using the process schematic that drops a bomb of the invention.It is shown in figure within a predetermined period of time, guided missile is from initial bit Beginning is set, under self gravity, the gentle Dynamic Synthesis effect of ejection force, the posture changing of guided missile in launch process.
Finally it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although Present invention has been described in detail with reference to the aforementioned embodiments, and those skilled in the art should understand that;It still can be right The technical solution that foregoing embodiments are recorded is modified or equivalent replacement of some of the technical features;And these It modifies or replaces, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.

Claims (10)

1. a kind of parallel embedded grids methodology of achievable overlay region rapid Optimum characterized by comprising
The initial interpolation boundary of the nested grid of one group parts of label, a group parts include at least one moving component and extremely A few stationary parts;The information of flow of the nested grid characterization all parts;
The edge that initial interpolation boundary is pushed away to all parts is obtained positioned at the intermediate region interpolation boundary of the adjacent component, So that the overlapping region of the nested grid reaches minimum;
Determine that interpolation relationship of the interpolation boundary in the grid of the adjacent component, the interpolation relationship are indicated in the phase Relative position of the interpolation boundary point of A component in the adjacent component in the grid of B component in adjacent component;
According to interpolation relationship, the information of flow in the grid of the B component is passed to the interpolation boundary point of the A component.
2. the method according to claim 1, wherein the initial interpolation of the nested grid of one group parts of the label Boundary, comprising:
The auxiliary cartesian grid of ambient surface grids is created for each component, marks the cartesian grid unit of inclusion surface grids For edge cells;
Judge whether the grid of A component intersects with the auxiliary cartesian grid of other component;The A component is to appoint in a group parts One component;
If so, judging whether the grid of the A component intersects with the edge cells of other component;
If so, whether the object plane grid cell for judging that the grid lines of A component and the edge cells of other component include intersects;
If so, the mesh point on the grid lines of label A component is initial interpolation boundary point.
3. according to the method described in claim 2, it is characterized in that, the side that initial interpolation boundary is pushed away to all parts Edge is obtained positioned at the intermediate region interpolation boundary of the adjacent component, comprising:
Wall surface distance is calculated for the mesh point of each component;
Chained list is promoted for the initial interpolation boundary creation of A component, and sets the hole face reference point of each chained list point as itself point, hole face Distance is 0;The hole face is the non-streaming site that initial interpolation boundary surrounds;
Promote A component propulsion chained list, it is more each promote chained list point wall surface distance and hole identity distance from;
If hole identity distance, from being less than with a distance from wall surface, marking chained list point is hole point, flow field point adjacent thereto is marked to insert labeled as new It is worth boundary point;
The propulsion chained list is deleted, new propulsion chained list is established;Above-mentioned forward operation is repeated, until A component does not generate new interpolation Until point.
4. according to the method described in claim 3, it is characterized in that, the determination interpolation boundary is in the adjacent component Interpolation relationship in grid, comprising:
The outer boundary mesh point of marker motion component is interpolation boundary point.
Judge interpolation boundary point in A component whether the mesh overlay with other component;
If so, finding out the donor element of interpolation boundary point by intersection algorithm, and corresponding interpolation coefficient is calculated, institute State the grid cell that donor element is the B component where interpolation boundary point.
5. method according to claim 1-4, which is characterized in that the method is real by multiple task parallelisms It is existing.
6. a kind of parallel nested grid device of achievable overlay region rapid Optimum characterized by comprising
Mark module, the initial interpolation boundary of the nested grid for one group parts of label, a group parts include at least one A moving component and at least one stationary parts;The information of flow of the nested grid characterization all parts;
Module is pushed away, for initial interpolation boundary to be pushed away to the edge of all parts, is obtained positioned at the centre of the adjacent component REGION INTERPOLATION boundary, so that the overlapping region of the nested grid reaches minimum;
Determining module, for determining that interpolation relationship of the interpolation boundary in the grid of the adjacent component, the interpolation are closed System indicates opposite position of the interpolation boundary point of the A component in the adjacent component in the adjacent component in the grid of B component It sets;
Transfer module, for according to interpolation relationship, the information of flow in the grid of the B component to be passed to the A component Interpolation boundary point.
7. the apparatus according to claim 1, which is characterized in that the mark module includes:
Creating unit, for creating the auxiliary cartesian grid of ambient surface grids for each component, label inclusion surface grids Cartesian grid unit is edge cells;
Judging unit, for judging whether the grid of A component intersects with the auxiliary cartesian grid of other component;The A component For any of group parts component;If so, judge the A component grid whether the edge cells phase with other component It hands over;If so, whether the object plane grid cell for judging that the grid lines of A component and the edge cells of other component include intersects;If so, Then marking the mesh point on the grid lines of A component is initial interpolation boundary point.
8. device according to claim 7, which is characterized in that the module that pushes away includes:
Computing unit, for calculating wall surface distance for the mesh point of each component;
Table unit is created, promotes chained list for the initial interpolation boundary creation for A component, and set the hole face reference of each chained list point Point is itself point, and hole identity distance is from being 0;The hole face is the non-streaming site that interpolation boundary surrounds;
Comparing unit is promoted, for promoting the propulsion chained list of A component, more each wall surface distance and hole identity distance for promoting chained list point From;
Hole dot element is marked, if marking chained list point is hole point for hole identity distance from being less than with a distance from wall surface, is marked adjacent thereto Flow field point is labeled as new interpolation boundary point;
Unit is deleted, for deleting the propulsion chained list, establishes new propulsion chained list;Above-mentioned forward operation is repeated, until A component does not have Have until generating new interpolation point.
9. device according to claim 8, which is characterized in that the determining module includes:
Demarcation of boundary dot element, the outer boundary mesh point for marker motion component are interpolation boundary point;
Judge overlapped elements, for judge the interpolation boundary point in A component whether the mesh overlay with other component;
Query unit for if so, find out the donor element of interpolation boundary point by intersection algorithm, and calculates corresponding Interpolation coefficient, the donor element be interpolation boundary point where B component grid cell.
10. a kind of computer-readable storage medium, instruction is stored on the storage medium, which is characterized in that described instruction The step of any one of claim 1-5 the method is realized when being executed by processor.
CN201910384986.XA 2019-05-09 2019-05-09 Parallel nested grid method and device capable of realizing rapid optimization of overlapping area Active CN110110446B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910384986.XA CN110110446B (en) 2019-05-09 2019-05-09 Parallel nested grid method and device capable of realizing rapid optimization of overlapping area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910384986.XA CN110110446B (en) 2019-05-09 2019-05-09 Parallel nested grid method and device capable of realizing rapid optimization of overlapping area

Publications (2)

Publication Number Publication Date
CN110110446A true CN110110446A (en) 2019-08-09
CN110110446B CN110110446B (en) 2023-01-13

Family

ID=67489075

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910384986.XA Active CN110110446B (en) 2019-05-09 2019-05-09 Parallel nested grid method and device capable of realizing rapid optimization of overlapping area

Country Status (1)

Country Link
CN (1) CN110110446B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111859825A (en) * 2020-07-28 2020-10-30 北京航空航天大学 Method and equipment for simulating unsteady non-pressure flow numerical value with arbitrary flow-solid interface
CN112100937A (en) * 2020-09-11 2020-12-18 南京航空航天大学 Efficient overlapped grid parallel assembly method based on wall surface distance
CN113095012A (en) * 2021-05-17 2021-07-09 华能陇东能源有限责任公司 Splicing and fusing method for numerical simulation calculation results of wind power plant flow field partitions
CN114119929A (en) * 2022-01-26 2022-03-01 南京航空航天大学 Object plane array propulsion parallel computing method in aircraft flow field analysis
CN114119882A (en) * 2022-01-25 2022-03-01 南京航空航天大学 Efficient nested grid host unit searching method in aircraft dynamic flow field analysis
CN115357849A (en) * 2022-10-24 2022-11-18 中国空气动力研究与发展中心计算空气动力研究所 Method and device for calculating wall surface distance under Cartesian grid

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101599104A (en) * 2009-07-16 2009-12-09 北京航空航天大学 A kind of analogy method of blade flutter boundary of aviation turbine engine
US20150260016A1 (en) * 2014-03-17 2015-09-17 Saudi Arabian Oil Company Modeling intersecting faults and complex wellbores in reservoir simulation
CN108459909A (en) * 2018-02-27 2018-08-28 北京临近空间飞行器系统工程研究所 A kind of Multi-bodies Separation mesh overlay method and system suitable for parallel processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101599104A (en) * 2009-07-16 2009-12-09 北京航空航天大学 A kind of analogy method of blade flutter boundary of aviation turbine engine
US20150260016A1 (en) * 2014-03-17 2015-09-17 Saudi Arabian Oil Company Modeling intersecting faults and complex wellbores in reservoir simulation
CN108459909A (en) * 2018-02-27 2018-08-28 北京临近空间飞行器系统工程研究所 A kind of Multi-bodies Separation mesh overlay method and system suitable for parallel processing

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111859825A (en) * 2020-07-28 2020-10-30 北京航空航天大学 Method and equipment for simulating unsteady non-pressure flow numerical value with arbitrary flow-solid interface
CN112100937A (en) * 2020-09-11 2020-12-18 南京航空航天大学 Efficient overlapped grid parallel assembly method based on wall surface distance
CN112100937B (en) * 2020-09-11 2021-05-28 南京航空航天大学 Efficient overlapped grid parallel assembly method based on wall surface distance
CN113095012A (en) * 2021-05-17 2021-07-09 华能陇东能源有限责任公司 Splicing and fusing method for numerical simulation calculation results of wind power plant flow field partitions
CN113095012B (en) * 2021-05-17 2022-11-29 华能陇东能源有限责任公司 Splicing and fusing method for numerical simulation calculation results of wind power plant flow field partitions
CN114119882A (en) * 2022-01-25 2022-03-01 南京航空航天大学 Efficient nested grid host unit searching method in aircraft dynamic flow field analysis
CN114119882B (en) * 2022-01-25 2022-07-12 南京航空航天大学 Efficient nested grid host unit searching method in aircraft dynamic flow field analysis
CN114119929A (en) * 2022-01-26 2022-03-01 南京航空航天大学 Object plane array propulsion parallel computing method in aircraft flow field analysis
CN114119929B (en) * 2022-01-26 2022-07-12 南京航空航天大学 Object plane array propulsion parallel computing method in aircraft flow field analysis
CN115357849A (en) * 2022-10-24 2022-11-18 中国空气动力研究与发展中心计算空气动力研究所 Method and device for calculating wall surface distance under Cartesian grid

Also Published As

Publication number Publication date
CN110110446B (en) 2023-01-13

Similar Documents

Publication Publication Date Title
CN110110446A (en) A kind of parallel embedded grids methodology and device of achievable overlay region rapid Optimum
Forsythe et al. Coupled flight simulator and CFD calculations of ship airwake using kestrel
CN105700549A (en) Unmanned plane multi-track planning method based on sequence ecological niche PSO (particle swarm optimization) algorithm
Kampolis et al. A multilevel approach to single-and multiobjective aerodynamic optimization
CN106529674A (en) Multiple-unmanned-aerial-vehicle cooperated multi-target distribution method
CN107832551A (en) A kind of modularization Architecture simulation system and method towards Space Equipment
CN108268042A (en) A kind of path planning algorithm based on improvement Visual Graph construction
CN111006693A (en) Intelligent aircraft track planning system and method thereof
CN104484500A (en) Air combat behavior modeling method based on fitting reinforcement learning
CN107977491A (en) The Aerodynamic Heating appraisal procedure in aircraft airvane gap in the case of a kind of unstable state
CN113156954A (en) Multi-agent cluster obstacle avoidance method based on reinforcement learning
CN109978286A (en) It is a kind of to be diversion thunderstorm Route planner based on the more aircrafts for improving ant group algorithm
CN108459505A (en) A kind of unconventional layout aircraft fast modeling method of suitable control Iterative Design
BR102019027269A2 (en) to optimize a parametric model of aircraft performance
CN109902359A (en) The housing construction optimum design method of unmanned plane race, Flying-wing
Tang et al. Multi-level hybridized optimization methods coupling local search deterministic and global search evolutionary algorithms
Dong et al. Trial input method and own-aircraft state prediction in autonomous air combat
Chen et al. Adaptive multiple task assignments for uavs using discrete particle swarm optimization
Ciarella et al. A Multi-fidelity, multi-disciplinary analysis and optimization framework for the design of morphing UAV wing
Gao et al. A constraint approximation assisted PSO for computationally expensive constrained problems
CN112800679A (en) Photoelectric equipment defense planning method, device, equipment and storage medium
CN116050304B (en) Intelligent fish flow field simulation control method, system, equipment and storage medium
Qiannan et al. Path planning of UAVs formation based on improved ant colony optimization algorithm
Lofthouse et al. Static and dynamic simulations of a generic UCAV geometry using the kestrel flow solver
Guzelbey et al. A Review of Aerodynamic Shape Optimization for a Missile

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