CN115416854B - 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
CN115416854B
CN115416854B CN202211381850.1A CN202211381850A CN115416854B CN 115416854 B CN115416854 B CN 115416854B CN 202211381850 A CN202211381850 A CN 202211381850A CN 115416854 B CN115416854 B CN 115416854B
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.)
Active
Application number
CN202211381850.1A
Other languages
Chinese (zh)
Other versions
CN115416854A (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 deicing prevention, 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 aircraft flight accidents, and icing on the leading edges of the wings and the empennage of the aircraft can cause serious flight accidents due to increased wing profile resistance, reduced lift force, reduced critical attack angle and deteriorated maneuverability and stability, so that the aircraft is widely concerned 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 and the electric conduction method are not widely applied because the differential pressure method and the electric conduction method are easily interfered by the outside, and the icing thickness information cannot be effectively given.
The most widely used products in the present foreign molded products are magnetostrictive icing detection systems of ROSEMOUNT company in America, flat-film icing sensors of VibroMeter company in Sweden, optical fiber icing sensors produced in British, 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 multi-point 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 photoelectric conversion, electric signal amplification, multi-way switches and other back-end signal acquisition and processing, and with the analysis point as high as tens of orders of magnitude and hundreds of orders of magnitude, the system integration has no engineering significance.
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 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 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 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
mass of impinging water per unit area, n is a local freezing coefficient,
Figure 150908DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure 715882DEST_PATH_IMAGE007
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 880016DEST_PATH_IMAGE008
is the incoming flow velocity.
Further, the heat flow of the ice surface at the initial stage of the freezing
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 a droplet,
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
for the convection heat exchange in the water film, T w the temperature of the liquid drops on the wall surface after impact,
Figure 857571DEST_PATH_IMAGE018
in order to obtain 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 water vapor of the incoming airThe mass concentration of (a) to (b),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 the water droplet collection coefficient at the stagnation point
Figure 954894DEST_PATH_IMAGE007
The number of the carbon atoms 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 closest to the stagnation point as the 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.
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 variation 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
And substituting the 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 piece, 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-wing, 20-temperature sensor array, 30-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 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 458426DEST_PATH_IMAGE001
(1),
Figure 637734DEST_PATH_IMAGE002
(2),
Figure 896546DEST_PATH_IMAGE023
(3),
wherein, the first and the second end of the pipe are connected with each other,
Figure 8859DEST_PATH_IMAGE004
for the heat flow of the ice surface at the initial stage of icing,
Figure 812867DEST_PATH_IMAGE005
mass of impinging water per unit area, n is a local freezing coefficient,
Figure 745051DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure 541974DEST_PATH_IMAGE007
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 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
Transfer of heat from the water film to the 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 the 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 a droplet,
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 first and the second end of the pipe are connected with each other,
Figure 526338DEST_PATH_IMAGE027
is the mass of ice formed per unit area,T w the wall temperature after the impact of the droplets was set to 0C,
Figure 79810DEST_PATH_IMAGE018
in order to achieve 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,to 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 mentioning that the heat transfer and mass transfer during the evaporation process on the water surface are similar, and the temperature distribution and the concentration distribution are mutually connected, so that the heat transfer coefficient can be realized through convectionObtaining the mass transfer coefficient m
Wherein the heat convection coefficienthCan pass through and flow velocity
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 the incoming flow
Figure 682885DEST_PATH_IMAGE018
Lower convective heat transfer coefficienthThen in the subsequent calculation, according to the incoming flow speed
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), the liquid water content LWC, the icing thickness H and the local freezing coefficient n of three unknown parameters 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. As known by persons skilled in the art, cloud and mist parameter values are not easy to obtain directly, and the method can be used for directly obtaining the liquid water content LWC on the basis of calculating the icing thickness, and is convenient and fast. 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, substituting the parameters into the above relationship, so as to directly obtain a corresponding MVD value, that is, the water droplet median diameter MVD in the cloud and mist parameters is obtained.
The MVD values for non-stationary points can also be obtained by the following method: firstly, determining an impact limit according to a local freezing coefficient n, and then determining a position 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
The 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 nearest to the stagnation point as an impact limit; when n < 0.95, the minimum temperature change 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 steep temperature rise and steep temperature fall positions of the acquisition points near the stagnation point on the time scale, taking the time difference as the icing time length, and calculating the maximum icing thickness according to the time length. The time point of the steep temperature rise is the time point of starting icing, the time point of the steep temperature 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 obtained through calculation based on the formula (2), the maximum icing thickness at the stagnation point can be obtained, and the maximum icing thickness of the object surface to be measured can also 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 a stagnation point position, and substitutes a calculation formula to calculate the icing thickness at the stagnation point position and the cloud and mist parameters based on at least two temperature values near the stagnation point position.
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 freezing detection method can self-adaptively find the stagnation position according to the detected temperature value, the temperature sensor does not need to consider the position to arrange when the periphery of the cylinder is arranged, and the arrangement can be carried out at will. 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 (7)

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 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 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 DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
wherein, in the step (A),
Figure DEST_PATH_IMAGE004
the heat flow of the ice surface at the initial stage of icing,
Figure DEST_PATH_IMAGE005
is the mass of impinging water per unit area, n is the local freezing coefficient,
Figure DEST_PATH_IMAGE006
is the density of the ice, and is,
Figure DEST_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 DEST_PATH_IMAGE009
is the incoming flow velocity;
heat flow of ice surface at initial stage of icing
Figure 607139DEST_PATH_IMAGE004
Calculating by numerical simulation:
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE011
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE014
Figure DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE016
the energy put into it for impacting the drop,
Figure DEST_PATH_IMAGE017
in order to evaporate the energy, the evaporator is provided with a heat exchanger,
Figure DEST_PATH_IMAGE018
in order to freeze the latent heat of phase change,
Figure DEST_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 DEST_PATH_IMAGE020
in order to achieve 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 DEST_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.
2. The icing detection method based on temperature measurement as claimed in claim 1, wherein in step S30, the liquid water content LWC is calculated at the same time, and the water drop median diameter MVD is calculated according to the following relation:
Figure DEST_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 DEST_PATH_IMAGE024
and
Figure 181951DEST_PATH_IMAGE008
in which the water droplet collection coefficient at the stagnation point
Figure 842739DEST_PATH_IMAGE008
The number of the carbon atoms is 1,
Figure DEST_PATH_IMAGE025
collected by a wind speed sensor.
3. The icing detection method based on temperature measurement according to the claim 1, 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 closest to the stagnation point as the 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.
4. 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; 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 283210DEST_PATH_IMAGE008
And will be
Figure 433569DEST_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 966181DEST_PATH_IMAGE023
(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 227398DEST_PATH_IMAGE024
and
Figure DEST_PATH_IMAGE026
the functional relationship of (a) to (b),
Figure 522376DEST_PATH_IMAGE025
collected by a wind speed sensor.
5. 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.
6. 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-5, 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.
7. The icing detection device based on temperature measurement as claimed in claim 6, 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 CN115416854A (en) 2022-12-02
CN115416854B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Citations (4)

* 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
US6347767B1 (en) * 1997-11-05 2002-02-19 Futuris As. Method of and apparatus for detection of ice accretion
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
CN114970308A (en) * 2021-12-30 2022-08-30 成都流体动力创新中心 Aircraft icing prediction method and system and computer program product

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
FR2914906B1 (en) * 2007-04-11 2009-10-30 Intertechnique Soc Par Actions METHOD AND DEVICE FOR DETECTING FROST AND / OR GIVING CONDITIONS ON AIRCRAFT IN FLIGHT
US7439877B1 (en) * 2007-05-18 2008-10-21 Philip Onni Jarvinen Total impedance and complex dielectric property ice detection system
US9625248B2 (en) * 2010-07-05 2017-04-18 Saab Ab Device and method for measuring ice thickness
US8462354B2 (en) * 2010-10-12 2013-06-11 Lumen International Inc. Aircraft icing detector
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
GB201404462D0 (en) * 2014-03-13 2014-04-30 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
EP3620379A1 (en) * 2018-09-10 2020-03-11 INTA, Instituto Nacional de Technica 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
CN113420503B (en) * 2021-06-11 2022-08-12 中国民航大学 Icing detection method, system, terminal and application based on temperature sequence similarity measurement
CN114180072B (en) * 2022-02-16 2022-04-12 中国空气动力研究与发展中心低速空气动力研究所 Icing thickness detection method
CN114676560B (en) * 2022-03-07 2023-04-11 中国空气动力研究与发展中心低速空气动力研究所 Cloud and mist parameter calculation and ice shape prediction method and medium based on optical fiber icing detection

Patent Citations (4)

* 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
US6347767B1 (en) * 1997-11-05 2002-02-19 Futuris As. Method of and apparatus for detection of ice accretion
CN114970308A (en) * 2021-12-30 2022-08-30 成都流体动力创新中心 Aircraft icing prediction method and system and computer program product
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

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于SH波的旋翼结冰冰层厚度检测机理研究;张雨杏;《中国优秀硕士学位论文全文数据库工程科技II辑》;20190215(第2期);全文 *

Also Published As

Publication number Publication date
CN115416854A (en) 2022-12-02

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
CN110879127A (en) Coupling type icing detection device and method based on airfoil
CN114676560B (en) Cloud and mist parameter calculation and ice shape prediction method and medium based on optical fiber icing detection
CN104180980B (en) External wind turbine freezing detector and wind turbine freezing detecting method
CN111396269B (en) Multi-time-step unsteady icing calculation method and system and storage medium
CN113420503B (en) Icing detection method, system, terminal and application based on temperature sequence similarity measurement
CN113899663B (en) MVD (mechanical vapor recompression) measuring method for complex icing cloud mist in large icing wind tunnel
CN109556551B (en) Icing thickness monitoring method based on interface temperature
CN105548285A (en) Method for measuring thin ice by adopting finger-inserting micro-structural capacitive thin ice sensor
CN110127061A (en) The method for improving icing detector detection bar susceptibility
CN201697932U (en) Capacitance type anemoclinograph of wind turbine generator set
CN206876374U (en) A kind of super-cooling waterdrop ice detection
CN110567357A (en) Dynamic strain piezoelectric ceramic icing detection sensor
CN113899784B (en) Method for measuring content of complex icing cloud mist water in large icing wind tunnel
CN105160599A (en) Power line icing assessment method based on filtering and support vector machine
CN103901505B (en) A kind of cloud-top height detection method based on wet bulb effect and device
CN109598033B (en) Icing parameter calculation method based on rotating multiple conductors
Lv et al. Pitot tube-based icing detection: Effect of ice blocking on pressure
CN112550723B (en) Super-cooled large water drop icing probe and detector
CN115060374A (en) Infrared thermopile sensor capable of improving absorption efficiency and MEMS process manufacturing method thereof
CN111104754B (en) Icing parameter calculation method based on particle swarm optimization
CN115649455B (en) Method and device for judging icing based on electrothermal deicing signal
CN115248022B (en) Icing thickness detection method for power transmission line
CN117761264A (en) liquid water content detector and detection method based on total temperature measurement technology

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