CN106202674B - A kind of two-dimension numerical calculation method of non-axis symmetry element - Google Patents
A kind of two-dimension numerical calculation method of non-axis symmetry element Download PDFInfo
- Publication number
- CN106202674B CN106202674B CN201610515885.8A CN201610515885A CN106202674B CN 106202674 B CN106202674 B CN 106202674B CN 201610515885 A CN201610515885 A CN 201610515885A CN 106202674 B CN106202674 B CN 106202674B
- Authority
- CN
- China
- Prior art keywords
- axis symmetry
- circumferential speed
- model
- symmetry element
- physical model
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of two-dimension numerical calculation method of non-axis symmetry element characterized by comprising simplifies the physical model according to the physical model of non-axis symmetry element;Two-dimentional computation model is established according to the type of the physical model, determines the range of computational domain, and according to the computational domain grid division;Circumferential speed is established than model according to the type of the physical model, determines the range of fluid domain, and calculates the circumferential speed ratio and empty revolving speed of the fluid domainAccording to the grid and the empty revolving speedCalculate element entrance port velocity vector Vi, element entrance port velocity vector V is calculated according to grid and empty revolving speediIf ViLast value V with iti‑1Difference be less than given threshold V0, then the two-dimension numerical calculation is completed, and otherwise recalculates circumferential speed ratio.The present invention solves the problems, such as that non-axis symmetry element loses element in Two-dimensional numerical simulation and causes error to the effect of fluid, improves the accuracy and practicability of air system two-dimensional analog.
Description
Technical field
The present invention relates to air system Two-dimensional numerical simulation method more particularly to a kind of two Dimension Numerical Values of non-axis symmetry element
Calculation method.
Background technique
Element is abundant in the air system flow path of gas turbine, and flow path is complicated, can directly carry out complete machine air system three
The computer prices valuableness difficulty of dimension value is larger, and the calculating analysis means for carrying out flow path level at present are mainly one-dimensional flow path meter
The two-dimensional/three-dimensional for calculating cooperation part calculates, and the high-efficient but simplified processing of one-dimensional calculating is more, not to the reflection of internal flow details
Foot, and the details that local two-dimensional/three-dimensional simulation obtains is enriched but computational efficiency is low, and cannot reflect the air system overall situation, these
The scope of application of method is restricted, and the reference value of calculated result is limited.It is owed since the technology of processing asymmetric element exists
It lacks, limits two-dimentional CFD (Computational Fluid Dynamics, Fluid Mechanics Computation) technology in gas turbine air
The application of system regions.
There are the element largely rotated or structure in air system, this class component or structure are by movement to flowing through therebetween
Air works.Have from structure completely about the symmetrical element of engine shaft, such as disk chamber, central aeration pipe;Also have
About the asymmetric element of shaft, such as with engine axis out-of-alignment hole, nozzle, gap.Axial symmetry element is pressed into the faces of grade
It is smaller that product principle is reduced to error when two dimensional model, but non-axis symmetry element then loses element in Two-dimensional numerical simulation to fluid
Effect and circumferential flow information.Therefore one can be brought when simple two dimensionization numerical simulation is applied to engine air system
Series of problems uses preceding method for example, the flow-through element of rotation is in rotation state or element local environment is Qiang Xuanliu at this time
Simulation can lose a large amount of flowing information, such as rotating in a circumferential direction strongly for preswirl nozzle exit flow, when air-flow via hole air-flow with it is quiet
Interaction etc. between only/rotating hole.Large error will be brought by doing simple simplification to the asymmetric element of this kind of occasion.
Summary of the invention
The present invention provides a kind of two-dimension numerical calculation method of non-axis symmetry element, solves non-axis symmetry element in two-dimemsional number
The problem of losing element larger simulation error caused by the effect and circumferential flow information of fluid when value simulation, improves air
The accuracy and practicability of system two-dimensional analog.
In order to achieve the above object, the present invention the following technical schemes are provided:
A kind of two-dimension numerical calculation method of non-axis symmetry element, comprising the following steps:
Simplify the physical model according to the physical model of non-axis symmetry element;
Two-dimentional computation model is established according to the type of the physical model, determines the range of computational domain, and according to the meter
Calculate domain grid division;
Circumferential speed is established according to the type of the physical model and determines the range of fluid domain than model, and described in calculating
The circumferential speed ratio of fluid domain and empty revolving speedThe circumferential speed is than the circumferential speed before and after characterizing the non-axis symmetry element
Variable quantity, the void revolving speedCharacterize the effect of the non-axis symmetry element stream field;
According to the grid and the empty revolving speedCalculate element entrance port velocity vector ViIf the element entrance port speed
Spend vector Vi-Vi-1Difference be less than given threshold V0, then the two-dimension numerical calculation is completed, and is otherwise returned and is calculated the empty revolving speed
Preferably, described to include: according to the simplified physical model of the physical model of non-axis symmetry element
By etc. circulation areas principle simplify circulating area;
Remove the non-axis symmetry element corner little to flow effect, escape;
Simplify gasket, the locking plate, counterweight of non-axis symmetry element;
Merge the adjacent structure with same movement feature.
Preferably, described that two-dimentional computation model is established according to the type of the physical model, determine the range of computational domain, and
According to the computational domain grid division, comprising:
Determine the model treatment method of hole class or convex type;
Compromise to mutual conflicting non-axis symmetry element.
Preferably, the type according to the physical model establishes circumferential speed and compares model, comprising: hole class circumferential speed
Compare model than model and convex type circumferential speed.
Preferably, the circumferential speed variable quantity is circumferential speed ratio CV, wherein Cv=V2,θ/V1,θ, V2,θFor entrance gas
The average circumferential speed of stream, V1,θFor the average circumferential speed of exit air-flow.
Preferably, circumferential speed ratio CVFor
Wherein, l is the thickness in hole, and d is aperture, and α is inlet air flow and axially bored line angle, and β is radial angle, γ critical angle
Degree, A are the area in hole, AcrossArea, Sr are led directly to for the circulation in hole1、Sr2For swirl ratio.
Preferably, the empty revolving speed includes the empty revolving speed of hole class component and the empty revolving speed of raised class component;
The empty rotational speed omega of the hole class component1=(Cv-1)V1,θ/ rx, wherein V1,θFor import circumferential speed, r2For outlet
Radial position;
It is described protrusion class component empty revolving speed be
Wherein, V1,θFor protruding component import circumferential speed, r is protruding component radius of turn, and d is that protruding component equivalent is straight
Diameter, ω ' are the revolving speed of rotation system.
The present invention provides a kind of two-dimension numerical calculation method of non-axis symmetry element, empty by the way that non-axis symmetry element is arranged
The method of rotation obtains an empty revolving speed, keeps the effect of non-axis exchange element stream field identical with Real Flow Field.Solution axis pair by no means
Element is claimed to lose element very big simulation error caused by the effect and circumferential flow information of fluid in Two-dimensional numerical simulation
The problem of, improve the accuracy and practicability of air system two-dimensional analog.
Detailed description of the invention
In order to illustrate more clearly of specific embodiments of the present invention, attached drawing needed in the embodiment will be made below
Simply introduce.
Fig. 1: being a kind of two-dimension numerical calculation method flow diagram of non-axis symmetry element provided by the invention;
Fig. 2: being a kind of two-dimentional CFD analysis flow chart diagram based on empty spinning solution provided in an embodiment of the present invention;
Fig. 3: pore model schematic diagram provided in an embodiment of the present invention;
Fig. 4: being a kind of dish structure provided in an embodiment of the present invention and its air flow path schematic diagram;
Fig. 5: being dish structure protrusion streamline contrast schematic diagram provided in an embodiment of the present invention;
Fig. 6: being dish structure pressure contour contrast schematic diagram provided in an embodiment of the present invention;
Fig. 7: being dish structure circumferential speed isopleth contrast schematic diagram provided in an embodiment of the present invention.
Specific embodiment
The scheme of embodiment in order to enable those skilled in the art to better understand the present invention with reference to the accompanying drawing and is implemented
Mode is described in further detail the embodiment of the present invention.
There is the effect for ignoring its circumferential flow to asymmetric element for current two dimension CFD technology, is easily making two-dimemsional number
Value simulation causes very big error with Real Flow Field, and the present invention provides a kind of two-dimension numerical calculation method of non-axis symmetry element,
Method by non-axis symmetry element being arranged empty rotation obtains an empty revolving speed, makes the effect of non-axis exchange element stream field
It is identical with Real Flow Field.It solves non-axis symmetry element and loses effect and peripheral flow of the element to fluid in Two-dimensional numerical simulation
Caused by dynamic information very big simulation error the problem of, improve the accuracy and practicability of air system two-dimensional analog.
As shown in Figure 1, be the two-dimension numerical calculation method flow diagram of non-symmetric element of one kind provided by the invention, including
Following steps:
S1: the physical model is simplified according to the physical model of non-axis symmetry element;
S2: two-dimentional computation model is established according to the type of the physical model, determines the range of computational domain, and according to described
Computational domain grid division;
S3: circumferential speed is established according to the type of the physical model and determines the range of fluid domain than model, and calculates institute
State the circumferential speed ratio and empty revolving speed of fluid domainThe circumferential speed is than the circumferential speed before and after characterizing the non-axis symmetry element
Spend variable quantity, the void revolving speedCharacterize the effect of the non-axis symmetry element stream field;
S4: according to the grid and the empty revolving speedCalculate element entrance port velocity vector ViIf the element into
Mouth velocity vector ViLast value V with iti-1Difference be less than given threshold V0, then the two-dimension numerical calculation is completed, otherwise
Recalculate the empty revolving speed
Specifically, as shown in Fig. 2, for a kind of two-dimentional CFD analysis based on empty spinning solution provided in an embodiment of the present invention
Flow chart, comprising the following steps:
S11: the input of input this process of physical model is the physical model and boundary condition of research object, and calculating purpose is
Obtain the information of flow such as circumferential speed distribution, the chamber pressure distribution of research object;
S12: simplified model does model cleaning and simplifies, remove convection current to the CAD diagram paper or physical model of research object
It is dynamic to influence little corner, escape etc., simplify gasket, locking plate, counterweight etc., merges adjacent or with same movement feature knot
Structure obtains simplified model;
S13: the flow path of identification asymmetric element analysis physical model confirms the structural parameters and fortune of each element in model
Dynamic information, distinguishes symmetric element and non-axis symmetry element, identifies the type and motion morphology of non-axis symmetry element;
S14: two dimensionization establishes the two dimensionization principle and method of each non-axis symmetry element after S13;According to demand and
The range of computational domain is determined to the influence of calculated result;
S15: establishing two dimensional model and do planarization process to non-axis symmetry element according to the principle determined after S14,
It mutual conflicting non-axis symmetry element should should influence do to compromise according to respective in the process, and then incorporate entire flow path computational domain,
Establish two-dimentional computation model;Grid dividing is done to computational domain by the requirement of numerical simulation;
S16: grid, which is imported, imports solver for grid after S15, and boundary condition is arranged, completes according to a conventional method for the first time
It calculates;
S17: circumferential speed is established than model and identifies non-axis symmetry element final finishing identified non-axis symmetry member in S13
Part type prejudges flow regime to each non-axis symmetry element, establishes circumferential speed accordingly and compare model;
S18: empty revolving speed is calculated after S17, each non-axis symmetry element entrance port velocity component is read, according to the speed of import
Component calculates circumferential speed ratio, then calculates the element in the empty revolving speed of the state;Empty revolving speed is assigned to each non-axis symmetry member
Part is calculated in solver to meeting convergence criterion;
S19: the non-axis symmetry element entrance port speed that epicycle CFD is calculated and the difference that upper wheel calculates are checked, if difference is less than
Designated value stops calculating.It is post-processed.Otherwise S17 step is returned.
Further, described to include: according to the simplified physical model of the physical model of non-axis symmetry element
Step 1: the removal non-axis symmetry element corner little to flow effect, escape;
Step 2: simplify gasket, the locking plate, counterweight of non-axis symmetry element;
Step 3: merge the adjacent structure with same movement feature.
Meanwhile it is described according to the two-dimentional computation model of the type of physical model foundation, determine the range of computational domain, and root
According to the computational domain grid division, comprising: determine the model treatment method of hole class or convex type and to mutual conflicting non-axis pair
Element is claimed to compromise.
Further, the type according to the physical model establishes circumferential speed and compares model, comprising: hole class circumferential speed
Compare model than model and convex type circumferential speed.
Circumferential speed variable quantity is by circumferential speed ratio CVIt characterizes, wherein Cv=V2,θ/V1,θ, V2,θFor being averaged for entrance air-flow
Circumferential speed, V1,θFor the average circumferential speed of exit air-flow.
Circumferential speed ratio CVFor
Wherein,
L is the thickness in hole, and d is aperture, and α is inlet air flow and axially bored line angle, and β is radial angle, γ critical angle, A
For the area in hole, AcrossArea, Sr are led directly to for the circulation in hole1、Sr2For swirl ratio.It can refer to shown in Fig. 3, implement for the present invention
The pore model schematic diagram that example provides.
Further, the empty revolving speed includes the empty revolving speed of hole class component and the empty revolving speed of raised class component;The hole class member
The empty rotational speed omega of part1=(Cv-1)V1,θ/r2, wherein V1,θFor import circumferential speed, r2To export radial position;The convex type
The empty revolving speed of element, wherein V1,θFor protruding component import circumferential speed, r is protruding component rotation
Turn radius, d is protruding component equivalent diameter, and ω ' is the revolving speed of rotation system.
In practical applications, as shown in figure 4, being a kind of dish structure provided in an embodiment of the present invention and its air flow path signal
Figure, the cavity of the dish structure is air flow circuit, and air-flow enters disk chamber from below, is flowed out from hole, disk is intracavitary bolt, disk 1,2
In rotation status, wherein bolt and hole are two non-axis symmetry elements in model.It is outside that air-flow enters disk chamber rear portion
Flowed out directly through hole, another part needs to flow through the space between bolt, and after flowed out through hole, i.e., the region essence between bolt
Upper is also airflow channel, and the two-dimension numerical calculation method through the asymmetric shaft element can get such as Fig. 5, Fig. 6, effect shown in Fig. 7
Fruit, Fig. 5 are dish structure protrusion streamline contrast schematic diagram provided in an embodiment of the present invention, and Fig. 6 is disk provided in an embodiment of the present invention
Structural pressure isopleth contrast schematic diagram, from Fig. 6 as it can be seen that conventional method calculate as the result is shown the low radius pressure loss compared with
Greatly, chamber pressure in top has a certain recovery, and is calculated using herein described method, and low radius loss is smaller, and top chamber pressure increases with radius
Increase greatly, consistent with practical flow phenomenon, pressure contour caused by the intracavitary stirring as protrusion of disk under is intensive, pressure
Gradient is bigger than conventional method.Fig. 7 is dish structure circumferential speed isopleth contrast schematic diagram provided in an embodiment of the present invention, from
As it can be seen that circumferential speed of the air-flow before and after rotating hole varies less in the analog result of conventional method, isopleth is in longitudinal row
Column, are then observed from the calculation method result of the application, and circumferential speed has significant change before and after air-flow passes through rotating hole, equivalent
Line morphology and conventional method result are dramatically different.The circumferential speed point of protruding component position is also obtained using empty spinning solution
Cloth.
As it can be seen that the present invention provides a kind of two-dimension numerical calculation method of non-axis symmetry element, by non-axis symmetry element
The method of empty rotation is set, an empty revolving speed is obtained, keeps the effect of non-axis exchange element stream field identical with Real Flow Field.It solves
Non-axis symmetry element loses element very big mould caused by the effect and circumferential flow information of fluid in Two-dimensional numerical simulation
The problem of quasi- error, improves the accuracy and practicability of air system two-dimensional analog.
Structure, feature and effect of the invention, the above institute is described in detail according to diagrammatically shown embodiment above
Only presently preferred embodiments of the present invention is stated, but the present invention does not limit the scope of implementation as shown in the drawings, it is all according to structure of the invention
Think made change or equivalent example modified to equivalent change, when not going beyond the spirit of the description and the drawings,
It should all be within the scope of the present invention.
Claims (4)
1. a kind of two-dimension numerical calculation method of non-axis symmetry element, which comprises the following steps:
Simplify the physical model according to the physical model of non-axis symmetry element;
Two-dimentional computation model is established according to the type of the physical model, determines the range of computational domain, and according to the computational domain
Grid division;
Circumferential speed is established according to the type of the physical model and determines the range of fluid domain than model, and calculates the fluid
The circumferential speed ratio in domain and empty revolving speed, the circumferential speed change than the circumferential speed before and after characterizing the non-axis symmetry element
Amount, the void revolving speed characterize the effect of the non-axis symmetry element stream field;
Element entrance port velocity vector V is calculated according to the grid and the empty revolving speediIf the element entrance port velocity vector
ViLast value V with iti-1Difference be less than given threshold V0, then the two-dimension numerical calculation is completed, and otherwise recalculates void
Revolving speed.
2. the two-dimension numerical calculation method of non-axis symmetry element according to claim 1, which is characterized in that described according to institute
The type for stating physical model establishes two-dimentional computation model, determines the range of computational domain, and according to the computational domain grid division, packet
It includes:
Establish the model of hole class or convex type;
Compromise to mutual conflicting non-axis symmetry element.
3. the two-dimension numerical calculation method of non-axis symmetry element according to claim 1, which is characterized in that described according to institute
The type for stating physical model establishes circumferential speed and compares model, comprising: hole class circumferential speed is than model and convex type circumferential speed ratio
Model.
4. the two-dimension numerical calculation method of non-axis symmetry element according to any one of claims 1 to 3, which is characterized in that
The circumferential speed variable quantity is circumferential speed ratio CV, wherein Cv=V2,θ/V1,θ, V2,θFor the average circumferential speed of entrance air-flow
Degree, V1,θFor the average circumferential speed of exit air-flow.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610515885.8A CN106202674B (en) | 2016-07-01 | 2016-07-01 | A kind of two-dimension numerical calculation method of non-axis symmetry element |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610515885.8A CN106202674B (en) | 2016-07-01 | 2016-07-01 | A kind of two-dimension numerical calculation method of non-axis symmetry element |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106202674A CN106202674A (en) | 2016-12-07 |
CN106202674B true CN106202674B (en) | 2019-06-18 |
Family
ID=57465528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610515885.8A Active CN106202674B (en) | 2016-07-01 | 2016-07-01 | A kind of two-dimension numerical calculation method of non-axis symmetry element |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106202674B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104346499A (en) * | 2014-11-19 | 2015-02-11 | 上海交通大学 | Multi-fan turbine engine design method based on computer platform |
CN105378735A (en) * | 2013-07-09 | 2016-03-02 | 斯奈克玛 | Method for modeling non-axisymmetric surface |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110072772A1 (en) * | 2008-05-22 | 2011-03-31 | Enertechnix, Inc | Skimmer for Concentrating an Aerosol and Uses Thereof |
-
2016
- 2016-07-01 CN CN201610515885.8A patent/CN106202674B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105378735A (en) * | 2013-07-09 | 2016-03-02 | 斯奈克玛 | Method for modeling non-axisymmetric surface |
CN104346499A (en) * | 2014-11-19 | 2015-02-11 | 上海交通大学 | Multi-fan turbine engine design method based on computer platform |
Non-Patent Citations (3)
Title |
---|
发动机进气系统不均匀性的三维数值模拟;王晗等;《小型内燃机与摩托车》;20070630;第36卷(第3期);第41-44页 |
微型涡轮的旋转盘腔内流场二维与三维数值模拟的比较;杨小利等;《机械设计与制造》;20050831(第8期);第69-71页 |
接收孔元件的周向速度系数模型与验证;张建超等;《重庆理工大学学报》;20140731;第28卷(第7期);第43-49页 |
Also Published As
Publication number | Publication date |
---|---|
CN106202674A (en) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guha et al. | Secondary motion in three-dimensional branching networks | |
Minemura et al. | A theoretical study on air bubble motion in a centrifugal pump impeller | |
CN104361625B (en) | A kind of cloud data compaction algorithm retained with border based on ray principle | |
CN106202674B (en) | A kind of two-dimension numerical calculation method of non-axis symmetry element | |
CN103995938B (en) | A kind of compressor Middle casing end wall profile formative method | |
Ahasan et al. | Effect of angle-of-attacks on deterministic lateral displacement (DLD) with symmetric airfoil pillars | |
CN106762811A (en) | The pipeline compressor model of discharge coefficient 0.0495 level and method for designing impeller | |
Machado et al. | Local Extraction of Bifurcation Lines. | |
Ji et al. | Computer 3D vision-aided full-3D optimization of a centrifugal impeller | |
CN109815590A (en) | A kind of multi stage axial flow compressor three dimendional blade formative method and blade based on petiolarea boundary-layer | |
CN106794391A (en) | For the new extraction cells of centrifugal distribution chromatograph, and the centrifugal distribution chromatograph comprising such extraction cells | |
CN106704256A (en) | Model level of pipeline compressor with discharge coefficient being 0.0322 and design method of impeller thereof | |
CN107066686A (en) | A kind of axial-flow pump impeller Hydraulic Optimizing Design method based on genetic algorithm | |
CN113393462B (en) | Focal region segmentation method and device, electronic equipment and storage medium | |
CN113361028B (en) | Two-dimensional design method of volute | |
CN114692338B (en) | Comprehensive optimization design method for cavitation and efficiency of low-temperature centrifugal pump | |
Song et al. | Numerical study of a bio‐centrifugal blood pump with straight impeller blade profiles | |
Song et al. | Numerical study of a centrifugal blood pump with different impeller profiles | |
Han et al. | Flow Control of Radial Inlet Chamber and Downstream Effects on a Centrifugal Compressor Stage | |
CN108561331A (en) | 0.0264 light substance high energy head compressor model grade of discharge coefficient and design method | |
Dix et al. | Vortex simulations on a 3-sphere | |
CN107906047A (en) | 0.0472 light substance high energy head compressor model level of discharge coefficient and design method | |
CN107906050A (en) | 0.014 light substance high energy head compressor model level of discharge coefficient and design method | |
Akgöl et al. | Influence of Geometric Parameters on the Average Outlet Velocity of the Bladeless Fan | |
CN215113825U (en) | Heat exchange air equalizing system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |