CN116644258A - Anti-icing heat load calculation method and system - Google Patents
Anti-icing heat load calculation method and system Download PDFInfo
- Publication number
- CN116644258A CN116644258A CN202310883924.XA CN202310883924A CN116644258A CN 116644258 A CN116644258 A CN 116644258A CN 202310883924 A CN202310883924 A CN 202310883924A CN 116644258 A CN116644258 A CN 116644258A
- Authority
- CN
- China
- Prior art keywords
- icing
- grid
- heat
- calculating
- thermal load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 45
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 60
- 239000007787 solid Substances 0.000 claims description 58
- 238000001704 evaporation Methods 0.000 claims description 26
- 230000008020 evaporation Effects 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 20
- 230000000704 physical effect Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 7
- 238000004134 energy conservation Methods 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000032823 cell division Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/002—Thermal testing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64D—EQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENT OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
- B64D15/00—De-icing or preventing icing on exterior surfaces of aircraft
- B64D15/12—De-icing or preventing icing on exterior surfaces of aircraft by electric heating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Air Conditioning Control Device (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The application relates to an anti-icing heat load calculation method and a system, wherein the calculation method comprises the following steps: based on the divided wall grid cells, for any divided grid cell, acquiring the respective outer surface temperatures of the grid cell and the adjacent grid cells; calculating the anti-icing heat load of the outer surfaces of each grid unit and the adjacent grid units; calculating the respective center temperatures of the grid cells and the adjacent grid cells; the inner surface anti-icing thermal load of the grid cell is calculated and taken as the anti-icing thermal load of the grid cell. The outer surface temperature is converted into the center temperature, so that the anti-icing heat load of the outer surface is calculated based on internal heat conduction, and the anti-icing heat load of the inner surface is used as the anti-icing heat load, so that the calculated anti-icing heat load is closer to the load required by a real anti-icing system.
Description
Technical Field
The application relates to the technical field of aviation anti-icing, in particular to an anti-icing heat load calculation method and system capable of considering solid internal heat conduction.
Background
When an aircraft passes through a cloud layer containing supercooled water drops, icing phenomenon can occur on the windward surface of the aircraft, and the icing phenomenon can seriously damage the flight safety of the aircraft. It is therefore often necessary to design and install an anti-icing system at critical locations on the aircraft. The thermal anti-icing system has the advantages of high efficiency, good reliability and the like, and is mainly an anti-icing means adopted on aircraft and aeroengines at present. During the design process of the thermal anti-icing system, the anti-icing thermal load is calculated firstly, and then the designed anti-icing system is further checked and optimized. In the existing anti-icing heat load calculation method, the anti-icing heat load of the outer surface (the surface contacted with air) of the solid wall surface is calculated, and the obtained anti-icing heat load has the problem of larger error with the actual anti-icing heat load.
The foregoing is provided merely for the purpose of facilitating understanding of the technical solutions of the present application and is not intended to represent an admission that the foregoing is prior art.
Disclosure of Invention
The application aims to solve the technical problem of providing an anti-icing heat load calculation method and an anti-icing heat load calculation system, which have the characteristic that the calculated anti-icing heat load is closer to the actual anti-icing heat load.
In a first aspect, an embodiment provides an anti-icing thermal load calculation method, where grid cell division is performed on a solid wall object, where the solid wall object is divided into N wall grid cells, where N is a natural number greater than 1; for any one of the divided grid cells, the anti-icing thermal load calculation method includes:
acquiring the respective outer surface temperatures of the grid cells and the adjacent grid cells;
calculating the anti-icing heat load of the outer surfaces of each grid unit and the adjacent grid units based on the latent heat absorbed by evaporation of water drops, the heat of convection heat exchange between a solid wall object and air and the energy of impinging water drops;
calculating respective center temperatures of the grid cells and adjacent grid cells thereof based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads;
calculating an inner surface anti-icing thermal load of the grid cells based on the respective outer surface anti-icing thermal load and the respective center temperature;
and taking the anti-icing thermal load of the inner surface of the grid unit as the anti-icing thermal load of the grid unit.
In one embodiment, the anti-icing thermal load comprises a dry anti-icing thermal load and/or a wet anti-icing thermal load.
In one embodiment, the acquiring the respective outer surface temperatures of the grid cell and the adjacent grid cells includes:
for dry anti-icing, the formula is basedCalculating the respective outer surface temperatures of the grid cells and adjacent grid cells, wherein T s Is the outer surface temperature; />For evaporating quality, < >>；/>Impact mass for water droplets corresponding to the grid cells; b= 346.1; />，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is the constant, P 0 T is the atmospheric pressure of the environment 0 H is the convection heat exchange coefficient of the corresponding grid unit for the external environment temperature; h and->The method comprises the steps of performing simulation calculation on an air flow field and a water drop field outside a solid wall object;
and/or the number of the groups of groups,
for wet anti-icing, the grid cells and their adjacent grid cells each have an outer surface temperatureIs a constant.
In one embodiment, the calculating the anti-icing thermal load of the outer surface of each of the grid cells and the adjacent grid cells based on the latent heat absorbed by the evaporation of the water drops, the heat of convection heat exchange between the solid wall object and the air, and the energy of impinging the water drops comprises:
for dry anti-icing, the formula is basedCalculating the anti-icing heat load of the outer surfaces of each grid unit and each adjacent grid unit;
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, and is a constant; />Is the convection heat of solid wall object and air, +.>；For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flying speed;
for wet anti-icing, the formula is basedCalculating the anti-icing heat load of the outer surfaces of each grid unit and each adjacent grid unit;
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, is constant, < ->For evaporating quality, < >>，，/>，，/>，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is constant, P 0 T is the atmospheric pressure of the environment 0 H is the convection heat exchange coefficient of the corresponding grid unit for the external environment temperature; />Is the convection heat of solid wall object and air, +.>；/>For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flying speed>Impact mass for water droplets corresponding to the grid cells; h and->The method is obtained by performing simulation calculation on an air flow field and a water drop field outside the solid wall object.
In one embodiment, the calculating the respective center temperatures of the grid cells and the adjacent grid cells based on the respective outer surface temperatures and the respective outer surface anti-icing thermal loads includes:
for dry and/or wet anti-icing, the formula is basedCalculating the temperature T of the center of the nth grid cell n Wherein N is more than or equal to 1 and less than or equal to N; />Anti-icing thermal load for the outer surface of the nth grid cell; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; t (T) sn An outer surface temperature for the nth grid cell;
based on the formula and />Calculating a temperature T of a center corresponding to a grid cell adjacent to the grid cell n-1 and Tn+1 Wherein N is more than or equal to 1 and less than or equal to N; />An ice and heat load is prevented for the outer surface of the n-1 th grid cell adjacent to the n-th cell; />An outer surface of the (n+1) th grid cell adjacent to the (n) th cell is protected from an ice and heat load; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; /> and />Is the temperature of the outer surface corresponding to the grid cell adjacent to the grid cell.
In one embodiment, the calculating the grid cell inner surface anti-icing thermal load based on the respective outer surface anti-icing thermal load and the respective center temperature comprises:
for dry and/or wet anti-icing, based on energy conservation, according to the formulaCalculating an anti-icing thermal load of the inner surface of the grid cell;
wherein ,，/>，/>for the distance from the center of the nth grid cell to the center of the n+1th grid cell, +.>Is the distance from the center of the nth grid cell to the center of the n-1 th grid cell.
In a second aspect, an embodiment provides an anti-icing thermal load computing system comprising:
the grid dividing unit is used for dividing the grid of the solid wall object into N wall grid units, wherein N is a natural number greater than 1;
any grid cell anti-icing thermal load calculation unit, including,
an outer surface temperature obtaining subunit, configured to obtain respective outer surface temperatures of the grid unit and adjacent grid units;
the outer surface anti-icing heat load calculation subunit is used for calculating the outer surface anti-icing heat load of each grid unit and the adjacent grid units based on the latent heat absorbed by water drop evaporation, the convective heat exchange heat of the solid wall object and air and the energy of impinging the water drop;
a center temperature calculating subunit, configured to calculate respective center temperatures of the grid cells and adjacent grid cells thereof based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads;
an inner surface anti-icing thermal load calculation subunit for calculating an inner surface anti-icing thermal load of the grid unit based on the respective outer surface anti-icing thermal load and the respective center temperature;
and the anti-icing heat load acquisition unit is used for acquiring the anti-icing heat load of the inner surface of the grid unit as the anti-icing heat load of the grid unit.
In a third aspect, an embodiment provides a computer readable storage medium having stored therein a program capable of being loaded by a processor and executing any one of the above-described anti-icing thermal load calculation methods.
The beneficial effects of the application are as follows:
according to the method, the internal heat conduction of the solid wall object is considered, and the external temperature is converted into the central temperature, so that the anti-icing heat load of the external surface is calculated based on the internal heat conduction, the anti-icing heat load of the internal surface is used as the anti-icing heat load, the difficulty in acquiring the anti-icing heat load of the internal surface of the solid wall object is solved, and meanwhile, the calculated anti-icing heat load is more similar to the load required by a real anti-icing system.
Drawings
FIG. 1 is a schematic diagram of an anti-icing thermal load calculation principle according to an embodiment of the present application;
FIG. 2 is a schematic flow chart of an anti-icing thermal load calculation method according to an embodiment of the present application;
fig. 3 is a schematic diagram of an embodiment of the present application based on a grid cell anti-icing thermal load calculation principle.
Detailed Description
The application will be described in further detail below with reference to the drawings by means of specific embodiments. Wherein like elements in different embodiments are numbered alike in association. In the following embodiments, numerous specific details are set forth in order to provide a better understanding of the present application. However, one skilled in the art will readily recognize that some of the features may be omitted, or replaced by other elements, materials, or methods in different situations. In some instances, related operations of the present application have not been shown or described in the specification in order to avoid obscuring the core portions of the present application, and may be unnecessary to persons skilled in the art from a detailed description of the related operations, which may be presented in the description and general knowledge of one skilled in the art.
Furthermore, the described features, operations, or characteristics of the description may be combined in any suitable manner in various embodiments. Also, various steps or acts in the method descriptions may be interchanged or modified in a manner apparent to those of ordinary skill in the art. Thus, the various orders in the description and drawings are for clarity of description of only certain embodiments, and are not meant to be required orders unless otherwise indicated.
For convenience of explanation of the inventive concept of the present application, the following will briefly explain the aviation technology.
When an aircraft passes through a cloud layer containing supercooled water drops, icing phenomenon can occur on the windward surface of the aircraft, and the icing phenomenon can seriously damage the flight safety of the aircraft. It is therefore often necessary to design and install an anti-icing system at critical locations on the aircraft. The thermal anti-icing system has the advantages of high efficiency, good reliability and the like, and is mainly an anti-icing means adopted on aircraft and aeroengines at present. During the design process of the thermal anti-icing system, the anti-icing thermal load is calculated firstly, and then the designed anti-icing system is further checked and optimized. In the conventional method for calculating the anti-icing heat load, the anti-icing heat load of the outer surface (the surface in contact with the air) of the solid wall is calculated, and the actual anti-icing heat load is usually applied to the inner surface of the solid wall. When heat is transferred from the inner surface to the outer surface of the solid, heat conduction and heat transfer can occur in the solid, so that the anti-icing heat load obtained by heat conduction calculation in the solid wall surface is not considered to have errors compared with the load required by a real anti-icing system.
In view of this, the present application provides a method for calculating an anti-icing thermal load, please refer to fig. 1, in which the internal heat conduction of a solid wall object is considered, and the external surface anti-icing thermal load is converted into the central temperature, so that the internal surface anti-icing thermal load of the solid wall object is calculated based on the internal heat conduction, and the anti-icing thermal load of the internal surface is taken as the anti-icing thermal load, which solves the difficulty of obtaining the internal surface anti-icing thermal load of the solid wall object, and also makes the calculated anti-icing thermal load closer to the load required by a real anti-icing system. The following describes an example of an anti-icing load calculation for an aircraft in flight. The anti-icing heat load calculation can comprise dry anti-icing heat load calculation and wet anti-icing heat load calculation, the anti-icing heat load method can be only used for dry anti-icing heat load calculation, can be only used for wet anti-icing heat load calculation, can be also used for dry anti-icing heat load calculation and wet anti-icing heat load calculation, water collected on the surface of the dry anti-icing instruction can be completely evaporated, and when wet anti-icing, anti-icing heat can only ensure that icing phenomenon does not occur in a protection area, but liquid water still exists, and the anti-icing state of the anti-icing surface can be selected for calculation. The following description will be given separately.
Prior to performing the method, the prior art method (such as CFD method) is used to determine the design conditions (including the atmospheric pressure P of the environment 0 At an external ambient temperature T 0 Air inflow speed (or flying speed) U 0 Water droplet diameter MVD, and liquid water content LWC) to develop a supercooled water droplet flow field simulation of outside air,obtaining the convective heat transfer coefficient h and the water drop breakdown quality of each position of the solid wall surfaceThe anti-icing load calculation method of the present application is then performed.
In one embodiment, in the method for calculating the anti-icing load, the solid wall object (the outer wall surface of the aircraft) is divided into N wall grid cells, where N is a natural number greater than 1. The anti-icing heat load of each grid unit can be calculated, so that the anti-icing heat load of the whole aircraft outer wall surface can be obtained. Referring to fig. 2, for any one of the divided grid cells, the anti-icing thermal load calculation method includes:
step S10, acquiring the respective outer surface temperatures of the grid cells and the adjacent grid cells.
For dry anti-icing, in one embodiment, a method of acquiring the respective outer surface temperatures of a grid cell and its neighboring grid cells includes:
based on the formulaThe respective outer surface temperatures of the grid cells and their neighbors are calculated. Wherein T is s Is the outer surface temperature; />For evaporating quality, < >>；/>Impact mass for water droplets corresponding to the grid cells; b= 346.1; />，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is the constant, P 0 For the atmospheric pressure of the environment, the pressure can be determined byThe height of the environment or the detection result, T 0 The external environment temperature can be obtained through detection, and h is the convection heat exchange coefficient of the corresponding grid unit; h and->The method is obtained by performing simulation calculation on an air flow field and a water drop field outside the solid wall object.
Thus, we can obtain a dry anti-icing with an outside surface temperature T for any one grid cell sn And the outer surface temperature T of the adjacent grid unit of any grid unit s(n-1) and Ts(n+1) 。
For wet anti-icing, in one embodiment, the temperatures of the outer surfaces of each of the grid cells and their adjacent grid cells may each be a constant, which is set as the case may be, and which is typically greater than 0 o C is generally 2 to tens of degrees centigrade, and can be larger, and the larger the C is, the safer the C is, but the larger the energy consumption is.
Step S20, calculating the anti-icing thermal load of the outer surfaces of each grid cell and the adjacent grid cells based on the latent heat absorbed by the evaporation of the water droplets, the heat of convection heat exchange between the solid wall object and the air, and the energy of the impinging water droplets.
For dry anti-icing, in one embodiment, the formula is basedThe outer surface anti-icing heat load of each grid cell and each adjacent grid cell is calculated.
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, and is a constant; />Is the convection heat of solid wall object and air, +.>；For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flight speed.
Thus, we can obtain a dry anti-icing for any grid cell with an anti-icing thermal load on its outer surface ofAnd the outer surface of the adjacent grid cell of the arbitrary grid cell is anti-icing heat load +.> and />。
For wet anti-icing, in one embodiment, the formula is basedThe outer surface anti-icing heat load of each grid cell and each adjacent grid cell is calculated.
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, is constant, < ->For evaporation quality, < > in wet anti-icing->，，/>，，/>，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is constant, P 0 T is the atmospheric pressure of the environment 0 H is the convection heat exchange coefficient of the corresponding grid unit for the external environment temperature; />Is the convection heat of solid wall object and air, +.>；/>For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flying speed>Impact mass for water droplets corresponding to the grid cells; h and->The method is obtained by performing simulation calculation on an air flow field and a water drop field outside the solid wall object.
Thus, we can obtain a dry anti-icing for any grid cell with an anti-icing thermal load on its outer surface ofAnd the outer surface of the adjacent grid cell of the arbitrary grid cell is anti-icing heat load +.> and />。
Step S30, calculating the respective center temperatures of the grid cells and the adjacent grid cells based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads.
For dry anti-icing, please refer to fig. 3, in one embodiment, based on the formulaCalculating the temperature T of the center of the nth grid cell n Wherein N is more than or equal to 1 and less than or equal to N; />Anti-icing thermal load for the outer surface of the nth grid cell; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; t (T) sn Is the outer surface temperature of the nth grid cell.
Based on the formula and />Calculating a temperature T of a center corresponding to a grid cell adjacent to the grid cell n-1 and Tn+1 Wherein N is more than or equal to 1 and less than or equal to N; />An ice and heat load is prevented for the outer surface of the n-1 th grid cell adjacent to the n-th cell; />An outer surface of the (n+1) th grid cell adjacent to the (n) th cell is protected from an ice and heat load; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; /> and />Is the temperature of the outer surface corresponding to the grid cell adjacent to the grid cell.
For wet anti-icing, in one embodiment, as such, based on the formulaCalculating the temperature T of the center of the nth grid cell n Wherein N is more than or equal to 1 and less than or equal to N; />Anti-icing thermal load for the outer surface of the nth grid cell; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; t (T) sn Is the outer surface temperature of the nth grid cell.
Based on the formula and />Calculating a temperature T of a center corresponding to a grid cell adjacent to the grid cell n-1 and Tn+1 Wherein N is more than or equal to 1 and less than or equal to N; />An ice and heat load is prevented for the outer surface of the n-1 th grid cell adjacent to the n-th cell; />An outer surface of the (n+1) th grid cell adjacent to the (n) th cell is protected from an ice and heat load; d is the thickness of the solid wall object; k is a solid wall objectIs a physical property constant; /> and />Is the temperature of the outer surface corresponding to the grid cell adjacent to the grid cell.
In FIG. 3, taking the current grid cell to be calculated as a #2 grid cell as an example, the outside surface temperature T of the #2 grid cell is calculated s2 And an outer surface anti-icing thermal loadI.e. according to the formula->Obtaining the center temperature T of the #2 grid cell 2 Calculating the outer surface temperature of each of #1 and #3 grid cells adjacent to the #2 grid cell>Andan external surface anti-icing thermal load-> and />Then it can be according to the formula +.> and />Obtaining the respective center temperatures T of the #1 and #3 grid cells 1 and T3 . With this, the anti-icing thermal load of the inner surface of the #2 grid cell can be further obtained.
Step S40, calculating the inner surface anti-icing thermal load of the grid cell based on the respective outer surface anti-icing thermal load and the respective center temperature.
For dry anti-icing, in one embodiment, based on energy conservation, the formula is followedThe inner surface anti-icing thermal load of the grid cell is calculated.
wherein ,，/>，/>for the distance from the center of the nth grid cell to the center of the n+1th grid cell, +.>Is the distance from the center of the nth grid cell to the center of the n-1 th grid cell.
For wet anti-icing, in one embodiment, as such, based on energy conservation, according to the formulaThe inner surface anti-icing thermal load of the grid cell is calculated.
wherein ,，/>，/>for the distance from the center of the nth grid cell to the center of the n+1th grid cell, +.>Is the distance from the center of the nth grid cell to the center of the n-1 th grid cell.
And S50, taking the anti-icing heat load of the inner surface of the grid unit as the anti-icing heat load of the grid unit.
The calculated dry and/or wet ice-proof inner surface ice-proof thermal load is outputted as the dry and/or wet ice-proof thermal load of the grid cell.
Based on the above, considering the internal heat conduction of the solid wall object, the external surface anti-icing heat load is calculated based on the internal heat conduction by converting the external surface temperature into the central temperature, and the anti-icing heat load of the internal surface is used as the anti-icing heat load, so that the difficulty of acquiring the anti-icing heat load of the internal surface of the solid wall object is solved, and the calculated anti-icing heat load is also more similar to the load required by a real anti-icing system.
It should be noted that the temperature in the formula of the above embodiment of the present application is calculated using kelvin, and it will be understood that those skilled in the art may make equivalent transformations if calculating using fahrenheit or celsius instead.
In one embodiment of the application, an anti-icing load system is provided comprising:
and the grid dividing unit is used for dividing the grid of the solid wall object into N wall grid units, wherein N is a natural number greater than 1.
The anti-icing thermal load calculation unit for any grid cell, including,
and the outer surface temperature acquisition subunit is used for acquiring the outer surface temperature of each grid unit and the adjacent grid units thereof.
And the outer surface anti-icing heat load calculating subunit is used for calculating the anti-icing heat load of each outer surface of the grid unit and the adjacent grid units based on the latent heat absorbed by the evaporation of the water drops, the convective heat exchange heat between the solid wall object and the air and the energy of the impinging water drops.
And the central temperature calculating subunit is used for calculating the respective central temperatures of the grid units and the adjacent grid units based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads.
And the inner surface anti-icing heat load calculation subunit is used for calculating the inner surface anti-icing heat load of the grid unit based on the respective outer surface anti-icing heat load and the respective center temperature.
And the anti-icing heat load acquisition unit is used for acquiring the anti-icing heat load of the inner surface of the grid unit as the anti-icing heat load of the grid unit.
In one embodiment of the present application, a computer-readable storage medium is provided, in which a program is stored, the program being capable of being loaded by a processor and executing the above-described anti-icing thermal load calculation method.
Those skilled in the art will appreciate that all or part of the functions of the various methods in the above embodiments may be implemented by hardware, or may be implemented by a computer program. When all or part of the functions in the above embodiments are implemented by means of a computer program, the program may be stored in a computer readable storage medium, and the storage medium may include: read-only memory, random access memory, magnetic disk, optical disk, hard disk, etc., and the program is executed by a computer to realize the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor, all or part of the functions described above can be realized. In addition, when all or part of the functions in the above embodiments are implemented by means of a computer program, the program may be stored in a storage medium such as a server, another computer, a magnetic disk, an optical disk, a flash disk, or a removable hard disk, and the program in the above embodiments may be implemented by downloading or copying the program into a memory of a local device or updating a version of a system of the local device, and when the program in the memory is executed by a processor.
The foregoing description of the application has been presented for purposes of illustration and description, and is not intended to be limiting. Several simple deductions, modifications or substitutions may also be made by a person skilled in the art to which the application pertains, based on the idea of the application.
Claims (8)
1. The anti-icing heat load calculation method is characterized by comprising the steps of dividing a solid wall object into N wall grid units, wherein N is a natural number larger than 1; for any one of the divided grid cells, the anti-icing thermal load calculation method includes:
acquiring the respective outer surface temperatures of the grid cells and the adjacent grid cells;
calculating the anti-icing heat load of the outer surfaces of each grid unit and the adjacent grid units based on the latent heat absorbed by evaporation of water drops, the heat of convection heat exchange between a solid wall object and air and the energy of impinging water drops;
calculating respective center temperatures of the grid cells and adjacent grid cells thereof based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads;
calculating an inner surface anti-icing thermal load of the grid cells based on the respective outer surface anti-icing thermal load and the respective center temperature;
and taking the anti-icing thermal load of the inner surface of the grid unit as the anti-icing thermal load of the grid unit.
2. The anti-icing thermal load calculation method of claim 1 wherein said anti-icing thermal load comprises a dry anti-icing thermal load and/or a wet anti-icing thermal load.
3. The method of calculating an anti-icing thermal load as recited in claim 2 wherein said obtaining the respective outer surface temperatures of said grid cells and their neighbors comprises:
for dry anti-icing, the formula is basedCalculating the respective outer surface temperatures of the grid cells and adjacent grid cells, wherein T s Is the outer surface temperature; />For evaporating quality, < >>；/>Impact mass for water droplets corresponding to the grid cells; b= 346.1; />，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is the constant, P 0 T is the atmospheric pressure of the environment 0 H is the convection heat exchange coefficient of the corresponding grid unit for the external environment temperature; h and->The method comprises the steps of performing simulation calculation on an air flow field and a water drop field outside a solid wall object;
and/or the number of the groups of groups,
for wet anti-icing, the grid cells and their adjacent grid cells each have an outer surface temperatureIs a constant.
4. The method of calculating the anti-icing thermal load according to claim 3, wherein the calculating the anti-icing thermal load of the outer surface of each of the mesh units and the adjacent mesh units based on the heat of convection between the solid wall object and the air by the evaporation and absorption of latent heat of the water droplets and the energy of the impinging water droplets comprises:
for dry anti-icing, the formula is basedCalculating the anti-icing heat load of the outer surfaces of each grid unit and each adjacent grid unit;
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, and is a constant; />Is the convection heat of solid wall object and air, +.>；/>For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flying speed;
for wet anti-icing, the formula is basedCalculating the anti-icing heat load of the outer surfaces of each grid unit and each adjacent grid unit;
wherein ,an external surface is protected from ice and heat load; />Absorb latent heat for water drop evaporation, < >>Le is the latent heat of evaporation of the water droplets, is constant, < ->For evaporating quality, < >>，，/>，，/>，/>Is the constant of specific heat capacity of air, lewis is the Lewis number of air, is constant, P 0 T is the atmospheric pressure of the environment 0 H is the convection heat exchange coefficient of the corresponding grid unit for the external environment temperature; />Is the convection heat of solid wall object and air, +.>；/>For the energy of striking the water drops +.>，C pw Is the specific heat capacity of water, is a constant, U 0 Is the flying speed>Impact mass for water droplets corresponding to the grid cells; h and->All according to the pairAnd (5) performing simulation calculation on the air flow field and the water drop field outside the solid wall object.
5. The method of calculating the anti-icing thermal load according to claim 4, wherein calculating the respective center temperatures of the grid cells and the adjacent grid cells based on the respective outer surface temperatures and the respective outer surface anti-icing thermal loads comprises:
for dry and/or wet anti-icing, the formula is basedCalculating the temperature T of the center of the nth grid cell n Wherein N is more than or equal to 1 and less than or equal to N; />Anti-icing thermal load for the outer surface of the nth grid cell; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; t (T) sn An outer surface temperature for the nth grid cell;
based on the formula and />Calculating a temperature T of a center corresponding to a grid cell adjacent to the grid cell n-1 and Tn+1 Wherein N is more than or equal to 1 and less than or equal to N; />An ice and heat load is prevented for the outer surface of the n-1 th grid cell adjacent to the n-th cell; />An outer surface of the (n+1) th grid cell adjacent to the (n) th cell is protected from an ice and heat load; d is the thickness of the solid wall object; k is the heat conductivity coefficient of the solid wall object and is a physical property constant; /> and />Is the temperature of the outer surface corresponding to the grid cell adjacent to the grid cell.
6. The method of calculating an anti-icing thermal load as recited in claim 5 wherein said calculating an inner surface anti-icing thermal load of said grid cells based on said respective outer surface anti-icing thermal load and respective center temperature comprises:
for dry and/or wet anti-icing, based on energy conservation, according to the formulaCalculating an anti-icing thermal load of the inner surface of the grid cell;
wherein ,，/>，/>for the distance from the center of the nth grid cell to the center of the n+1th grid cell, +.>Is the distance from the center of the nth grid cell to the center of the n-1 th grid cell.
7. An anti-icing thermal load computing system comprising:
the grid dividing unit is used for dividing the grid of the solid wall object into N wall grid units, wherein N is a natural number greater than 1;
any grid cell anti-icing thermal load calculation unit, including,
an outer surface temperature obtaining subunit, configured to obtain respective outer surface temperatures of the grid unit and adjacent grid units;
the outer surface anti-icing heat load calculation subunit is used for calculating the outer surface anti-icing heat load of each grid unit and the adjacent grid units based on the latent heat absorbed by water drop evaporation, the convective heat exchange heat of the solid wall object and air and the energy of impinging the water drop;
a center temperature calculating subunit, configured to calculate respective center temperatures of the grid cells and adjacent grid cells thereof based on the respective outer surface temperatures and the respective outer surface anti-icing heat loads;
an inner surface anti-icing thermal load calculation subunit for calculating an inner surface anti-icing thermal load of the grid unit based on the respective outer surface anti-icing thermal load and the respective center temperature;
and the anti-icing heat load acquisition unit is used for acquiring the anti-icing heat load of the inner surface of the grid unit as the anti-icing heat load of the grid unit.
8. A computer-readable storage medium, wherein a program is stored in the medium, the program being capable of being loaded by a processor and executing the anti-icing heat load calculation method according to one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310883924.XA CN116644258B (en) | 2023-07-19 | 2023-07-19 | Anti-icing heat load calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310883924.XA CN116644258B (en) | 2023-07-19 | 2023-07-19 | Anti-icing heat load calculation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116644258A true CN116644258A (en) | 2023-08-25 |
| CN116644258B CN116644258B (en) | 2023-12-22 |
Family
ID=87615663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310883924.XA Active CN116644258B (en) | 2023-07-19 | 2023-07-19 | Anti-icing heat load calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116644258B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118607416A (en) * | 2024-08-05 | 2024-09-06 | 中国空气动力研究与发展中心低速空气动力研究所 | Design method and medium for anti-icing heat flow density of electric heating anti-icing dry surface |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08247590A (en) * | 1995-03-08 | 1996-09-27 | Toshiba Corp | Ice heat storage apparatus |
| JPH0949650A (en) * | 1995-08-04 | 1997-02-18 | Shinryo Corp | Thawing system for ice heat storage |
| CN105035355A (en) * | 2015-08-12 | 2015-11-11 | 中国航空工业集团公司西安飞机设计研究所 | Thermal load power simulation method and thermal load simulation system |
| CN107792377A (en) * | 2016-09-02 | 2018-03-13 | 埃姆普里萨有限公司 | Aircraft air-conditioning system air stream adjusts |
| CN109696290A (en) * | 2019-02-18 | 2019-04-30 | 中国空气动力研究与发展中心低速空气动力研究所 | A kind of wind energy conversion system wing panel Anti-icing Heat Loads demand measuring system |
| CN112989727A (en) * | 2021-05-10 | 2021-06-18 | 中国空气动力研究与发展中心低速空气动力研究所 | Wall surface temperature simulation method of anti-icing system |
| CN113959593A (en) * | 2021-10-22 | 2022-01-21 | 中国航发沈阳发动机研究所 | Method for solving surface temperature of anti-icing component |
| CN114413536A (en) * | 2021-12-15 | 2022-04-29 | 合肥通用机械研究院有限公司 | Heat load calculation method, detector, refrigeration control method and system for ice rink |
| CN114940266A (en) * | 2021-12-31 | 2022-08-26 | 成都流体动力创新中心 | Method and system for predicting skin surface temperature capable of being maintained by complex anti-icing cavity |
| CN115587506A (en) * | 2022-12-09 | 2023-01-10 | 四川大学 | Design method of electric heating ice preventing and removing system |
| CN115817822A (en) * | 2023-02-09 | 2023-03-21 | 中国空气动力研究与发展中心低速空气动力研究所 | Heat load distribution design method of electric heating anti-icing system |
-
2023
- 2023-07-19 CN CN202310883924.XA patent/CN116644258B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08247590A (en) * | 1995-03-08 | 1996-09-27 | Toshiba Corp | Ice heat storage apparatus |
| JPH0949650A (en) * | 1995-08-04 | 1997-02-18 | Shinryo Corp | Thawing system for ice heat storage |
| CN105035355A (en) * | 2015-08-12 | 2015-11-11 | 中国航空工业集团公司西安飞机设计研究所 | Thermal load power simulation method and thermal load simulation system |
| CN107792377A (en) * | 2016-09-02 | 2018-03-13 | 埃姆普里萨有限公司 | Aircraft air-conditioning system air stream adjusts |
| CN109696290A (en) * | 2019-02-18 | 2019-04-30 | 中国空气动力研究与发展中心低速空气动力研究所 | A kind of wind energy conversion system wing panel Anti-icing Heat Loads demand measuring system |
| CN112989727A (en) * | 2021-05-10 | 2021-06-18 | 中国空气动力研究与发展中心低速空气动力研究所 | Wall surface temperature simulation method of anti-icing system |
| CN113959593A (en) * | 2021-10-22 | 2022-01-21 | 中国航发沈阳发动机研究所 | Method for solving surface temperature of anti-icing component |
| CN114413536A (en) * | 2021-12-15 | 2022-04-29 | 合肥通用机械研究院有限公司 | Heat load calculation method, detector, refrigeration control method and system for ice rink |
| CN114940266A (en) * | 2021-12-31 | 2022-08-26 | 成都流体动力创新中心 | Method and system for predicting skin surface temperature capable of being maintained by complex anti-icing cavity |
| CN115587506A (en) * | 2022-12-09 | 2023-01-10 | 四川大学 | Design method of electric heating ice preventing and removing system |
| CN115817822A (en) * | 2023-02-09 | 2023-03-21 | 中国空气动力研究与发展中心低速空气动力研究所 | Heat load distribution design method of electric heating anti-icing system |
Non-Patent Citations (5)
| Title |
|---|
| NINGLI CHEN 等: "A mathematical model based on unstructured mesh for ice accretion", 《AIP ADVANCES》, pages 125149 - 1 * |
| NINGLI CHEN 等: "A three-dimensional mathematical model for simulating ice accretion on helicopter rotors", 《PHYSICS OF FLUIDS》, pages 1 - 12 * |
| RICHARD HANN 等: "Experimental Heal Loads for Electrothermal anti-icing and De-icing on UAVs", 《AEROSPACE》, pages 1 - 15 * |
| 朱光亚: "电加热防冰部件加热功率的分布特性研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, pages 031 - 212 * |
| 肖春华: "飞机电热除冰过程的传热特性及其影响研究", 《中国博士学位论文全文数据库 工程科技II辑》, pages 031 - 3 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118607416A (en) * | 2024-08-05 | 2024-09-06 | 中国空气动力研究与发展中心低速空气动力研究所 | Design method and medium for anti-icing heat flow density of electric heating anti-icing dry surface |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116644258B (en) | 2023-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yirtici et al. | Predictions of ice formations on wind turbine blades and power production losses due to icing | |
| Wu et al. | Study on the effect of shape-stabilized phase change materials on spacecraft thermal control in extreme thermal environment | |
| CN116644258B (en) | Anti-icing heat load calculation method and system | |
| Shen et al. | Unsteady simulation of aircraft electro-thermal deicing process with temperature-based method | |
| Pourbagian et al. | Constrained problem formulations for power optimization of aircraft electro-thermal anti-icing systems | |
| Pourbagian et al. | Surrogate-based optimization of electrothermal wing anti-icing systems | |
| Jiang et al. | Model for ice wet growth on composite insulator and its experimental validation | |
| Gori et al. | Modeling in-flight ice accretion under uncertain conditions | |
| CN115817822B (en) | Heat load distribution design method of electric heating anti-icing system | |
| Zhu et al. | Numerical simulation and experimental verification of the airfoil electrothermal deicing system performance | |
| Fouladi et al. | Quasi-steady modeling of ice accretion on a helicopter fuselage in forward flight | |
| Bu et al. | Numerical analysis of a swept wing hot air ice protection system | |
| Zhang et al. | Structural Parameter Sensitivity Analysis of an Aircraft Anti‐Icing Cavity Based on Thermal Efficiency | |
| Gallia et al. | Robust Optimization of a Thermal Anti-Ice Protection System in Uncertain Cloud Conditions | |
| CN118607416B (en) | Design method and medium for anti-icing heat flow density of electric heating anti-icing dry surface | |
| Barzi | Multidisciplinary optimization of in-flight electro-thermal ice protection systems | |
| Ibrahim et al. | Transient atmospheric ice accretion on wind turbine blades | |
| Wang et al. | Numerical Simulation of Ice Crystal Supercooled Droplet Mixed Phase Icing Based on the Improved Messinger Model | |
| Ma et al. | A novel technique for predicting the thermal behavior of stratospheric balloon | |
| CN116611266B (en) | Unsteady state electric heating anti-icing system simulation method | |
| Pourbagian et al. | Power and design optimization of electro-thermal anti-icing systems via FENSAP-ICE | |
| Zhang et al. | Measured moisture accumulation in aircraft walls during simulated commercial flights | |
| US20220297819A1 (en) | Hybrid surface designs with passive diffusion-limited deposition capabilities | |
| Szilder et al. | Numerical simulation of ice accretion influenced by anti-icing heating system | |
| Coppinger | Icing in the Electric Age |
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 |