• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Numerical Simulation of Trace Water Icing Characteristics on Airfoil Under Cryogenic Conditions

    2023-05-21 02:37:22,*,,,

    ,*,,,

    1.School of Energy and Power Engineering,Xi’an Jiaotong University,Xi’an 710049,P.R.China;

    2.Facility Design and Instrumentation Institute,China Aerodynamics Research and Development Center,Mianyang,621000,P.R.China

    Abstract: Compared with conventional wind tunnels,cryogenic wind tunnels can supply testing airflow with higher Reynolds numbers.However,the residual trace water vapor in the cryogenic wind tunnel may affect the accuracy of aerodynamic data at low temperatures.Due to the difficulty of experimental methods,the characteristics of trace water icing still need to be revealed.Based on the Eulerian method,a multiphase flow model of wet air,subcooled droplets and ice crystals is established in this paper.The icing process caused by subcooled droplets and ice crystals condensed from water vapor on the airfoil under the condition of cryogenic wind tunnels is simulated.The results show good consistency with experimental data.The effects of the total water content,melt ratio,airflow temperature,pressure,Mach number and other parameters on icing process are analyzed,and the effects of the airfoil icing on aerodynamic characteristics are also quantified,which may provide a theoretical reference for the actual operation of cryogenic wind tunnels.

    Key words:cryogenic wind tunnel;trace water;airfoil icing;numerical simulation

    0 Introduction

    The wind tunnel is an indispensable equipment to study the aerodynamic characteristics of complex viscous flow field around the test article.In the thriving aerospace industry,the research and development of aircraft and space vehicles put forward urgent requirements for the capability of wind tunnels to operate at high Reynolds numbers.It is difficult for conventional wind tunnels to realize a further increase of the Reynolds number.Problems such as power consumption and structural strength may be caused by increasing the tunnel size or operation power.At present,it is considered that the most effective way is to cool down the gas medium,reduce the viscosity and increase the density,thus to achieve much a higher Reynolds number and a lower power consumption.Liquid nitrogen is sprayed into the tunnel to cool down the gas,and the internal temperature ranges from 110 to 320 K[1].The water content of 1×10-6corresponds to the dew point of 197 K.When the wind tunnel is running at the lowest operating temperature of 110 K,the remaining water in the airflow will still exist in the form of ice crystals or supercooled droplets,which may result in frost on the test article or equipment.In order to achieve an ultra-low dew point and prevent water vapor from condensation or desublimation into droplets or ice crystals,the water content in the airflow should be reduced to less than 1×10-6before experiments,leading to much more stringent operating conditions of the wind tunnel.Therefore,mastering the icing mechanism of water vapor in the cryogenic wind tunnels and quantifying the frost formation could provide a technical support for optimizing the operating process of cryogenic wind tunnels,ensuring the accuracy of aerodynamic data and preventing damage from instruments.

    The mechanism of airfoil icing has been studied through the ground test or the numerical simulation.The ground test requires the support of icing wind tunnel system and has the disadvantages of difficulties in simulating the meteorology conditions in the atmosphere,high cost,etc.Recently,the numerical simulation has become an important method for its low cost in manpower,material resources and time.

    In experimental studies,experiments on icing under mixed-phase or glaciated conditions have been carried out with the help of wind tunnels.In order to study the effect of icing conditions on the performance of thermal ice protection systems and ice crystals’ behavior on the airfoil,Al-Khalil et al.[2]conducted experiments of mixed-phase and glaciated icing at different temperatures and water contents.They divided the icing conditions into glaze condition and rime condition according to geometric features and texture of ice accretion.They also found that ice crystal erosion could significantly weaken the ice accretion process.Mason et al.[3]experimentally simulated the warm air or cold ice conditions of glaciated icing in jet engines.The results obtained from operations under different airflow temperatures illustrated that ice was only formed on the airfoil when liquid water existed,i.e.under mixed-phase conditions.Currie et al.[4]experimentally quantified the mixed-phase sticking efficiency for ice crystal accretion.The results showed that the sticking efficiency was highly correlated with the ratio of liquid water content(LWC)to total water content(TWC)in the airflow,reaching a maximum value when the melt ratio fell in a specific range.The reason was that sufficient ice crystals were needed to cool down the airfoil surface below the freezing point,and liquid water should be present on the airfoil for ice crystals to effectively adhere to the surface.In addition,the erosion effects of ice crystals will narrow down the range.Baumert et al.[5]carried out experimental investigations into mixed-phase icing at the Icing Wind Tunnel Braunschweig.According to the data,the ambient temperature influenced the icing process to a great extent,and the sticking efficiency of ice particles was dependent on the melt ratio,which agreed with the previous literature[4].

    For the numerical simulation,several models have been proposed to predict the physical behaviors in the icing process.Villedieu et al.[6]testified the capabilities of IGLOO2D to simulate the glaciated and mixed-phase ice accretion.The software considers some phenomena ignored before,such as heat transfer between air and particles,evaporation and melting,particles’ behavior when impinging dry and wet surfaces,etc.The authors simplified the movement and melting process of ice crystals into three stages and created governing equations based on the thermodynamic characteristics of the stages.They defined the particle sphericity to represent the irregularity of the particle shape,and improved the drag coefficient model and the heat transfer model.Considering the effect of erosion caused by ice crystals,Trontin et al.[7]improved the ice crystal sticking efficiency model,and corrected the empirical parameters according to the experimental findings of Currie et al.[4]Then,the numerical simulations of supercooled large droplets,ice crystals and mixing conditions were carried out.The proposed model could well predict the ice shape under a wide range of icing conditions,and could reproduce the distinctive feature of icing plateau phenomenon.Bu et al.[8]extended the Messinger’s model[9]to the mixed-phase icing applications by analyzing the heat and mass transfer on the airfoil under mixed-phase conditions.A sticking efficiency correlation of ice particles was also added to the classical model to take into account the sticking effects of ice particles.As for numerical analyses into the mechanism,Qi et al.[10]studied the influence of temperature and content of water vapor on the desublimation process under constant Mach and Reynolds number conditions.The results showed that the maximum nucleation rate and maximum ice mass fraction increase with the temperature and water vapor content in the freestream,and the water vapor content had a more significant effect on it.The authors also categorized the desublimation process of water vapor into two patterns.The difference was whether the water vapor desublimated or not in the upstream of the airfoil.

    Previous studies have been focused on airfoil icing in the atmosphere.In order to solve the “moisture contamination” problem in cryogenic wind tunnels,it is also necessary to investigate the icing process of trace water under cryogenic conditions.In this study,the icing characteristics of trace water on airfoil in cryogenic wind tunnels are discussed by means of numerical simulation.The validated model is used to explore the icing process considering the effects of total water content,melt ratio,airflow temperature and pressure,Mach number,surface roughness and icing duration.

    1 Numerical Method

    1.1 Ice accretion model

    On the airfoil surface,a liquid film may be formed due to impinging droplets.The droplet may freeze immediately after impingement,or part of them may freeze,the other part of the water gradually freezes while flowing to nearby surfaces.The film may run back,forced by the airflow shear stress,centrifugal force,or gravity.Based on the classical Messinger model,the liquid film icing process is modeled[11].Fig.1 shows the heat and mass transfer phenomena modeled when solving for ice accretion and runback.

    Fig.1 Water film icing model on a solid surface

    Based on the operating parameters and surface conditions,part of the film may freeze,evaporate or sublimate.The height of the liquid film can be determined at all grid points on the solid surface,and the velocity is a function of the coordinatesx=(x1,x2) on the surface andynormal to the surface.The problem is simplified by introducing a linear profile for the film velocityVf(x,y) normal to the wall,with zero velocity imposed at the wall

    whereτa,wall(x) is the air shear stress,which is the main driving force for the film;andμfthe film dynamic viscosity.By averaging the film velocity across the film thickness,a mean velocity can be derived as

    wherehfis the film height.

    Two PDEs are solved on all solid surfaces.

    The first one is the mass conservation equation

    whereρf,,represent the film density,evaporative mass flux and instantaneous mass accumulation of ice in the control volume,respectively.The three terms on the right-hand side correspond to the mass transfer by water droplet impingement,evaporation and ice accretion.The second partial differential equation(PDE)is the energy conservation equation

    wherecfandcscorrespond to the specific heat capacity of water film and ice,respectively;are the equilibrium temperature and the recovery temperature,respectively;LevapandLfusionthe latent heat of evaporation and fusion,respectively;Qanti-icingis the anti-icing heat flux.The first three terms on the right-hand side correspond to the heat transfer generated by the impinging supercooled water droplets,evaporation and ice accretion,respectively.The last three terms are the radiative,convective and conductive heat fluxes.

    Compatibility relations are needed to solve the equations set

    Eq.(5)ensures that there is no liquid film when the equilibrium temperature is below the freezing point,and that no ice is formed if the film is above the freezing point.

    1.2 Ice roughness model

    Ice roughness will affect the flow pattern and heat transfer on the surface,which may change the morphology of the ice accretion layer[12].Experimental evidence shows that there are ice beads existing on the solid surface[13].They will keep stationary temporarily to grow or coalesce into larger beads.

    The roughness distribution is calculated by a beading model[14].It is assumed that approximately 85% of the local grid facet is covered by the ice beads,and the equivalent sand-grain roughness height is equal to the bead height.The local critical bead height is calculated using the shear stress and pressure gradient data from the airflow solution.When the bead height reaches the critical value,in wet zones,the bead will stop growing and run back through the grid cells.In dry zones,the bead will completely freeze and stop growing.

    The critical bead height is computed from the force analysis.Neglecting the internal pressure,a bead is under the action of the gravitational force component parallel to the surfaceFg,the component of surface tension along the moving axisFσ,the pressure drag forceFp,and the friction drag forceFτ.The balance of the forces is given by

    Once the balance is disturbed,the bead starts moving,and the bead height at this time is defined as the critical height.Every term in Eq.(6)is expressed as the polynomial function of the bead heighthb,and the smallest positive root of the equation is the local critical bead heighthb,cr.

    The total wetted areaAweton the grid facet is expressed as

    whereAbandNbrepresent the wetted area of a single bead and the number of beads,respectively.Nbis calculated by

    whereMbeadlayer,Vb,ρbare the total mass of bead layer,the volume and density of a bead,respectively.The initial value ofAwetis zero,and it will grow larger with time.WhenAwetreaches the maximum ratio 85%,hbwill increase due to the continuously impinging droplets and the coalescence of small beads,andNbwill decrease.Nbwill be recalculated by Eq.(8).Ifhb≥hb,cr,the beads in that grid are assumed to stop growing.Ifhb<hb,cr,the bead height is calculated withMbeadlayer,Nbandρb.

    1.3 Calculation flow and algorithm

    The governing equations above are discretized by the finite element method and solved in the software Ansys Fluent Icing.The turbulence model is Spalart-Allmaras,which is designed specifically for aeronautics and aerospace applications involving wall-bounded flows.It is shown to generate good results for boundary layers subjected to adverse pressure gradients.An approach based on unidirectional coupling Eulerian two-phase method is developed to compute the airflow and the movement of droplets and ice crystals.The temperature and velocity of droplets and crystals are assumed to be equal to that of the airflow.Due to the minuteness and sparsity of the particles,it is considered that they do not affect the macroscopic flow field.

    The icing process is calculated with a multiple quasi-steady approach,i.e.multi-shot,as shown in Fig.2.A multi-shot run executes four solvers in sequence.First,the impingement of droplets and ice crystals is computed in the module DROP3D.Second,the results obtained such as droplet collection coefficient and crystal sticking fraction will be imported to the next module ICE3D to perform ice accretion calculation.Thrid,the ice thickness and surface roughness data are used for mesh update and flow calculation in the ensuing shot.The entire process is repeated until the icing duration reaches the preset value.This approach is designed to simulate the evolving flow field as accurately as possible.Therefore,the number of shots depends on the influence of ice on the flow field.

    Fig.2 Process of multi-shot

    1.4 Computational domain and grid generation

    According to the profile data of NASA SC(2)-0610 airfoil,a two-dimensional computational domain with a pressure-far-field boundary condition is established,as shown in Fig.3.C-shaped structured mesh is chosen in the meshing.The chord length of the airfoil is 0.1 m.The mesh is refined near the airfoil,as shown in Fig.4.Non-slip boundary condition is adopted on the airfoil surface.

    Fig.3 Profile of airfoil and computational domain

    Fig.4 Mesh details

    1.5 Parameter setup

    Due to the extremely low temperature in the cryogenic tunnel,the water content is proved to be supersaturated for most working conditions in this study.The fractional pressure of water vapor is assumed to be equal to the saturated vapor pressure.Since the airflow velocity is high,the stagnation process is assumed to be adiabatic.After analyzing the thermodynamic process which vapor and ice crystals undergo,it is assumed that the vapor would not condense or desublimate and ice crystals would not melt in all working conditions in this paper.Hereby,only the icing process of droplets and ice crystals is carried out in the simulation.

    It is generally considered that ice crystals will not frost on dry and cold surfaces without liquid water[3].However,frost on the airfoil observed during the operation of the cryogenic wind tunnel at about 200 K indicates the possible existence of liquid water in the airflow under such extreme conditions.Manka et al.[15]found that a small fraction of the nanosized droplets remained unfrozen far below 200 K in an experimental determination of quenched droplets.Sellberg et al.[16]examined the structure factors of water in micron-scale droplets cooled to about 227 K.Even under such conditions,a few droplets(100—3 600)remained liquid over a time span of a few milliseconds.In addition,most of the temperature range of cryogenic wind tunnels falls in the region of “Supercooled liquid water”,“No man’s land” or “Ultra-viscous liquid” in the phase diagram of non-crystalline water[17],as shown in Fig.5.Therefore,it can be inferred that deep supercooled liquid water may still exist in the cryogenic wind tunnel.As a preliminary attempt,the simulation method of mixed-phase icing can be adopted to study the icing process in the cryogenic wind tunnel,assuming that LWC≈0.03(TWC-VWC)and ice crystal content(ICC)=TWC-VWC-LWC in this paper,where VWC represents the vapor water content.The ratio 0.03 is derived from Sellberg’s experimental data.

    Fig.5 Phase diagram of non-crystalline water[17]

    Parameters of the reference condition are shown in Table 1,where the unit of TWC “1×10-6”represents the mass fraction of water in air.In each experiment,only one parameter is variable,while other parameters remain unchanged.

    Table 1 Parameters of the reference condition

    2 Results and Discussion

    2.1 Grid independence validation

    Meshes of different grid density are used to calculate the aerodynamic data of the original NACA0012 airfoil.As can be seen from Fig.6,the mesh with 370 000 grids is suitable for calculation with drag coefficient variation less than 1%.

    Fig.6 Grid independence validation

    In the icing calculation,different numbers of shots are carried out to ensure that it has little effect on the results.As shown in Fig.7,Run 19 is taken as an example,the ice shape of 2 shots is approximate to that of 3 shots,and the mass deviation of the accumulated ice is about 1%.Therefore,it can be considered that the calculation of 2 shots is appropriate.

    Fig.7 Results of different number of shots (Run 19)

    2.2 Model validation

    An icing experiment is conducted on NACA0012 airfoil with chord length of 91.44 cm in the Cox &Co.Icing Tunnel in Ref.[2].According to two groups of parameters in Ref.[2],i.e.Run 10 and Run 19,as shown in Table 2,the same parame-ters are used in the software to simulate the icing process.

    As for Run 10 and Run 19,numerical results of 3 shots and 2 shots are used to compare with the experimental data respectively.In Fig.8,it is shown that the calculated results are in good agreement with the experimental data,with average deviations of 22% and 15%,respectively,indicating that the numerical model has a relative high accuracy.

    Fig.8 Comparison of ice shape

    2.3 Effects of TWC

    Figs.9(a,b)show the profiles of ice accretion layers under different TWC.For TWC between 0.1×10-6to 10×10-6,the type of the ice accreted is rime ice.This is characterized by regular ice shape and ice accumulation at the leading edge,which has little disturbance on the aerodynamic characteristics of the airfoil.When TWC≥50×10-6,the ice shape becomes irregular,the type changes from rime ice to mixed ice with obvious concave and protruding geometric features,and the boundary of the ice extends to the middle part of the airfoil.The characteristic of mixed ice is that,due to insufficient heat transfer near the leading edge,glaze ice with a higher density will occur.It is not easy to be eroded,so the profile is smooth.In the middle of the airfoil where heat transfer is more sufficient,the loose-textured and low-density rime ice will occur,which is easy to be eroded by ice crystals.Therefore,the local profile is not smooth on microscopic level.In addition,there is ice near the trailing edge in each working condition.The local ice thickness under 200×10-6condition is about 0.2 mm.

    In case of TWC=200×10-6,there is an untypical ice accretion in a mushroom shape with a sharp edge in the middle of the accretion.This accretion can be regarded as the cone-shaped product of continuous condensation of droplets on the basis of hemispherical ice accretion with TWC=50×10-6—100×10-6.According to the erosion model established by Trontin et al.[7],it is assumed that the erosion mainly depends on the tangential velocity of the impinging ice crystals.Therefore,more serious the erosion appears when the ice accretion is further away from the center line,so the final ice accretion is in the shape of a cone.Currie et al.[18]conducted an ice crystal icing experiment on a hemisphere,observing the ice accretion of a similar shape.It is considered that less liquid water remains in the accretion layer under lower temperatures.The harder the ice accretion is,the more likely the sharp edge is to resist erosion caused by ice crystals.Thus,the mushroom-shaped base can be formed.Because of gravity,the ice accretion layer slightly droops.

    Fig.9(c)shows a linear dependency between the total mass of ice and TWC.Apparently,the reason lies in the increasing quantity of impinging droplets and crystals.

    Fig.9 Ice thickness distribution and total mass of ice at different TWC

    As shown in Fig.10,with the increase of TWC,the drag coefficient of the icing airfoil increases monotonically.The lift coefficient increases first and then decreases.The lift coefficient reaches the maximum at TWC=50×10-6.When TWC=0.1×10-6—10×10-6,the lift and drag coefficients are not much different from the clean airfoil because the rime ice on the leading edge remains streamlined.On the one hand,for TWC≥50×10-6,due to the formation of hemispherical glaze ice at the leading edge,the bluff and irregularly shaped body enlarges the area of flow separation.Plus,the roughness of ice increases with TWC.As a result,there is a steep rise in the drag coefficient due to the two factors mentioned earlier.The lift coefficient at TWC≥50×10-6is significantly higher than other working conditions,for there is a large low-pressure area above the upper surface of the airfoil.The airfoil can obtain greater lift force under the effect of the pressure difference between the upper and lower surfaces of the airfoil.

    Fig.10 Lift and drag coefficients at different TWCs

    According to Fig.11 and Fig.12,there is an area where the gauge pressure is less than 0.The static pressure here is lower than the far-field pressure,while the static temperature is higher.This indicates that the eddies convert the kinetic energy and the pressure energy into the internal energy,reducing both the local velocity and static pressure.

    Fig.11 Pressure coefficient contour at different TWCs

    Fig.12 Static temperature contour (TWC=200×10-6)

    2.4 Effects of melt ratio LWC/TWC

    Different from the mixed-phase icing at relatively high temperatures,where an icing plateau exists,with the increase of melt ratio,the mass of ice increases linearly and monotonically,as shown in Fig.13.As the surface temperature of airfoil is extremely low,about 149 K,the droplets quickly freeze after impingement,and the ice crystals are sparsely dispersed in air,resulting in no ice crystal adhesion on the airfoil.The ice crystals almost contribute nothing to ice accretion.All the ice on the airfoil is generated by the solidification of subcooled droplets,so the icing severity is monotonically correlated with LWC.

    Fig.13 Ice thickness distribution and total mass of ice at different LWC/TWC

    In Fig.14,only in cases of LWC/TWC≤0.25,the lift coefficient of the icing airfoil is close to that of the clean airfoil,and in other cases,the lift and drag coefficients show significant deviations.

    Fig.14 Lift and drag coefficients at different LWC/TWC

    2.5 Effects of airflow total temperature T0

    For TWC of 1×10-6and airflow pressure of 200 kPa,the dew point is approximately 204 K.In order to explore the influence of temperature near the dew point on icing,four test conditions between 200—206 K are additionally investigated.

    As the ambient temperature decreases,the runback phenomenon attenuates.Within the temperature range of cryogenic wind tunnels,there is scarcely any runback water after droplets impact the airfoil.In addition,the LWC is too low to form a liquid film.Therefore,the type of ice is rime ice,as shown in Fig.15(a).As TWC is kept unchanged and the total temperature of airflowT0rises,the content of saturated water vapor increases accordingly.The content of droplets and ice crystals decreases.Nevertheless,water vapor does not sublimate under all working conditions in this experiment,so the amount of accumulated ice is reduced,as shown in Fig.15(b).Near the dew point,LWC decreases sharply with the increasing temperature,and the total mass of ice shares the same trend.AtT0=260 K,all water exists in the form of vapor,so the total mass of ice is 0.

    Fig.15 Ice thickness distribution and total mass of ice at different airflow total temperatures

    As depicted in Fig.16,owing to the small thickness and streamlined shape of the accretion layer,the aerodynamic data of the contaminated airfoil is similar to that of the clean airfoil on the whole.In cases ofT0=200 K andT0=204 K,the lift and drag coefficients of the airfoil are significantly reduced.

    Fig.16 Lift and drag coefficients at different airflow total temperatures

    2.6 Effects of airflow total pressure p0

    As can be seen from Fig.17,a rime type ice accretion layer occurs on the leading edge.TWC increases with ascending ofp0,thus the thickness and total mass of ice are positively correlated withp0.The lift and drag coefficients of the contaminated airfoils are lower than those of the clean airfoil.Asp0increases,the deviations also increase gradually,as shown in Fig.18.

    Fig.17 Ice thickness distribution and total mass of ice at different airflow total pressures

    Fig.18 Lift and drag coefficients at different airflow total pressures

    2.7 Effects of Mach number

    With the increase of Mach number,the mass and thickness of ice increase,and the increment attenuates gradually,as shown in Fig.19.

    Fig.19 Ice thickness distribution and total mass of ice at different Mach numbers

    The increase of Mach number results in two opposite effects on icing.In higher speed range,more liquid water will be accumulated on the airfoil in the same time span.This leads to more ice accretion.At the same time,the kinetic energy of ice crystals will increase,deepening the extent of erosion.The kinetic energy dissipation caused by impingement will also melt part of the accumulated ice.Under the shear stress of airflow,the droplets are more likely to flow towards the rear part of the airfoil,resulting in the expansion of the wetting area inYdirection.According to Fig.20,the differences of lift and drag coefficients between contaminated and clean airfoils decrease from positive to negative values with the growing Mach number.In cases ofMa=0.16 and 0.8,the two coefficients deviate evidently from those of the clean airfoil.

    Fig.20 Lift and drag coefficients at different Mach numbers

    2.8 Effects of airfoil surface roughness

    It is recognized that the surface roughness has a significant effect on the frost thickness caused by water vapor desublimation[19].Within the present range of the airfoil surface roughness,it can hardly affect the amount of ice on the airfoil,as shown in Fig.21.Comparing the results,it can be speculated that the water vapor in the air really does not desublimate under such airflow condition.As illustrated in Fig.22,under different working conditions,the surface roughness distributions of ice accretion layer are generally the same and higher than that of the clean airfoil.Hence,the differences of the aerodynamic parameters between icing and clean airfoils are basically unchanged.

    Fig.22 Lift and drag coefficients at different surface roughnesses

    2.9 Effects of icing duration

    With time elapses,there will be more ice accumulation on the surface,as shown in Fig.23.The thickness and mass of ice at the leading edge generally increase linearly with the icing duration extended from 1 h to 5 h.As displayed in Fig.24,the lift and drag coefficients increase accordingly as well.

    Fig.23 Ice thickness distribution and total mass of ice at different icing durations (1×10-6)

    Fig.24 Lift and drag coefficients at different icing durations

    3 Conclusions

    The trace water icing process on the airfoil in cryogenic wind tunnels is simulated based on Messinger’s thermodynamic model.Based on this,the effects of water content parameters,airflow parameters,airfoil surface characteristics,icing duration on ice accretion are investigated.Main conclusions can be drawn as follows:

    (1)Based on the Eulerian method,a multiphase flow model is established to predict the ice formation of trace water on the airfoil under cryogenic conditions.Comparing with the experimental results in the literature,the validated numerical model is suitable to calculate the ice accretion layer.

    (2)TWC and LWC/TWC are the predominant factors affecting the amount of ice on the airfoil.With the increase of TWC,the amount of ice increases linearly.Under high TWC conditions,ice accretion layer at the leading edge will seriously affect the aerodynamic performance of the airfoil.The amount of ice increases monotonically with LWC/TWC.Under corresponding conditions,ice crystals contribute little to the ice accretion process.

    (3)When the airflow total temperature increases,the mass of ice decreases,and the mass changes sharply around the dew point.With the increase of airflow total pressure,the amount of ice generally increases linearly.With the ascending of Mach number,the amount of ice increases,but the increment gets smaller.

    (4)The surface roughness of the airfoil has little effect on the total amount of ice.The amount of ice increases linearly with icing duration.

    日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 电影成人av| 国产男女超爽视频在线观看| 性色av一级| 少妇的丰满在线观看| 男女啪啪激烈高潮av片| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 看免费成人av毛片| 国产精品久久久久久精品古装| 日日摸夜夜添夜夜爱| 日日撸夜夜添| 永久免费av网站大全| kizo精华| 飞空精品影院首页| 一本久久精品| 777米奇影视久久| 免费人妻精品一区二区三区视频| 国产欧美亚洲国产| 下体分泌物呈黄色| 欧美日韩亚洲国产一区二区在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 成人毛片60女人毛片免费| 一个人免费看片子| 午夜激情av网站| av天堂久久9| 少妇的丰满在线观看| 国产成人精品久久久久久| 少妇的逼水好多| 欧美亚洲日本最大视频资源| 免费在线观看视频国产中文字幕亚洲 | 色94色欧美一区二区| 日韩制服丝袜自拍偷拍| 侵犯人妻中文字幕一二三四区| 亚洲国产色片| 久久久国产一区二区| 最近最新中文字幕大全免费视频 | 日本-黄色视频高清免费观看| a级片在线免费高清观看视频| 亚洲精品中文字幕在线视频| 国产免费视频播放在线视频| 啦啦啦在线观看免费高清www| 久久精品亚洲av国产电影网| 亚洲成人av在线免费| 亚洲精华国产精华液的使用体验| 韩国精品一区二区三区| 曰老女人黄片| 最近中文字幕高清免费大全6| 亚洲av福利一区| 国产 一区精品| 国产免费福利视频在线观看| 捣出白浆h1v1| 亚洲三级黄色毛片| 三上悠亚av全集在线观看| 国产精品久久久久久精品电影小说| 黄色怎么调成土黄色| 日本av免费视频播放| 亚洲视频免费观看视频| 午夜精品国产一区二区电影| 亚洲av综合色区一区| 免费观看性生交大片5| 日韩一卡2卡3卡4卡2021年| 视频在线观看一区二区三区| 国产高清不卡午夜福利| 亚洲欧美成人精品一区二区| 一级爰片在线观看| 久热久热在线精品观看| 国产av国产精品国产| 又粗又硬又长又爽又黄的视频| 欧美bdsm另类| 亚洲成色77777| 蜜桃在线观看..| 免费黄频网站在线观看国产| 麻豆乱淫一区二区| 亚洲av免费高清在线观看| 天堂8中文在线网| av天堂久久9| 美女国产视频在线观看| 97在线视频观看| 欧美国产精品va在线观看不卡| 欧美日韩亚洲高清精品| 蜜桃国产av成人99| 国产精品熟女久久久久浪| 一二三四中文在线观看免费高清| 免费av中文字幕在线| 中国三级夫妇交换| 91aial.com中文字幕在线观看| 精品福利永久在线观看| 黑丝袜美女国产一区| 成年人免费黄色播放视频| 国产av一区二区精品久久| 亚洲国产毛片av蜜桃av| 久久久久久人妻| 韩国高清视频一区二区三区| 亚洲精品日本国产第一区| 国产成人a∨麻豆精品| 亚洲国产精品一区三区| 一级毛片电影观看| 街头女战士在线观看网站| 女人久久www免费人成看片| 亚洲精品美女久久久久99蜜臀 | www.熟女人妻精品国产| 国精品久久久久久国模美| 成人国语在线视频| 亚洲av成人精品一二三区| 久久99精品国语久久久| 午夜福利在线观看免费完整高清在| 九草在线视频观看| 日韩一区二区视频免费看| 免费黄网站久久成人精品| 婷婷色av中文字幕| 最近的中文字幕免费完整| 超碰成人久久| 久久午夜福利片| 天堂8中文在线网| 又粗又硬又长又爽又黄的视频| 性少妇av在线| 美女国产视频在线观看| 国产女主播在线喷水免费视频网站| 九草在线视频观看| 午夜免费男女啪啪视频观看| 欧美在线黄色| 波多野结衣av一区二区av| 亚洲av在线观看美女高潮| 久久精品国产亚洲av天美| 女人高潮潮喷娇喘18禁视频| 在现免费观看毛片| 26uuu在线亚洲综合色| 亚洲中文av在线| 亚洲欧美成人精品一区二区| 久久精品久久久久久久性| 少妇人妻精品综合一区二区| 日本欧美国产在线视频| 国产免费视频播放在线视频| 91精品国产国语对白视频| 一本大道久久a久久精品| 在线观看美女被高潮喷水网站| 久久久精品94久久精品| 黑人猛操日本美女一级片| 一边亲一边摸免费视频| 亚洲综合色网址| 中文精品一卡2卡3卡4更新| 亚洲成国产人片在线观看| 亚洲av男天堂| 美女国产视频在线观看| 亚洲国产av影院在线观看| 只有这里有精品99| 男的添女的下面高潮视频| 香蕉丝袜av| 伊人亚洲综合成人网| 老司机影院毛片| 国产一区二区三区综合在线观看| 熟女av电影| 欧美人与善性xxx| 亚洲久久久国产精品| 亚洲成色77777| 午夜福利在线免费观看网站| 久久久久久久久久人人人人人人| 久久久国产精品麻豆| 欧美激情高清一区二区三区 | 亚洲图色成人| 伊人久久国产一区二区| 欧美日韩av久久| a级毛片黄视频| 成年人午夜在线观看视频| 色婷婷久久久亚洲欧美| 亚洲国产色片| 国产精品久久久久久精品电影小说| 97人妻天天添夜夜摸| 九色亚洲精品在线播放| 一边摸一边做爽爽视频免费| 建设人人有责人人尽责人人享有的| 免费看av在线观看网站| 国产精品蜜桃在线观看| 香蕉精品网在线| 日韩制服丝袜自拍偷拍| 91成人精品电影| 亚洲第一av免费看| 青春草视频在线免费观看| 久久久精品免费免费高清| 国产精品不卡视频一区二区| 99国产精品免费福利视频| 女人久久www免费人成看片| 久久国产精品大桥未久av| 久久人妻熟女aⅴ| 男的添女的下面高潮视频| 精品酒店卫生间| 成人亚洲欧美一区二区av| 国产精品二区激情视频| 18+在线观看网站| 久热久热在线精品观看| 一区二区三区乱码不卡18| 色婷婷久久久亚洲欧美| 男男h啪啪无遮挡| 国产精品嫩草影院av在线观看| 亚洲三区欧美一区| 国产片特级美女逼逼视频| 午夜激情av网站| 激情视频va一区二区三区| 美女视频免费永久观看网站| 99久久人妻综合| 午夜影院在线不卡| 午夜福利乱码中文字幕| 在线观看免费视频网站a站| 国产精品不卡视频一区二区| 91精品三级在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久国产精品大桥未久av| 18禁国产床啪视频网站| 大香蕉久久网| 国产女主播在线喷水免费视频网站| 九草在线视频观看| 免费在线观看完整版高清| 亚洲国产日韩一区二区| 看十八女毛片水多多多| 国产女主播在线喷水免费视频网站| 国产精品一国产av| 色婷婷av一区二区三区视频| 免费在线观看黄色视频的| 欧美另类一区| 一区二区三区乱码不卡18| 男男h啪啪无遮挡| 亚洲综合精品二区| www日本在线高清视频| 色婷婷久久久亚洲欧美| 最近2019中文字幕mv第一页| 久热这里只有精品99| 少妇精品久久久久久久| 国产日韩欧美亚洲二区| 日本欧美视频一区| 国产成人一区二区在线| 久久久久久久久久人人人人人人| 成人毛片a级毛片在线播放| 久久精品亚洲av国产电影网| 国产一区二区在线观看av| 黑人猛操日本美女一级片| 亚洲精品日本国产第一区| 99久国产av精品国产电影| 九色亚洲精品在线播放| 18禁动态无遮挡网站| 美女午夜性视频免费| 欧美+日韩+精品| 国产高清不卡午夜福利| 1024视频免费在线观看| 亚洲欧美成人综合另类久久久| 亚洲av福利一区| 午夜福利影视在线免费观看| 一级黄片播放器| av在线app专区| 国产成人午夜福利电影在线观看| 欧美日韩综合久久久久久| 中文字幕人妻熟女乱码| 爱豆传媒免费全集在线观看| 激情五月婷婷亚洲| 久久午夜福利片| 激情视频va一区二区三区| 一个人免费看片子| 精品少妇黑人巨大在线播放| 欧美精品一区二区大全| 免费高清在线观看视频在线观看| 成人国产麻豆网| 国产一区二区在线观看av| av又黄又爽大尺度在线免费看| 日韩成人av中文字幕在线观看| 亚洲精品国产av成人精品| 韩国av在线不卡| 国产在线一区二区三区精| 少妇被粗大的猛进出69影院| 久久国产精品大桥未久av| 国产极品天堂在线| 一级毛片 在线播放| 制服人妻中文乱码| 国语对白做爰xxxⅹ性视频网站| 国产精品国产三级国产专区5o| 91成人精品电影| 99re6热这里在线精品视频| 午夜福利在线观看免费完整高清在| 亚洲美女黄色视频免费看| 亚洲精品乱久久久久久| 久久 成人 亚洲| 美女国产高潮福利片在线看| 亚洲精品国产av蜜桃| 国产精品香港三级国产av潘金莲 | 日韩三级伦理在线观看| 久久久久久久久免费视频了| 免费在线观看视频国产中文字幕亚洲 | 黄片无遮挡物在线观看| 亚洲av成人精品一二三区| 一级毛片电影观看| av电影中文网址| 999精品在线视频| 两个人免费观看高清视频| 老司机影院毛片| 国产精品一国产av| √禁漫天堂资源中文www| 99九九在线精品视频| 999精品在线视频| 午夜免费观看性视频| 国产av码专区亚洲av| 国产熟女午夜一区二区三区| 一边亲一边摸免费视频| 天天躁日日躁夜夜躁夜夜| 日韩精品免费视频一区二区三区| videosex国产| 麻豆av在线久日| 18禁裸乳无遮挡动漫免费视频| av在线播放精品| 麻豆精品久久久久久蜜桃| 青春草视频在线免费观看| 国产男女超爽视频在线观看| 成人黄色视频免费在线看| 国产成人精品婷婷| 亚洲国产最新在线播放| 男人舔女人的私密视频| 免费在线观看完整版高清| 欧美日韩一区二区视频在线观看视频在线| 看非洲黑人一级黄片| 91久久精品国产一区二区三区| av女优亚洲男人天堂| 国产欧美日韩一区二区三区在线| 精品视频人人做人人爽| 你懂的网址亚洲精品在线观看| 日韩不卡一区二区三区视频在线| 久久久久久久精品精品| 一区二区三区激情视频| 1024香蕉在线观看| 婷婷色麻豆天堂久久| 亚洲欧洲日产国产| 亚洲欧美清纯卡通| 麻豆乱淫一区二区| 成年女人在线观看亚洲视频| 寂寞人妻少妇视频99o| 中文精品一卡2卡3卡4更新| 热re99久久国产66热| 丝袜脚勾引网站| 免费少妇av软件| 欧美人与善性xxx| 久久久久久久久久久久大奶| 国产精品三级大全| 少妇 在线观看| 波多野结衣av一区二区av| 久久av网站| 一级,二级,三级黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜美足系列| 纯流量卡能插随身wifi吗| 亚洲男人天堂网一区| 免费人妻精品一区二区三区视频| www.熟女人妻精品国产| 777久久人妻少妇嫩草av网站| 日韩视频在线欧美| 在线看a的网站| 又粗又硬又长又爽又黄的视频| av国产久精品久网站免费入址| 欧美日韩亚洲高清精品| 欧美日韩综合久久久久久| 又粗又硬又长又爽又黄的视频| 九色亚洲精品在线播放| 精品国产一区二区三区四区第35| 肉色欧美久久久久久久蜜桃| 99热网站在线观看| 99久久综合免费| 国产乱人偷精品视频| 精品国产乱码久久久久久男人| 少妇的丰满在线观看| 国产片特级美女逼逼视频| 午夜福利乱码中文字幕| 亚洲av中文av极速乱| 国产精品一区二区在线不卡| 精品国产超薄肉色丝袜足j| 在线精品无人区一区二区三| 欧美精品一区二区免费开放| 日日摸夜夜添夜夜爱| 亚洲三区欧美一区| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 国产日韩一区二区三区精品不卡| 亚洲三区欧美一区| 人妻一区二区av| 91aial.com中文字幕在线观看| 在线观看美女被高潮喷水网站| 亚洲成人一二三区av| 看免费成人av毛片| 又大又黄又爽视频免费| 91精品国产国语对白视频| 久久人人97超碰香蕉20202| 亚洲成色77777| 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 亚洲三级黄色毛片| 母亲3免费完整高清在线观看 | 精品第一国产精品| 国产无遮挡羞羞视频在线观看| 亚洲精品美女久久av网站| 美女国产视频在线观看| 国产老妇伦熟女老妇高清| 亚洲男人天堂网一区| 国产精品成人在线| 日本猛色少妇xxxxx猛交久久| 日韩 亚洲 欧美在线| 国产又爽黄色视频| 高清黄色对白视频在线免费看| 亚洲欧洲日产国产| 麻豆精品久久久久久蜜桃| 亚洲精品美女久久久久99蜜臀 | 在线免费观看不下载黄p国产| 日韩在线高清观看一区二区三区| 大香蕉久久成人网| 日本av手机在线免费观看| 婷婷成人精品国产| 美女高潮到喷水免费观看| 国产成人免费观看mmmm| 国产成人aa在线观看| 美女午夜性视频免费| 18在线观看网站| 久久97久久精品| av线在线观看网站| 亚洲精品aⅴ在线观看| av天堂久久9| 久久人人97超碰香蕉20202| 国产免费福利视频在线观看| 久久久久人妻精品一区果冻| 久久久久久人妻| 寂寞人妻少妇视频99o| 中国国产av一级| 亚洲av综合色区一区| 五月开心婷婷网| 一本色道久久久久久精品综合| 亚洲五月色婷婷综合| 大香蕉久久网| 99久久中文字幕三级久久日本| 曰老女人黄片| 亚洲国产毛片av蜜桃av| 成人毛片a级毛片在线播放| 亚洲国产精品一区三区| 乱人伦中国视频| 一级毛片黄色毛片免费观看视频| 国产精品人妻久久久影院| 秋霞在线观看毛片| 久久久久久久久久人人人人人人| 一边亲一边摸免费视频| 欧美 日韩 精品 国产| 亚洲av电影在线进入| 久久久国产一区二区| 麻豆精品久久久久久蜜桃| 啦啦啦啦在线视频资源| 性高湖久久久久久久久免费观看| 国产一区亚洲一区在线观看| 人妻系列 视频| 国产白丝娇喘喷水9色精品| 精品久久蜜臀av无| 一级黄片播放器| 成人亚洲精品一区在线观看| 亚洲精品国产一区二区精华液| 蜜桃国产av成人99| 久久久久久久久久人人人人人人| 日本猛色少妇xxxxx猛交久久| 亚洲成人一二三区av| 香蕉丝袜av| 久久影院123| 国产亚洲精品第一综合不卡| 91aial.com中文字幕在线观看| av国产久精品久网站免费入址| 久久精品aⅴ一区二区三区四区 | 2018国产大陆天天弄谢| 狠狠精品人妻久久久久久综合| 国产午夜精品一二区理论片| 最黄视频免费看| 男人爽女人下面视频在线观看| www.av在线官网国产| xxx大片免费视频| 蜜桃国产av成人99| 午夜福利在线观看免费完整高清在| 精品第一国产精品| 在线观看人妻少妇| 国产精品偷伦视频观看了| 亚洲欧美清纯卡通| 捣出白浆h1v1| 欧美中文综合在线视频| 韩国精品一区二区三区| 最新的欧美精品一区二区| 一本久久精品| 捣出白浆h1v1| 2018国产大陆天天弄谢| 亚洲美女搞黄在线观看| 亚洲精品日韩在线中文字幕| 欧美人与性动交α欧美软件| 久久久久国产网址| 熟女少妇亚洲综合色aaa.| 久久久久久久国产电影| 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| 蜜桃在线观看..| 国产精品亚洲av一区麻豆 | 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 尾随美女入室| 成人国产麻豆网| 一本久久精品| 天天影视国产精品| 边亲边吃奶的免费视频| 欧美少妇被猛烈插入视频| 国产综合精华液| 老熟女久久久| 亚洲精品视频女| 日日摸夜夜添夜夜爱| av视频免费观看在线观看| 国产精品一国产av| 亚洲欧美清纯卡通| 成年人免费黄色播放视频| 一边亲一边摸免费视频| 伊人久久大香线蕉亚洲五| 日韩一本色道免费dvd| 男男h啪啪无遮挡| 国产成人一区二区在线| 国产精品人妻久久久影院| 熟妇人妻不卡中文字幕| 男女啪啪激烈高潮av片| 在线观看免费高清a一片| 午夜福利一区二区在线看| 色婷婷av一区二区三区视频| 黄色 视频免费看| 永久免费av网站大全| 亚洲国产欧美网| 国产精品不卡视频一区二区| 一本色道久久久久久精品综合| 国产精品久久久久久精品古装| 成人毛片a级毛片在线播放| 日日爽夜夜爽网站| 18+在线观看网站| 可以免费在线观看a视频的电影网站 | 国产白丝娇喘喷水9色精品| 中文字幕色久视频| 亚洲美女搞黄在线观看| 两性夫妻黄色片| 水蜜桃什么品种好| 99国产综合亚洲精品| 精品第一国产精品| 老女人水多毛片| 国产精品秋霞免费鲁丝片| 久久久久久免费高清国产稀缺| 母亲3免费完整高清在线观看 | 午夜精品国产一区二区电影| 亚洲av国产av综合av卡| 色婷婷久久久亚洲欧美| 久久女婷五月综合色啪小说| 高清不卡的av网站| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| av在线老鸭窝| 免费播放大片免费观看视频在线观看| 一边摸一边做爽爽视频免费| 午夜福利在线观看免费完整高清在| 卡戴珊不雅视频在线播放| 亚洲欧美成人精品一区二区| 亚洲内射少妇av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲美女视频黄频| 999精品在线视频| 制服诱惑二区| 久久久久久久国产电影| 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 午夜福利在线免费观看网站| av国产精品久久久久影院| 久久亚洲国产成人精品v| 国产精品三级大全| 久久婷婷青草| 国产亚洲欧美精品永久| 精品亚洲成国产av| 人人妻人人澡人人看| 亚洲精品中文字幕在线视频| 成人毛片a级毛片在线播放| 成年人午夜在线观看视频| 欧美激情高清一区二区三区 | 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 国产一区亚洲一区在线观看| 国产野战对白在线观看| 日韩一卡2卡3卡4卡2021年| 女的被弄到高潮叫床怎么办| 如日韩欧美国产精品一区二区三区| 少妇熟女欧美另类| 亚洲精品,欧美精品| 韩国av在线不卡| 久久99精品国语久久久| 人妻系列 视频| 亚洲精品,欧美精品| 丰满饥渴人妻一区二区三| 亚洲人成电影观看| 少妇人妻精品综合一区二区| 亚洲精品自拍成人| 中文字幕精品免费在线观看视频| 精品国产露脸久久av麻豆| 性色avwww在线观看| 99国产综合亚洲精品| av国产久精品久网站免费入址| 大香蕉久久网| av卡一久久| 99精国产麻豆久久婷婷| 久久国内精品自在自线图片| 肉色欧美久久久久久久蜜桃| 欧美国产精品va在线观看不卡| 制服丝袜香蕉在线| 国产精品嫩草影院av在线观看| 久久人妻熟女aⅴ| 国产精品嫩草影院av在线观看| 纵有疾风起免费观看全集完整版| 亚洲精品日本国产第一区| 亚洲一码二码三码区别大吗| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人 |