CN114048552B - Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal - Google Patents

Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal Download PDF

Info

Publication number
CN114048552B
CN114048552B CN202210030530.5A CN202210030530A CN114048552B CN 114048552 B CN114048552 B CN 114048552B CN 202210030530 A CN202210030530 A CN 202210030530A CN 114048552 B CN114048552 B CN 114048552B
Authority
CN
China
Prior art keywords
mass flow
rotor
large water
water drops
calculating
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
CN202210030530.5A
Other languages
Chinese (zh)
Other versions
CN114048552A (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 CN202210030530.5A priority Critical patent/CN114048552B/en
Publication of CN114048552A publication Critical patent/CN114048552A/en
Application granted granted Critical
Publication of CN114048552B publication Critical patent/CN114048552B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Measuring Volume Flow (AREA)

Abstract

The invention discloses a rotor wing surface large water drop mass flow calculation method considering secondary impact and a terminal, belonging to the technical field of numerical simulation, wherein the method comprises the following steps: calculating the mass flow of the secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily
Figure 429399DEST_PATH_IMAGE001
(ii) a According to the mass flow of the secondary large water drops
Figure 252998DEST_PATH_IMAGE002
Combined with the mass flow of large water drops on the surface of the rotor under the condition of no large water drop splashing temp And mass loss rate of water drops on rotor surface under condition of large water drop splashingf m Calculating the final large drop mass on the rotor surfaceFlow rate. In the invention, in the calculation of the mass flow of the large water drops on the surface of the rotor wing, the final mass flow change of the large water drops caused by the secondary impact of the splashed supercooled large water drops on the surface of the rotor wing is considered, so that more accurate final mass flow calculation of the large water drops is realized, the numerical simulation of the icing of the rotor wing is closer to the physical reality, the numerical simulation precision of the icing of the rotor wing is improved, and the flight safety of an aircraft is ensured.

Description

Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
Technical Field
The invention relates to the technical field of numerical simulation, in particular to a rotor wing surface large water drop mass flow calculation method considering secondary impact and a terminal.
Background
The helicopter is widely applied to military operations, civil rescue and other aspects, and along with the wider field expansion of the use of the helicopter, the helicopter is required to have all-weather flight capability, and the probability of the helicopter encountering the dangerous condition of rotor icing is increased. The rotor freezes can change the rotor appearance, increase rotor weight, reduce rotor rotation rate, leads to rotor lift to descend, and the freezing that drops on the rotor simultaneously probably causes the striking harm to organism, weapon system, electronic equipment etc. increases flight uncertainty. Therefore, the research on rotor icing is of great significance to the improvement of flight safety and flight efficiency.
The range of the flight height of the helicopter is wide, freezing weather conditions such as freezing rain and frigid hair rain often occur, and the particle size of water drops in the air is Large and is called super-cooled Large water drops (SLD, the average particle size is larger than 50 μm). Compared with small water drops, supercooled large water drops widely exist in cloud mist, have poor pneumatic following performance and strong dynamic characteristics, cause more complex icing process and cause more serious damage, and increase the difficulty of SLD (scale-reduction potential) icing numerical value prediction and test evaluation.
The mass flow of water drops on the surface of the blade is a key intermediate quantity of icing calculation, and the accurate calculation of the mass flow directly influences the prediction accuracy of the rotor icing shape. The icing of the rotor wing is different from the icing of the fixed wing, and power needs to be given in the test, so that the test is limited too much and is more difficult, and numerical simulation becomes a primary means for predicting the mass flow of water drops on the surface of the rotor wing. In the rotor flow field, receive the rotatory influence of paddle, rotatory flow characteristic of washing down will appear in big water droplet, and big water droplet strikes in the paddle leading edge, can take place to splash in most times, leads to the reduction of paddle water droplet mass flow. At present, the study on the splashing characteristic of the large water drops on the surface of the rotor is less, and the processing mode of the large water drop mass flow is not clear, so that the study on the calculation method of the large water drop mass flow on the surface of the rotor needs to be further carried out.
Disclosure of Invention
The invention aims to solve the problem that the mass flow of large water drops on the surface of a rotor wing cannot be accurately calculated in the prior art, and provides a method and a terminal for calculating the mass flow of the large water drops on the surface of the rotor wing by considering secondary impact.
The purpose of the invention is realized by the following technical scheme: the rotor wing surface large water drop mass flow calculation method considering the secondary impact specifically comprises the following steps:
calculating the mass flow of the secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily
Figure 875797DEST_PATH_IMAGE001
According to the mass flow of the secondary large water drops
Figure 976608DEST_PATH_IMAGE002
Combined with the mass flow of large water drops on the surface of the rotor under the condition of no large water drop splashing temp And mass loss rate of water drops on rotor surface under condition of large water drop splashingf m Calculating the final large drop mass flow on the rotor surface
In one example, the calculation of the second time of the water drop on the rotor surface under the condition that the water drop impacts the rotor secondarilyMass flow rate
Figure 159328DEST_PATH_IMAGE002
The method comprises the following substeps:
calculating the mass flow total of secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily imp
Calculating the mass flow total of the primary large water drops entering the flow field in
Calculating the proportion of the mass flow total of the secondary large water drops to the mass flow total of the primary large water dropsλ
According to the ratio of lambda to the mass flow of large water drops temp Calculating the mass flow of the secondary large water drop
Figure 108041DEST_PATH_IMAGE003
In one example, the secondary large drop mass flow total imp The calculation formula of (2) is as follows:
Figure 237671DEST_PATH_IMAGE004
wherein (A), (B), (C), (D), (C), (B), (C)x,y,z) Grid points representing a calculated profile surface of the rotor; ω represents the set of grid points of the rotor's calculated profile surface.
In one example, the primary large drop mass flow total in The calculation formula of (2) is as follows:
Figure 442256DEST_PATH_IMAGE005
wherein,nrepresenting the blade rotational speed;rrepresenting the radius of the paddle disk;
Figure 61457DEST_PATH_IMAGE006
representing the angle between the water droplet and the impact surface.
In one example, the secondary large droplet mass flowMeasurement of
Figure 469435DEST_PATH_IMAGE007
The calculation formula of (2) is as follows:
Figure 199494DEST_PATH_IMAGE008
wherein (A), (B), (C), (D), (C), (B), (C)x,y,z) Representing grid points of the calculated profile surface of the rotor.
In one example, the final large drop mass flow rateThe calculation formula of (2) is as follows:
Figure 760050DEST_PATH_IMAGE009
wherein (A), (B), (C), (D), (C), (B), (C)x,y,z) Grid points representing a calculated profile surface of the rotor;λthe ratio of the total mass flow of the secondary large water drops to the total mass flow of the primary large water drops is shown.
In one example, the method further comprises the steps of:
calculating the mass flow of large water drops on the surface of the rotor under the condition of no splashing temp
Calculating the mass loss rate of water drops on the surface of the rotor wing under the splashing conditionf m
In one example, the calculation calculates the mass flow of large water droplets on the rotor surface without splashing temp The method comprises the following substeps:
generating rotor calculation profilesΓGrid ofΩ
Computational gridΩThe air flow field information of (a);
calculating water drop flow field information according to the air flow field information;
according to a gridΩAnd calculating the mass flow of the big water drop according to the water drop flow field information temp
In one example, the calculation calculates a mass loss rate of water droplets on the surface of the rotor under splash conditionsf m IncludedThe following substeps:
calculating impact parameters of water drops impacting rotor surfaceK
Based on the minimum mass loss rate, according to the impact parametersKCalculating the mass loss rate according to the water drop flow field informationf m
It should be further noted that the technical features corresponding to the above examples can be combined with each other or replaced to form a new technical solution.
The present invention also includes a storage medium having stored thereon computer instructions operable to perform the steps of the method for calculating a mass flow of large water droplets on a rotor surface that accounts for secondary impacts formed in accordance with any one or more of the above-described example compositions.
The invention also includes a terminal comprising a memory having stored thereon computer instructions executable on the processor, and a processor that when executed performs the steps of the method of calculating a mass flow of large water droplets from a rotor surface that accounts for secondary impacts formed by any one or more of the example compositions.
Compared with the prior art, the invention has the beneficial effects that:
in the invention, in the calculation of the mass flow of the large water drops on the surface of the rotor wing, the final mass flow change of the large water drops caused by the secondary impact of the splashed supercooled large water drops on the surface of the rotor wing is considered, so that more accurate final mass flow calculation of the large water drops is realized, the numerical simulation of the icing of the rotor wing is closer to the physical reality, the numerical simulation precision of the icing of the rotor wing is improved, and the flight safety of an aircraft is ensured.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the invention without limiting the invention.
FIG. 1 is a flow chart of a method in an example of the invention;
fig. 2 is a flow chart of a method of a preferred embodiment of the present invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it should be noted that directions or positional relationships indicated by "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", and the like are directions or positional relationships described based on the drawings, and are only for convenience of description and simplification of description, and do not indicate or imply that the device or element referred to must have a specific orientation, be configured and operated in a specific orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first" and "second" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the description of the present invention, it should be noted that, unless otherwise explicitly stated or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
In addition, the technical features involved in the different embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
The invention aims to solve the problems that the mass flow change of large water drops caused by the secondary impact of water drops on the surface of a rotor wing is not considered in the prior art, so that the final mass flow calculation accuracy of the large water drops on the surface of the rotor wing is not high, the icing numerical simulation precision of the rotor wing is low, and accurate theoretical support cannot be provided for the safe flight of an aircraft.
In one example, a method for calculating the mass flow of large water drops on the surface of a rotor by considering secondary impact, as shown in fig. 1, specifically comprises the following steps:
s1: calculating the mass flow of the secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily
Figure 550152DEST_PATH_IMAGE010
(ii) a Wherein the mass flow rate of large water droplets represents the mass of large water droplets passing through the rotor surface per unit area per unit time. In the present application, the large water droplets (water droplets) are specifically supercooled large water droplets SLD having a particle size of more than 50 μm. Further, the water droplet secondary impact rotor specifically is: after water drops collide with the surface of the rotor wing and splash, water drops on the upper wing surface of the rotor wing blade enter a rotary flow field and collide with the rotor blade again, secondary water drops are collected on the rotor wing blade, and the final mass flow of the large water drops on the surface of the rotor wing can be influenced by the mass flow (secondary mass flow of the large water drops) of the large water drops entering the flow field
S2: according to the mass flow of the secondary large water drops
Figure 711006DEST_PATH_IMAGE007
Combined with the mass flow of large water drops on the surface of the rotor under the condition of no large water drop splashing temp And mass loss rate of water drops on rotor surface under condition of large water drop splashingf m Calculating the final large drop mass flow on the rotor surface. The mass loss refers to a mass loss caused by splashing of large supercooled water droplets, and correspondingly, the mass loss rate represents a ratio of a water droplet loss amount to a collection amount on the surface of the rotor under a splashing condition.
This example has considered the final big water droplet mass flow change that splashes water droplet secondary striking rotor surface and arouse in the big water droplet mass flow on rotor surface calculates to this realizes more accurate final big water droplet mass flow and calculates, makes the frozen numerical simulation of rotor more press close to the physics reality, has improved the frozen numerical simulation precision of rotor, guarantees the flight safety of aircraft with this.
In one example, the mass flow of the second largest water droplet on the rotor surface is calculated under the condition that the water droplet impacts the rotor secondarily
Figure 979176DEST_PATH_IMAGE002
The method comprises the following substeps:
s11: calculating the total mass flow of the primary large water drops on the surface of the rotor imp At the moment, no large water drops splash on the surface of the rotor wing;
s12: calculating the mass flow total of the primary large water drops entering the flow field in (ii) a Wherein the total mass flow of the primary large water drops in Namely the mass flow of the water drops entering the flow field for the first time.
S13: calculating the proportion of the mass flow total of the secondary large water drops to the mass flow total of the primary large water dropsλ
S14: according to the ratio of lambda to the mass flow of large water drops temp Calculating the mass flow of the secondary large water drop
Figure 158354DEST_PATH_IMAGE003
In one example, the secondary large drop mass flow total imp The calculation formula of (2) is as follows:
Figure 384936DEST_PATH_IMAGE011
wherein (A), (B), (C), (D), (C), (B), (C)x,y,z) Representing calculated profile of rotorΓA grid point of the surface; omega denotes the calculated profile of the rotorΓA set of grid points of the surface.
In one example, the total mass flow of one large drop in The calculation formula of (2) is as follows:
Figure 501927DEST_PATH_IMAGE012
wherein,nrepresenting the blade rotational speed;rrepresenting the radius of the paddle disk;
Figure 573789DEST_PATH_IMAGE006
representing the angle between the water droplet and the impact surface.
In one example, the ratioλThe calculation formula of (2) is as follows:
Figure 374517DEST_PATH_IMAGE013
in one example, the secondary large drop mass flow rate
Figure 240842DEST_PATH_IMAGE007
The calculation formula of (2) is as follows:
Figure 376288DEST_PATH_IMAGE014
the mass flow of water droplets is now considered to be comparable to the collection efficiency at the first impact.
In one example, the final large drop mass flowThe calculation formula of (2) is as follows:
Figure 986261DEST_PATH_IMAGE015
the final large drop mass flow rate in this exampleWith large water drop mass flow temp Mass loss ratef m Mass flow of secondary large water drop
Figure 140030DEST_PATH_IMAGE016
Correlation, i.e. final large drop mass flowThe mass flow of the large water drops splashed out of the surface of the rotor wing is deducted by the mass flow of the large water drops without considering the splashing condition, and then the mass flow of the secondary large water drops is added.
In an example, the method further comprises the steps of:
s01: calculating the mass flow of large water drops on the surface of the rotor under the condition of no splashing temp (ii) a Wherein, the big water droplet of supercooling can not splash on the surface of the rotor wing after impacting without splashing.
S02: calculating the mass loss rate of water drops on the surface of the rotor wing under the splashing conditionf m . Wherein splashing is the splashing that occurs after a water droplet hits the rotor, leading to a certain water droplet quality loss.
In one example, the mass flow of large water droplets on the rotor surface without splashing is calculated temp The method comprises the following substeps:
s011: generating rotor calculation profilesΓGrid ofΩ(ii) a Wherein the rotor calculates the profileΓI.e. vector models, meshes of rotorsΩAnd the vector model is a basic calculation composition unit containing the vector model and the calculation domain. In particular, a gridΩThe generation specifically comprises: using mesh generation software to compute the shapeΓAs inputs, the computational topology and the grid Ω are generated. The grid generation software includes Gridgen, Pointwise, GridStar, etc., and the GridStar grid generation software used in this embodiment generates the calculation shapeΓGrid omega.
S012: computational gridΩThe air flow field information of (a); specifically, air flow field calculation software or program is employed to calculate the profileΓGrid ofΩAnd setting information such as an air flow field calculation method, boundary conditions, calculation conditions and the like as input, and calculating to obtain air flow field information P corresponding to the grids.
S013: calculating water drop flow field information according to the air flow field information; specifically, water drop flow field calculation software or program is adopted, the air flow field P is used as input, information such as a water drop flow field calculation method, boundary conditions and calculation conditions is set, and the shape is calculatedΓThe surface of the grid is set as a wall suction boundary condition, and finally, a grid is obtained through calculationΩRecording the corresponding water drop flow field information as W; more specifically, the information on the water drop flow field includes the average particle size of the water drop and waterDrop density, normal velocity at the drop airfoil, drop surface tension coefficient, drop dynamic viscosity, drop incidence frequency, calculated profileΓWater content, water droplet velocity, incoming flow water content, far field water droplet velocity, etc.
S014: according to a gridΩAnd calculating the mass flow of the big water drop according to the water drop flow field information temp . In particular, large drop mass flow temp Can be calculated by calculating the profileΓSurface grid points (x,y,z) Water droplet volume fraction ofαAnd velocityvThe calculation is obtained, and the specific calculation formula is as follows:
Figure 177256DEST_PATH_IMAGE017
wherein,ρwhich is the density of the water droplets,
Figure 799999DEST_PATH_IMAGE018
is in the shape ofΓNormal to the surface.
In one example, the mass loss rate of water droplets on the rotor surface under splash conditions is calculatedf m The method comprises the following substeps:
s021: calculating impact parameters of water drops impacting rotor surfaceK(ii) a Specifically, given the spatter criterion, the impact parameterKThe calculation formula of (2) is as follows:
Figure 213663DEST_PATH_IMAGE019
wherein,ρrepresents the water drop density;drepresents an average particle diameter of water droplets;vrepresents the normal velocity at the water droplet airfoil;σrepresents the surface tension coefficient of water drops;μrepresents the dynamic viscosity of water drops;θrepresenting the angle between the water droplet and the collision surface;
Figure 474136DEST_PATH_IMAGE020
which represents the frequency of incidence of the water droplets,
Figure 151105DEST_PATH_IMAGE021
s021: based on the minimum mass loss rate, according to the impact parametersKCalculating the mass loss rate according to the water drop flow field informationf m . In this example, the minimum mass loss rate is preferably 0.2, and the calculation formula of the mass loss rate is:
Figure DEST_PATH_IMAGE022
further, the minimum mass loss rate is a constant and may also be a function of the incoming flow conditions with respect to average particle size, angle of incidence, impact velocity, etc., in this case:
a=a (d-a 2 ) 2 +a 3
wherein,arepresents the minimum mass loss rate;a 1 ,a 2 ,a 3 coefficients representing positive real numbers;dthe average particle size of water droplets is shown. As a preference, it is preferred that,a 1 =9.92×10-6a 2 =50;a 3 =0.12。
the above examples are combined to obtain the preferred example of the present invention, as shown in fig. 2, which specifically includes the following steps:
s1': calculating the mass flow of large water drops on the surface of the rotor under the condition of no splashing temp
S2': calculating the mass loss rate of water drops on the surface of the rotor wing under the splashing conditionf m
S3': calculating the mass flow of the secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily
Figure 916936DEST_PATH_IMAGE010
S4': according to the mass flow of the secondary large water drops
Figure 291548DEST_PATH_IMAGE007
Combined with the mass flow of large water drops on the surface of the rotor under the condition of no large water drop splashing temp And mass loss rate of water drops on rotor surface under condition of large water drop splashingf m Calculating the final large drop mass flow on the rotor surface
The invention provides a rotor wing surface large water drop mass flow calculation method considering secondary impact under a large water drop splashing condition. After a rotor wing air flow field and a rotor wing water drop flow field are solved, water drop impact and collection information on the surface of a blade is obtained, then the water drop mass loss rate caused by splashing is calculated based on a splashing judgment criterion, a secondary impact and collection amount calculation method of the splashed water drops is further provided, and finally a rotor wing surface water drop mass flow calculation result is obtained. Through introducing the thought of secondary striking to make the computational process of rotor big water droplet collection rate more close to physics reality, the computational result is more accurate.
The present embodiment provides a storage medium, which has the same inventive concept as the rotor surface large water droplet mass flow rate calculation method considering secondary impact formed by any one or a combination of the above examples, and on which computer instructions are stored, and the computer instructions are executed when executed to perform the steps of the rotor surface large water droplet mass flow rate calculation method considering secondary impact formed by any one or a combination of the above examples.
Based on such understanding, the technical solution of the present embodiment or parts of the technical solution may be essentially implemented in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The present embodiment also provides a terminal, which has the same inventive concept as the rotor surface large water droplet mass flow calculation method considering secondary impact formed by any one or more of the above examples in combination, and includes a memory and a processor, where the memory stores computer instructions executable on the processor, and the processor executes the computer instructions to perform the steps of the rotor surface large water droplet mass flow calculation method considering secondary impact formed by any one or more of the above examples in combination. The processor may be a single or multi-core central processing unit or a specific integrated circuit, or one or more integrated circuits configured to implement the present invention.
Each functional unit in the embodiments provided by the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The above detailed description is for the purpose of describing the invention in detail, and it should not be construed that the detailed description is limited to the description, and it will be apparent to those skilled in the art that various modifications and substitutions can be made without departing from the spirit of the invention.

Claims (8)

1. The rotor wing surface large water drop mass flow calculation method considering the secondary impact is characterized by comprising the following steps of: the method comprises the following steps:
calculating the mass flow of the secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily
Figure 790966DEST_PATH_IMAGE001
The mass flow of the secondary large water drops
Figure 500296DEST_PATH_IMAGE001
The calculation formula of (2) is as follows:
Figure 559257DEST_PATH_IMAGE002
according to the mass flow of the secondary large water drops
Figure 848287DEST_PATH_IMAGE001
Combined with the mass flow of large water drops on the surface of the rotor under the condition of no large water drop splashing temp And mass loss rate of water drops on rotor surface under condition of large water drop splashingf m Calculating the final large drop mass flow on the rotor surface
The final mass flow of large water dropletsThe calculation formula of (2) is as follows:
Figure 731929DEST_PATH_IMAGE003
wherein (A), (B), (C), (D), (C), (B), (C)x,y,z) Grid points representing a calculated profile surface of the rotor;λthe ratio of the total mass flow of the secondary large water drops to the total mass flow of the primary large water drops is shown.
2. The rotor surface large water droplet mass flow calculation method considering secondary impact according to claim 1, characterized in that: calculate the secondary mass flow of big water droplet on rotor surface under the water droplet secondary striking rotor condition
Figure 674478DEST_PATH_IMAGE001
The method comprises the following substeps:
calculating the mass flow total of secondary large water drops on the surface of the rotor under the condition that the water drops impact the rotor secondarily imp
Calculating the mass flow total of the primary large water drops entering the flow field in
Calculating secondary mass flow of large water dropsRatio of total amount to total mass flow of large water dropsλ
According to the ratio of lambda to the mass flow of large water drops temp Calculating the mass flow of the secondary large water drop
Figure 112412DEST_PATH_IMAGE001
3. The rotor surface large water drop mass flow calculation method considering secondary impact according to claim 2, characterized in that: the total mass flow of the secondary large water drops imp The calculation formula of (2) is as follows:
Figure 736292DEST_PATH_IMAGE004
where ω represents the set of grid points of the rotor's calculated profile surface.
4. The rotor surface large water drop mass flow calculation method considering secondary impact according to claim 2, characterized in that: the total mass flow of the primary large water drops in The calculation formula of (2) is as follows:
Figure 208861DEST_PATH_IMAGE005
wherein,nrepresenting the blade rotational speed;rrepresenting the radius of the paddle disk;
Figure 994415DEST_PATH_IMAGE006
representing the angle between the water droplet and the impact surface.
5. The rotor surface large water droplet mass flow calculation method considering secondary impact according to claim 1, characterized in that: the method further comprises the steps of:
calculating the strip without large water drop splashLarge water drop mass flow on the surface of the lower rotor temp
Calculating the mass loss rate of water drops on the surface of the rotor wing under the splashing conditionf m
6. The rotor surface large water droplet mass flow calculation method considering secondary impact according to claim 5, characterized in that: calculating the mass flow of the large water drops on the surface of the rotor wing under the condition that the large water drops do not splash temp The method comprises the following substeps:
generating rotor calculation profilesΓGrid ofΩ
Computational gridΩThe air flow field information of (a);
calculating water drop flow field information according to the air flow field information;
according to a gridΩAnd calculating the mass flow of the big water drop according to the water drop flow field information temp
7. The rotor surface large water droplet mass flow calculation method considering secondary impact according to claim 6, characterized in that: calculating the mass loss rate of water drops on the surface of the rotor wing under the splashing conditionf m The method comprises the following substeps:
calculating impact parameters of water drops impacting rotor surfaceK
Based on the minimum mass loss rate, according to the impact parametersKCalculating the mass loss rate according to the water drop flow field informationf m
8. A terminal comprising a memory and a processor, the memory having stored thereon computer instructions executable on the processor, the terminal comprising: the processor, when executing the computer instructions, performs the steps of the method of calculating the mass flow of large water droplets on the surface of a rotor taking into account secondary impacts as claimed in any one of claims 1 to 7.
CN202210030530.5A 2022-01-12 2022-01-12 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal Active CN114048552B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210030530.5A CN114048552B (en) 2022-01-12 2022-01-12 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210030530.5A CN114048552B (en) 2022-01-12 2022-01-12 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal

Publications (2)

Publication Number Publication Date
CN114048552A CN114048552A (en) 2022-02-15
CN114048552B true CN114048552B (en) 2022-03-29

Family

ID=80196297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210030530.5A Active CN114048552B (en) 2022-01-12 2022-01-12 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal

Country Status (1)

Country Link
CN (1) CN114048552B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2807522A1 (en) * 2000-04-07 2001-10-12 Aerospatiale Matra Airbus DEVICE FOR DETERMINING THE VALUES OF AT LEAST ONE PARAMETER OF PARTICLES, IN PARTICULAR WATER DROPLETS
CN102682144A (en) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 Flight icing numerical value simulation method of helicopter rotor wing
CA2775292A1 (en) * 2011-04-29 2012-10-29 Rosemount Aerospace Inc. Apparatus and method for detecting aircraft icing conditions
CN206876374U (en) * 2017-06-07 2018-01-12 中国空气动力研究与发展中心低速空气动力研究所 A kind of super-cooling waterdrop ice detection
CN112345196A (en) * 2021-01-11 2021-02-09 中国空气动力研究与发展中心低速空气动力研究所 Super-cooled large water drop splash simulation test device and method
CN112906140A (en) * 2021-04-28 2021-06-04 中国空气动力研究与发展中心计算空气动力研究所 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate
CN113255067A (en) * 2021-07-16 2021-08-13 中国空气动力研究与发展中心计算空气动力研究所 Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104268399B (en) * 2014-09-24 2017-03-15 空气动力学国家重点实验室 The computational methods of model parameter in icing wind tunnel test under the conditions of big water droplet are subcooled
CN112985347B (en) * 2021-05-11 2021-08-03 中国空气动力研究与发展中心低速空气动力研究所 Method for calculating roughness of icing surface of airplane

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2807522A1 (en) * 2000-04-07 2001-10-12 Aerospatiale Matra Airbus DEVICE FOR DETERMINING THE VALUES OF AT LEAST ONE PARAMETER OF PARTICLES, IN PARTICULAR WATER DROPLETS
CA2775292A1 (en) * 2011-04-29 2012-10-29 Rosemount Aerospace Inc. Apparatus and method for detecting aircraft icing conditions
CN102682144A (en) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 Flight icing numerical value simulation method of helicopter rotor wing
CN206876374U (en) * 2017-06-07 2018-01-12 中国空气动力研究与发展中心低速空气动力研究所 A kind of super-cooling waterdrop ice detection
CN112345196A (en) * 2021-01-11 2021-02-09 中国空气动力研究与发展中心低速空气动力研究所 Super-cooled large water drop splash simulation test device and method
CN112906140A (en) * 2021-04-28 2021-06-04 中国空气动力研究与发展中心计算空气动力研究所 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate
CN113255067A (en) * 2021-07-16 2021-08-13 中国空气动力研究与发展中心计算空气动力研究所 Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
Effect of surface wettability on impact-freezing of supercooled large water droplet;WeiliangKong 等;《Experimental Thermal and Fluid Science》;20220101;第130卷;1-15 *
Effects of surface roughness on splashing characteristics of large droplets with digital inline holographic imaging;QiaoWang 等;《Cold Regions Science and Technology》;20211130;第191卷;1-14 *
关于飞机结冰的水滴撞击特性计算与结冰相似准则研究;王梓旭;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20091015(第10期);C031-14 *
变形破碎特性对SLD结冰过程影响;桑为民 等;《航空动力学报》;20170731;第32卷(第7期);1537-1544 *
结冰条件下过冷水滴撞击特性及热平衡分析;闵现花;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20101015(第10期);C031-10 *
过冷大水滴变形与破碎的影响因素;李维浩 等;《航空学报》;20181225;第39卷(第12期);1-8 *
过冷大水滴条件下机翼结冰数值仿真;吴佩佩 等;《计算机仿真》;20140930;第31卷(第9期);51-55 *
飞行防冰与结冰中的流动传热计算方法研究;朱剑鋆;《中国优秀博硕士学位论文全文数据库(博士)》;20200315(第03期);A009-6 *

Also Published As

Publication number Publication date
CN114048552A (en) 2022-02-15

Similar Documents

Publication Publication Date Title
CN113255067B (en) Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal
US20040155151A1 (en) Morphogenetic modelling of in-flight icing
Cao et al. Numerical Simulation of Supercooled Large Droplet Icing Phenomenon: A Review: Y. Cao, M. Xin
CN112906140B (en) Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate
Yeong et al. Aerodynamic optimization of micro aerial vehicle
CN114386273B (en) Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal
CN111396269A (en) Multi-time-step unsteady icing calculation method and system and storage medium
CN111388908A (en) Aviation fire extinguishing parameter acquisition method and system
EP2215572B1 (en) Computer-aided method for predicting particles capture by a surface of a moving object
CN115526078A (en) A Calculation Method of Water Droplet Freezing Coefficient on Wind Turbine Blades
CN114048552B (en) Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
WO2019186151A1 (en) Methods and apparatus for simulating liquid collection on aerodynamic components
GUFFOND et al. Overview of icing research at ONERA
CN114863039B (en) A kind of water droplet initial release area determination method and computer readable storage medium
Rutkowski et al. Numerical study of droplet splashing and re-impingement
CN117010204A (en) Wing icing calculation method and system, storage medium, and electronic equipment
CN115828421A (en) Helicopter noise early warning method
CN111079296B (en) Aircraft component and aircraft flight load assessment method
CN107577851B (en) A kind of buffeting special efficacy modeling method based on acceleration power spectral density
Wright et al. Computational simulation of large droplet icing
CN119378443A (en) A wind turbine blade icing risk assessment and early warning method and system
CN115098825A (en) Water drop release method and collection rate acquisition method for on-way encryption and medium
Hann Numerical Simulation of In-Flight Icing of Unmanned Aerial Vehicles
Lennie et al. Development of ice throw model for wind turbine simulation software QBlade
CN114036727A (en) Terminal guidance-based fast planning method for glide trajectory

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