CN114386273A - Rotor surface large water drop collection rate calculation method considering secondary impact and terminal - Google Patents
Rotor surface large water drop collection rate calculation method considering secondary impact and terminal Download PDFInfo
- Publication number
- CN114386273A CN114386273A CN202210030651.XA CN202210030651A CN114386273A CN 114386273 A CN114386273 A CN 114386273A CN 202210030651 A CN202210030651 A CN 202210030651A CN 114386273 A CN114386273 A CN 114386273A
- Authority
- CN
- China
- Prior art keywords
- rotor
- large water
- collection rate
- temp
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a rotor wing surface large water drop collection rate calculation method considering secondary impact and a terminal, belonging to the technical field of numerical simulation, wherein the method comprises the following steps: calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact on rotor wing'temp(ii) a According to the secondary large water droplet collection rate beta'tempCombined with the collection rate beta of large water drops on the surface of the rotor under the condition that the large water drops do not splashtempAnd the mass loss rate f of water drops on the surface of the rotor under the condition of large water drop splashingmThe final large water droplet collection rate beta of the rotor surface was calculated. In the calculation of the collection rate of the large water drops on the surface of the rotor wing, the final large water drop collection rate change caused by secondary impact of the splashed supercooled large water drops on the surface of the rotor wing is considered, so that a more accurate final large water drop collection rate meter is realizedAnd the numerical simulation of the rotor icing is closer to the physical reality, and the numerical simulation precision of the rotor icing is improved, so that the flight safety of the aircraft is ensured.
Description
Technical Field
The invention relates to the technical field of numerical simulation, in particular to a rotor wing surface large water drop collection rate 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 large water drop collection rate is a key intermediate quantity of icing calculation, the accurate calculation of the large water drop collection rate directly influences the results of ice shape prediction and aerodynamic influence, and the accurate calculation method for the rotor wing large water drop collection rate can be established to improve the accuracy of icing numerical simulation. 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 the first means for predicting the collection rate of large water drops on the surface of the rotor wing. In the rotor flow field, receive centrifugal force to influence, rotatory flow characteristic of washing down will appear in big water droplet, and big water droplet strikes in fast rotating's paddle leading edge, can take place to splash in most times, leads to the reduction of paddle water droplet collection volume. At present, the study on the splashing characteristics of the large water drops on the surface of the rotor is less, and the processing mode of the large water drop collection rate is not clear, so that the study on the calculation method of the large water drop collection rate on the surface of the rotor needs to be further carried out.
Disclosure of Invention
The invention aims to solve the problem that the collection rate of large water drops on the surface of a rotor cannot be accurately calculated in the prior art, and provides a calculation method and a terminal for the collection rate of the large water drops on the surface of the rotor, which take secondary impact into consideration.
The purpose of the invention is realized by the following technical scheme: the method for calculating the collection rate of the large water drops on the surface of the rotor wing considering the secondary impact specifically comprises the following steps:
calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact on rotor wing'temp；
According to the secondary large water droplet collection rate beta'tempCombined with the collection rate beta of large water drops on the surface of the rotor under the condition that the large water drops do not splashtempAnd the mass loss rate f of water drops on the surface of the rotor under the condition of large water drop splashingmThe final large water droplet collection rate beta of the rotor surface was calculated.
In one example, the calculation of the secondary large water droplet collection rate β 'of the rotor surface under the condition of large water droplet secondary impact on the rotor is performed'tempThe method comprises the following substeps:
calculating the total quantity beta of the secondary large water drop collection rate on the surface of the rotor under the condition that the large water drops impact the rotor secondarilyimp；
Calculating the ratio lambda of the total collection rate of the secondary water drops to the total collection rate of the primary large water drops;
according to the proportion of lambda and the collection rate of large water drops betatempCalculating secondary big water drop collection rate beta'temp。
In one example, the total secondary large droplet collection rate βimpIs calculated byThe formula is as follows:
wherein (x, y, z) represents grid points of the calculated contoured surface of the rotor; ω represents the set of grid points of the rotor's calculated profile surface.
In one example, the ratio λ is calculated as:
wherein (x, y, z) represents grid points of the rotor's computational profile surface;
the efficiency of hitting a water droplet on the surface of the rotor is shown as 1, and the total sum is carried out;
a set of grid points representing the secondary impact of a large water droplet on the calculated profile surface of the rotor.
In one example, the secondary large water droplet collection rate β'tempThe calculation formula of (2) is as follows:
β′temp(x,y,z)＝λ·βtemp(x,y,z)
where (x, y, z) represents grid points of the calculated contoured surface of the rotor.
In one example, the final large water droplet collection rate β is calculated by the formula:
β(x,y,z)＝βtemp(x,y,z)-fm(x,y,z)·βtemp(x,y,z)+λ·βtemp(x,y,z)＝(1-fm(x,y,z)+λ)·βtemp(x,y,z)
wherein (x, y, z) represents grid points of the calculated contoured surface of the rotor; λ represents the ratio of the total amount of the second water droplet collection to the total amount of the first large droplet collection rate.
In one example, the method further comprises the steps of:
calculating the collection rate beta of large water drops on the surface of the rotor under the condition of no splashingtemp；
Calculating the mass loss rate f of water drops on the surface of the rotor wing under the splashing conditionm。
In one example, the calculation calculates a large water droplet collection rate β for the rotor surface without spattertempThe method comprises the following substeps:
generating a grid omega of a rotor wing calculation shape gamma;
calculating the information of the air flow field of the grid omega;
calculating water drop flow field information according to the air flow field information;
calculating the collection rate beta of the big water drops according to the grid omega and the water drop flow field informationtemp。
In one example, the calculation calculates the mass loss rate f of water droplets on the surface of the rotor under splashing conditionsmThe method comprises the following substeps:
calculating an impact parameter K of water drops impacting the surface of the rotor wing;
calculating the mass loss rate f according to the impact parameter K and the water drop flow field information on the basis of considering the minimum mass loss ratem。
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 collection rate of large water droplets on a rotor surface considering secondary impacts formed from any one or more of the compositions of the examples described above.
The invention also includes a terminal comprising a memory and a processor, the memory having stored thereon computer instructions executable on the processor, the processor executing the computer instructions to perform the steps of the method for calculating a collection rate of large water droplets on a rotor surface considering secondary impacts formed from any one or more of the example compositions.
Compared with the prior art, the invention has the beneficial effects that:
in the calculation of the collection rate of the large water drops on the surface of the rotor wing, the final large water drop collection rate change caused by secondary impact of the splashed supercooled large water drops on the surface of the rotor wing is considered, so that the more accurate final large water drop collection rate calculation 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 change of the collection rate of large water drops caused by the secondary impact of the large water drops on the surface of a rotor wing is not considered in the prior art, so that the calculation accuracy of the final collection rate of the large water drops on the surface of the rotor wing is not high, the numerical simulation precision of the icing 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 a collection rate of large water drops on a rotor surface considering secondary impact, as shown in fig. 1, specifically includes the following steps:
s1: calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact on rotor wing'temp(ii) a The large water drop collection rate, namely the local large water drop collection coefficient, represents the ratio of the actually collected water amount on a certain local area of an object plane, namely the rotor surface of the present application, to the maximum value of the water amount possibly collected on the rotor surface, and represents the water drop impact range on the rotor surface and the water amount distribution in the impact area, and is the most important water drop impact characteristic parameter. 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 rotor is specifically hit to big water droplet secondary: after the water droplet strikes the rotor surface and splashes, the water droplet on the upper wing surface of the rotor blade enters the rotating flow field and collides with the rotor blade again, so that the secondary water droplet collection phenomenon occurs on the rotor blade, and the final large water droplet collection rate beta on the rotor surface can be influenced by the secondary large water droplet collection rate.
S2: according to the secondary large water droplet collection rate beta'tempCombined with the collection rate beta of large water drops on the surface of the rotor under the condition that the large water drops do not splashtempAnd the mass loss rate f of water drops on the surface of the rotor under the condition of large water drop splashingmThe final large water droplet collection rate beta of the rotor surface was calculated. The mass loss represents 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 collection rate change that splashes big water droplet secondary striking rotor surface and arouse in the big water droplet collection rate calculation on rotor surface to this realizes more accurate final big water droplet collection rate calculation, makes the numerical simulation that the rotor freezes more press close to the physics reality, has improved the numerical simulation precision that the rotor freezes, with this flight safety of guaranteeing the aircraft.
In one example, a secondary large water droplet collection rate β 'of the rotor surface under large water droplet secondary impact rotor conditions is calculated'tempThe method comprises the following substeps:
s11: calculating the total quantity beta of the secondary large water drop collection rate on the surface of the rotor under the condition that the large water drops impact the rotor secondarilyimp；
S12: calculating the ratio lambda of the total collection rate of the secondary water drops to the total collection rate of the primary large water drops;
s13: according to the proportion of lambda and the collection rate of large water drops betatempCalculating secondary big water drop collection rate beta'temp。
In one example, the total secondary large droplet collection rate βimpThe calculation formula of (2) is as follows:
wherein (x, y, z) represent grid points of the calculated profile Γ surface of the rotor; ω represents the set of grid points of the rotor's calculated profile Γ surface.
In one example, the ratio λ is calculated as:
wherein the content of the first and second substances,
representing the set of grid points where the calculated profile Γ surface of the rotor is impacted twice by a large water drop.
In an example, secondary large droplet collection rate β'tempThe calculation formula of (2) is as follows:
β′temp(x,y,z)＝λ·βtemp(x,y,z)
the water droplet collection efficiency is now considered to be comparable to the first impact collection efficiency.
In one example, the final large droplet collection rate β is calculated as:
β(x,y,z)＝βtemp(x,y,z)-fm(x,y,z)·βtemp(x,y,z)+λ·βtemp(x,y,z)＝(1-fm(x,y,z)+λ)·βtemp(x,y,z)
in the present example, the final large water droplet collection rate beta and the large water droplet collection rate betatempMass loss rate fmAnd secondary large water droplet collection rate beta'tempA correlation is made, i.e. the final large water droplet collection rate β is the large water droplet collection rate without considering the splash minus the large water droplet collection rate of the splashed rotor surface plus the secondary large water droplet collection rate.
In an example, the method further comprises the steps of:
s01: calculating the collection rate beta of large water drops on the surface of the rotor under the condition of no splashingtemp(ii) a Wherein, the non-splash condition means that the super-cooled large water drops do not splash after impacting the surface of the rotor wing, and the collection rate of the large water drops is 100 percent of the numerical calculation value.
S02: calculating the mass loss rate f of water drops on the surface of the rotor wing under the splashing conditionm. 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 magnitude of the rotor surface under no-spatter conditions is calculatedWater droplet collection rate betatempThe method comprises the following substeps:
s011: generating a grid omega of a rotor wing calculation shape gamma; the rotor wing calculates the shape gamma, namely a vector model of the rotor wing, and the grid omega is a basic calculation composition unit comprising the vector model and a calculation domain. Specifically, the grid Ω generation specifically includes: and (3) adopting grid generation software to take the calculation shape gamma as input so as to generate a calculation topology and a grid omega. The grid generation software includes GRIDGEN, POINTWISE, GRIDSTAR, and the like, and the GRIDSTAR grid generation software employed in the present embodiment generates the grid Ω of the calculation shape Γ.
S012: calculating the information of the air flow field of the grid omega; specifically, air flow field calculation software or a program is adopted, a grid omega of a calculation shape gamma is used as input, information such as an air flow field calculation method, boundary conditions and calculation conditions is set, and air flow field information P corresponding to the grid is obtained through calculation.
S013: calculating water drop flow field information according to the air flow field information; specifically, water drop flow field calculation software or a program is adopted, an air flow field p is used as input, information such as a water drop flow field calculation method, boundary conditions, calculation conditions and the like is set, the surface of the calculation shape gamma is set as a wall surface suction boundary condition, and water drop flow field information corresponding to the grid omega is obtained through calculation and is marked as W; more specifically, the water droplet flow field information includes a water droplet average particle diameter, a water droplet density, a normal velocity at a water droplet airfoil, a water droplet surface tension coefficient, a water droplet dynamic viscosity, a water droplet incidence frequency, a water content of a calculated profile Γ, a water droplet velocity, an incoming flow water content, a far-field water droplet velocity, and the like.
S014: calculating the collection rate beta of the big water drops according to the grid omega and the water drop flow field informationtemp. In particular, the large water droplet collection rate βtempThe volume fraction alpha and the velocity v of water drops at the grid points (x, y, z) of the surface of the profile gamma can be calculated, and the volume fraction alpha of incoming flow of the far field∞And tip speed v∞The calculation is obtained, and the specific calculation formula is as follows:
where n is the normal to the surface of the profile Γ.
In one example, the mass loss rate f of water droplets on the rotor surface under splash conditions is calculatedmThe method comprises the following substeps:
s021: calculating an impact parameter K of water drops impacting the surface of the rotor wing; specifically, given the spatter criterion, the formula for the impact parameter K is:
where ρ represents a water droplet density; d represents an average particle diameter of water droplets; v represents the normal velocity at the water droplet airfoil; σ represents a water droplet surface tension coefficient; μ represents the dynamic viscosity of water droplets; θ represents the angle between the water droplet and the collision surface; Λ represents the water drop incident frequency, Λ ═ 1.5 α∞ 13。
S021: calculating the mass loss rate f according to the impact parameter K and the water drop flow field information on the basis of considering the minimum mass loss ratem. In this example, the minimum mass loss rate is preferably 0.2, and the calculation formula of the mass loss rate is:
fm(x,y,z)＝max{0.7(1-sinθ)[1-e-0.0092(K(x,y,z)-200)],0.2}
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＝a1·(d-a2)2+a3
wherein a represents a minimum mass loss rate; a is1,a2,a3Coefficients representing positive real numbers; d represents the average particle diameter of water droplets. Preferably, a1＝9.92×10-6；a2＝50；a3＝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 no splash conditionLarge water drop collection rate beta on lower rotor surfacetemp；
S2': calculating the mass loss rate f of water drops on the surface of the rotor wing under the splashing conditionm；
S3': calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact on rotor wing'temp；
S4': according to the secondary large water droplet collection rate beta'tempCombined with the collection rate beta of large water drops on the surface of the rotor under the condition that the large water drops do not splashtempAnd the mass loss rate f of water drops on the surface of the rotor under the condition of large water drop splashingmThe final large water droplet collection rate beta of the rotor surface was calculated.
The invention provides a method for calculating the water drop collection rate of the surface of a rotor wing by considering secondary impact under the condition of large water drop splashing. 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 calculation method of the splashed water drops is further provided, and a rotor wing surface large water drop collection rate result is finally 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 collection rate calculation method considering secondary impact formed by any one or a combination of the above examples, and on which computer instructions are stored, which when executed, perform the steps of the rotor surface large water droplet collection 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 method for calculating a collection rate of large water droplets on a rotor surface 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 method for calculating a collection rate of large water droplets on a rotor surface 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 (10)
1. The rotor wing surface large water drop collection rate calculation method considering the secondary impact is characterized by comprising the following steps of: the method comprises the following steps:
calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact on rotor wing'temp；
According to the secondary large water droplet collection rate beta'tempCombined with the collection rate beta of large water drops on the surface of the rotor under the condition that the large water drops do not splashtempAnd large water droplet splashing stripMass loss rate f of water drops on surface of lower rotormThe final large water droplet collection rate beta of the rotor surface was calculated.
2. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 1, wherein: calculating secondary large water drop collection rate beta 'of rotor surface under condition of large water drop secondary impact rotor wing'tempThe method comprises the following substeps:
calculating the total quantity beta of the secondary large water drop collection rate on the surface of the rotor under the condition that the large water drops impact the rotor secondarilyimp；
Calculating the ratio lambda of the total collection rate of the secondary water drops to the total collection rate of the primary large water drops;
according to the proportion of lambda and the collection rate of large water drops betatempCalculating secondary big water drop collection rate beta'temp。
3. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 2, characterized in that: the total secondary large water drop collection rate betaimpThe calculation formula of (2) is as follows:
wherein (x, y, z) represents grid points of the calculated contoured surface of the rotor; ω represents the set of grid points of the rotor's calculated profile surface.
4. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 2, characterized in that: the calculation formula of the ratio lambda is as follows:
wherein (x, y, z) represents grid points of the rotor's computational profile surface;
the efficiency of hitting a water droplet on the surface of the rotor is shown as 1, and the total sum is carried out;
a set of grid points representing the secondary impact of a large water droplet on the calculated profile surface of the rotor.
5. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 2, characterized in that: the secondary large water droplet collection rate beta'tempThe calculation formula of (2) is as follows:
β′temp(x,y,z)＝λ·βtemp(x,y,z)
where (x, y, z) represents grid points of the calculated contoured surface of the rotor.
6. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 1, wherein: the final formula for calculating the collection rate beta of the large water drops is as follows:
β(x,y,z)＝βtemp(x,y,z)-fm(x,y,z)·βtemp(x,y,z)+λ·βtemp(x,y,z)
＝(1-fm(x,y,z)+λ)·βtemp(x,y,z)
wherein (x, y, z) represents grid points of the calculated contoured surface of the rotor; λ represents the ratio of the total amount of the second water droplet collection to the total amount of the first large droplet collection rate.
7. The rotor surface large water droplet collection rate calculation method considering the secondary impact according to claim 1, wherein: the method further comprises the steps of:
calculating the collection rate beta of the large water drops on the surface of the rotor wing under the condition of not generating the splashing of the large water dropstemp；
Calculating the mass loss rate f of water drops on the surface of the rotor wing under the splashing conditionm。
8. The rotor surface large water droplet collection rate calculation method considering secondary impact according to claim 7, wherein: the calculation of the large water droplet collection rate beta of the rotor surface under the condition of no water droplet splashingtempThe method comprises the following substeps:
generating a grid omega of a rotor wing calculation shape gamma;
calculating the information of the air flow field of the grid omega;
calculating water drop flow field information according to the air flow field information;
calculating the collection rate beta of the big water drops according to the grid omega and the water drop flow field informationtemp。
9. The rotor surface large water droplet collection rate calculation method considering secondary impact according to claim 8, wherein: calculating the mass loss rate f of water drops on the surface of the rotor wing under the splashing conditionmThe method comprises the following substeps:
calculating an impact parameter K of water drops impacting the surface of the rotor wing;
calculating the mass loss rate f according to the impact parameter K and the water drop flow field information on the basis of considering the minimum mass loss ratem。
10. 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 a rotor surface large water droplet collection rate considering secondary impacts according to any one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210030651.XA CN114386273B (en) | 2022-01-12 | 2022-01-12 | Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210030651.XA CN114386273B (en) | 2022-01-12 | 2022-01-12 | Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114386273A true CN114386273A (en) | 2022-04-22 |
| CN114386273B CN114386273B (en) | 2023-06-02 |
Family
ID=81201478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210030651.XA Active CN114386273B (en) | 2022-01-12 | 2022-01-12 | Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114386273B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116401834A (en) * | 2023-03-17 | 2023-07-07 | 中国民航大学 | Method and device for calculating impact characteristics of water drops of airplane, equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
-
2022
- 2022-01-12 CN CN202210030651.XA patent/CN114386273B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 (3)
| Title |
|---|
| 周峰;张淼;黄炜;: "二元翼型结冰数值模拟研究", 民用飞机设计与研究 * |
| 桑为民;贾韫泽;鲁天;: "水滴撞击飞溅效应对过冷大水滴结冰影响研究", 西北工业大学学报 * |
| 申晓斌;张志强;林贵平;穆作栋;卜雪琴;: "旋转部件复杂表面水滴撞击计算", 空气动力学学报 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116401834A (en) * | 2023-03-17 | 2023-07-07 | 中国民航大学 | Method and device for calculating impact characteristics of water drops of airplane, equipment and storage medium |
| CN116401834B (en) * | 2023-03-17 | 2024-02-13 | 中国民航大学 | Method and device for calculating impact characteristics of water drops of airplane, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114386273B (en) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113255067B (en) | Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal | |
| CN112906140B (en) | Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate | |
| Wright et al. | Semi-empirical modelling of SLD physics | |
| US6868721B2 (en) | Morphogenetic modelling of in-flight icing | |
| Cao et al. | Numerical simulation of supercooled large droplet icing phenomenon: a review | |
| CN103903481B (en) | The threshold value and envelope curve design method of a kind of ground proximity warning system | |
| CN114386273A (en) | Rotor surface large water drop collection rate calculation method considering secondary impact and terminal | |
| CN112678189B (en) | Improved icing sensor installation position determining method | |
| CN109558650A (en) | The analysis method that lifting airscrew icing influences rotary wing performance | |
| CN111396269A (en) | Multi-time-step unsteady icing calculation method and system and storage medium | |
| CN108090246A (en) | A kind of anti-ice components Calculation Method of Temperature Field | |
| Hashino et al. | The Spectral Ice Habit Prediction System (SHIPS). Part IV: Box model simulations of the habit-dependent aggregation process | |
| EP2215572B1 (en) | Computer-aided method for predicting particles capture by a surface of a moving object | |
| CN114048552B (en) | Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal | |
| Rutkowski et al. | Numerical study of droplet splashing and re-impingement | |
| GUFFOND et al. | Overview of icing research at ONERA | |
| CN111079296B (en) | Aircraft component and aircraft flight load assessment method | |
| WO2019186151A1 (en) | Methods and apparatus for simulating liquid collection on aerodynamic components | |
| Wright et al. | Computational simulation of large droplet icing | |
| Lennie et al. | Development of ice throw model for wind turbine simulation software QBlade | |
| CN108090257B (en) | Method for evaluating probability of collision of flying birds with wind shield over oar | |
| Isobe et al. | Numerical simulation of ice accretion on the rotor blade of a jet engine considering splash and bounce | |
| Li et al. | A numerical study on rime ice accretion characteristics for wind turbine blades | |
| CN117875224B (en) | Non-equal-division continuous multi-particle-size water drop coupling Euler method | |
| CN111516898A (en) | Rotor blade icing aerodynamic characteristic evaluation method, electronic product and storage device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |