CN112507572A - Optimized matching evaluation method for icing condition of large supercooled water drops - Google Patents

Optimized matching evaluation method for icing condition of large supercooled water drops Download PDF

Info

Publication number
CN112507572A
CN112507572A CN202110167936.3A CN202110167936A CN112507572A CN 112507572 A CN112507572 A CN 112507572A CN 202110167936 A CN202110167936 A CN 202110167936A CN 112507572 A CN112507572 A CN 112507572A
Authority
CN
China
Prior art keywords
droplet
diameter
characteristic
interval
measured
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
CN202110167936.3A
Other languages
Chinese (zh)
Other versions
CN112507572B (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.)
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Low Speed 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 Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202110167936.3A priority Critical patent/CN112507572B/en
Publication of CN112507572A publication Critical patent/CN112507572A/en
Application granted granted Critical
Publication of CN112507572B publication Critical patent/CN112507572B/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/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

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)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention is suitable for the technical field of wind tunnel tests, and provides an optimal matching evaluation method for the icing condition of supercooled large water drops, which comprises the following steps: selecting a plurality of characteristic liquid drop points from the standard liquid drop distribution; dividing the standard droplet distribution into a plurality of characteristic droplet diameter intervals; measuring to obtain measured liquid drop distribution, and calculating the accumulated volume fraction of the liquid drop diameter of each characteristic liquid drop point in the measured liquid drop distribution; dividing the distribution of the measured liquid drops into a plurality of measured liquid drop diameter intervals; calculating a match bias factorF(ii) a Calculating the optimal threshold value of the matching deviationF c (ii) a According to the matching deviation factorFSum match deviation optimum thresholdF c And determining the matching degree between the measured droplet distribution and the standard droplet distribution. In the invention, the optimal threshold value of the matching deviation is calculated in addition to the matching deviation factor, so that the quality of the matching effect can be further measuredAnd the evaluation is more objective.

Description

Optimized matching evaluation method for icing condition of large supercooled water drops
Technical Field
The invention belongs to the technical field of wind tunnel tests, and particularly relates to an optimization matching evaluation method for an icing condition of large supercooled water drops.
Background
When an airplane flies in a cloud layer, supercooled water drops (namely liquid water drops with the temperature lower than the freezing point) in the cloud layer continuously impact the windward side of the airplane, so that the icing phenomenon of the surface of the airplane is caused. Aircraft icing is widespread in flight practice and poses a serious threat to flight safety.
Supercooled large water droplets generally refer to supercooled water droplets with a droplet diameter of more than 100 μm, and in view of the serious damage of the icing of the supercooled large water droplets to the flight safety of an airplane, the federal aviation administration has officially released appendix O of 25 parts of federal aviation regulations in 2014, given supercooled large water droplet icing meteorological conditions, and explicitly proposed droplet distribution characteristics specific to the supercooled large water droplet icing meteorological conditions, as shown in fig. 1, four typical standard droplet distribution diagrams are shown, which include fine freezing rain with MVD <40 μm, fine freezing rain with MVD >40 μm, freezing rain with MVD <40 μm, freezing rain with MVD >40 μm, wherein the abscissa in fig. 1 is the droplet diameter, the ordinate is the cumulative volume fraction, and the cumulative volume fraction is defined as: a ratio of droplet volume to total droplet volume that is less than a set droplet diameter; in addition, MVD is expressed as volume median diameter.
The icing wind tunnel is an important ground test device for developing airplane icing research and verifying an airplane component ice prevention and removal system, and plays an important role in airplane icing airworthiness examination. The comprehensive test simulation of the icing meteorological condition of the supercooled large water drops is the test simulation capability of the key development of the conventional icing wind tunnel, wherein the liquid drop distribution characteristic of the icing meteorological condition of the supercooled large water drops is the main test simulation content.
When the icing wind tunnel test is performed, the droplet distribution simulated by the icing wind tunnel needs to be matched with the droplet distribution under the icing condition of the supercooled large droplets, for example, in a certain icing wind tunnel test, the freezing rain with the MVD of less than 40 μm needs to be simulated, and then, the droplet distribution simulated by the icing wind tunnel is expected to be as close as possible to the droplet distribution of the freezing rain with the MVD of less than 40 μm, specifically, the droplet distribution simulated by the icing wind tunnel is obtained through measurement or simulation, the droplet distribution of the freezing rain with the MVD of less than 40 μm is given by fig. 1 (namely, belongs to a standard distribution, or is a test standard), and when the measured or simulated droplet distribution simulated by the icing wind tunnel is closer to the droplet distribution of the freezing rain with the MVD of less than 40 μm given in fig. 1, the higher matching degree is proved.
As in the related document 1 (fiscal, SLD simulation method in icing wind tunnel and experimental verification study thereof), fig. 7 shows a comparison of the cumulative volume distribution of the SLD icing conditions determined by numerical simulation and airworthiness standard.
However, in the related document 1, when the degree of coincidence between the numerical simulation result and the airworthiness standard drop spectrum curve is judged, a qualitative visual inspection judgment method is adopted, the method is obviously affected by artificial subjective factors, the quality of the matching degree is difficult to objectively and comprehensively evaluate, the quality of the simulation capability of the icing wind tunnel supercooling large water drop test cannot be objectively evaluated, and the development of the simulation capability of the icing wind tunnel supercooling large water drop test is not facilitated.
In the aspect of other research institutions, main icing wind tunnels such as an IRT icing wind tunnel of the American NASA Glenn center, an NRC AIWT icing wind tunnel of the national research institute of Canada, an Olympic RTA icing wind tunnel and the like all develop the size test research of droplets under the condition of ice formation of large supercooled water droplets, and the aerodynamic research and development center of China also develops preliminary test research, but the methods all stay in a qualitative visual inspection judgment stage and lack a quantitative and objective matching evaluation method for the droplet size test under the condition of ice formation of large supercooled water droplets.
In summary, in the prior art, when the matching evaluation is performed on the droplet size test under the condition of the icing of the supercooled large droplets, a qualitative visual judgment method is adopted, and a quantitative and objective method is not formed at home and abroad aiming at the problem of the matching evaluation of the droplet size test under the condition of the icing of the supercooled large droplets.
Disclosure of Invention
The invention aims to provide an optimal matching evaluation method for the icing condition of supercooled large water drops, and aims to solve the technical problem that only qualitative evaluation can be performed in the prior art.
The invention provides an optimal matching evaluation method for an icing condition of supercooled large water drops, which comprises the following steps of:
step S10: selecting a plurality of characteristic liquid drop points from the standard liquid drop distribution;
step S20: dividing the standard droplet distribution into a plurality of characteristic droplet diameter intervals;
step S30: measuring to obtain measured liquid drop distribution, and calculating the accumulated volume fraction of the liquid drop diameter of each characteristic liquid drop point in the measured liquid drop distribution;
step S40: dividing the distribution of the measured liquid drops into a plurality of measured liquid drop diameter intervals;
step S50: calculating a match bias factorF；
Step S60: calculating the optimal threshold value of the matching deviationF c ；
Step S70: according to the matching deviation factorFSum match deviation optimum thresholdF c And determining the matching degree between the measured droplet distribution and the standard droplet distribution.
Alternatively, in step S10, the serial number of the characteristic droplet point is recorded asi，1≤i≤N，NTo have a characteristic drop point of maximum drop diameter, firstiThe diameter of a droplet at a characteristic droplet point is expressed asD v,i Of 1 atiThe cumulative volume fraction of each characteristic drop point is expressed asV c,i The cumulative volume fraction of the characteristic drop point includes at least four values of 0.1, 0.5, 0.9 and 1.
Alternatively, in step S20, the serial number of the characteristic droplet diameter section is described asjNumber of characteristic droplet diameter intervaljThe serial number of the characteristic liquid drop point with larger liquid drop diameter in the characteristic liquid drop diameter intervaliCorrespond tojThe range of droplet diameters for each characteristic droplet diameter interval is denoted as ΔD v，j Of 1 atjThe volume fraction of the characteristic droplet diameter interval is expressed asV j (ii) a When in usejIn the case of =1, the droplet diameter range Δ of the first characteristic droplet diameter intervalD v，1 =[D 0 ～D v，1 ]Volume fraction of the first characteristic droplet diameter intervalV 1 =V c,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range delta of characteristic droplet diameter intervalD v，j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of characteristic droplet diameter intervalV j =V c，j -V c，j-1 (ii) a Wherein,D 0 is the minimum droplet diameter in a standard droplet distribution.
Optionally, the step S60 includes the following steps:
step S61: constructing an optimal matching cumulative distribution functionV c,opt (D)：
Figure 913828DEST_PATH_IMAGE001
Wherein the standard droplet distribution is formed by a standard first spray and a standard second spray together,Das a function of the diameter of the droplets,βis the ratio of the optimal total droplet volume content of the standard first spray and the standard second spray,MVD 1 is the optimal volume median diameter for the standard first spray,q 1 the optimal droplet distribution width control factor for a standard first spray,MVD 2 is the optimal volume median diameter for the standard second spray,q 2 the optimal droplet distribution width control factor for the standard second spray,Ran empirical cumulative distribution function for a typical artificial spray;
step S62: first, thejThe optimal matching drop volume fraction of the characteristic drop diameter interval is expressed asV opt，j When is coming into contact withj(ii) the optimally matched drop volume fraction of the first characteristic drop diameter interval when =1V opt，1 =V c,opt (D v，1 )(ii) a When the content is less than or equal to 2j≤NWhen it comes tojOptimally matched drop volume fraction for individual characteristic drop diameter intervalsV opt，j =V c,opt (D v，j )-V c,opt (D v，j-1 )(ii) a Interval of standard distribution [ 2 ]D v，N ～∞]Of the optimum matching drop volume fractionV opt，N+1 =1-V c,opt (D v，N )；
Step S63: calculating the optimal threshold value of the matching deviation according to the following formulaF c ：
Figure 541119DEST_PATH_IMAGE002
Wherein, V*The volume fraction is limited for the maximum droplet diameter,f j is as followsjThe volume fraction relative deviation weighting coefficient of the liquid drop in the characteristic liquid drop diameter interval,f N+1 is a standard distribution intervalD v，N ～∞]The drop volume fraction relative deviation weighting factor.
Alternatively, in step S30, the serial number of the measured droplet point is recorded ask，1≤k≤M，MTo measure the drop point with the largest drop diameter, firstkThe droplet diameter of each measured droplet point is expressed asD vm,k Of 1 atkThe cumulative volume fraction of each measured drop point is expressed asV cm,k Of 1 atiThe cumulative volume fraction of the droplet diameter at each characteristic droplet point in the measured droplet distribution is expressed asV cms,i ，V cms,i Is calculated as follows:
when in useD v,i <D vm,1 When the temperature of the water is higher than the set temperature,
Figure 580619DEST_PATH_IMAGE003
；
when in useD vm,k <D v,i <D vm,k+1 And 1 is less than or equal tokWhen the temperature of the water is higher than the set temperature,
Figure 445807DEST_PATH_IMAGE004
；
when in useD vm,M <D v,i When the temperature of the water is higher than the set temperature,V cms,i =1。
optionally, in the step S40, when the stepNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoNThe number of each measured droplet diameter interval is also recorded asjMeasuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 。
Optionally in the second placeNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 132134DEST_PATH_IMAGE005
wherein,f j is as followsjCharacteristic droplet diameter intervalThe drop volume fraction relative deviation weighting factor.
Optionally, in the step S40, when the stepNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoN+1The number of the measured liquid drop diameter interval is recorded asjAndN+1measuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atN +1The range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，N+1 Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j Of 1 atN+1The volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，N+1 (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 (ii) a First, theN+1Droplet diameter range Delta of interval of measured droplet diameterD vm，N+1 =[D v,N ～D vm,M ]Of 1 atN+1Volume fraction of interval of measured droplet diameterV ms，N+1 =1-V cms，N 。
Optionally in the second placeNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMIndividual surveyDroplet diameter of measuring droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 715562DEST_PATH_IMAGE006
。
optionally, in step S70, the matching degree between the measured droplet distribution and the standard droplet distribution is divided into four grades of good, medium and poor matching deviation factors corresponding to the four grades of good, medium and poorFSum match deviation optimum thresholdF c The relationship of (a) to (b) is as follows:
and (3) excellent: 0 is less than or equal toF≤F c ；
Good:F c <F≤20%；
the method comprises the following steps: 20 percent of<F≤40%；
Difference: 40 percent of<F。
Compared with the prior art, the invention has the technical effects that:
1. in the invention, the optimal threshold value of the matching deviation is calculated besides the matching deviation factor, so that the quality of the matching effect can be further measured, and the evaluation is more objective;
2. in the invention, an optimal matching cumulative distribution function is constructed and used for expressing the bimodal distribution characteristics of the standard liquid drop distribution, so that the calculated optimal threshold value of the matching deviation can objectively evaluate the matching degree between the measured liquid drop distribution and the standard liquid drop distribution;
3. in the invention, the matching degree between the measured droplet distribution and the standard droplet distribution is evaluated through the matching deviation factor and the matching deviation optimal threshold, and the method is a specific quantitative method, not a qualitative visual judgment method in the prior art, and can be objectively evaluated;
4. in the invention, the standard liquid drop distribution is divided into a plurality of characteristic liquid drop diameter intervals, and the measured liquid drop distribution is divided into a plurality of measured liquid drop diameter intervals, so that the partition matching degree of the measured liquid drop distribution and the standard liquid drop distribution between the intervals is calculated by considering the matching deviation factor, and the overall matching degree between the measured liquid drop distribution and the standard liquid drop distribution can be objectively evaluated.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments of the present invention or in the description of the prior art will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a graph of four typical standard drop profiles;
FIG. 2 is a flowchart of a method for matching and evaluating icing conditions of supercooled large water droplets according to an embodiment of the present invention;
FIG. 3 is a simplified diagram of a standard droplet distribution in accordance with a first embodiment of the present invention;
FIG. 4 is a measured droplet profile according to a first embodiment of the present invention;
FIG. 5 is a combination diagram of the first embodiment of the present invention;
FIG. 6 is a first schematic diagram illustrating the division of the measured droplet diameter intervals of the measured droplet distribution map according to the first embodiment of the present invention;
FIG. 7 is a combination diagram of the first embodiment of the present invention;
FIG. 8 is a second schematic view of the division of the measured droplet diameter interval of the measured droplet distribution map in the first embodiment of the present invention;
FIG. 9 is a diagram illustrating an optimal matching evaluation method for the icing condition of supercooled large water droplets according to a second embodiment of the present invention;
FIG. 10 is a standard droplet distribution diagram in this experimental example;
FIG. 11 is a measured droplet distribution diagram in this experimental example;
FIG. 12 is a combination view of the experimental example in which the measurement droplet point is hidden.
Detailed Description
Aspects of the present invention will be described more fully hereinafter with reference to the accompanying drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to any specific structure or function presented throughout this disclosure. Rather, these aspects are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art. Based on the teachings herein one skilled in the art should appreciate that the scope of the present invention is intended to encompass any aspect disclosed herein, whether alone or in combination with any other aspect of the invention to accomplish any aspect disclosed herein. For example, it may be implemented using any number of the apparatus or performing methods set forth herein. In addition, the scope of the present invention is intended to cover apparatuses or methods implemented with other structure, functionality, or structure and functionality in addition to the various aspects of the invention set forth herein. It is to be understood that any aspect disclosed herein may be embodied by one or more elements of a claim.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. The terms "comprises," "comprising," and the like, as used herein, specify the presence of stated features, steps, operations, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, or components.
All terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art unless otherwise defined. It is noted that the terms used herein should be interpreted as having a meaning that is consistent with the context of this specification and should not be interpreted in an idealized or overly formal sense.
The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments.
Example one
As shown in fig. 2, an embodiment of the present invention provides a matching evaluation method for icing conditions of supercooled large water droplets, which includes the following steps:
step S10: selecting a plurality of characteristic liquid drop points from the standard liquid drop distribution;
as shown in the simplified diagram of the standard drop distribution of fig. 3, the solid line curve in fig. 3 represents one of the four typical standard drop distribution diagrams in fig. 1; the abscissa in fig. 3 is the droplet diameter and the ordinate is the cumulative volume fraction.
The circle center in fig. 3 is the selected characteristic drop point.
Step S20: dividing the standard droplet distribution into a plurality of characteristic droplet diameter intervals;
with continued reference to fig. 3, the regions between the imaginary straight lines in fig. 3, each passing through a characteristic drop point and perpendicular to the horizontal axis representing the drop diameter, and the region between the leftmost imaginary straight line and the real straight line, collectively form a plurality of characteristic drop diameter intervals. The star point in fig. 3 represents the drop point with the smallest drop diameter in the standard drop distribution, through which the solid straight line in fig. 3 passes, and is perpendicular to the horizontal axis representing the drop diameter.
It should be noted that the characteristic droplet diameter section is substantially a closed figure, and the closed figure is defined by a standard droplet distribution curve, an imaginary straight line passing through the characteristic droplet point, and a horizontal axis, or defined by a standard droplet distribution curve, a real straight line passing through a star point, and a horizontal axis. Hereinafter, the measurement of the droplet diameter interval is similarly expressed, and the details are not described below.
Step S30: measuring to obtain measured liquid drop distribution, and calculating the accumulated volume fraction of the liquid drop diameter of each characteristic liquid drop point in the measured liquid drop distribution;
as shown in fig. 4, which is a distribution diagram of the measured droplet, the abscissa in fig. 4 is the diameter of the droplet, and the ordinate is the cumulative volume fraction, and the solid circle in fig. 4 is the measured droplet point.
The simplified graph of the standard drop distribution in fig. 3 and the measured drop distribution graph in fig. 4 are each placed in a graph to form a combined graph, as shown in fig. 5 for combined graph one, with the abscissa of fig. 5 being the drop diameter and the ordinate being the cumulative volume fraction.
Similar to fig. 3, each characteristic droplet point is taken as a virtual straight line along a direction perpendicular to the horizontal axis, an intersection point is formed between the virtual straight line and a virtual curve representing the distribution diagram of the measured droplet, the intersection point is represented by a hollow frame, the abscissa of the hollow frame is the droplet diameter of the characteristic droplet point, and the ordinate of the hollow frame is the accumulated volume fraction of the droplet diameter of the characteristic droplet point in the distribution of the measured droplet.
Step S40: dividing the distribution of the measured liquid drops into a plurality of measured liquid drop diameter intervals;
for convenience of description, the measured droplet distribution diagram in fig. 5 is extracted, as shown in fig. 6, a first measured droplet diameter interval of the measured droplet distribution diagram is divided into schematic diagrams, in fig. 6, the abscissa is the droplet diameter, the ordinate is the cumulative volume fraction, and a plurality of measured droplet diameter intervals are formed by a region between imaginary straight lines, a region between an imaginary straight line at the leftmost side and a real straight line at the leftmost side, and a region between an imaginary straight line at the rightmost side and a real straight line at the rightmost side, wherein each imaginary straight line passes through a characteristic droplet point, also passes through a hollow frame in fig. 5 and 6, and is perpendicular to a horizontal axis representing the droplet diameter, the real straight line at the rightmost side in fig. 6 passes through the measured droplet point having the largest measured droplet diameter, and the real straight line at the leftmost side in fig. 6 passes through a star point.
Step S50: a matching deviation factor is calculated.
After the matching deviation factor is calculated, the matching result of the icing condition of the supercooled large water drops can be evaluated, and specifically, the matching result of the standard liquid drop distribution and the measured liquid drop distribution can be obtained; the larger the matching deviation factor, the worse the matching degree of the standard droplet distribution and the measured droplet distribution, and the smaller the matching deviation factor, the better the matching degree of the standard droplet distribution and the measured droplet distribution.
In the invention, the matching degree between the measured droplet distribution and the standard droplet distribution is evaluated by matching deviation factors, and the method is a specific quantitative method, but not a qualitative visual judgment method in the prior art.
Further, in step S10, the serial number of the characteristic droplet point is described asi，1≤i≤N，NTo have a characteristic drop point of maximum drop diameter, firstiThe diameter of a droplet at a characteristic droplet point is expressed asD v,i Of 1 atiThe cumulative volume fraction of each characteristic drop point is expressed asV c,i The cumulative volume fraction of the characteristic drop point includes at least four values of 0.1, 0.5, 0.9 and 1.
As in fig. 3, the1The diameter of a droplet at a characteristic droplet point is expressed asD v,1 Of 1 at1The cumulative volume fraction of each characteristic drop point is expressed asV c,1 Of 1 atNThe diameter of a droplet at a characteristic droplet point is expressed asD v,N Of 1 atNThe cumulative volume fraction of each characteristic drop point is expressed asV c,N . First, theNThe diameter of the droplet of one characteristic droplet point is larger than that of the other characteristic droplet pointsNCumulative volume fraction of individual characteristic drop pointV c,N =1。
Further, in step S20, the serial numbers of the characteristic droplet diameter sections are described asjNumber of characteristic droplet diameter intervaljThe serial number of the characteristic liquid drop point with larger liquid drop diameter in the characteristic liquid drop diameter intervaliCorrespond tojThe range of droplet diameters for each characteristic droplet diameter interval is denoted as ΔD v，j Of 1 atjThe volume fraction of the characteristic droplet diameter interval is expressed asV j (ii) a When in usejIn the case of =1, the droplet diameter range Δ of the first characteristic droplet diameter intervalD v，1 =[D 0 ～D v，1 ]Volume fraction of the first characteristic droplet diameter intervalV 1 =V c,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range delta of characteristic droplet diameter intervalD v，j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of characteristic droplet diameter intervalV j =V c，j -V c，j-1 (ii) a Wherein,D 0 is the minimum droplet diameter in a standard droplet distribution.
With continued reference to FIG. 3, whenjWhere =1, the first characteristic droplet diameter interval is defined by a real straight line and a leftmost imaginary straight line, and in this case,i=j(ii) a When the content is less than or equal to 2j≤NWhen it comes tojVirtual straight lines at two sides of the characteristic liquid drop diameter interval respectively pass through the firsti-1Characteristic drop point andithe number of the characteristic liquid drop diameter intervaljThe serial number of the characteristic drop point with larger drop diameter in the interval of the characteristic drop diameter is set asiCorrespondingly, at this time,i=j。
further, in step S30, the number of the measured droplet point is described ask，1≤k≤M，MTo measure the drop point with the largest drop diameter, firstkThe droplet diameter of each measured droplet point is expressed asD vm,k Of 1 atkThe cumulative volume fraction of each measured drop point is expressed asV cm,k Of 1 atiThe cumulative volume fraction of the droplet diameter at each characteristic droplet point in the measured droplet distribution is expressed asV cms,i ，V cms,i Is calculated as follows:
when in useD v,i <D vm,1 When the temperature of the water is higher than the set temperature,
Figure 496437DEST_PATH_IMAGE007
；
when in useD vm,k <D v,i <D vm,k+1 And 1 is less than or equal tokWhen the temperature of the water is higher than the set temperature,
Figure 58874DEST_PATH_IMAGE004
；
when in useD vm,M <D v,i When the temperature of the water is higher than the set temperature,V cms,i =1。
first, theMThe droplet diameter of one measuring droplet point is larger than that of the other measuring droplet points,D vm,M is the firstMDroplet diameter of individual characteristic droplet point, cumulative volume fraction of maximum droplet diameter measured droplet pointV cm,M And a firstNCumulative volume fraction of individual characteristic drop pointV c,N The same applies to 1.
Further, in the step S40, when the first step is performedNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoNThe number of each measured droplet diameter interval is also recorded asjMeasuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 ；
FIG. 7 shows a second combination diagram, in which the abscissa of FIG. 7 is the droplet diameter and the ordinate is the cumulative volume fraction, and the second combination diagram of FIG. 7NDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M Extracting the measured droplet distribution diagram in fig. 7, as shown in fig. 8, dividing the measured droplet diameter interval of the measured droplet distribution diagram into two schematic diagrams, in fig. 8, the abscissa is the dropletThe diameters, the ordinate of which is the cumulative volume fraction, the area between the imaginary straight lines, each passing through a characteristic drop point, also the open box in fig. 7 and 8, and the area between the leftmost imaginary straight line and the leftmost real straight line, which passes through a star point, in fig. 7 and 8, together form a plurality of measured drop diameter intervals.
Further, in the second placeNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 368632DEST_PATH_IMAGE005
wherein,f j is as followsjThe volume fraction of droplets within a characteristic droplet diameter interval is weighted by the relative deviation.
Further, when it comes toNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoN+1The number of the measured liquid drop diameter interval is recorded asjAndN+1measuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atN+1The range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，N+1 Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j Of 1 atN+1The volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，N+1 (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 (ii) a First, theN+1Droplet diameter range Delta of interval of measured droplet diameterD vm，N+1 =[D v,N ～D vm,M ]Of 1 atN+1Volume fraction of interval of measured droplet diameterV ms，N+1 =1-V cms，N 。
FIG. 5 is a view showing a combination of the first and second drawings in FIG. 5NDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M And the measured droplet distribution diagram in fig. 5 is extracted, as shown in fig. 6, the measured droplet diameter intervals of the measured droplet distribution diagram are divided into a schematic diagram i, in fig. 6, a plurality of measured droplet diameter intervals are jointly formed by a region between the imaginary straight lines, a region between the leftmost imaginary straight line and the leftmost real straight line, and a region between the rightmost imaginary straight line and the rightmost real straight line, wherein each imaginary straight line passes through a characteristic droplet point, also passes through a hollow frame in fig. 5 and 6, and is perpendicular to a horizontal axis representing the droplet diameter, the leftmost real straight line in fig. 5 and 6 passes through a star point, and the rightmost real straight line in fig. 5 and 6 passes through a second real straight lineMOne measures the drop point.
Further, in the second placeNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 314723DEST_PATH_IMAGE008
wherein, V*The volume fraction is limited for the maximum droplet diameter,f j is as followsjThe volume fraction relative deviation weighting coefficient of the liquid drop in the characteristic liquid drop diameter interval,f N+1 is a standard distribution intervalD v，N ～∞]The drop volume fraction relative deviation weighting factor.
In the invention, the standard liquid drop distribution is divided into a plurality of characteristic liquid drop diameter intervals, and the measured liquid drop distribution is divided into a plurality of measured liquid drop diameter intervals, so that the partition matching degree of the measured liquid drop distribution and the standard liquid drop distribution between the intervals is calculated by considering the matching deviation factor, and the overall matching degree between the measured liquid drop distribution and the standard liquid drop distribution can be objectively evaluated.
Example two
As shown in fig. 9, the second embodiment of the present invention provides an optimal matching evaluation method for icing conditions of supercooled large water droplets, which includes the following steps:
step S10: selecting a plurality of characteristic liquid drop points from the standard liquid drop distribution;
step S20: dividing the standard droplet distribution into a plurality of characteristic droplet diameter intervals;
step S30: measuring to obtain measured liquid drop distribution, and calculating the accumulated volume fraction of the liquid drop diameter of each characteristic liquid drop point in the measured liquid drop distribution;
step S40: dividing the distribution of the measured liquid drops into a plurality of measured liquid drop diameter intervals;
step S50: calculating a match bias factorF；
Step S60: calculating the optimal threshold value of the matching deviationF c ；
Step S70: according to the matching deviation factorFSum match deviation optimum thresholdF c Determining the matching degree between the measured droplet distribution and the standard droplet distribution;
wherein, step S10～Step S50 is the same as in embodiment one, and therefore, regarding stepS10～The detailed steps of step S50 are not described herein.
Further, step S60 includes the following steps:
step S61: constructing an optimal matching cumulative distribution functionV c,opt (D)：
Figure 151485DEST_PATH_IMAGE001
Wherein the standard droplet distribution is formed by a standard first spray and a standard second spray together,Das a function of the diameter of the droplets,βis the ratio of the optimal total droplet volume content of the standard first spray and the standard second spray,MVD 1 is the optimal volume median diameter for the standard first spray,q 1 the optimal droplet distribution width control factor for a standard first spray,MVD 2 is the optimal volume median diameter for the standard second spray,q 2 the optimal droplet distribution width control factor for the standard second spray,Ran empirical cumulative distribution function for a typical artificial spray;
the main purpose of this step is to construct the above formula considering that the droplet size distribution of the supercooled large droplet icing condition given in appendix O in section 25 of airworthiness regulations has a significant bimodal distribution characteristic, requiring two artificial sprays of different droplet distribution morphologies to match.
Step S62: first, thejThe optimal matching drop volume fraction of the characteristic drop diameter interval is expressed asV opt，j When is coming into contact withj(ii) the optimally matched drop volume fraction of the first characteristic drop diameter interval when =1V opt，1 =V c,opt (D v，1 )(ii) a When the content is less than or equal to 2j≤NWhen it comes tojOptimally matched drop volume fraction for individual characteristic drop diameter intervalsV opt，j =V c,opt (D v，j )-V c,opt (D v，j-1 )(ii) a Interval of standard distribution [ 2 ]D v，N ～∞]Of the optimum matching drop volume fractionV opt，N+1 =1-V c,opt (D v，N )；
Step S63: calculating the optimal threshold value of the matching deviation according to the following formulaF c ：
Figure 771691DEST_PATH_IMAGE002
。
In step S70, the matching degree between the measured droplet distribution and the standard droplet distribution can be divided into a plurality of grades, such as four grades of excellent, good, medium and poor, and the matching deviation factors corresponding to the gradesFSum match deviation optimum thresholdF c The relationship of (a) to (b) is as follows:
and (3) excellent: 0 is less than or equal toF≤F c ；
Good:F c <F≤20%；
the method comprises the following steps: 20 percent of<F≤40%；
Difference: 40 percent of<F。
The second embodiment of the invention calculates the optimal threshold value of the matching deviation on the basis of the first embodimentF c And the quality of the matching effect can be further measured.
Examples of the experiments
The present invention is further illustrated below by an experimental example, in which the standard droplet distribution is that of a fine freezing rain having an MVD <40 μm.
As shown in FIG. 10, the abscissa is the droplet diameter and the ordinate is the cumulative volume fraction at MVD<Selection on standard droplet distribution of 40 μm frozen drizzleN=10 characteristic droplet points, the droplet diameters of the 10 characteristic droplet points being 15, 36, 106, 197, 247, 276, 309, 354, 390, 473 in order from small to large, and the cumulative volume fractions of the 10 characteristic droplet points being 0.1, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 0.95, 0.975, 1 in order from small to large.
The characteristic droplet diameter interval is [1 ] in sequence～15]、[15～36]、[36～106]、[106～197]、[197～247]、[247～276]、[276～309]、[309～354]、[354～390]、[390～473]
The volume fraction of the characteristic droplet diameter interval is 0.1, 0.2, 0.1, 0.05, 0.025 and 0.025 in sequence.
The drop distribution was measured as shown in FIG. 11, with the abscissa of the drop diameter and the ordinate of the cumulative volume fraction, totalingMAnd =60 measured droplet points, wherein the droplet diameter of the measured droplet points is in the range of 1-2500 μm, and the cumulative volume fractions of the droplet diameters of the characteristic droplet points in the measured droplet distribution are 0.09, 0.288, 0.494, 0.709, 0.786, 0.822, 0.857, 0.896, 0.92 and 0.961 in sequence.
First, theNDroplet diameter of characteristic droplet pointD v,N =473 μm, thMMeasuring the droplet diameter of the droplet pointD vm,M =2500 μm, it is therefore necessary to divide the measured droplet distribution intoN+1The volume fractions of the measured droplet diameter intervals are calculated to be 0.09, 0.198, 0.207, 0.215, 0.077, 0.036, 0.035, 0.039, 0.024, 0.04 and 0.039 in sequence.
f j Sequentially takes the values of 0.1, 0.2, 0.1, 0.05, 0.025 and 0.025,f N+1 the value is 0.025.
Calculated matching deviation factorF=14.93%。
A combined plot hiding measured droplet points is shown in fig. 12 with droplet diameter on the abscissa and cumulative volume fraction on the ordinate.
Further, the following experiment can be continued:
the empirical cumulative distribution function of a typical artificial spray is chosen asRosin-RammlerDistribution function, then, optimally matching cumulative distribution functionV c,opt (D)Then it is:
Figure 442232DEST_PATH_IMAGE009
wherein,β=0.58，MVD 1 =24μm，q 1 =1.8，MVD 2 =192μm，q 2 =1.9；
through calculation, the optimal threshold value of the deviation is matchedF c =3.76%
Due to the fact thatFc<F<20%, so the degree of match between the measured droplet distribution and the standard droplet distribution is a "good" rating.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (10)

1. An optimal matching evaluation method for an icing condition of supercooled large water drops is characterized by comprising the following steps:
step S10: selecting a plurality of characteristic liquid drop points from the standard liquid drop distribution;
step S20: dividing the standard droplet distribution into a plurality of characteristic droplet diameter intervals;
step S30: measuring to obtain measured liquid drop distribution, and calculating the accumulated volume fraction of the liquid drop diameter of each characteristic liquid drop point in the measured liquid drop distribution;
step S40: dividing the distribution of the measured liquid drops into a plurality of measured liquid drop diameter intervals;
step S50: calculating a match bias factorF；
Step S60: calculating the optimal threshold value of the matching deviationF c ；
Step S70: according to the matching deviation factorFSum match deviation optimum thresholdF c And determining the matching degree between the measured droplet distribution and the standard droplet distribution.
2. The method for optimizing matching evaluation of supercooled large droplet icing condition according to claim 1, wherein said step S10 is performedThe serial number of the characteristic drop point is recorded asi，1≤i≤N，NTo have a characteristic drop point of maximum drop diameter, firstiThe diameter of a droplet at a characteristic droplet point is expressed asD v,i Of 1 atiThe cumulative volume fraction of each characteristic drop point is expressed asV c,i The cumulative volume fraction of the characteristic drop point includes at least four values of 0.1, 0.5, 0.9 and 1.
3. The method as claimed in claim 2, wherein in step S20, the serial number of the characteristic droplet diameter interval is recorded asjNumber of characteristic droplet diameter intervaljThe serial number of the characteristic liquid drop point with larger liquid drop diameter in the characteristic liquid drop diameter intervaliCorrespond tojThe range of droplet diameters for each characteristic droplet diameter interval is denoted as ΔD v，j Of 1 atjThe volume fraction of the characteristic droplet diameter interval is expressed asV j (ii) a When in usejIn the case of =1, the droplet diameter range Δ of the first characteristic droplet diameter intervalD v，1 =[D 0 ～D v，1 ]Volume fraction of the first characteristic droplet diameter intervalV 1 =V c,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range delta of characteristic droplet diameter intervalD v，j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of characteristic droplet diameter intervalV j =V c，j -V c，j-1 (ii) a Wherein,D 0 is the minimum droplet diameter in a standard droplet distribution.
4. The method for optimizing matching evaluation of supercooled large droplet icing condition according to claim 3, wherein said step S60 includes the steps of:
step S61: constructing an optimal matching cumulative distribution functionV c,opt (D)：
Figure 916954DEST_PATH_IMAGE001
Wherein the standard droplet distribution is formed by a standard first spray and a standard second spray together,Das a function of the diameter of the droplets,βis the ratio of the optimal total droplet volume content of the standard first spray and the standard second spray,MVD 1 is the optimal volume median diameter for the standard first spray,q 1 the optimal droplet distribution width control factor for a standard first spray,MVD 2 is the optimal volume median diameter for the standard second spray,q 2 the optimal droplet distribution width control factor for the standard second spray,Ran empirical cumulative distribution function for a typical artificial spray;
step S62: first, thejThe optimal matching drop volume fraction of the characteristic drop diameter interval is expressed asV opt，j When is coming into contact withj(ii) the optimally matched drop volume fraction of the first characteristic drop diameter interval when =1V opt，1 =V c,opt (D v，1 )(ii) a When the content is less than or equal to 2j≤NWhen it comes tojOptimally matched drop volume fraction for individual characteristic drop diameter intervalsV opt，j =V c,opt (D v，j )-V c,opt (D v，j-1 )(ii) a Interval of standard distribution [ 2 ]D v，N ～∞]Of the optimum matching drop volume fractionV opt，N+1 =1-V c,opt (D v，N )；
Step S63: calculating the optimal threshold value of the matching deviation according to the following formulaF c ：
Figure 305210DEST_PATH_IMAGE002
Wherein, V*The volume fraction is limited for the maximum droplet diameter,f j is as followsjThe volume fraction relative deviation weighting coefficient of the liquid drop in the characteristic liquid drop diameter interval,f N+1 is a standard distribution intervalD v，N ～∞]The drop volume fraction relative deviation weighting factor.
5. The method as claimed in claim 4, wherein in step S30, the serial number of the measured dropping point is recorded ask，1≤k≤M，MTo measure the drop point with the largest drop diameter, firstkThe droplet diameter of each measured droplet point is expressed asD vm,k Of 1 atkThe cumulative volume fraction of each measured drop point is expressed asV cm,k Of 1 atiThe cumulative volume fraction of the droplet diameter at each characteristic droplet point in the measured droplet distribution is expressed asV cms,i ，V cms,i Is calculated as follows:
when in useD v,i <D vm,1 When the temperature of the water is higher than the set temperature,
Figure 265076DEST_PATH_IMAGE003
；
when in useD vm,k <D v,i <D vm,k+1 And 1 is less than or equal tokWhen the temperature of the water is higher than the set temperature,
Figure 711712DEST_PATH_IMAGE004
；
when in useD vm,M <D v,i When the temperature of the water is higher than the set temperature,V cms,i =1。
6. the method of claim 5, wherein in step S40, when the step S40 is executedNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoNThe number of each measured droplet diameter interval is also recorded asjMeasuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 。
7. The method of claim 6, wherein the evaluation of the optimum matching of supercooled large droplet icing conditions is performed in the first stepNDroplet diameter of characteristic droplet pointD v,N Is greater than or equal toMMeasuring the droplet diameter of the droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 84924DEST_PATH_IMAGE005
wherein,f j is as followsjThe volume fraction of droplets within a characteristic droplet diameter interval is weighted by the relative deviation.
8. The method of claim 5, wherein in step S40, when the step S40 is executedNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M While dividing the measured droplet distribution intoN+1The number of the measured liquid drop diameter interval is recorded asjAndN+ 1measuring the number of intervals of droplet diameterjAnd the serial number of the characteristic liquid drop point with larger liquid drop diameter in the interval of measuring the liquid drop diameteriCorrespond tojThe range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，j Of 1 atN+1The range of droplet diameters for each interval of measured droplet diameters is denoted ΔD vm，N+1 Of 1 atjThe volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，j Of 1 atN +1The volume fraction of the interval of the diameter of each measured droplet is expressed asV ms，N+1 (ii) a When in usej(ii) droplet diameter range Delta of the first measured droplet diameter interval when =1D vm,1 =[D 0 ～D v，1 ]First measuring the volume fraction of the droplet diameter intervalV ms,1 =V cms,1 (ii) a When the content is less than or equal to 2j≤NWhen it comes tojDroplet diameter range Delta of interval of measured droplet diameterD vm,j =[D v，j-1 ～D v，j ]Of 1 atjVolume fraction of interval of measured droplet diameterV ms，j =V cms，j -V cms，j-1 (ii) a First, theN+1Droplet diameter range Delta of interval of measured droplet diameterD vm，N+1 =[D v,N ～D vm,M ]Of 1 atN+1Volume fraction of interval of measured droplet diameterV ms，N+1 =1-V cms，N 。
9. A process as claimed in claim 8The optimal matching evaluation method for the icing condition of supercooled large water drops is characterized in thatNDroplet diameter of characteristic droplet pointD v,N Is less than that ofMMeasuring the droplet diameter of the droplet pointD vm,M In the case of (3), the step S50 calculates a matching deviation factor F by the following formula:
Figure 11292DEST_PATH_IMAGE006
。
10. the method of claim 1, wherein in step S70, the matching degree between the measured droplet distribution and the standard droplet distribution is divided into four grades of good, medium and poor, and the matching deviation factors corresponding to the four grades of good, medium and poorFSum match deviation optimum thresholdF c The relationship of (a) to (b) is as follows:
and (3) excellent: 0 is less than or equal toF≤F c ；
Good:F c <F≤20%；
the method comprises the following steps: 20 percent of<F≤40%；
Difference: 40 percent of<F。
CN202110167936.3A 2021-02-07 2021-02-07 Optimized matching evaluation method for icing condition of large supercooled water drops Active CN112507572B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110167936.3A CN112507572B (en) 2021-02-07 2021-02-07 Optimized matching evaluation method for icing condition of large supercooled water drops

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110167936.3A CN112507572B (en) 2021-02-07 2021-02-07 Optimized matching evaluation method for icing condition of large supercooled water drops

Publications (2)

Publication Number Publication Date
CN112507572A true CN112507572A (en) 2021-03-16
CN112507572B CN112507572B (en) 2021-04-30

Family

ID=74952748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110167936.3A Active CN112507572B (en) 2021-02-07 2021-02-07 Optimized matching evaluation method for icing condition of large supercooled water drops

Country Status (1)

Country Link
CN (1) CN112507572B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113252281A (en) * 2021-06-02 2021-08-13 中国空气动力研究与发展中心低速空气动力研究所 Method for reconstructing size distribution of icing cloud droplets
CN118603488A (en) * 2024-08-01 2024-09-06 中国空气动力研究与发展中心低速空气动力研究所 Water drop impact characteristic capturing method based on water sensing paper and coloring agent

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120274938A1 (en) * 2011-04-29 2012-11-01 Rosemount Aerospace Inc. Apparatus and method for detecting aircraft icing conditions
CN103910062A (en) * 2013-01-02 2014-07-09 波音公司 Automated Water Drop Measurement and Ice Detection System
CN104268399A (en) * 2014-09-24 2015-01-07 空气动力学国家重点实验室 Computing method of model parameters in icing wind tunnel experiment under supercooled large droplet condition
CN105992941A (en) * 2013-10-24 2016-10-05 密执安州立大学董事会 Ice and supercooled water detection system
CN107106925A (en) * 2014-10-18 2017-08-29 艾拉索尔戴纳米克斯有限公司 Lasting supersaturation for the condensation growth of particle
CN111563354A (en) * 2020-06-10 2020-08-21 中国空气动力研究与发展中心低速空气动力研究所 Icing wind tunnel test similarity conversion method based on numerical simulation
CN111738481A (en) * 2020-04-01 2020-10-02 南京航空航天大学 Airplane icing meteorological parameter MVD prediction method based on BP neural network

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120274938A1 (en) * 2011-04-29 2012-11-01 Rosemount Aerospace Inc. Apparatus and method for detecting aircraft icing conditions
CN103910062A (en) * 2013-01-02 2014-07-09 波音公司 Automated Water Drop Measurement and Ice Detection System
CN105992941A (en) * 2013-10-24 2016-10-05 密执安州立大学董事会 Ice and supercooled water detection system
CN104268399A (en) * 2014-09-24 2015-01-07 空气动力学国家重点实验室 Computing method of model parameters in icing wind tunnel experiment under supercooled large droplet condition
CN107106925A (en) * 2014-10-18 2017-08-29 艾拉索尔戴纳米克斯有限公司 Lasting supersaturation for the condensation growth of particle
CN111738481A (en) * 2020-04-01 2020-10-02 南京航空航天大学 Airplane icing meteorological parameter MVD prediction method based on BP neural network
CN111563354A (en) * 2020-06-10 2020-08-21 中国空气动力研究与发展中心低速空气动力研究所 Icing wind tunnel test similarity conversion method based on numerical simulation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JOSHUA D.BLAKE ET AL.: "Effects of Surface Characteristics and Droplet Diameter on the Freezing of Supercooled Water Droplets Impacting a Cooled Substrate", 《TRANSACTIONS OF JAPANESE SOCIETY FOR MEDICAL AND BIOLOGICAL ENGINEERING》 *
王桥: "过冷大水滴动力学特性的温度影响实验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113252281A (en) * 2021-06-02 2021-08-13 中国空气动力研究与发展中心低速空气动力研究所 Method for reconstructing size distribution of icing cloud droplets
CN113252281B (en) * 2021-06-02 2021-09-21 中国空气动力研究与发展中心低速空气动力研究所 Method for reconstructing size distribution of icing cloud droplets
CN118603488A (en) * 2024-08-01 2024-09-06 中国空气动力研究与发展中心低速空气动力研究所 Water drop impact characteristic capturing method based on water sensing paper and coloring agent

Also Published As

Publication number Publication date
CN112507572B (en) 2021-04-30

Similar Documents

Publication Publication Date Title
CN112507572B (en) Optimized matching evaluation method for icing condition of large supercooled water drops
Korolev et al. A new look at the environmental conditions favorable to secondary ice production
CN109214675B (en) Power transmission line icing risk assessment method
Li et al. Aircraft measurements of the vertical distribution and activation property of aerosol particles over the Loess Plateau in China
Mao et al. Linkage between the Arctic Oscillation and winter extreme precipitation over central-southern China
CN113793021B (en) Three-dimensional power transmission corridor mountain fire monitoring device and method
CN106570780B (en) Power transmission line galloping early warning method based on grey correlation theory
Ren et al. Change in climatic extremes over mainland China based on an integrated extreme climate index
Zheng et al. Precipitation microphysical processes in the inner rainband of tropical cyclone Kajiki (2019) over the South China Sea revealed by polarimetric radar
Takahashi et al. Projected future change of MJO and its extratropical teleconnection in East Asia during the northern winter simulated in IPCC AR4 models
Chenchen et al. Improvement in the forecasting of heavy rainfall over South China in the DSAEF_LTP model by introducing the intensity of the tropical cyclone
Karalis et al. Effects of secondary ice processes on a stratocumulus to cumulus transition during a cold-air outbreak
Leroy et al. HAIC/HIWC field campaigns-Specific findings on ice crystals characteristics in high ice water content cloud regions
CN113252281B (en) Method for reconstructing size distribution of icing cloud droplets
Hagen et al. Effect of canopy leaf distribution on sand transport and abrasion energy
CN108983322A (en) A kind of prediction technique of winter rainfall snowfall
Liu et al. The effects of ice nucleation on the microphysical processes and precipitation for a heavy rainfall event in Beijing
Zhao et al. Reconstruction of conceptual prediction model for the three rainfall patterns in the summer of eastern China under global warming
Kucera et al. Features of the Weather Modification Assesment Project in Southwest Region of Saudi Arabia
Zhu et al. High-rise building group regional fire risk assessment model based on AHP
Makino et al. Evaluation of Lateral Attenuation for Aircraft Takeoff-roll Noise by Multi-point Measurement
Jones How was the MRI model applied
Delobbe et al. Simulation of marine stratocumulus: Effect of precipitation parameterization and sensitivity to droplet number concentration
DE102010013607B4 (en) Method for measuring the energy content of storms
CN118211298A (en) Method and system for evaluating safety level of stand in windy weather

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