CN114444417B - Vortex flow form judging method suitable for curve vortex shaft and storage medium - Google Patents

Vortex flow form judging method suitable for curve vortex shaft and storage medium Download PDF

Info

Publication number
CN114444417B
CN114444417B CN202210373310.2A CN202210373310A CN114444417B CN 114444417 B CN114444417 B CN 114444417B CN 202210373310 A CN202210373310 A CN 202210373310A CN 114444417 B CN114444417 B CN 114444417B
Authority
CN
China
Prior art keywords
vortex
shaft
axis
flow
characteristic value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210373310.2A
Other languages
Chinese (zh)
Other versions
CN114444417A (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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202210373310.2A priority Critical patent/CN114444417B/en
Publication of CN114444417A publication Critical patent/CN114444417A/en
Application granted granted Critical
Publication of CN114444417B publication Critical patent/CN114444417B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

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

The invention is suitable for the technical field of fluid mechanics, provides a vortex flow form judging method suitable for a curved vortex shaft and a storage medium, establishes an orthogonal curve coordinate system taking the vortex shaft of a vortex as a coordinate axis, and firstly calculates the vortex shaft of the vortex
Figure 729144DEST_PATH_IMAGE001
Characteristic value of directionλ(q 1 )When the characteristic value is greater than 0, on the section vertical to the vortex shaft direction, the vortex flow form is spiral inward; when the characteristic value is less than 0, the vortex flow pattern is spiral outward on the cross section perpendicular to the vortex axis direction. The invention provides a novel method for judging the flow form of vortex, which is calculated based on a vortex shaft which is closer to a real curve and can effectively judge the flow characteristics of the vortex.

Description

Vortex flow form judging method suitable for curve vortex shaft and storage medium
Technical Field
The invention relates to the technical field of fluid mechanics, in particular to a method for judging a vortex motion form suitable for a curved vortex shaft and a storage medium.
Background
Vortices are very common phenomena in nature, such as tornadoes in nature, vortexes formed in water flow, and the like, and are tendons in movement of fluid media. Swirling motion is also common and important in the fluid-mechanics-related industries, for example, a high-speed vehicle can form a vortex around it, and if a rear-driving, lighter-weight vehicle becomes caught in it, a runaway can occur; in an aircraft, for example, vortex flow formed by fluid motion is important for generating and stably controlling the lift force of the aircraft. The research on the characteristics of the vortex motion is carried out, and the method is very beneficial to the development of a control means of the vortex motion.
In the field of fluid mechanics industry, wind tunnel/water tunnel experiments or numerical simulation means are generally adopted to develop flow property and flow control research. In the optimization of many problems, vortex flow control is an important part of the optimization. Therefore, accurate analysis of the swirl flow characteristics is helpful for researching relevant phenomena, mechanisms and rules of swirl motion, and finally an effective flow control method is developed.
In the past research, the vortex scroll is always considered to be a straight line, but in the practical problem, the vortex scroll is curved, so that the research result of researching the vortex scroll as a straight line has larger errors, and a judgment method closer to the real state of the vortex is urgently needed to be developed, so that the flow characteristic of the vortex is judged, and the accuracy of the analysis of the motion characteristic of the vortex in the experiment and numerical simulation is improved.
Disclosure of Invention
The invention aims to provide a method for judging the flow form of a vortex, which is suitable for a curved vortex shaft and can accurately judge the flow form of the vortex.
S10. calculating vortex shaft
Figure 797376DEST_PATH_IMAGE001
Characteristic value of directionλ(q 1 ),
Figure 172993DEST_PATH_IMAGE002
Wherein the content of the first and second substances,
Figure 270307DEST_PATH_IMAGE001
is a vortex shaft of vortex, and an orthogonal curve coordinate system is established based on the vortex shaft
Figure 946139DEST_PATH_IMAGE003
Figure 246539DEST_PATH_IMAGE004
Figure 843874DEST_PATH_IMAGE005
And
Figure 756598DEST_PATH_IMAGE006
are respectively an edge
Figure 286936DEST_PATH_IMAGE001
Figure 774549DEST_PATH_IMAGE007
And
Figure 639606DEST_PATH_IMAGE008
the speed of the direction of the beam is,
Figure 605288DEST_PATH_IMAGE009
Figure 740865DEST_PATH_IMAGE010
and
Figure 399380DEST_PATH_IMAGE011
is the Lame coefficient, M1Is that
Figure 768044DEST_PATH_IMAGE001
The local Mach number on the vortex shaft, p is a pressure value, and rho is the density of the fluid medium;
s20, judging the vortex flow form based on the characteristic value of the vortex shaft direction:
when in use
Figure 786685DEST_PATH_IMAGE012
In the direction of the vortex axis
Figure 26036DEST_PATH_IMAGE001
On the vertical section, the vortex flow shape is spiral inward;
when in use
Figure 855452DEST_PATH_IMAGE013
In the direction of the vortex axis
Figure 196565DEST_PATH_IMAGE001
On the vertical section, the vortex flow pattern is spiral outward.
Further, the Lame coefficient is as follows:
Figure 504050DEST_PATH_IMAGE014
wherein
Figure 129066DEST_PATH_IMAGE015
Are coordinate values of the orthonormal coordinate system.
Further, when the direction of the vortex axis is
Figure 378651DEST_PATH_IMAGE007
Or
Figure 190749DEST_PATH_IMAGE008
In the direction, the direction of the vortex axis is calculated in the following way
Figure 567504DEST_PATH_IMAGE007
Or
Figure 532180DEST_PATH_IMAGE008
Characteristic value of directionp(q 2 )Orp(q 3 )
Figure 703398DEST_PATH_IMAGE016
Figure 252060DEST_PATH_IMAGE017
Based onλ(q 2 )Orλ(q 3 )The vortex flow pattern is determined in step S20.
The present invention further includes a storage medium storing a computer program of the method for discriminating a vortex motion pattern of a curved vortex shaft as described above.
Compared with the prior art, the method for judging the vortex motion form of the curved vortex shaft has the following beneficial effects that:
the method for judging the vortex flow characteristics sets the vortex shaft of the vortex to be a curve, is closer to the real situation, and can effectively judge the flow characteristics of the vortex by calculating based on the real curve-shaped vortex shaft through simulation calculation, thereby providing a new method for experiments and numerical simulation calculation analysis.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments of the present invention or in the description of the prior art will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a method for determining a vortex motion pattern of a curved vortex shaft according to an embodiment of the present invention;
FIG. 2 is a schematic representation of an orthogonal curvilinear coordinate system of the vortex of the present invention;
FIG. 3 is a schematic view of a cavity calculation model according to an embodiment of the present invention;
FIG. 4 is a diagram illustrating the calculation results of an exemplary embodiment of the present invention;
FIG. 5 shows a cross-sectional spiral eddy current pattern at Y = -0.99D according to the present invention;
FIG. 6 shows a cross-sectional spiral eddy current pattern at Y = -0.9D according to the present invention;
fig. 7 shows a cross-sectional spiral vortex line pattern at Y = -0.8D according to the present invention.
Detailed Description
The following description provides many different embodiments, or examples, for implementing different features of the invention. The particular examples set forth below are illustrative only and are not intended to be limiting.
A method for judging the vortex motion mode of a curved vortex shaft is disclosed, as shown in figure 1, and comprises the following steps:
s10. calculating vortex shaft
Figure 166927DEST_PATH_IMAGE001
Characteristic value of directionλ(q 1 )
Figure 500956DEST_PATH_IMAGE002
Wherein the content of the first and second substances,
Figure 593808DEST_PATH_IMAGE001
is a vortex shaft of vortex, and an orthogonal curve coordinate system is established based on the vortex shaft
Figure 646078DEST_PATH_IMAGE003
As shown in fig. 2, in this example,
Figure 99056DEST_PATH_IMAGE004
Figure 536859DEST_PATH_IMAGE005
and
Figure 49880DEST_PATH_IMAGE006
are respectively an edge
Figure 589446DEST_PATH_IMAGE001
Figure 331268DEST_PATH_IMAGE007
And
Figure 639890DEST_PATH_IMAGE008
the speed of the direction of the beam is,
Figure 323812DEST_PATH_IMAGE009
Figure 68783DEST_PATH_IMAGE010
and
Figure 129143DEST_PATH_IMAGE011
is the Lame coefficient, M1Is that
Figure 292271DEST_PATH_IMAGE001
The local Mach number on the vortex shaft, p is a pressure value, and rho is the density of the fluid medium;
the calculation method of the Lamei coefficient comprises the following steps:
Figure 897827DEST_PATH_IMAGE014
wherein
Figure 146406DEST_PATH_IMAGE015
Is a standard orthogonal coordinate system
Figure 744877DEST_PATH_IMAGE018
The coordinate values of (2).
S20, judging the vortex flow form based on the characteristic value of the vortex shaft direction:
when the temperature is higher than the set temperature
Figure 746200DEST_PATH_IMAGE012
In the direction of the vortex axis
Figure 303083DEST_PATH_IMAGE001
On the vertical section, the vortex flow shape is spiral inward;
when in use
Figure 773379DEST_PATH_IMAGE013
In the direction of the vortex axis
Figure 406834DEST_PATH_IMAGE001
On the vertical section, the vortex flow pattern is spiral outward.
The above is q in an orthogonal curve coordinate system1The shaft is a vortex shaft of the vortex for calculation and judgment, and those skilled in the art can understand that in the actual operation process, the vortex shaft direction of the vortex can be always determined as q of the orthogonal curve coordinate system1The shaft is judged by adopting the method; of course, the same concept can be used to define the vortex axis of the vortex in q of the orthogonal curve coordinate system2Axis or q3Shaft:
when the vortex axis direction of the vortex is q of an orthogonal curve coordinate system2When the shaft is driven,
Figure 278975DEST_PATH_IMAGE016
wherein M is2Is that
Figure 990448DEST_PATH_IMAGE007
Local mach number on the scroll;
when in use
Figure 682461DEST_PATH_IMAGE019
In the direction q of the vortex axis2On the vertical section, the vortex flow shape is spiral inward;
when in use
Figure 153893DEST_PATH_IMAGE020
In the direction q of the vortex axis2On the vertical section, the vortex flow pattern is spiral outward.
When the vortex axis direction of the vortex is q of an orthogonal curve coordinate system3When the shaft is driven,
Figure 631273DEST_PATH_IMAGE021
wherein M is3Is that
Figure 998801DEST_PATH_IMAGE022
Local mach number on the scroll;
when in use
Figure 443688DEST_PATH_IMAGE023
In the direction q with the vortex axis3On the vertical section, the vortex flow form is spiral inward;
when in use
Figure 436921DEST_PATH_IMAGE024
In the cross section perpendicular to the scroll direction q3, the spiral flow pattern is spiral outward.
Example of calculation
The invention applies the above discrimination method to the cavity flow simulation calculation, the calculation model is shown in figure 3, the length L, the width W and the depth D of the cavity are respectively 200mm, 66.67mm and 33.33mm, and the depth is setDegree D, flow initiation conditions: the incoming flow Mach number is 0.9, and the Reynolds number is 5.68 multiplied by 105The gas temperature is 247K, the input condition of the inflow of fig. 3, in which the flow is developed, which in the initial condition brings about a movement of the gas flow inside the cavity.
Calculation results as shown in fig. 4, it can be seen that two vortices are generated at the bottom of the cavity rising from the bottom, and it can be seen that the vortex axes of the two vortices are curved.
The discriminant function of the method is used for analyzing the vortex state. Taking the right vortex in FIG. 4 as an example, respectively
Figure 18075DEST_PATH_IMAGE025
Figure 572816DEST_PATH_IMAGE026
And
Figure 239420DEST_PATH_IMAGE027
three horizontal planes are taken, as shown in fig. 5, 6 and 7. The calculation data is substituted into the calculation formula of the application to calculate and find
Figure 787076DEST_PATH_IMAGE025
A is>0, according to the judgment rule, the vortex spirals inwards and is consistent with the result in the figure 5; in that
Figure 3163DEST_PATH_IMAGE026
A is<0, according to the judgment rule, the vortex spirals outwards and is consistent with the result in the figure 6; in that
Figure 712493DEST_PATH_IMAGE027
A is<0, according to the judgment rule, the vortex spirals outwards, which is consistent with the result in fig. 7. The results show that the method can effectively judge the motion state of the vortex with the vortex shaft being a curve.
The present invention further provides a storage medium storing a computer program of the method for discriminating a vortex motion pattern of a curved vortex shaft, including but not limited to a floppy disk, an optical disk, a DVD, a hard disk, a flash Memory, a U disk, a CF card, an SD card, an MMC card, an SM card, a Memory Stick (Memory Stick), an XD card, etc.
The above description is intended to be illustrative of the preferred embodiment of the present invention and should not be taken as limiting the invention, but rather, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.

Claims (4)

1. A method for judging the vortex motion form of curved vortex shaft is characterized by that,
s10. calculating vortex shaft
Figure 568783DEST_PATH_IMAGE001
Characteristic value of directionλ(q 1 ),
Figure 951835DEST_PATH_IMAGE002
Wherein the content of the first and second substances,
Figure 200414DEST_PATH_IMAGE001
is a vortex shaft of vortex, and an orthogonal curve coordinate system is established based on the vortex shaft
Figure 798886DEST_PATH_IMAGE003
Figure 82099DEST_PATH_IMAGE004
Figure 170141DEST_PATH_IMAGE005
Are respectively an edge
Figure 640437DEST_PATH_IMAGE001
And
Figure 308178DEST_PATH_IMAGE006
the speed of the direction of the beam is,
Figure 445899DEST_PATH_IMAGE007
Figure 908104DEST_PATH_IMAGE008
and
Figure 928013DEST_PATH_IMAGE009
is the Lame coefficient, M1Is that
Figure 868287DEST_PATH_IMAGE001
The local Mach number on the vortex shaft, p is a pressure value, and rho is the density of the fluid medium;
s20, judging the vortex flow form based on the characteristic value of the vortex shaft direction:
when in use
Figure 860514DEST_PATH_IMAGE010
In the direction of the vortex axis
Figure 555937DEST_PATH_IMAGE001
On the vertical section, the vortex flow shape is spiral inward;
when in use
Figure 266404DEST_PATH_IMAGE011
In the direction of the vortex axis
Figure 744790DEST_PATH_IMAGE001
On the vertical section, the vortex flow pattern is spiral outward.
2. The method as claimed in claim 1, wherein the Lame coefficient is:
Figure 857103DEST_PATH_IMAGE012
wherein
Figure 661111DEST_PATH_IMAGE013
Are coordinate values of the orthonormal coordinate system.
3. The method according to claim 1 or 2, wherein the vortex flow pattern is determined when the vortex axis direction is set to be a direction of a curved vortex axis
Figure 655611DEST_PATH_IMAGE014
Or
Figure 937688DEST_PATH_IMAGE006
In the direction, the direction of the vortex axis is calculated in the following way
Figure 904507DEST_PATH_IMAGE014
Or
Figure 941733DEST_PATH_IMAGE006
Characteristic value of directionλ(q 2 )Orλ(q 3 )
Figure 358284DEST_PATH_IMAGE015
Wherein, M2Is that
Figure 444051DEST_PATH_IMAGE014
Local mach number on the scroll;
Figure 62114DEST_PATH_IMAGE016
is a rim
Figure 473504DEST_PATH_IMAGE014
The speed of the direction;
Figure 114701DEST_PATH_IMAGE017
wherein, M3Is that
Figure 4160DEST_PATH_IMAGE006
Local mach number on the scroll;
based onλ(q 2 )Orλ(q 3 )The vortex flow pattern is judged in step S20.
4. A storage medium storing a computer program for a method of discriminating a vortex motion pattern applied to a curved vortex shaft according to any one of claims 1 to 3.
CN202210373310.2A 2022-04-11 2022-04-11 Vortex flow form judging method suitable for curve vortex shaft and storage medium Active CN114444417B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210373310.2A CN114444417B (en) 2022-04-11 2022-04-11 Vortex flow form judging method suitable for curve vortex shaft and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210373310.2A CN114444417B (en) 2022-04-11 2022-04-11 Vortex flow form judging method suitable for curve vortex shaft and storage medium

Publications (2)

Publication Number Publication Date
CN114444417A CN114444417A (en) 2022-05-06
CN114444417B true CN114444417B (en) 2022-06-10

Family

ID=81360523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210373310.2A Active CN114444417B (en) 2022-04-11 2022-04-11 Vortex flow form judging method suitable for curve vortex shaft and storage medium

Country Status (1)

Country Link
CN (1) CN114444417B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012309A (en) * 2010-07-01 2011-04-13 北京航空航天大学 Vortex break point position judging method based on two-dimensional velocity field
DE102011112121A1 (en) * 2011-09-02 2013-03-07 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method for determining fluid flow of data, involves calculating spatial course of swirl axis of vortex flow in dependence on vortex flow induced flow measurement data and measuring position of vortex-flow-induced flow measurement data

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109859311B (en) * 2019-01-29 2021-05-18 河海大学 Simulation method based on Liutex-Omega vortex recognition theory

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012309A (en) * 2010-07-01 2011-04-13 北京航空航天大学 Vortex break point position judging method based on two-dimensional velocity field
DE102011112121A1 (en) * 2011-09-02 2013-03-07 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method for determining fluid flow of data, involves calculating spatial course of swirl axis of vortex flow in dependence on vortex flow induced flow measurement data and measuring position of vortex-flow-induced flow measurement data

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Optimization Calculation of Vortex Type Granulation Devices: Application of Software Products and Computer Modeling;Artem Artyukhov et al;《ceur-ws.org》;20201230;第1-12页 *
三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态;张涵信等;《空气动力学学报》;20120815;第第30卷卷(第04期);第421-430+436页 *
三维定常分离流和涡运动的定性分析研究;张涵信 等;《空气动力学学报》;19920401;第10卷(第1期);第8-20页 *
截面涡结构的描述方法及其在叶栅流场中的应用;王松涛等;《哈尔滨工业大学学报》;20020228;第第34卷卷(第01期);第128-131页 *
末级涡轮气动设计方法及非轴对称端壁造型技术研究;管继伟;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20070615(第6期);第C031-79页 *
超声三速角翼背风区旋涡运动的数值模拟;张树海,张涵信,朱国林;《空气动力学学报》;19961228;第第14卷卷(第04期);第430-435页 *
跨、超声速三角翼背风区旋涡运动的数值模拟;张树海等;《空气动力学学报》;19970328;第第15卷卷(第01期);第121-129页 *

Also Published As

Publication number Publication date
CN114444417A (en) 2022-05-06

Similar Documents

Publication Publication Date Title
Asnaghi et al. Large eddy simulations of cavitating tip vortex flows
Liu et al. Investigation of vortical structures and turbulence characteristics in corner separation in a linear compressor cascade using DDES
Gbadebo et al. Three-dimensional separations in axial compressors
Yin et al. Study on the air core formation of a gas–liquid separator
Zambonini et al. Corner separation dynamics in a linear compressor cascade
Juan et al. The impact of casing groove location on stall margin and tip clearance flow in a low-speed axial compressor
Karasu et al. Performance assessment of transition models for three-dimensional flow over NACA4412 wings at low Reynolds numbers
Sonoda et al. A study of advanced high-loaded transonic turbine airfoils
Li et al. Efficient design method for applying vortex generators in turbomachinery
Liu et al. Design of directional probes for high-frequency turbine measurements
Celic et al. Comparison of eddy-viscosity turbulence models in flows with adverse pressure gradient
Wang et al. Investigation into the phenomenon of flow deviation in the S-shaped discharge passage of a slanted axial-flow pumping system
Tucker et al. Hybrid les approach for practical turbomachinery flows—part ii: Further applications
Maynard et al. Unsteady structure of compressor tip leakage flows
Leishman et al. Effects of inlet ramp surfaces on the aerodynamic behavior of bleed hole and bleed slot off-take configurations
Kohri et al. Prediction of the performance of the engine cooling fan with CFD simulation
Bi et al. Numerical study of tip leakage vortex around a NACA0009 hydrofoil
Singh et al. Excitation of shear layer due to surface roughness near the leading edge: An experiment
CN114444417B (en) Vortex flow form judging method suitable for curve vortex shaft and storage medium
Talamelli et al. Linear instability characteristics of incompressible coaxial jets
Dawkins et al. The Unsteady Topology of Corner Separations
Burdett et al. Analysis of averaging methods for nonuniform total pressure fields
Ghosh et al. Topology optimization and experimental validation of an additively manufactured U-bend channel
Zhu et al. Experimental and numerical investigation of the flow field in the tip region of an axial ventilation fan
Asnaghi et al. Implicit large eddy simulation of tip vortex on an elliptical foil

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