CN115416854A - Icing detection device and icing detection method based on temperature measurement - Google Patents

Icing detection device and icing detection method based on temperature measurement Download PDF

Info

Publication number
CN115416854A
CN115416854A CN202211381850.1A CN202211381850A CN115416854A CN 115416854 A CN115416854 A CN 115416854A CN 202211381850 A CN202211381850 A CN 202211381850A CN 115416854 A CN115416854 A CN 115416854A
Authority
CN
China
Prior art keywords
temperature
icing
water
point
detection method
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
CN202211381850.1A
Other languages
Chinese (zh)
Other versions
CN115416854B (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 CN202211381850.1A priority Critical patent/CN115416854B/en
Publication of CN115416854A publication Critical patent/CN115416854A/en
Application granted granted Critical
Publication of CN115416854B publication Critical patent/CN115416854B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64DEQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENTS OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
    • B64D15/00De-icing or preventing icing on exterior surfaces of aircraft
    • B64D15/20Means for detecting icing or initiating de-icing

Abstract

The invention relates to the technical field of deicing prevention and provides an icing detection device and an icing detection method based on temperature measurement. Moreover, by adopting the detection device and the detection method, the cloud and mist parameters can be calculated while the icing thickness is detected: the water content LWC and the water drop median diameter MVD solve the technical problem of difficult cloud parameter detection. The detection device can be made into a universal part, and is equivalent to a novel icing sensor.

Description

Icing detection device and icing detection method based on temperature measurement
Technical Field
The invention relates to the technical field of ice prevention and deicing, in particular to an icing detection device and an icing detection method based on temperature measurement.
Background
Icing is one of the main causes of flight accidents of an aircraft, and icing on the leading edges of wings and empennages of the aircraft can cause increase of wing section resistance, reduction of lift force, reduction of critical angle of attack, and deterioration of maneuverability and stability, thus causing serious flight accidents, and therefore, the icing is widely paid attention and researched by people.
The existing airplane icing detection technology develops various measurement methods such as an optical method, a thermal method, an electrical method, a mechanical method, a waveguide method and the like, and some old icing detection methods such as a visual method, an eye-shielding method, an X-ray method and the like are basically eliminated at present due to the larger defects of the old icing detection methods; the differential pressure method, the electric conduction method and other technologies are not widely applied because of the external interference and the failure to effectively provide the information of the thickness of the ice.
In the current foreign molded products, the most applied are magnetostrictive icing detection systems of American ROSEMOUNT company, flat-film icing sensors of Sweden VibroMeter company, optical fiber icing sensors produced in the UK, and the like.
However, the existing icing sensors are designed for single-point ice measurement, and if the multi-point ice thickness needs to be measured, a solution scheme of installing a large number of icing sensors is needed, so that the scale and the cost of a detection system are rapidly increased, and the difficulty of arrangement and maintenance is brought along. Because the optical fiber sensor is relatively cheap and the system is easy to integrate, most of the current multipoint icing detection sensors adopt an optical fiber scheme. However, the integration and signal processing of large-scale optical fiber detectors require a large number of independent hardware sets for back-end signal acquisition and processing, such as photoelectric conversion, electric signal amplification, multi-way switches and the like, and the integration of such systems has no engineering significance with the analysis point as high as tens of orders or hundreds of orders.
Therefore, how to simplify the icing detection device to detect the accumulated ice in a large area and reduce the system scale and cost is an urgent technical problem to be solved in the field.
Disclosure of Invention
In order to solve the defects in the prior art, the application provides an icing detection device and an icing detection method based on temperature measurement, a temperature sensor array is adopted for temperature detection, and icing thickness and icing cloud and mist parameters are calculated based on detected temperature values.
The invention provides an icing detection method based on temperature measurement, which is characterized in that a temperature sensor array is arranged on the surface of an object to be detected, the position of each temperature sensor corresponds to an acquisition point, and the detection method comprises the following steps:
s10, acquiring the temperature value of each acquisition point according to the time sequence, and obtaining the change trend of the temperature of each acquisition point along with the time;
s20, obtaining the change trend of the temperature at each acquisition point along with time according to the temperature value of each acquisition point; comparing the temperature change trend of each acquisition point, and determining the acquisition point corresponding to the maximum temperature rise value as a standing point position;
s30, selecting temperature values T acquired by at least two temperature sensors near the stagnation position, and substituting the temperature values T into the following formula to calculate the icing thickness H of the stagnation position:
Figure 530571DEST_PATH_IMAGE001
Figure 796467DEST_PATH_IMAGE002
Figure 874145DEST_PATH_IMAGE003
wherein the content of the first and second substances,
Figure 500167DEST_PATH_IMAGE004
for the heat flow of the ice surface at the initial stage of icing,
Figure 448531DEST_PATH_IMAGE005
is the mass of impinging water per unit area, n is the local freezing coefficient,
Figure 150908DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure 715882DEST_PATH_IMAGE007
for water droplet collection rate, LWC is the liquid water content,C w is the specific heat of the water, and the specific heat of the water,
Figure 880016DEST_PATH_IMAGE008
is the incoming flow velocity.
Further, the heat flow of the ice surface at the initial stage of icing
Figure 214045DEST_PATH_IMAGE004
Calculating by numerical simulation:
Figure 821744DEST_PATH_IMAGE009
Figure 139593DEST_PATH_IMAGE010
Figure 107418DEST_PATH_IMAGE011
Figure 561533DEST_PATH_IMAGE012
Figure 340133DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 879699DEST_PATH_IMAGE014
the energy put into it for impacting the drop,
Figure 385635DEST_PATH_IMAGE015
in order to evaporate the energy, the evaporator is provided with a heat exchanger,
Figure 694257DEST_PATH_IMAGE016
in order to freeze the latent heat of phase change,
Figure 643758DEST_PATH_IMAGE017
is waterThe convection heat exchange is carried out in the membrane, T w the temperature of the liquid drops on the wall surface after impact,
Figure 857571DEST_PATH_IMAGE018
in order to achieve the temperature of the incoming flow,V w is the speed of the liquid drops and is,L is the latent heat of solidification in the liquid state, m in order to be the mass transfer coefficient,is a convective heat transfer coefficient, C s And
Figure 183510DEST_PATH_IMAGE019
respectively the mass concentration of the water vapor near the water surface and the mass concentration of the water vapor of the incoming air,L e is the latent heat of evaporation per unit mass, t s And t 1 Respectively the wall temperature and the ambient temperature outside the boundary layer.
Further, in step S30, the liquid water content LWC is calculated at the same time, and the median diameter MVD of the water drop is calculated according to the following relation:
Figure 799168DEST_PATH_IMAGE020
(ii) a Wherein f is the temperature of the object surface to be measured obtained in the process of simulating the flow field of air and water drops by adopting an Euler-Euler method at 0 ℃,
Figure 919571DEST_PATH_IMAGE021
Figure 89521DEST_PATH_IMAGE022
and with
Figure 671681DEST_PATH_IMAGE007
In which water droplet collection coefficient at stagnation point
Figure 954894DEST_PATH_IMAGE007
Is a number of 1, and the number of the main chain is 1,
Figure 511778DEST_PATH_IMAGE022
collected by a wind speed sensor.
Further, in step S30, a local freezing coefficient n is calculated at the same time; when n is more than or equal to 0.95, selecting the position with the minimum temperature change value nearest to the stagnation point as an impact limit; and when n is less than 0.95, selecting the position with the minimum temperature change value nearest to the stagnation point as an impact limit in the preset time at the initial stage of icing.
Further, in step S30, the liquid water content LWC is calculated at the same time; determining the water drop collection coefficient of the position to be measured according to the change relation of the water drop collection coefficient between the stagnation point and the impact limit
Figure 982073DEST_PATH_IMAGE007
And will be
Figure 633503DEST_PATH_IMAGE007
Substituting LWC into the following relational expression to calculate the median diameter MVD of the water drop at the position to be measured:
Figure 974486DEST_PATH_IMAGE020
(ii) a Wherein, the first and the second end of the pipe are connected with each other,when the temperature of an object surface to be measured obtained in the process of simulating the flow field of air and water drops by adopting an Euler-Euler method is 0 ℃,
Figure 702271DEST_PATH_IMAGE021
Figure 174709DEST_PATH_IMAGE022
and
Figure 646142DEST_PATH_IMAGE007
the functional relationship of (a) to (b),
Figure 841631DEST_PATH_IMAGE022
collected by a wind speed sensor.
And further, calculating the time difference of the steep temperature rise and steep temperature fall positions of the acquisition point near the stagnation point on a time scale, taking the time difference as the icing time length, and calculating the maximum icing thickness according to the time length.
The invention provides an icing detection device based on temperature measurement, which is used for executing the icing detection method based on temperature measurement and comprises a temperature sensor array, a temperature acquisition module and a calculation module;
the temperature sensor array is arranged on the surface of an object to be measured or arranged around the circumference of the cylinder;
the acquisition module executes step S10;
the calculation module performs steps S20-S30.
Further, when the temperature sensor array is arranged, the distance between the temperature sensor array and the temperature sensor array needs to satisfy the following requirements: there is no thermal contact between the temperature sensors.
Compared with the prior art, the icing detection device and the icing detection method based on temperature measurement at least have the following beneficial effects:
(1) According to the icing detection device and method, the conventional icing sensor is replaced by combining the temperature sensor array with proper theoretical calculation, so that the size of the detection device can be reduced to a certain extent, the system cost is reduced, the monitoring area is increased, and the system reliability is improved;
(2) According to the freezing detection device and method, the cloud and mist parameters can be obtained while the freezing thickness is calculated, and the complexity of independently erecting a cloud and mist parameter detection system is reduced;
(3) The utility model provides a detection device freezes can make a general part, is equivalent to an integrated detector that freezes, a neotype detector that freezes.
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 schematic flow chart of a method for detecting icing based on temperature measurement according to embodiment 1 of the present application;
FIG. 2 is a schematic view of the arrangement of temperature sensors in embodiment 1 of the present application;
FIG. 3 is a schematic structural diagram of an icing detection device based on temperature measurement according to embodiment 2 of the present application;
fig. 4 is a schematic structural diagram of an icing detection device based on temperature measurement according to embodiment 3 of the present application.
In the figure, 10 is an airfoil, 20 is a temperature sensor array, and 30 is a cylinder.
Detailed Description
The following description provides many different embodiments, or examples, for implementing different features of the invention. The particular examples set forth below are illustrative only and are not intended to be limiting.
Example 1
A freezing detection method based on temperature measurement does not need to be provided with a plurality of freezing sensors, only needs to be provided with the temperature sensors, the temperature sensors can select thermocouples with small volume and small heat capacity, the whole volume is much smaller than that of the freezing sensors, and the system scale and the cost are reduced more favorably.
As shown in fig. 1, a temperature sensor array is arranged on a surface of an object to be detected, as shown in fig. 2, a temperature sensor array 20 is arranged on a surface of a wing 10, and a position of each temperature sensor corresponds to an acquisition point, where the number of temperature sensors in the temperature sensor array is not specifically limited, and the detection method includes the following steps:
s10, collecting the temperature value of each collection point according to the time sequence;
s20, obtaining the change trend of the temperature at each acquisition point along with the time according to the temperature value of each acquisition point; comparing the temperature change trend of each acquisition point, and determining the acquisition point corresponding to the maximum temperature rise value as a standing point position;
s30, selecting temperature values T collected by at least two temperature sensors near the stagnation position, substituting the temperature values T into the following formula to calculate the icing thickness H of the stagnation position:
Figure 458426DEST_PATH_IMAGE001
(1),
Figure 637734DEST_PATH_IMAGE002
(2),
Figure 896546DEST_PATH_IMAGE023
(3),
wherein the content of the first and second substances,
Figure 8859DEST_PATH_IMAGE004
for the heat flow of the ice surface at the initial stage of icing,
Figure 812867DEST_PATH_IMAGE005
is the mass of impinging water per unit area, n is the local freezing coefficient,
Figure 745051DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure 541974DEST_PATH_IMAGE007
for water droplet collection rate, LWC is the liquid water content,C w is the specific heat of the water, and the specific heat of the water,
Figure 508793DEST_PATH_IMAGE008
is the incoming flow velocity.
Wherein, an Euler-Euler method is adopted to simulate the flow field of air and water drops, and the heat flow of the ice surface at the initial stage of icing is calculated through numerical simulation
Figure 483703DEST_PATH_IMAGE004
Heat transfer from water film to ice layer
Figure 168762DEST_PATH_IMAGE024
And the heat transfer from the ice layer to the surface of the structure is 0, the convection heat transfer in the thin water film on the surface of the ice layer is neglected, and the temperature on a freezing interface is assumed to be 0 ℃.
Heat flow from ice surface at initial stage of icing
Figure 769376DEST_PATH_IMAGE004
Can be approximated as:
Figure 590702DEST_PATH_IMAGE009
(4)
wherein, the first and the second end of the pipe are connected with each other,
Figure 189042DEST_PATH_IMAGE014
the energy put into it for impacting the drop,
Figure 33501DEST_PATH_IMAGE015
in order to evaporate the energy, the evaporator is provided with a heat exchanger,
Figure 437807DEST_PATH_IMAGE016
in order to freeze the latent heat of phase change,
Figure 113639DEST_PATH_IMAGE017
convection heat exchange is carried out inside the water film;
Figure 695930DEST_PATH_IMAGE010
(5),
Figure 824423DEST_PATH_IMAGE025
(6),
Figure 993540DEST_PATH_IMAGE012
(7),
Figure 507567DEST_PATH_IMAGE026
(8),
wherein the content of the first and second substances,
Figure 526338DEST_PATH_IMAGE027
is the mass of ice formed per unit area,T w the wall temperature after the impact of the droplet (set to 0 ℃),
Figure 79810DEST_PATH_IMAGE018
in order to obtain the temperature of the incoming flow,V w is the droplet velocity, equal to the incoming flow velocity
Figure 560339DEST_PATH_IMAGE008
L Is the latent heat of solidification in the liquid state, m in order to be the mass transfer coefficient,is a convective heat transfer coefficient, C s And
Figure 679605DEST_PATH_IMAGE019
respectively the mass concentration of the water vapor near the water surface and the mass concentration of the water vapor of the incoming air,L e is the latent heat of evaporation per unit mass, t s And t 1 Respectively the wall temperature and the ambient temperature outside the boundary layer. It is worth noting that the heat transfer and mass transfer during the water surface evaporation process are similar, and the temperature distribution and concentration distribution are interconnected, so that the heat transfer coefficient can be changed by convectionObtaining the mass transfer coefficient m
Wherein the heat convection coefficienthCan pass through and flow speed
Figure 587387DEST_PATH_IMAGE008
And the temperature of incoming flow
Figure 690472DEST_PATH_IMAGE018
Is calculated from the function of (a). Namely, an Euler-Euler method is adopted to simulate the flow field of air and water drops to obtain different incoming flow velocities
Figure 725424DEST_PATH_IMAGE008
And the temperature of incoming flow
Figure 682885DEST_PATH_IMAGE018
Lower convective heat transfer coefficienthThen, in subsequent calculations, based on the incoming flow velocity
Figure 246721DEST_PATH_IMAGE008
And the temperature of incoming flow
Figure 758474DEST_PATH_IMAGE018
Determination of heat transfer coefficient of flowh
Substituting formulas (5) - (8) into formula (4), substituting (4) into (1), substituting (3) into (2), substituting at least two temperature values into formula (1), and simultaneously establishing formulas (1) and (2), three unknown parameters of liquid water content LWC, icing thickness H and local freezing coefficient n can be calculated.
Wherein, the icing thickness H is the icing thickness at the stagnation point. In addition, cloud and mist parameter values can be obtained through calculation. The technical personnel in the field know that the cloud and mist parameter value is not easy to obtain directly, and the method can directly obtain the LWC (liquid water content) on the basis of calculating the icing thickness, and is convenient. And preferably, another cloud parameter value MVD can be further calculated by adopting the following method.
Firstly, simulating air and water drop flow field by using Euler-Euler method to obtain the wall surface temperature of the specific shape structure of the object surface to be measured, which is 0 ℃, different MVDs, different LWCs and different incoming flow speeds
Figure 580806DEST_PATH_IMAGE008
Lower water droplet collection rate
Figure 127193DEST_PATH_IMAGE007
To obtain
Figure 861931DEST_PATH_IMAGE007
And the velocity of incoming flow
Figure 392139DEST_PATH_IMAGE008
The relationship between MVD and LWC;
secondly, on the basis of the LWC calculated in the previous step, the incoming flow speed
Figure 768893DEST_PATH_IMAGE008
Measured with a wind speed sensor and stationary
Figure 982837DEST_PATH_IMAGE007
And =1, by substituting these parameters into the above relationship, the corresponding MVD value can be directly obtained, that is, the median diameter MVD of the water droplets in the cloud and mist parameters is obtained.
The MVD value at a non-stationary point can also be obtained by the following method: firstly, determining an impact limit according to a local freezing coefficient n, and then determining a local freezing coefficient n according to a distance between a stagnation point and the impact limit
Figure 668902DEST_PATH_IMAGE007
Determination of change relationship
Figure 437138DEST_PATH_IMAGE007
Value is further based on
Figure 617584DEST_PATH_IMAGE020
An MVD value is obtained. Specifically, the method comprises the following steps:
when n is more than or equal to 0.95, selecting the position with the minimum temperature change value closest to the stagnation point as the impact limit; when n < 0.95, the temperature change minimum position closest to the stagnation point is selected as the impact limit within a predetermined time (for example, 10 s) in the initial stage of icing.
At a stagnation point
Figure 669722DEST_PATH_IMAGE007
=1, impact limit position
Figure 995530DEST_PATH_IMAGE007
=0, from a stagnation position to an impact limit position
Figure 516641DEST_PATH_IMAGE007
The value is exponentially reduced, and the position relation between the position of the point to be measured and the stagnation point and the impact limit can be obtained
Figure 235198DEST_PATH_IMAGE007
Value of will
Figure 938581DEST_PATH_IMAGE007
Value, calculated LWC value substitution
Figure 451602DEST_PATH_IMAGE020
And obtaining the MVD value of the point to be measured.
Preferably, the maximum icing thickness can also be calculated from the temperature profile at the stagnation point:
and calculating the time difference of the positions of the acquisition points near the stagnation point, which are subjected to the steep temperature rise and the steep temperature fall on the time scale, wherein the time difference is used as the icing time length, and the maximum icing thickness is calculated according to the icing time length. The time point of temperature steep rise is the time point of starting icing, the time point of temperature steep drop is the point of ending icing, the time difference between the time point of starting icing and the time point of ending icing is icing time length, the icing speed and the icing time length are calculated based on the formula (2), and the maximum icing thickness at the stationary point and the maximum icing thickness of the object surface to be measured can be obtained.
Thus, in this embodiment 1, the icing thickness and the maximum icing thickness can be calculated based on the temperature values acquired by the temperature sensor array, and the cloud and mist parameters which are not easily obtained by the conventional method are also calculated: LWC and MVD; impact limit positions can also be obtained.
The method is very suitable for detecting the icing of the aircraft, mainly because the change of the flight attitude of the aircraft can cause the change of the location of the stagnation point in the flight process, the method does not need to calculate the location of the stagnation point in advance, but searches the location of the stagnation point again in the icing calculation process for calculation. Therefore, the temperature sensor array can be randomly arranged when arranged, the conventional stagnation position does not need to be considered, and the sensor is more convenient to mount. And only the temperature sensor is installed, and more integrated modules are not provided, so that the volume is small, and the cost is low.
Example 2
The embodiment provides an icing detecting device based on temperature measurement, which is used for executing the icing detecting method based on temperature measurement as described in embodiment 1, and as shown in fig. 3, the icing detecting device comprises a temperature sensor array, a temperature acquisition module and a calculation module;
the temperature sensor array is arranged on the surface of an object to be measured, as shown in figure 2; thermocouples with small volume and small heat capacity are preferred as temperature sensors. There is no thermal contact between the temperature sensors, i.e. the distance between the temperature sensors in the array of temperature sensors is such that there is no thermal contact between them. Preferably, the surface of the sensor is covered with a weather-resistant film, and a heat insulating layer is arranged between the surface of the object to be measured and the sensor.
The acquisition module executes step S10, that is, acquires each acquisition point (i.e., temperature value of each temperature sensor) according to a time sequence;
the calculation module executes steps S20 to S30, that is, processes the acquired temperature values, determines the stagnation position, and substitutes a calculation formula based on at least two temperature values near the stagnation position to calculate the icing thickness at the stagnation position and the cloud and mist parameters.
Example 3
Embodiment 3 also provides an icing detection device based on temperature measurement for executing an icing detection method based on temperature measurement as described in embodiment 1. The difference between the embodiment and the embodiment 2 is that in the embodiment 2, the temperature sensor is arranged on the object surface to be measured to measure the icing condition of the object surface to be measured, and in the embodiment, the temperature sensor is arranged on the periphery of the cylinder 30 to be made into an independent icing detection device, and the icing detection device is installed near the object surface to be measured to be used, which is equivalent to an icing detector.
Specifically, the icing detection device comprises a temperature sensor array, a temperature acquisition module and a calculation module;
the temperature sensor array is arranged on the periphery of the cylinder 30, as shown in fig. 4; because the icing detection method can self-adaptively find the stagnation position according to the detected temperature value, the temperature sensor does not need to be arranged at any position when the periphery of the cylinder is arranged, and the icing detection method can be used for randomly arranging the temperature sensor. Preferably, the temperature sensor array is arranged at least one circle on the periphery of the cylinder, so that the temperature sensor array can be freely arranged on the surface of the object to be measured when the device is used, and the inflow direction is not required to be considered. Compared with the existing icing detection device, the icing detection device has the advantages of low manufacturing cost, long service life and good performance stability.
The above description is intended to be illustrative of the preferred embodiment of the present invention and should not be taken as limiting the invention, but rather, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.

Claims (8)

1. An icing detection method based on temperature measurement is characterized in that a temperature sensor array is arranged on an object surface to be detected, the position of each temperature sensor corresponds to an acquisition point, and the detection method comprises the following steps:
s10, collecting the temperature value of each collection point according to the time sequence, and obtaining the change trend of the temperature on each collection point along with the time;
s20, obtaining the change trend of the temperature at each acquisition point along with time according to the temperature value of each acquisition point; comparing the temperature change trend of each acquisition point, and determining the acquisition point corresponding to the maximum temperature rise value as a stationary point position;
s30, selecting temperature values T collected by at least two temperature sensors near the stagnation position, substituting the temperature values T into the following formula to calculate the icing thickness H of the stagnation position:
Figure 800244DEST_PATH_IMAGE001
Figure 389488DEST_PATH_IMAGE002
Figure 700383DEST_PATH_IMAGE003
wherein the content of the first and second substances,
Figure 751385DEST_PATH_IMAGE004
for the heat flow of the ice surface at the initial stage of icing,
Figure 96916DEST_PATH_IMAGE005
is the mass of impinging water per unit area, n is the local freezing coefficient,
Figure 653799DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure 327357DEST_PATH_IMAGE008
LWC is the liquid water content for water droplet collection efficiency,C w is the specific heat of the water, and the specific heat of the water,
Figure 57415DEST_PATH_IMAGE009
is the incoming flow velocity.
2. Icing detection method based on temperature measurement according to claim 1, characterized in that the heat flow of the ice surface during the initial stage of icing is
Figure 382086DEST_PATH_IMAGE004
Calculating by numerical simulation:
Figure 172188DEST_PATH_IMAGE010
Figure 333042DEST_PATH_IMAGE011
Figure 601212DEST_PATH_IMAGE012
Figure 780390DEST_PATH_IMAGE014
Figure 944655DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 717439DEST_PATH_IMAGE016
the energy put into it for impacting the drop,
Figure 399087DEST_PATH_IMAGE017
in order to evaporate the energy, the evaporator is provided with a heat exchanger,
Figure 573716DEST_PATH_IMAGE018
in order to freeze the latent heat of phase change,
Figure 830254DEST_PATH_IMAGE019
for the convection heat exchange in the water film,T w the temperature of the liquid drops on the wall surface after impact,
Figure 90334DEST_PATH_IMAGE020
in order to obtain the temperature of the incoming flow,V w in order to be the velocity of the liquid droplets,L is the latent heat of solidification in the liquid state, m in order to be the mass transfer coefficient,to convective heat transfer coefficient, C s And
Figure 575673DEST_PATH_IMAGE021
respectively the mass concentration of the water vapor near the water surface and the mass concentration of the water vapor of the incoming air,L e is the latent heat of evaporation per unit mass, t s And t 1 Respectively the wall temperature and the ambient temperature outside the boundary layer.
3. The icing detection method based on temperature measurement according to claim 2, characterized in that in step S30, the liquid water content LWC is calculated and the water drop median diameter MVD is calculated according to the following relation:
Figure 604809DEST_PATH_IMAGE023
(ii) a Wherein f is the temperature of the object surface to be measured obtained in the process of simulating the flow field of air and water drops by adopting an Euler-Euler method at 0 ℃,
Figure 110877DEST_PATH_IMAGE024
and with
Figure 717308DEST_PATH_IMAGE008
In which water droplet collection coefficient at stagnation point
Figure 130971DEST_PATH_IMAGE008
The number of the carbon atoms is 1,
Figure 155559DEST_PATH_IMAGE025
collected by a wind speed sensor.
4. The icing detection method based on temperature measurement according to the claim 2, characterized in that in step S30, the local freezing coefficient n is calculated simultaneously; when n is more than or equal to 0.95, selecting the position with the minimum temperature change value nearest to the stagnation point as an impact limit; and when n is less than 0.95, selecting the position with the minimum temperature change value closest to the stagnation point as the impact limit in the preset time at the initial stage of icing.
5. The icing detection method based on temperature measurement according to claim 3, wherein in step S30, the liquid water content LWC is calculated; determining the water drop collection coefficient of the position to be measured according to the variation relation of the water drop collection coefficient between the stagnation point and the impact limit
Figure 629266DEST_PATH_IMAGE008
And will be
Figure 722993DEST_PATH_IMAGE008
Substituting LWC into the following relational expression to calculate the median diameter MVD of the water drop at the position to be measured:
Figure 674768DEST_PATH_IMAGE023
(ii) a Wherein the content of the first and second substances,when the temperature of an object surface to be measured obtained in the process of simulating the flow field of air and water drops by adopting an Euler-Euler method is 0 ℃,
Figure 881759DEST_PATH_IMAGE024
and
Figure DEST_PATH_IMAGE026
the functional relationship of (a) to (b),
Figure 932891DEST_PATH_IMAGE025
collected by a wind speed sensor.
6. The icing detection method based on temperature measurement is characterized in that the time difference of the positions of the steep temperature rise and the steep temperature fall of the acquisition point near the stagnation point on the time scale is calculated and used as the icing time length, and the maximum icing thickness is calculated according to the icing time length.
7. An icing detection device based on temperature measurement, which is used for executing the icing detection method based on temperature measurement according to any one of claims 1-6, and comprises a temperature sensor array, a temperature acquisition module and a calculation module;
the temperature sensor array is arranged on the surface of an object to be measured or arranged around the circumference of the cylinder;
the acquisition module executes step S10;
the calculation module performs steps S20-S30.
8. The icing detection device based on temperature measurement as claimed in claim 7, wherein the distance of the temperature sensor array is required to be as follows: there is no thermal contact between the temperature sensors.
CN202211381850.1A 2022-11-07 2022-11-07 Icing detection device and icing detection method based on temperature measurement Active CN115416854B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211381850.1A CN115416854B (en) 2022-11-07 2022-11-07 Icing detection device and icing detection method based on temperature measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211381850.1A CN115416854B (en) 2022-11-07 2022-11-07 Icing detection device and icing detection method based on temperature measurement

Publications (2)

Publication Number Publication Date
CN115416854A true CN115416854A (en) 2022-12-02
CN115416854B CN115416854B (en) 2023-01-24

Family

ID=84207418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211381850.1A Active CN115416854B (en) 2022-11-07 2022-11-07 Icing detection device and icing detection method based on temperature measurement

Country Status (1)

Country Link
CN (1) CN115416854B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117168331A (en) * 2023-11-02 2023-12-05 山西锦烁生物医药科技有限公司 Real-time detection method for thickness of ice layer of natural ice rink based on optical fiber sensor

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB575257A (en) * 1942-10-24 1946-02-11 Honeywell Brown Ltd Electrical apparatus for measuring the thickness of ice on a surface
US5005015A (en) * 1989-08-07 1991-04-02 General Electric Company Ice detection system
US6239601B1 (en) * 1996-03-20 2001-05-29 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Thickness measurement device for ice, or ice mixed with water or other liquid
US6328467B1 (en) * 1999-05-07 2001-12-11 University Of Tennessee Research Corp. Method and apparatus for detecting ice or frost deposition
US6347767B1 (en) * 1997-11-05 2002-02-19 Futuris As. Method of and apparatus for detection of ice accretion
US7439877B1 (en) * 2007-05-18 2008-10-21 Philip Onni Jarvinen Total impedance and complex dielectric property ice detection system
US20100116940A1 (en) * 2007-04-11 2010-05-13 Nicolas Picco Method and device for detecting rime and/or rime conditions on a flying aircraft
US20120085868A1 (en) * 2010-10-12 2012-04-12 Lumen International Inc. Aircraft icing detector
US20130238282A1 (en) * 2010-07-05 2013-09-12 Eduardo Figueroa-Karlström Device and method for measuring ice thickness
CN103983376A (en) * 2014-05-14 2014-08-13 太原理工大学 Automatic online monitoring device for temperature gradient and thickness of ice layer
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
US20150260501A1 (en) * 2014-03-13 2015-09-17 Rolls-Royce Plc Ice thickness measurement sensor
CN111738481A (en) * 2020-04-01 2020-10-02 南京航空航天大学 Airplane icing meteorological parameter MVD prediction method based on BP neural network
CN113420503A (en) * 2021-06-11 2021-09-21 中国民航大学 Icing detection method, system, terminal and application based on temperature sequence similarity measurement
US20210316869A1 (en) * 2018-09-10 2021-10-14 Instituto Nacional De Tecnica Aeroespacial System and Method for Detecting Ice Formation on a Body
RU2763473C1 (en) * 2021-05-26 2021-12-29 Виктор Александрович Бараусов Method and device for determining ice thickness on the working surface of ice-formation sensor
CN114180072A (en) * 2022-02-16 2022-03-15 中国空气动力研究与发展中心低速空气动力研究所 Icing thickness detection method
CN114516403A (en) * 2022-03-07 2022-05-20 中国空气动力研究与发展中心低速空气动力研究所 Power control method of electric heating ice prevention and removal system based on optical fiber icing detection
CN114676560A (en) * 2022-03-07 2022-06-28 中国空气动力研究与发展中心低速空气动力研究所 Cloud and mist parameter calculation and ice shape prediction method and medium based on optical fiber icing detection
CN114970308A (en) * 2021-12-30 2022-08-30 成都流体动力创新中心 Aircraft icing prediction method and system and computer program product

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB575257A (en) * 1942-10-24 1946-02-11 Honeywell Brown Ltd Electrical apparatus for measuring the thickness of ice on a surface
US5005015A (en) * 1989-08-07 1991-04-02 General Electric Company Ice detection system
US6239601B1 (en) * 1996-03-20 2001-05-29 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Thickness measurement device for ice, or ice mixed with water or other liquid
US6347767B1 (en) * 1997-11-05 2002-02-19 Futuris As. Method of and apparatus for detection of ice accretion
US6328467B1 (en) * 1999-05-07 2001-12-11 University Of Tennessee Research Corp. Method and apparatus for detecting ice or frost deposition
US20100116940A1 (en) * 2007-04-11 2010-05-13 Nicolas Picco Method and device for detecting rime and/or rime conditions on a flying aircraft
US7439877B1 (en) * 2007-05-18 2008-10-21 Philip Onni Jarvinen Total impedance and complex dielectric property ice detection system
US20130238282A1 (en) * 2010-07-05 2013-09-12 Eduardo Figueroa-Karlström Device and method for measuring ice thickness
US20120085868A1 (en) * 2010-10-12 2012-04-12 Lumen International Inc. Aircraft icing detector
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
US20150260501A1 (en) * 2014-03-13 2015-09-17 Rolls-Royce Plc Ice thickness measurement sensor
CN103983376A (en) * 2014-05-14 2014-08-13 太原理工大学 Automatic online monitoring device for temperature gradient and thickness of ice layer
US20210316869A1 (en) * 2018-09-10 2021-10-14 Instituto Nacional De Tecnica Aeroespacial System and Method for Detecting Ice Formation on a Body
CN111738481A (en) * 2020-04-01 2020-10-02 南京航空航天大学 Airplane icing meteorological parameter MVD prediction method based on BP neural network
RU2763473C1 (en) * 2021-05-26 2021-12-29 Виктор Александрович Бараусов Method and device for determining ice thickness on the working surface of ice-formation sensor
CN113420503A (en) * 2021-06-11 2021-09-21 中国民航大学 Icing detection method, system, terminal and application based on temperature sequence similarity measurement
CN114970308A (en) * 2021-12-30 2022-08-30 成都流体动力创新中心 Aircraft icing prediction method and system and computer program product
CN114180072A (en) * 2022-02-16 2022-03-15 中国空气动力研究与发展中心低速空气动力研究所 Icing thickness detection method
CN114516403A (en) * 2022-03-07 2022-05-20 中国空气动力研究与发展中心低速空气动力研究所 Power control method of electric heating ice prevention and removal system based on optical fiber icing detection
CN114676560A (en) * 2022-03-07 2022-06-28 中国空气动力研究与发展中心低速空气动力研究所 Cloud and mist parameter calculation and ice shape prediction method and medium based on optical fiber icing detection

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张雨杏: "基于SH波的旋翼结冰冰层厚度检测机理研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
於萧萧等: "基于水滴运动的机翼结冰数值模拟", 《江苏航空》 *
祁相莹: "基于传热原理的结冰条件探测方法及仿真研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 *
胡站伟: "结冰风洞云雾参数变化规律初探", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117168331A (en) * 2023-11-02 2023-12-05 山西锦烁生物医药科技有限公司 Real-time detection method for thickness of ice layer of natural ice rink based on optical fiber sensor
CN117168331B (en) * 2023-11-02 2024-01-02 山西锦烁生物医药科技有限公司 Real-time detection method for thickness of ice layer of natural ice rink based on optical fiber sensor

Also Published As

Publication number Publication date
CN115416854B (en) 2023-01-24

Similar Documents

Publication Publication Date Title
CN107132278B (en) Multi-cylinder array icing detection method
CN115416854B (en) Icing detection device and icing detection method based on temperature measurement
US8408871B2 (en) Method and apparatus for measuring air flow condition at a wind turbine blade
CN110879127A (en) Coupling type icing detection device and method based on airfoil
CN104180980B (en) External wind turbine freezing detector and wind turbine freezing detecting method
CN114676560B (en) Cloud and mist parameter calculation and ice shape prediction method and medium based on optical fiber icing detection
CN102706389A (en) System and method for predicting ice accretion on surface of aircraft through rolling forecast within short time
CN107607579B (en) Wet cooling tower thermal performance online monitoring system and method
CN113420503B (en) Icing detection method, system, terminal and application based on temperature sequence similarity measurement
CN109556551B (en) Icing thickness monitoring method based on interface temperature
CN113899663B (en) MVD (mechanical vapor recompression) measuring method for complex icing cloud mist in large icing wind tunnel
CN110127061A (en) The method for improving icing detector detection bar susceptibility
CN206876374U (en) A kind of super-cooling waterdrop ice detection
CN115507734A (en) Icing three-dimensional imaging method based on piezoelectric array
CN113899784B (en) Method for measuring content of complex icing cloud mist water in large icing wind tunnel
CN112572809B (en) Hybrid icing detection method suitable for unmanned aerial vehicle platform
Gui et al. Atmospheric icing process measurement utilizing impedance spectroscopy and thin film structure
CN109598033B (en) Icing parameter calculation method based on rotating multiple conductors
CN114877837B (en) Icing thickness measuring method and system based on far-field sound pressure level change
CN111173687A (en) On-line monitoring device and method for crack damage of wind power fan blade
CN110715750A (en) Anti-icing and deicing method for total temperature sensor
CN112550723B (en) Super-cooled large water drop icing probe and detector
CN115649455B (en) Method and device for judging icing based on electrothermal deicing signal
CN111216898A (en) Grating optical fiber icing detector
CN106645276B (en) A kind of measurement method of ice coating wire icicle surface Local Condensing Heat Transfer Coefficients

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