[go: up one dir, main page]

CN112906140A - Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate - Google Patents

Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate Download PDF

Info

Publication number
CN112906140A
CN112906140A CN202110462830.6A CN202110462830A CN112906140A CN 112906140 A CN112906140 A CN 112906140A CN 202110462830 A CN202110462830 A CN 202110462830A CN 112906140 A CN112906140 A CN 112906140A
Authority
CN
China
Prior art keywords
water droplet
water
mass loss
splashing
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
Application number
CN202110462830.6A
Other languages
Chinese (zh)
Other versions
CN112906140B (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 CN202110462830.6A priority Critical patent/CN112906140B/en
Publication of CN112906140A publication Critical patent/CN112906140A/en
Application granted granted Critical
Publication of CN112906140B publication Critical patent/CN112906140B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for calculating a water drop collection rate based on large water drop splashing and minimum mass loss rate, belonging to the technical field of aerospace, and comprising the following steps: determining a calculated profile of a simulated aircraft
Figure 943867DEST_PATH_IMAGE001
Grid omega and water drop flow field information; calculating the calculated profile without spatter
Figure 41398DEST_PATH_IMAGE001
First water droplet collection rate
Figure 423838DEST_PATH_IMAGE002
(ii) a Calculating the mass loss rate f caused by splashing on the basis of considering the minimum mass loss ratem(ii) a At the first water droplet collection rate
Figure 10939DEST_PATH_IMAGE002
On the basis of which the mass loss is removed to obtain a calculated profile
Figure 377198DEST_PATH_IMAGE001
Final water droplet collection rate
Figure DEST_PATH_IMAGE003
. According to the invention, on the basis of considering the minimum quality loss rate, the quality loss rate caused by splashing is calculated, namely the quality loss caused by splashing is subjected to artificial minimum limitation, so that the problem of little or even no quality loss in the splashing condition is effectively avoided, the influence of the splashing effect on the collection rate is improved to a certain extent, the numerical simulation is closer to the physical reality, and the prediction accuracy is high.

Description

Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate
Technical Field
The invention relates to the technical field of aerospace, in particular to a water droplet collection rate calculation method based on large water droplet splashing and minimum mass loss rate.
Background
The icing of the airplane directly affects the aerodynamic performance of the airplane, can cause serious accidents of airplane damage and death, and has important significance for developing the icing research of the airplane to improve the flight safety and the flight efficiency. Through intensive research for more than half a century, the freezing of the water drops obtains the theoretical and application result of the system. However, the research on the icing of Supercooled Large water droplets (SLD, with an average particle size of more than 50 μm) just starts, and a plurality of key problems need to be broken through and solved.
SLD icing was first discovered in the American eagle aviation ATR-72-212 flight accident, often in freezing weather conditions such as sleet and frigid rain. Compared with the small water drops, the SLD has larger size, so that the pneumatic following performance is poorer, the impact icing range is wider, and the icing quantity is more. Meanwhile, the SLD is used as a more obvious stress body, the dynamic behavior is more obvious, and the effects of deformation, breakage, rebound, splashing and the like are often generated, so that the icing process is more complicated, and the research on the problems of SLD icing test, numerical simulation, ice prevention/removal and the like is harder.
The SLD leaves the wall surface of the airplane partially or completely with certain kinetic energy, and finally influences the local water drop collection coefficient of the icing surface, namely SLD splashing is a main influence factor of the water drop collection rate, and the SLD splashing is easier to record and quantify and is a breakthrough of the current SLD icing research. Based on wind tunnel test data, various dynamic models related to the SLD appear, and certain results are obtained, however, when the conventional dynamic model is used for an icing wind tunnel test under the SLD condition, the collection rate of the large water droplet splash effect is still different from the corresponding collection rate of the water droplet under the real condition, and on the basis, research on a calculation method for the collection rate of the large water droplet splash effect needs to be further developed.
Disclosure of Invention
The invention aims to solve the problem that the collection rate of the large water droplet splashing effect in the prior art is not high enough in calculation accuracy, and provides a water droplet collection rate calculation method based on large water droplet splashing and the minimum mass loss rate.
The purpose of the invention is realized by the following technical scheme: a method of calculating the water droplet collection rate based on large droplet spatter and minimum mass loss rate, the method comprising calculating the profile of a simulated aircraft
Figure 471097DEST_PATH_IMAGE001
The surface water drop collection rate step specifically comprises the following steps:
determining a calculated profile of a simulated aircraft
Figure 151478DEST_PATH_IMAGE001
The grid omega and the water drop flow field information of the simulated aircraft;
calculating the calculation appearance under the condition of no splashing according to the grid omega and the water drop flow field information
Figure 31709DEST_PATH_IMAGE001
First water droplet collection rate
Figure 312386DEST_PATH_IMAGE002
；
On the basis of considering the minimum mass loss rate, calculating the mass loss rate f caused by splashing according to the grid omega and the water drop flow field informationm；
At the first water droplet collection rate
Figure 467424DEST_PATH_IMAGE002
On the basis of which the mass loss is removed to obtain a calculated profile
Figure 2311DEST_PATH_IMAGE001
Final water droplet collection rate
Figure 928810DEST_PATH_IMAGE003
。
As an option, the calculated profile of the simulated aircraft is determined
Figure 181936DEST_PATH_IMAGE001
The grid Ω specifically includes: using mesh generation software to compute the shape
Figure 717829DEST_PATH_IMAGE001
As input, a computational topology and mesh are generated.
As an option, the determining information of the water droplet flow field of the simulated aircraft specifically includes: calculating to obtain a calculated shape by adopting water drop flow field software or program and taking air flow field information as input
Figure 513746DEST_PATH_IMAGE001
The corresponding water drop flow field information of the grid.
As an option, the calculating of the air flow field information specifically includes: using air flow field calculation software or programs to calculate the profile
Figure 125993DEST_PATH_IMAGE001
The grids are used as input, and air flow field information corresponding to the grids is obtained through calculation.
As an option, the mass loss rate fmThe calculation formula of (2) is as follows:
fm=max{0.7（1-sin
Figure 617149DEST_PATH_IMAGE004
）[1-e-0.092(K-200)],
Figure 113989DEST_PATH_IMAGE005
}
wherein,
Figure 623468DEST_PATH_IMAGE004
is the angle between the water droplet and the collision surface; k is a water drop impact parameter;
Figure 390304DEST_PATH_IMAGE006
in order to minimize the rate of mass loss,
Figure 86865DEST_PATH_IMAGE006
are positive real numbers.
As an option, the minimum mass loss rate
Figure 200446DEST_PATH_IMAGE006
As a function of the incoming flow conditions with respect to any of average particle size, angle of incidence, and impact velocity.
As an option, the minimum mass loss rate
Figure 564431DEST_PATH_IMAGE006
The relationship with the average particle size of the water droplets is:
Figure 33327DEST_PATH_IMAGE006
=
Figure 686025DEST_PATH_IMAGE007
•
Figure 134455DEST_PATH_IMAGE008
+
Figure 618526DEST_PATH_IMAGE009
wherein d is the average particle size of water droplets;
Figure 884423DEST_PATH_IMAGE007
,
Figure 227679DEST_PATH_IMAGE010
,
Figure 712756DEST_PATH_IMAGE011
is a coefficient of a positive real number,
Figure 316913DEST_PATH_IMAGE007
=9.92×10-6，
Figure 753710DEST_PATH_IMAGE012
=50，
Figure 928471DEST_PATH_IMAGE009
=0.12。
as an option, the impact parameter calculation expression is:
Figure 108917DEST_PATH_IMAGE013
wherein,
Figure 833159DEST_PATH_IMAGE014
is the water drop density; d is the average particle size of the water droplets;
Figure 440858DEST_PATH_IMAGE015
is the normal velocity at the water droplet airfoil;
Figure 335870DEST_PATH_IMAGE016
is the water droplet surface tension coefficient;
Figure 444641DEST_PATH_IMAGE017
is the dynamic viscosity of water droplets;
Figure 898756DEST_PATH_IMAGE018
is the frequency of incidence of the water droplets,
Figure 552722DEST_PATH_IMAGE018
=1.5
Figure 561129DEST_PATH_IMAGE019
。
as an option, the first water droplet collection rate
Figure 208011DEST_PATH_IMAGE002
The calculation formula of (2) is as follows:
Figure 516633DEST_PATH_IMAGE020
•
Figure 840036DEST_PATH_IMAGE021
wherein,
Figure 866898DEST_PATH_IMAGE022
for calculating the shape
Figure 51891DEST_PATH_IMAGE001
Water content of (2); v represents the water drop velocity; n is a shape
Figure 480599DEST_PATH_IMAGE001
The normal vector of the surface;
Figure 476368DEST_PATH_IMAGE023
indicating the incoming flow moisture content;
Figure 724946DEST_PATH_IMAGE024
representing the far field water drop velocity.
As an option, the computing profile
Figure 713631DEST_PATH_IMAGE001
The surface droplet collection rate step further comprises determining as yes according to a given spatter determination criterionWhether the splashing happens or not, the splashing judgment criterion is specifically as follows: the impact parameter K is more than 200, and when the impact parameter is more than 200, the splash is judged to occur; the expression for the impact parameters is:
Figure 996845DEST_PATH_IMAGE013
wherein,
Figure 662050DEST_PATH_IMAGE014
is the water drop density; d is the average particle size of the water droplets;
Figure 397925DEST_PATH_IMAGE015
is the normal velocity at the water droplet airfoil;
Figure 190300DEST_PATH_IMAGE016
is the water droplet surface tension coefficient;
Figure 328021DEST_PATH_IMAGE017
is the dynamic viscosity of water droplets;
Figure 665592DEST_PATH_IMAGE018
is the frequency of incidence of the water droplets,
Figure 482239DEST_PATH_IMAGE018
=1.5
Figure 219250DEST_PATH_IMAGE019
；
Figure 585379DEST_PATH_IMAGE023
indicating the incoming flow moisture content;
Figure 218485DEST_PATH_IMAGE004
is the angle between the water droplet and the impact surface.
It should be further noted that the technical features corresponding to the above-mentioned method options can be combined with each other or replaced to form a new technical solution.
Compared with the prior art, the invention has the beneficial effects that:
according to the invention, on the basis of considering the minimum quality loss rate, the quality loss rate caused by splashing is calculated, namely the quality loss caused by splashing is subjected to artificial minimum limitation, so that the problem of little or even no quality loss when splashing occurs is effectively avoided, the influence of the splashing effect on the collection rate is improved to a certain extent, the numerical simulation is closer to the physical reality, and the prediction accuracy is high.
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 flowchart of a method of example 1 of the present invention;
FIG. 2 is a calculated profile of example 1 of the present invention
Figure 522428DEST_PATH_IMAGE001
A grid schematic of (a);
FIG. 3 is a two-dimensional airfoil of MS317 according to example 1 of the present invention;
FIG. 4 is a curve for calculating the collection rate of water droplets not splashed in example 1 of the present invention;
FIG. 5 is a curve of water droplet collection rate calculation obtained by the method of the present invention in example 1 of the present invention;
FIG. 6 is a graph comparing the results of the collection rate calculations of the method of the present invention in example 1 of the present invention and the prior art.
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.
According to the invention, an air flow field and a water drop flow field with calculated shapes are obtained through numerical simulation, then a splash judgment criterion is given, the mass loss rate is determined, meanwhile, the limitation of minimum mass loss is carried out in the loss rate calculation, so that certain mass loss is ensured after large water drops splash, and finally, a water drop collection rate calculation result considering the splash effect is obtained, so that the numerical simulation is closer to the physical reality.
Example 1
As shown in FIG. 1, in example 1, a method for calculating a water droplet collection rate based on a large water droplet splash and a minimum mass loss rate includes calculating a profile of a simulated aircraft
Figure 797551DEST_PATH_IMAGE001
The surface water drop collection rate step specifically comprises the following steps:
s01: determining a calculated profile of a simulated aircraft
Figure 519651DEST_PATH_IMAGE001
The grid omega and the water drop flow field information of the simulated aircraft; wherein the shape is calculated
Figure 448292DEST_PATH_IMAGE001
I.e. the profile of the aircraft on which the water drops impact by splashing, the invention is mainly directed to the wing profile of the aircraft, i.e. the calculated profile
Figure 380476DEST_PATH_IMAGE001
The vector model of the aircraft wing is a calculation shape of the MS317 airfoil profile based on the consideration of the splash and the minimum mass loss rate in the embodiment
Figure 591050DEST_PATH_IMAGE001
Surface water droplet collection rate, a two-dimensional MS317 airfoil plot as shown in fig. 2, where the abscissa represents the x-direction of the airfoil in mm; the ordinate represents the Y direction of the airfoil in mm.
S02: calculating the calculation appearance under the condition of no splashing according to the grid omega and the water drop flow field information
Figure 26711DEST_PATH_IMAGE001
First water droplet collection rate
Figure 391833DEST_PATH_IMAGE002
；
S03: on the basis of considering the minimum mass loss rate, calculating the mass loss rate f caused by splashing according to the grid omega and the water drop flow field informationm；
S04: at the first water droplet collection rate
Figure 811313DEST_PATH_IMAGE002
On the basis of which the mass loss is removed to obtain a calculated profile
Figure 772447DEST_PATH_IMAGE001
Final water droplet collection rate
Figure 452827DEST_PATH_IMAGE003
。
According to the invention, on the basis of considering the minimum quality loss rate, the quality loss rate caused by splashing is calculated, namely the quality loss caused by splashing is subjected to artificial minimum limitation, so that the problem of little or even no quality loss when splashing occurs is effectively avoided, the influence of the splashing effect on the collection rate is improved to a certain extent, the numerical simulation is closer to the physical reality, and the prediction accuracy is high.
Further, the calculated profile of the simulated aircraft is determined in step S01
Figure 864217DEST_PATH_IMAGE001
The grid Ω specifically includes:
using mesh generation software to compute the shape
Figure 144894DEST_PATH_IMAGE001
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
Figure 299932DEST_PATH_IMAGE001
As shown in fig. 3.
Further, the determination of the information of the water droplet flow field of the simulated aircraft in step S01 specifically includes: adopting water drop flow field software or program, using air flow field information P as input, and setting information of water drop flow field calculation method, boundary condition and calculation condition, etc., in this embodiment, the calculation appearance is
Figure 834819DEST_PATH_IMAGE001
Is set as a wall suction boundary condition, and then calculated to obtain a calculated profile
Figure 292476DEST_PATH_IMAGE001
The corresponding water droplet flow field information W of grid Ω. More specifically, the water drop flow field information comprises the average particle diameter of the water drop, the density of the water drop, the normal speed at the airfoil surface of the water drop, the surface tension coefficient of the water drop, the dynamic viscosity of the water drop, the incidence frequency of the water drop and the calculated shape
Figure 420969DEST_PATH_IMAGE001
Water content, water droplet velocity, incoming flow water content, far field water droplet velocity, etc.
Further, before the water droplet flow field information calculation step, an air flow field information calculation step is further included, and the method specifically includes:
using air flow field calculation software or programs to calculate the profile
Figure 973173DEST_PATH_IMAGE001
The grid omega is used as input, and information such as an air flow field calculation method, boundary conditions, calculation conditions and the like is set, in the embodiment, the attack angle of water drops is 0 degrees, and the speed of far-field water drops is
Figure 877413DEST_PATH_IMAGE025
And calculating to obtain the air flow field information P corresponding to the grids at 78 m/s.
Further, the mass loss rate fmThe calculation formula of (2) is as follows:
fm=max{0.7（1-sin
Figure 365026DEST_PATH_IMAGE026
）[1-e-0.092(K-200)],
Figure 371028DEST_PATH_IMAGE006
}
wherein,
Figure 477656DEST_PATH_IMAGE026
is the angle between the water droplet and the collision surface; k is a water drop impact parameter;
Figure 128080DEST_PATH_IMAGE006
in order to minimize the rate of mass loss,
Figure 176807DEST_PATH_IMAGE006
are positive real numbers. Preferably, the minimum mass loss rate is 0.2.
Further, minimum mass loss rate
Figure 14313DEST_PATH_IMAGE006
As a function of the incoming flow conditions with respect to any of average particle size, angle of incidence, and impact velocity. As an embodiment, minimum mass loss rate
Figure 423167DEST_PATH_IMAGE006
The relationship with the average particle size of the water droplets is:
Figure 52731DEST_PATH_IMAGE006
=
Figure 147726DEST_PATH_IMAGE027
•
Figure 347894DEST_PATH_IMAGE028
+
Figure 186538DEST_PATH_IMAGE029
wherein d is the average particle size of water droplets;
Figure 936188DEST_PATH_IMAGE027
,
Figure 202084DEST_PATH_IMAGE010
,
Figure 653663DEST_PATH_IMAGE011
is a coefficient of a positive real number,
Figure 948859DEST_PATH_IMAGE027
=9.92×10-6，
Figure 21857DEST_PATH_IMAGE030
=50，
Figure 193076DEST_PATH_IMAGE029
=0.12。
further, the impact parameter calculation expression is:
Figure 600792DEST_PATH_IMAGE013
wherein,
Figure 905872DEST_PATH_IMAGE014
is the water drop density; d is the average particle size of the water droplets;
Figure 239901DEST_PATH_IMAGE015
is the normal velocity at the water droplet airfoil;
Figure 519704DEST_PATH_IMAGE016
is the water droplet surface tension coefficient;
Figure 398405DEST_PATH_IMAGE017
is the dynamic viscosity of water droplets;
Figure 772754DEST_PATH_IMAGE018
is the frequency of incidence of the water droplets,
Figure 961290DEST_PATH_IMAGE018
=1.5
Figure 5469DEST_PATH_IMAGE019
(ii) a Is the incoming water content.
Further, the first water droplet collection rate
Figure 420401DEST_PATH_IMAGE031
The calculation formula of (2) is as follows:
Figure 208229DEST_PATH_IMAGE033
•
Figure 641484DEST_PATH_IMAGE035
wherein,
Figure 856565DEST_PATH_IMAGE022
for calculating the shape
Figure 883427DEST_PATH_IMAGE001
Water content of (2); v represents the water drop velocity; n is a shape
Figure 583267DEST_PATH_IMAGE001
The normal vector of the surface;
Figure 746395DEST_PATH_IMAGE023
indicating the incoming flow moisture content;
Figure 257011DEST_PATH_IMAGE024
the far-field water drop velocity is represented, the calculation curve of the water drop collection rate without splashing in the embodiment is shown in fig. 4, the abscissa of fig. 4 is the normalization result of the y-direction coordinate in fig. 2 (dimensionless chord length in fig. 2), wherein c is the chord length; the ordinate is the first water droplet collection rate
Figure 505590DEST_PATH_IMAGE031
。
Further, the final water droplet collection rate
Figure 245007DEST_PATH_IMAGE003
The calculation formula of (2) is as follows:
Figure 793800DEST_PATH_IMAGE003
=(1-fm)•
Figure 85104DEST_PATH_IMAGE002
. More specifically, in the present embodiment, the incoming flow water content
Figure 945612DEST_PATH_IMAGE023
Is 1kg/m3Far field water drop velocity
Figure 878933DEST_PATH_IMAGE024
At 78m/s, water droplet density
Figure 390555DEST_PATH_IMAGE014
Is 1000 kg/m3The average diameter d of the water drop is 92 mu m, and the surface tension coefficient of the water drop
Figure 852760DEST_PATH_IMAGE016
7.56, dynamic viscosity of water droplet
Figure 341511DEST_PATH_IMAGE017
Is 1.7921X 10-3Pas, frequency of water drop incidence
Figure 671998DEST_PATH_IMAGE036
In the present embodiment, the final water drop collection rate calculation curve is shown in fig. 5, and the abscissa of fig. 5 is the normalization result of the y-direction coordinate in fig. 2 (the dimensionless chord length in fig. 2), where c is the chord length; ordinate final water droplet collection rate
Figure 929804DEST_PATH_IMAGE003
As can be seen from fig. 4, the water droplet collecting rate calculated by the method of the present invention is smaller than the water droplet collecting rate when no splashing occurs, and is closer to the physical reality.
Further, before the step of calculating the quality loss rate caused by spatter in step S03, the method further includes:
judging whether the splash happens according to a given splash judging criterion, wherein the splash judging criterion is as follows: parameters of impactKAnd if the impact parameter is more than 200, judging that the splash occurs.
To further illustrate the technical effects of the present invention, the water droplet collection rate calculated by the method of the present invention, the water droplet collection rate obtained based on the no-spatter model, the water droplet collection rate obtained based on the Mundo model, and the actual test result are compared, for example, as shown in fig. 6, the abscissa of fig. 6 is the normalization result of the y-direction coordinate in fig. 2 (the dimensionless chord length in fig. 2), wherein c is the chord length; water droplet harvesting with final ordinateConcentration ratio
Figure 703856DEST_PATH_IMAGE003
It can be obviously seen that the No-splash model (No-SLD) does not consider the mass loss of broken large water drops, and the obtained results are all larger than the test results (Exp), and the water drop collection rate calculated by the method (deployed) under the premise of considering the splash of the large water drops and the minimum mass loss rate is closer to the test value in numerical simulation and higher in prediction precision compared with the existing Mundo model.
Example 2
The present embodiment provides a storage medium having the same inventive concept as embodiment 1, and having stored thereon computer instructions which, when executed, perform the steps of the droplet collection rate calculation method based on the large droplet splash and the minimum mass loss rate in embodiment 1.
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, which is stored in a storage medium and includes several instructions to enable 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 of 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.
Example 3
The present embodiment also provides a terminal, which has the same inventive concept as embodiment 1, 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 execute the steps of the method for calculating the water droplet collection rate based on the splashing of large water droplets and the minimum mass loss rate in embodiment 1. 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 method for calculating the water drop collection rate based on the large water drop splashing and the minimum mass loss rate is characterized by comprising the following steps of: the method includes calculating a profile of the simulated aircraft
Figure 760947DEST_PATH_IMAGE001
The surface water drop collection rate step specifically comprises the following steps:
determining a calculated profile of a simulated aircraft
Figure 302787DEST_PATH_IMAGE001
The grid omega and the water drop flow field information of the simulated aircraft;
calculating the calculation appearance under the condition of no splashing according to the grid omega and the water drop flow field information
Figure 677268DEST_PATH_IMAGE001
First water droplet collection rate
Figure 570137DEST_PATH_IMAGE002
；
On the basis of considering the minimum mass loss rate, calculating the mass loss rate f caused by splashing according to the grid omega and the water drop flow field informationm；
At the first water droplet collection rate
Figure 688266DEST_PATH_IMAGE002
On the basis of which the mass loss is removed to obtain a calculated profile
Figure 84612DEST_PATH_IMAGE001
Final water droplet collection rate
Figure 629994DEST_PATH_IMAGE003
。
2. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 1, wherein the method comprises the following steps: determining a calculated profile of a simulated aircraft
Figure 744581DEST_PATH_IMAGE001
The grid Ω specifically includes:
using mesh generation software to compute the shape
Figure 164936DEST_PATH_IMAGE001
As input, a computational topology and mesh are generated.
3. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 1, wherein the method comprises the following steps: the determining of the water droplet flow field information of the simulated aircraft specifically comprises:
calculating to obtain a calculated shape by adopting water drop flow field software or program and taking air flow field information as input
Figure 415788DEST_PATH_IMAGE001
The corresponding water drop flow field information of the grid.
4. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 3, wherein the method comprises the following steps: the calculation of the air flow field information specifically includes:
using air flow field calculation software or programsTo calculate the profile
Figure 132072DEST_PATH_IMAGE001
The grids are used as input, and air flow field information corresponding to the grids is obtained through calculation.
5. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 1, wherein the method comprises the following steps: the mass loss rate fmThe calculation formula of (2) is as follows:
fm=max{0.7（1-sin
Figure 468375DEST_PATH_IMAGE004
）[1-e-0.092(K-200)],
Figure 928306DEST_PATH_IMAGE005
}
wherein,
Figure 299245DEST_PATH_IMAGE004
is the angle between the water droplet and the collision surface; k is a water drop impact parameter;
Figure 186429DEST_PATH_IMAGE006
in order to minimize the rate of mass loss,
Figure 275608DEST_PATH_IMAGE006
are positive real numbers.
6. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 5, wherein the method comprises the following steps: the minimum mass loss rate
Figure 772186DEST_PATH_IMAGE006
As a function of the incoming flow conditions with respect to any of average particle size, angle of incidence, and impact velocity.
7. Root of herbaceous plantThe method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 5, wherein the method comprises the following steps: the minimum mass loss rate
Figure 732052DEST_PATH_IMAGE006
The relationship with the average particle size of the water droplets is:
Figure 180351DEST_PATH_IMAGE006
=
Figure 366613DEST_PATH_IMAGE007
•
Figure 292980DEST_PATH_IMAGE008
+
Figure 513877DEST_PATH_IMAGE009
wherein d is the average particle size of water droplets;
Figure 867498DEST_PATH_IMAGE007
,
Figure 541056DEST_PATH_IMAGE010
,
Figure 271115DEST_PATH_IMAGE011
is a coefficient of a positive real number,
Figure 110633DEST_PATH_IMAGE007
=9.92×10-6，
Figure 510521DEST_PATH_IMAGE012
=50，
Figure 61588DEST_PATH_IMAGE009
=0.12。
8. the method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 5, wherein the method comprises the following steps: the impact parameter calculation expression is as follows:
Figure 205125DEST_PATH_IMAGE013
wherein,
Figure 259668DEST_PATH_IMAGE014
is the water drop density; d is the average particle size of the water droplets;
Figure 96037DEST_PATH_IMAGE015
is the normal velocity at the water droplet airfoil;
Figure 868821DEST_PATH_IMAGE016
is the water droplet surface tension coefficient;
Figure 783425DEST_PATH_IMAGE017
is the dynamic viscosity of water droplets;
Figure 958055DEST_PATH_IMAGE018
is the frequency of incidence of the water droplets,
Figure 699746DEST_PATH_IMAGE018
=1.5
Figure 694247DEST_PATH_IMAGE019
。
9. the method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 1, wherein the method comprises the following steps: the first water droplet collection rate
Figure 304219DEST_PATH_IMAGE020
The calculation formula of (2) is as follows:
Figure 208722DEST_PATH_IMAGE021
•
Figure 511527DEST_PATH_IMAGE022
wherein,
Figure 868690DEST_PATH_IMAGE024
for calculating the shape
Figure 282354DEST_PATH_IMAGE001
Water content of (2); v represents the water drop velocity; n is a shape
Figure 805477DEST_PATH_IMAGE001
The normal vector of the surface;
Figure 13604DEST_PATH_IMAGE025
indicating the incoming flow moisture content;
Figure 858064DEST_PATH_IMAGE026
representing the far field water drop velocity.
10. The method for calculating the water droplet collection rate based on the large water droplet splashing and the minimum mass loss rate according to claim 1, wherein the method comprises the following steps: the calculated shape
Figure 809839DEST_PATH_IMAGE001
The surface water droplet collection rate step further comprises judging whether the splash occurs according to a given splash judgment criterion, wherein the splash judgment criterion is specifically as follows: the impact parameter K is more than 200, and when the impact parameter is more than 200, the splash is judged to occur; the expression for the impact parameters is:
Figure 688933DEST_PATH_IMAGE013
wherein,
Figure 333541DEST_PATH_IMAGE014
is the water drop density; d is the average particle size of the water droplets;
Figure 665297DEST_PATH_IMAGE015
is the normal velocity at the water droplet airfoil;
Figure 155184DEST_PATH_IMAGE016
is the water droplet surface tension coefficient;
Figure 121741DEST_PATH_IMAGE017
is the dynamic viscosity of water droplets;
Figure 937250DEST_PATH_IMAGE027
is the frequency of incidence of the water droplets,
Figure 756301DEST_PATH_IMAGE027
=1.5
Figure 315459DEST_PATH_IMAGE019
；
Figure 903566DEST_PATH_IMAGE025
indicating the incoming flow moisture content;
Figure 624397DEST_PATH_IMAGE028
is the angle between the water droplet and the impact surface.
CN202110462830.6A 2021-04-28 2021-04-28 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate Active CN112906140B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110462830.6A CN112906140B (en) 2021-04-28 2021-04-28 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110462830.6A CN112906140B (en) 2021-04-28 2021-04-28 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate

Publications (2)

Publication Number Publication Date
CN112906140A true CN112906140A (en) 2021-06-04
CN112906140B CN112906140B (en) 2021-08-03

Family

ID=76108943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110462830.6A Active CN112906140B (en) 2021-04-28 2021-04-28 Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate

Country Status (1)

Country Link
CN (1) CN112906140B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255067A (en) * 2021-07-16 2021-08-13 中国空气动力研究与发展中心计算空气动力研究所 Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal
CN114048552A (en) * 2022-01-12 2022-02-15 中国空气动力研究与发展中心计算空气动力研究所 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
CN114386273A (en) * 2022-01-12 2022-04-22 中国空气动力研究与发展中心计算空气动力研究所 Rotor surface large water drop collection rate calculation method considering secondary impact and terminal

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080032066A1 (en) * 2005-10-26 2008-02-07 Lars Stiblert Platforms, apparatuses, systems and methods for processing and analyzing substrates
CN104268399A (en) * 2014-09-24 2015-01-07 空气动力学国家重点实验室 Computing method of model parameters in icing wind tunnel experiment under supercooled large droplet condition
CN107092718A (en) * 2017-03-17 2017-08-25 中国人民解放军陆军航空兵学院 Method for numerical simulation during rainfall is met with a kind of aircraft flight
US20180033183A1 (en) * 2016-08-01 2018-02-01 Macrograph Co., Ltd. Splash particle simulation apparatus and method based on virtual particle
US20180293339A1 (en) * 2017-04-10 2018-10-11 Harris Corporation Systems and methods to achieve tailored fluid and energy interaction through surface texturing
CN111680256A (en) * 2020-06-16 2020-09-18 中国空气动力研究与发展中心计算空气动力研究所 Method for calculating icing density of airplane

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080032066A1 (en) * 2005-10-26 2008-02-07 Lars Stiblert Platforms, apparatuses, systems and methods for processing and analyzing substrates
CN104268399A (en) * 2014-09-24 2015-01-07 空气动力学国家重点实验室 Computing method of model parameters in icing wind tunnel experiment under supercooled large droplet condition
US20180033183A1 (en) * 2016-08-01 2018-02-01 Macrograph Co., Ltd. Splash particle simulation apparatus and method based on virtual particle
CN107092718A (en) * 2017-03-17 2017-08-25 中国人民解放军陆军航空兵学院 Method for numerical simulation during rainfall is met with a kind of aircraft flight
US20180293339A1 (en) * 2017-04-10 2018-10-11 Harris Corporation Systems and methods to achieve tailored fluid and energy interaction through surface texturing
CN111680256A (en) * 2020-06-16 2020-09-18 中国空气动力研究与发展中心计算空气动力研究所 Method for calculating icing density of airplane

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
S.C.TAN: ""A TENTATIVE MASS LOSS MODEL FOR SIMULATING WATER DROPLET SPLASH"", 《42ND AIAA AEROSPACE SCIENCES MEETING AND EXHIBIT》 *
ZHIHONG ZHOU等: ""Quantitative detection method for icing of horizontal‐axis wind turbines"", 《WIND ENERGY》 *
朱剑鋆: ""飞行防冰与结冰中的流动传热计算方法研究"", 《中国博士学位论文全文数据库》 *
白俊强等: "过冷大水滴情况下的积冰数值模拟(英文)", 《空气动力学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255067A (en) * 2021-07-16 2021-08-13 中国空气动力研究与发展中心计算空气动力研究所 Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal
CN114048552A (en) * 2022-01-12 2022-02-15 中国空气动力研究与发展中心计算空气动力研究所 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
CN114048552B (en) * 2022-01-12 2022-03-29 中国空气动力研究与发展中心计算空气动力研究所 Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
CN114386273A (en) * 2022-01-12 2022-04-22 中国空气动力研究与发展中心计算空气动力研究所 Rotor surface large water drop collection rate calculation method considering secondary impact and terminal
CN114386273B (en) * 2022-01-12 2023-06-02 中国空气动力研究与发展中心计算空气动力研究所 Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal

Also Published As

Publication number Publication date
CN112906140B (en) 2021-08-03

Similar Documents

Publication Publication Date Title
CN112906140B (en) Water drop collection rate calculation method based on large water drop splashing and minimum mass loss rate
CN113255067B (en) Rotor icing ice shape calculation method considering flapping effect, storage medium and terminal
US6868721B2 (en) Morphogenetic modelling of in-flight icing
CN108052772A (en) A kind of geometrical non-linearity static aeroelastic analysis method based on structure reduced-order model
CN111310381B (en) Three-dimensional water drop collection coefficient calculation method
Ozcer et al. FENSAP-ICE: Numerical prediction of ice roughness evolution, and its effects on ice shapes
CN106777783B (en) Method for predicting blade cracks of aircraft engine
CN111396269A (en) Multi-time-step unsteady icing calculation method and system and storage medium
Baeder et al. OVERTURNS Simulation of S-76 Rotor in Hover
WO2019186151A1 (en) Methods and apparatus for simulating liquid collection on aerodynamic components
CN110489939A (en) A kind of optimization method of Design of Propeller, optimization device and readable storage medium storing program for executing
CN114398844A (en) Steady-state anti-icing simulation method based on continuous water film flow
CN114386273B (en) Rotor wing surface large water drop collection rate calculation method considering secondary impact and terminal
CN111400826B (en) Method and system for predicting ice shape at any moment
Das et al. Ice shape prediction for turbofan rotating blades
CN115099048B (en) Icing numerical simulation method, system, equipment and medium based on Messinger model
Sagaga et al. Acoustic predictions for side-by-side rotor with ground effect
CN116611175A (en) Prediction method for free degree flutter of large aspect ratio aircraft body
CN115098825B (en) Method for releasing water drops through encryption along path, method for acquiring collection rate and medium
CN114048552B (en) Rotor wing surface large water drop mass flow calculation method considering secondary impact and terminal
CN117875224B (en) Non-equal-division continuous multi-particle-size water drop coupling Euler method
Isobe et al. Numerical simulation of ice accretion on the rotor blade of a jet engine considering splash and bounce
CN118194775B (en) Method and device for determining icing thickness of rotary negative pressure fan
CN110750836A (en) Aircraft pitching damping moment calculation method based on steady-state cone motion
Sathyabhama et al. Experimental and computational analysis of bio-inspired winglets for micro air vehicles

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