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

    Modeling and simulation of solvent behavior and temperature distribution within long stick propellants with large web thickness undergoing drying

    2023-11-11 04:06:58EnfaFuQianlingLiuYuLuanYaoZhuWeidongHeZhenggangXiao
    Defence Technology 2023年10期

    Enfa Fu, Qianling Liu, Yu Luan, Yao Zhu, Weidong He, Zhenggang Xiao

    Key Laboratory of Special Energy Materials, Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, People's Republic of China

    Keywords:Stick propellants Drying Large web thickness 3D numerical modeling Heat transfer Solvent behavior

    ABSTRACT Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin layer near the surface to block the free access of most solvent through for long stick propellants with large web thickness,which lead to lower drying efficiency and worse drying quality.This study aims to gain a comprehensive understanding of drying process and clarify the mechanism of the blocked layer near the propellant surface.A new three-dimensional coupled heat and mass transfer(3D-CHMT)model was successfully developed under transient conditions.The drying experiment results show that the 3DCHMT model could be applied to describe the drying process well since the relative error of the content of solvent between simulation and experiment values is only 5.5%.The solvent behavior simulation demonstrates that the mass transfer process can be divided into super-fast (SF) and subsequent minorfast(MF)stages,and the SF stage is vital to the prevention of the blocked layer against the free access for solvent molecules inside propellant grains.The effective solvent diffusion coefficient (Deff) of the propellant surface initially increases from 3.4 × 10-6 to 5.3 × 10-6 m2/s as the temperature increases, and then decreases to 4.1 × 10-8 m2/s at 60-100 min.The value of Deff of surface between 0-1.4 mm has a unique trend of change compared with other regions, and it is much lower than that of the internal at 100 min under simulation conditions.Meanwhile, the temperature of the propellant surface increases rapidly at the SF stage (0-100 min) and then very slowly thereafter.Both the evolution of Deff and temperature distribution demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and its thickness is about 1.4 mm.To mitigate the formation of blocked layer and improve its drying quality of finial propellant products effectively, it should be initially dried at lower drying temperature(30-40 °C)in 0-100 min and then dried at higher drying temperature (50-60 °C) to reduce drying time for later drying process in double base gun propellants.The present results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.

    1.Introduction

    Long stick propellants are widely used in artillery due to many advantages, such as simplified charge configuration and improved interior ballistic performance [1-3].The regular packing of long stick propellants may have higher loading densities than that of granular propellants since their packing is random.For long stick propellant with natural flow channels,the proper venting by partial cutting perpendicular to the axis of the stick can provide a rapid escape of burning gases and effectively reduce the pressure-wave during the burning process, thus improve the safety of propelling charge in service [1,3,4].

    Nomenclature Ci Concentration/(mol?m3)Cp Specific heat/(J?kg-1?K-1)Di Solvent diffusion coefficient/(m2?s-1)R Drying rate/(% per h)hm Mass transfer coefficient/(m?s-1)ht Heat transfer coefficient/(W?m-2?K-1)k Thermal conductivity/(W?m-1?K-1)K Evaporation rate constant/(s-1)mi Weight/g Me Equilibrium solvent content/%M0 Initial solvent content/%T Temperature/°C t Time/s v Velocity/(m?s-1)Volume/m3 dt Time interval/s x, y, z Coordinates SF Super-fast MF Minor-fast 3D-CHMT Three-dimensional coupled heat and mass transfer mt The instantaneous mass at drying moment t of propellant sample Ct The specific heat of stick propellant grains at drying moment t▽▽ Divergence Ni The internal solvent diffusion rate Ri The internal solvent evaporation rate μ The velocity vector P The mean relatively percentage deviation modulus V

    Drying is an indispensable production process of long stick propellants [5-7].It not only requires a mass of energy to remove excess solvent, but also causes changes in the physicochemical properties of the final product [8,9].During the drying process of propellant grains, inappropriate drying techniques may lead to poorer performance of propellants and deformation of propellant matrix shape[10,11].The drying should be conducted at a suitable temperature and time range.Urbanski [12] proposed that higher drying temperature at the beginning may cause the propellant grains to swell or crack, but no specific solution was given.And longer drying time may cause the unnecessary consumption of energy and lower drying efficiency.Especially, for long stick propellants with large web thickness[13],it is easy to form a blocked layer near the surface of propellant grains which can immensely prevent excess solvent inside the interior grain from transferring to the outside.The above problems can easily arose from improper determination of drying parameters most of which are selected empirically for long stick propellants with large web thickness during present drying techniques.

    So far,some experiment methods have been utilized to optimize important drying parameters of propellants [14-16], such as solvent content, drying temperature, air humidity and so on.Chen et al.[17]demonstrated that the content of residual solvents has a significant influence on the mechanical property of TEGDN propellants by using the auto-mass device.Cheng et al.[18] investigated the effect of the drying temperature on the safety performance of modified double-based propellants.The results show that the fast drying of propellant grain can be realized at 80°C by convective heat transfer.The above experimental methods need large amounts of experiments with high costs and longer time.Since the limitation of existing measurement techniques and hostile experimental conditions, mass and heat transfer of the drying process of propellant matrix are difficult to be directly observed in experiments generally, especially the measurement of solvent concentration variations and gradient [19].Meanwhile, few research has been conducted to investigate the form mechanism and prevention method of blocked layer near the propellant theoretically.Consequently,the accurate drying parameters of long stick propellants with large web thickness could not be determined easily.

    Thus,it is imperative to develop a suitable and scientific method to clarify the mechanism of block layer near the propellant surface and determine appropriate drying process parameters of long stick propellants with large web thickness theoretically.Although the numerical method has been well used in the determination of drying parameters and quality control of agriculture products such as vegetables and fruit [9,20], few studies on modeling and simulation of the drying process of long stick propellants with large web thickness have been carried out yet.It should be noted that stick propellants have some unique thermophysical properties different from above drying materials [21].So it is necessary to build a new mathematical model for long stick propellants with large web thickness.

    In this study, a new three-dimensional coupled heat and mass transfer (3D-CHMT) model with changing solvent diffusivity was developed for long stick propellants with large web thickness during the drying process.The drying experiment was conducted to validate the established model of long stick propellants with large web thickness.And the solvent behavior and temperature distribution within stick propellant grains were both investigated.The formation mechanism of the blocked layer near the propellant surface was further analyzed, and a facile method to inhibit the formation of blocked layer during drying process was established for long stick propellants with large web thickness.

    2.Materials and experimental investigations

    2.1.Preparation of long stick propellants with large web thickness

    The nitrocellulose-nitroglycerin paste was purchased from Sichuan Nitrocellulose Co.Ltd.The propellant samples were prepared by the semi-solvent method[22].First,200 g dried paste was added in a kneader.The mixed solvent of acetone and ethanol(about 70 mL)with the mass ratio of 1:1 was sprayed on the paste.Subsequently,a dough was obtained after gelatinizing the mixture for 2.5 h,and it was added into a mold with a 19-perforation die.A propellant strand with a web thickness of 2.3 mm and a perforation diameter of 0.5 mm was extruded from the mold.Finally, the propellant sample was dried at room temperature for 24 h.The propellant sample was usually cut with a feasible length,i.e less than 80 mm, for the closed bomb test and other experiments.In addition, the longer the length of the propellant sample is, the more computing resources and calculation time the simulation process needs.Therefore,in this study,the long stick propellant with large web thickness was cut with a length of 60 mm for further investigation.

    2.2.Drying experiment

    The drying experiment was conducted to determine the change of the content of solvent inside the stick propellant grains.Four drying temperatures (30°C, 40°C, 50°C and 60°C) were selected during the drying experiment.

    The instantaneous mass of the propellant sample mtwas measured by on-line auto mass measurement experiment setup developed by our group [17], and the experiment setup is illustrated in Fig.1.The mass of propellant sample used in the drying experiment is 10 g.

    The instantaneous content of solvent Ctinside the stick propellant grains is calculated by Eq.(1)

    where m0and mtis the initial mass and the instantaneous mass at drying moment t of propellant sample, Ctand C0are the instantaneous content of solvent at drying moment t and initial content of solvent of the propellant sample, respectively.

    3.Numerical model

    A new 3D-CHMT model was established to describe the drying process of long stick propellants with large web thickness.

    3.1.Assumptions

    Some assumptions used in the modeling of long stick propellants with large web thickness during the drying process are as follows:

    (1) Stick propellants are homogenous and isotropic materials.The shrinkage effect of propellant grains is neglected as the solvent content decreases.

    (2) The solvent within the stick propellant grains is transported by liquid flow-through,and the evaporation only takes place at the surface of the propellant grains.

    (3) The drying air in the interior perforation of the propellant grains is initially uniform before drying.

    (4) There is no heat generation within the stick propellant grains.Operating conditions of surface and interior perforation of the stick propellant grains are assumed alike.

    (5) The effect of radiation is neglected during drying process.The heat transfer between the drying air and propellant grains are the coupling of conduction, convection, and radiation.However, the radiation effect is extremely small and neglected in this study.

    Fig.1.Schematic diagram of automatic-acquisition drying experiment setup.

    3.2.Heat and mass transfer equations

    During the drying process,the drying air blows over propellant grains, which makes temperature of propellant surface increases due to the convective heat transfer.Meanwhile, the solvent inside stick propellant grains evaporates on account of increase of ambient temperature.Solvent molecule moves to propellant surface on the effect of convection and diffusion,then is carried away by drying air, which results in the gradual reduce of solvent content,as shown in Fig.2.The governing equation of 3D-CHMT model is established, and it includes heat and mass equations between drying air and stick propellant grains.

    During the drying process, the solvent molecule diffusion conforms to the second law of Fick diffusion[14],which is an unsteady process changing with time.The equations for mass transfer are shown in Eq.(2) and Eq.(3)

    where Ciis the instantaneous concentration of solvent diffusion,Riis the internal solvent evaporation rate,which represented the rate of solvent evaporation inside the propellant grains,Niis the internal solvent diffusion rate, Diis the solvent inter-diffusion coefficient,which is the physical quantity of solvent molecular diffusion ability that can accurately reflect the diffusion behavior of solvent molecule in the propellant matrix [23].

    The heat transfer in stick propellant grains is evaluated by Eq.(4)using Fourier's Law, and similar heat equations are provided by Bhargava et al.[24-26].

    where ρ is density of stick propellant grains, λ is thermal conductivity of stick propellant grains, Cpis the specific heat of stick propellant grains,Q is the amount of heat lost due to the evaporation.In this study, the solvent used in stick propellant grains is the mixture of acetone and ethanol.

    As heat is transferred from drying air to stick propellant grains,solvent molecule becomes excited by the addition of heat which increases the potential and kinetic of solvent molecule.The continuous addition of more heat causes solvent molecule to diffuse outside the propellant grains and evaporate as vapor.The solvent molecule carries heat away as they leave from stick propellant grains, thus a certain amount of heat is lost.The heat transfer for evaporation of solvent Q from propellant grains are evaluated using Eq.(5)

    Fig.2.Schematic of the drying process of stick propellant grains.

    where Wdis the total weight of stick propellant grains, Qstis the latent heat of the vaporization of solvent, Vtis the total volume of stick propellant grains, dCt/dt is the time differential of solvent.

    3.3.Initial and boundary conditions

    The initial and boundary conditions for heat transfer process are defined as: homogeneous temperature at the beginning of the drying process, convection on the surface of the propellant material.Similarly, the initial and boundary conditions for the mass transfer process are defined as: homogeneous solvent content at the beginning of drying process,convective condition on surface of propellant materials.A no-slip boundary condition is applied on stick propellant surface, and same boundary conditions have been used in previous research[27,28].

    Initial conditions

    where T0is the initial temperature, K.C0is the initial content of solvent inside the stick propellant grains, %.

    3.4.Numerical solution

    To analyze the drying process of long stick propellants with large web thickness clearly, the area between the end face A and face B was defined as the calculation region of this model,as shown in Fig.3.The face B is a hypothetical face, which is located in the middle of the propellant grains.The length from the end face A to face B is 30 mm.A divided grid of the 3D-CHMT model of the propellant grains is shown in Fig.4.The mesh,composed of 340000 quadrilateral elements,was created by the software on the physical field control condition.The time step of simulation process is 10000, and the time interval is 1 s.

    Some material parameters are determined experimentally to solve equations in this study, such as heat capacity of stick propellants at constant pressure (Cp), thermal conductivity (k) and effective solvent diffusivity (Deff) of stick propellants, as shown in Table 1.And the density (ρ) and latent heat of vaporization (λ) of stick propellants are instead taken from previous study [17,20].

    Fig.3.3D diagram of long stick propellant with large web thickness.

    Fig.4.Divided grid of 3D model of the stick propellant grains.

    This model was solved according to the flow chart illustrated in Fig.5.Heat Transfer in Solids and Transport of Diluted Species in COMSOL Multiphysics v4.5(COMSOL Inc.,Palo Alto,USA)were used to handle the heat and mass transfer of the drying process.Then,the time-dependent Backward Differentiation Formula solver was used with the Multifrontal Massively Parallel Sparse Direct Solver to perform the simulation with a time step.

    4.Results and discussion

    4.1.Experimental validation of established 3D-CHMT model

    To validate the established 3D-CMHT model, the average content of solvent of the stick propellant grains is estimated from numerical results and the predicated content of solvent is compared with the experimental values.A detailed comparison of the experimental results with predicted solvent content of the propellant sample at various drying temperatures is shown in Fig.6.The solvent loss process of the propellant sample consists of two parts:the curve of Period I(0-48 h)is very steep,and the curve of Period II (48-168 h) is gently inclined.In Period I, about 75%solvent mass of the propellant sample is removed.The solvent gradient within the propellant grains is relatively high at the beginning of drying, and it is decreased as the drying time increases.In Period II,the solvent loss reaches a stable stage,which is approximately below 2%.After 168 h, the content of residual solvent in the stick propellant grains has almost no obvious change.

    It is observed that the rate of solvent loss increases with drying temperature since the solvent diffusivity of the stick propellant grains depends on the temperature [17].The drying time to reach steady content of solvent is decreased as drying temperature T increases because this rise in drying temperature increases the solvent diffusion coefficient,that increases the diffusion rate of solvent from interior to propellant surface.Scilicet,the loss rate of solvent in the stick propellant grains is higher when the drying temperature T is higher, leading to shorter drying time.It shows that a higher temperature of drying air takes less drying time to reach a same residual content of solvent.For example,the time to reach 2%residual content of solvent is decreased by approximately 25%from 120 to 96 h as drying temperature increases from 30 to 40°C.The increase in T from 30 to 60°C reduces the drying time to reach the same residual content of solvent about 150% from 120 h to 48 h.

    It is observed in Fig.6 that the simulated results agree well with the experimental values.In this study,the capacity of the 3D-CHMTmodel to represent the experimental data is realized in the terms of the mean relatively percentage deviation modulus P[29,30]and it is calculated by Eq.(8)

    Table 1Physical parameters of stick propellants used in the simulation.

    Fig.5.Flow chart of the calculation strategy in the simulation.

    where N is the total number of compared samples,and eiand piare the experimental and predicated solvent content of samples at the same condition, respectively.

    The P-value reflects the accuracy and reliability of the mathematics model.Based on Eq.(9), the P between the predicated and experimental results is only 5.5%in this study.The results indicate that the simulation predicts the general trend of the experimental solvent content curve well, which confirms that the established model can be described the drying process of long stick propellants with large web thickness.It also shows that the established mathematics model of long stick propellants with large web thickness is reasonable, and the solvent and temperature distribution could be achieved based on that model.

    Fig.6.Comparison between experimental and theoretical content of solvent at various drying temperatures.

    4.2.Solvent behavior

    The solvent distribution in the long stick propellants with large web thickness during the drying process at drying temperature of 50°C is shown in Fig.7,with red shades indicating higher content of solvent and purple shades indicating the lower content of solvent.The stick propellant grains are with 20.0% initial content of solvent(mixed solvent of acetone and ethanol).At the initial drying stage, the free solvent molecules on the surface of the propellant sample are dried first and the propellant surface reaches lower content of solvent after a few minutes,resulting in solvent gradient inside the stick propellant grains.The solvent gradient is the rate at which solvent changes, or increases and decreases, between one region and another, and it is a major driving force for solvent transport.In turn, the gradient of solvent content inside the stick propellant grains continues to drive the migration of solvent from the interior to the surface.When the drying process proceeds further,the solvent in the interior of the propellant sample diffuses to the surface and is eliminated to surrounding drying air from the surface further.It is mentioned that unbound solvent is initially removed, and bound solvent takes time to transfer, but these mechanisms happen continuously.The content of solvent inside the propellant grains is gently decreased over time.Finally, the propellant surface achieves steady content of solvent after 100 min,and the content gradient of solvent in the surface is almost lower than 0.5%, indicating that the drying reaches to a steady stage.

    Fig.7.Surface of predicated solvent profile of stick propellant grains at drying temperature of 50 °C.

    To better investigate the solvent behavior of stick propellant grains, the slices are chosen at the height of 7.5,15 and 22.5 mm along the axial direction in the propellant grains,as shown in Fig.8.Fig.8 represents the content of solvent in the interior of the stick propellant grains.There is a significant solvent gradient in the 19 perforations of the stick propellant grains at 0-60 min,and they are all similar with that of the surface in the propellant samples.Due to the existence of perforations in the propellant samples,the solvent can be fast transferred by the evaporation on the perforation surface.But the heat contact area of the perforations with drying air is relatively smaller compared with the whole propellant grains.The results indicate that the perforation in the stick propellants has a slight positive effect on the process of solvent transfer.

    Meanwhile, an obvious difference appears at 0-200 min.Specifically,the content of solvent within the propellant grains has an obvious gradient in three different slices at 0-100 min while it has a more uniform distribution after 100 min, showing the early drying process is divided into two periods: super-fast (SF) stage(0-100 min)and subsequent minor-fast(MF)stage(100-200 min).The blocked layer near the propellant surface formed by over-fast drying prevents the residual solvent from transferring to the outside,and the over-fast drying usually takes place in the SF stage[17].Hence, the solvent behavior in the SF stage should be restrained to inhibit the formation of the blocked layer near the propellant surface during the drying process.The aim of this work is used to analyze the solvent behavior inside the stick propellant grains,which can guide the optimization of drying process of long stick propellants with large web thickness.

    Fig.8.Various internal slices of predicated solvent profile of the stick propellant grains at drying temperature of 50 °C.

    Fig.9.Diagram of profile in the stick propellant grains.

    To further investigate the formation mechanism of the blocked layer near the propellant surface(the profile of long stick propellant see Fig.9), the effective solvent (mixed solvent of acetone and ethanol)diffusion coefficient Deffof the propellant grains at drying temperature of 50°C and 30°C was calculated by the simulation results based on Eq.(3),as shown in Fig.10 and Fig.11.The Deffis the rate at which a diffusing substance is transported between opposite faces of a unit cube of a system when there is unit content difference between them, and it is related to the temperature and physical structure of propellant sample[31,32].At the beginning of drying process, the initial value of Deffof the propellant sample is about 3.4×10-6m2/s.During the initial drying process(0-10 min),the value of Deffis rapidly increased over time and its increasing rate of the surface is more than that of the interior.In this period,the solvent is fast transferred from the propellant grains into the drying air.At the following drying stage.At 60-100 min,the value of Deffwithin propellant grains all decreases over time and its decreasing rate in the surface is more than that of the interior.At the ending drying stage (after 100 min), the value of Deffof the surface decreases to only 8.1 × 10-8m2/s but the interior is 2.7 × 10-7m2/s, indicating that the propellant surface is severely dried and has formed a hard skin layer to block the free access of most solvent through.Compared Fig.11 with Fig.10, the similar phenomena do not appear after 100 min and little blocked layer near the propellant surface is observed for propellants dried under 30°C, showing initial lower drying temperature can effective inhabit the formation of blocked layer for long stick propellants with large web thickness.As a result,the solvent is being held at the internal and only few solvent molecules could diffuse out to the surface due to the formation of the blocked layer near the propellant surface.The above results show that the general effective solvent diffusion coefficient of the propellant grains is enormously decreased after the middle drying stage (60-100 min).That is the main reason why the residual solvent within long stick propellants with large web thickness could not be completely removed during the later drying process.And this further indicates that the blocked layer near the propellant surface plays an important role in the mass transfer process and it prevents more residual solvent from transferring to the outside of the propellant grains.

    Fig.10.Effective solvent diffusion coefficient contour in section of stick propellant grains at drying temperature of 50 °C.

    In this study, the solvent diffusion behavior inside the stick propellant grains conforms to the Vrentas-Duda model[33],and it is largely affected by temperature and physical structure of propellant grains.The blocked layer near the propellant surface could limit the movement of polymer chains(the main component of the stick propellant)and result in the reduce of its free volume,leading to the difficult diffusion of solvent molecules in the propellant grains.So the value-Deffof the propellant surface is decreased as the blocked layer is formed.

    Fig.11.Effective solvent diffusion coefficient contour in section of stick propellant grains at drying temperature of 30 °C.

    Fig.12.Effective solvent diffusion coefficient distribution on the centerline of face B at drying temperature of 50 °C and 100 min.

    Fig.12 shows the effective solvent diffusion coefficient (Deff)distribution on the centerline of face B at drying temperature of 50°C and 100 min, and all data points of Deffare chosen at the equidistance(0.5 mm)along the centerline.It is obviously observed in Fig.12 that the value-Deffat the internal has vibrational distribution along the centerline between the perforations.Combined with solvent transfer results, the solvent content at the surface is much lower than that of the internal at 100 min, and the solvent content is first increased and then decreased along the centerline from the perforation to adjoining perforation.The nitrocellulose chains, the main component of long stick propellants, have lots of hydroxyl groups and may be linked by hydrogen bonds between the chains at the lower solvent content,resulting in the significant reduction of free volume and ultimately decreasing the solvent diffusion coefficient of the propellant surface and inner perforation surface.So the value-Deffat the propellant has similar change rule along the centerline between the perforations at 100 min compared with the gradient distribution of solvent content, indicating that the solvent content has a slight positive effect on the value-Deffof the propellant grains at 100 min.As shown in Fig.10,the effective solvent diffusion coefficient of propellant surface between 0-1.4 mm has unique rate of change compared with other region,and it is much less than that of the internal, further indicating the blocked layer near the propellant surface has been formed and its thickness is about 1.4 mm.

    Fig.13.Effective solvent diffusion coefficient distribution on the centerline of face B at drying temperature of 30 °C and 100 min.

    Fig.13 shows the effective solvent diffusion coefficient distribution on the centerline of face B at drying temperature of 30°C and 100 min.As can be seen from Fig.13,the value-Deffat the stick propellant grains has vibrational distribution along the centerline between the perforations.Compared Fig.13 with Fig.12,the unique trend of distribution of value-Deffat the propellant surface between 0-1.4 mm disappears as the drying temperature decreases from 50 to 30°C.The results show that the formation of blocked layer has been effective mitigated at drying temperature of 30°C in 0-100 min.These results could be applied to improve the drying quality of final propellant products and prevent the formation of the blocked layer near the propellant surface.For example,the long stick propellants with large web thickness are initial dried at lower drying temperature (30-40°C) at 0-100 min and then dried at higher drying temperature (50-60°C) to reduce drying time for later drying process.For the case of this experiment,it can effective inhabit the formation of blocked layer and improve drying quality for long stick propellant with large web thickness based on this method.

    Based on this method (using different drying temperature combination)in this study,the stick propellant grains are first dried at 30°C in 0-100 min and then dried at 60°C after 100 min,and the evolution of solvent content within the propellant grains is shown in Fig.14.Compared with that propellant samples dried at 60°C,the drying time to achieve same solvent content(2%)is reduced by 15% and the residual content of solvent for finally propellant product is reduced by 20% under different drying temperature combination condition, showing the drying quality and drying efficiency are obviously improved under this condition for long stick propellants with large web thickness.

    Fig.14.Comparison between content of solvent for propellants dried at 60 °C and different drying temperature combination.

    4.3.Temperature distribution

    The simulation of temperature distribution is conducted to investigate the heat transfer in the stick propellant grains during the drying process.Fig.15 displays the evolution of average temperature of stick propellant grains under simulation condition.As shown in Fig.15,the increase of temperature at 60°C is faster than that at lower drying temperatures due to the larger temperature difference between stick propellant grains and drying air.As drying temperature increases,heat transfer via convection from drying air to propellant surface and subsequently by the means of conduction to the center of stick propellant grains occurs and causes the increase of temperature of the propellant sample.The slope of simulated temperature curves is very steep and the slope of the curve at 60°C is higher than that at other temperatures, showing that there is a rapid thermal energy transfer process.Moreover,the rate of heat transfer is reduced over time,which is reflected by the slope going to zero.

    Fig.15.Simulated data showing the average temperature of stick propellant grains at various drying temperatures.

    The temperature contour inside the propellant grains is estimated for all drying temperatures(30,40,50,and 60°C),and they have similar trend of change.In this paper, the temperature contours at drying temperature of 30°C and 50°C for further investigation are respectively displayed in Fig.16 and Fig.17.

    Fig.16 clearly shows that the surface of the propellant grains gets heated initially in 10 min because of convection heat transfer between the propellant and the drying air.The heat is conducted into the interior region of stick propellant grains as the interior one of stick propellant grains has a low temperature region compared with the surface region.At the early drying stage, the heat consumed by evaporation exceeds that transferred from propellant surface under lower drying temperature condition and evaporation rate in the internal is larger, resulting in a slightly dropping of internal temperatures [34].When drying process proceeds further,heat transferred from propellant surface is greater than that consumed by evaporation, and internal temperatures start to increase.The propellant grains finally approach the drying temperature and the surface temperature difference is no longer noticeable.The temperature gradient between the surface and interior of the propellant grains are 1.2°C at 10 min and 1.4°C at 30 min.But the temperature gradient decreases to 1.1,1.0 and 0.1°C at 60,100 and 200 min, respectively.It means that the propellant samples will reach their steady state when the drying time is increased till 200 min.

    Fig.17 shows the temperature profiles at the drying temperature of 50°C and at different drying times such as 0,10,30,60,100,and 200 min.The nature of Fig.17 is the same as explained in Fig.16.The temperature gradient between the surface and the center of propellant grains is 3.8,2.5,2.0,1.6 and 0.1°C at drying times of 10,30,60, 100 and 200 min, respectively.It shows that the propellant temperature increases when drying time increases, and it reaches its steady state when the drying time is 200 min approximately.

    Fig.18 shows the temperature distribution on the centerline of face B at different time points and 50°C.When drying process proceeds further, the temperature difference ΔT (at the same time difference)is decreased over time,and the ΔT at the first 100 min is obviously more than that of later drying stage,showing the transferring heat is decreased over time.The results indicate that the heat transfer between the propellant grains and drying air is not uniform in the time dimension during the drying process.

    Fig.16.Three internal slices of predicated temperature profile of the stick propellant grains at drying temperature of 30 °C.

    Fig.17.Three internal slices of predicated temperature profile of the stick propellant grains at drying temperature of 50 °C.

    Fig.18.Temperature distribution on the centerline of face B at different time points and 50 °C.

    Fig.19.Average temperature of long stick propellant grains and effective solvent diffusion coefficient of the propellant surface at drying temperature of 30 °C and 50 °C.

    Fig.19 shows the average temperature of long stick propellant grains and effective solvent diffusion coefficient of the propellant surface at drying temperature of 30°C and 50°C, respectively.According to Vrentas-Duda model, the higher the temperature is,the faster solvent molecules diffuse, and the larger the effective solvent diffusion coefficient is.Combined with heat transfer and effective solvent diffusion coefficient results, after 100 min, the temperature of the long stick propellant grains is gently increased until the steady state is reached,the value of Deffof the propellant surface continues to increase over time at drying temperature of 30°C but decrease at drying temperature of 50°C, as shown in Fig.19.Both these results demonstrate the blocked layer near the surface has been formed in the time period of approximately 0-100 min under initial higher drying temperature and it results in the reduction of the effective solvent diffusion coefficient of propellant surface.

    For further elucidate its mechanism theoretically,it is deserved to evaluate the solvent behavior,Deffand temperature distribution inside propellant grains under different drying temperature combination.However,the numerical simulation of the combination of different temperature drying process is very complicated and is difficult to converge based on present numerical methods in this study.Moreover, comparison with the actual temperature, all values of the stick propellant grains by established model are predictive due to the limitation of present measure technology.Further efforts need focus on these issues to analyze experimentally the heat transfer effect of different drying conditions on the drying process of long stick propellants with large web thickness.

    5.Conclusions

    This study contributes valuable knowledge for the clarification of the formation mechanism of blocked layer near the propellant surface and drying behavior of long stick propellants with large web thickness.In this paper, a new 3D-CHMT model was successfully developed and it can be used to investigate the solvent behavior and temperature distribution within the stick propellant grains.

    The solvent behavior shows that the predicted change of the content of solvent during the drying process is decreased over time,and the solvent transfer can be divided into SF and MF stages,and the SF stage can be restrained to inhibit the formation of the blocked layer near the propellant surface.The value of Deffof the propellant surface first increases from 3.4 × 10-6m2/s to 5.3 × 10-6m2/s with the increase of temperature, and then decreases to 4.1 × 10-8m2/s at 60-100 min.The value of Deffof the propellant surface between 0-1.4 mm has a unique trend of change compared with the other regions,and it is much lower than that of the internal at 100 min.Meanwhile,the temperature of propellant surface increases rapidly at the SF stage(0-100 min)and then very slowly thereafter.Both the evolution of Deff and temperature distribution results demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and the thickness of the formed blocked layer is about 1.4 mm.To avoid the formation of blocked layer and improve drying quality of finial propellant products effectively, it would be better to adjust the drying parameters.For example, for long stick propellants with large web thickness, it should be initially dried at lower drying temperature (30-40°C) in 0-100 min and then dried at higher drying temperature(50-60°C)to reduce drying time for later drying process.The results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    The research was supported by the National Natural Science Foundation of China (Grant No.22075146).

    国语自产精品视频在线第100页| 久久国产精品影院| 久久久国产成人精品二区| 国产一区二区三区在线臀色熟女| 国产精品影院久久| 欧美+日韩+精品| 国产精品电影一区二区三区| 精品人妻一区二区三区麻豆 | 国产精品日韩av在线免费观看| 在线天堂最新版资源| 中文字幕久久专区| 亚洲av一区综合| av在线观看视频网站免费| 欧美激情国产日韩精品一区| 中文字幕免费在线视频6| 91久久精品国产一区二区成人| 中文字幕人成人乱码亚洲影| 小说图片视频综合网站| 久久久久久久亚洲中文字幕 | 18+在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| а√天堂www在线а√下载| 成年女人看的毛片在线观看| 真实男女啪啪啪动态图| 露出奶头的视频| 亚洲乱码一区二区免费版| 老熟妇乱子伦视频在线观看| 观看免费一级毛片| 国产在线男女| 国产精品国产高清国产av| 国产午夜福利久久久久久| 久久久久亚洲av毛片大全| 午夜福利18| a级一级毛片免费在线观看| 精品无人区乱码1区二区| 日本 欧美在线| 久久天躁狠狠躁夜夜2o2o| 男女下面进入的视频免费午夜| 亚洲专区国产一区二区| 成人欧美大片| 日本黄色片子视频| 97热精品久久久久久| 国产三级黄色录像| 精品久久久久久久久亚洲 | 午夜免费男女啪啪视频观看 | av天堂在线播放| 婷婷亚洲欧美| 国产精品电影一区二区三区| 久久草成人影院| 97超级碰碰碰精品色视频在线观看| 中文字幕熟女人妻在线| 亚洲欧美日韩无卡精品| 黄色配什么色好看| 中文字幕av在线有码专区| 精品日产1卡2卡| 69人妻影院| 国产精品一及| 亚洲 国产 在线| 国产乱人伦免费视频| 免费一级毛片在线播放高清视频| 成人亚洲精品av一区二区| 国产淫片久久久久久久久 | 精品一区二区三区视频在线| 日韩高清综合在线| 国产综合懂色| 亚洲人与动物交配视频| 久久久精品欧美日韩精品| 中文字幕人成人乱码亚洲影| 日本成人三级电影网站| 免费观看精品视频网站| 国产亚洲欧美98| 久久久久久大精品| 自拍偷自拍亚洲精品老妇| 国产成+人综合+亚洲专区| 99热这里只有是精品在线观看 | 噜噜噜噜噜久久久久久91| 婷婷精品国产亚洲av在线| 色5月婷婷丁香| 婷婷丁香在线五月| 久久久色成人| 国产成年人精品一区二区| 久久国产精品影院| 国产 一区 欧美 日韩| 亚洲欧美激情综合另类| 深爱激情五月婷婷| 综合色av麻豆| 女生性感内裤真人,穿戴方法视频| av国产免费在线观看| 国产亚洲精品av在线| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人av| 国产极品精品免费视频能看的| 久久午夜福利片| 欧美午夜高清在线| 色哟哟哟哟哟哟| 国内精品久久久久精免费| 国产亚洲欧美98| 偷拍熟女少妇极品色| 国产精品av视频在线免费观看| 久久国产精品影院| 欧美精品啪啪一区二区三区| 9191精品国产免费久久| 国产成人啪精品午夜网站| a在线观看视频网站| 一个人免费在线观看的高清视频| 日本黄大片高清| 国产精品久久久久久久久免 | avwww免费| 亚洲 欧美 日韩 在线 免费| 日韩国内少妇激情av| 内地一区二区视频在线| 亚洲性夜色夜夜综合| 国产精品久久久久久精品电影| av中文乱码字幕在线| 老熟妇乱子伦视频在线观看| 国产aⅴ精品一区二区三区波| 欧美日韩福利视频一区二区| 嫩草影院入口| 日日干狠狠操夜夜爽| 亚洲av.av天堂| 成人毛片a级毛片在线播放| 久久精品国产亚洲av天美| av女优亚洲男人天堂| 最后的刺客免费高清国语| 热99在线观看视频| 成年版毛片免费区| 中文字幕av在线有码专区| 日本三级黄在线观看| 非洲黑人性xxxx精品又粗又长| 国产淫片久久久久久久久 | 国产一区二区三区视频了| 好看av亚洲va欧美ⅴa在| 高潮久久久久久久久久久不卡| 老司机午夜十八禁免费视频| 亚洲aⅴ乱码一区二区在线播放| 他把我摸到了高潮在线观看| ponron亚洲| 老司机深夜福利视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 免费在线观看亚洲国产| 日本a在线网址| 日韩欧美一区二区三区在线观看| 欧美+亚洲+日韩+国产| 国语自产精品视频在线第100页| 久久欧美精品欧美久久欧美| 男女那种视频在线观看| 午夜激情福利司机影院| 日本 欧美在线| 欧美一区二区国产精品久久精品| 亚洲精品日韩av片在线观看| 欧美日韩综合久久久久久 | 人人妻人人澡欧美一区二区| 我的老师免费观看完整版| 亚洲欧美日韩高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 精品人妻熟女av久视频| 18禁裸乳无遮挡免费网站照片| 日本免费a在线| 久久亚洲精品不卡| 国产三级中文精品| 激情在线观看视频在线高清| av在线老鸭窝| 麻豆av噜噜一区二区三区| 在现免费观看毛片| 久久久久久久精品吃奶| 舔av片在线| 精品国产三级普通话版| 人人妻人人澡欧美一区二区| 波多野结衣高清无吗| 99久久九九国产精品国产免费| 女人被狂操c到高潮| 99久久精品国产亚洲精品| 成人午夜高清在线视频| 欧美黑人欧美精品刺激| 成人特级黄色片久久久久久久| 中文字幕高清在线视频| 免费一级毛片在线播放高清视频| av在线老鸭窝| 精品人妻熟女av久视频| 性色avwww在线观看| 狠狠狠狠99中文字幕| 最近最新免费中文字幕在线| a级毛片a级免费在线| 91麻豆精品激情在线观看国产| 精品一区二区免费观看| 久久婷婷人人爽人人干人人爱| 国产高潮美女av| 久久香蕉精品热| 757午夜福利合集在线观看| 国产av麻豆久久久久久久| av黄色大香蕉| 亚洲av电影不卡..在线观看| 嫩草影院入口| 欧美激情在线99| 长腿黑丝高跟| 波多野结衣高清无吗| 欧美精品啪啪一区二区三区| 久久国产精品人妻蜜桃| 极品教师在线视频| 国产美女午夜福利| 国产单亲对白刺激| 18禁黄网站禁片午夜丰满| 国产精品久久久久久久久免 | 青草久久国产| 国产一区二区在线av高清观看| 久久精品国产清高在天天线| 国产午夜精品久久久久久一区二区三区 | 亚洲成人精品中文字幕电影| 欧美+亚洲+日韩+国产| 精品福利观看| 婷婷精品国产亚洲av| 免费av不卡在线播放| 内射极品少妇av片p| 色5月婷婷丁香| 欧美绝顶高潮抽搐喷水| 久久久久久久久中文| 亚洲在线观看片| 国产大屁股一区二区在线视频| 国产69精品久久久久777片| 亚洲国产色片| 99riav亚洲国产免费| 成人午夜高清在线视频| 日韩欧美免费精品| 淫秽高清视频在线观看| www.999成人在线观看| 欧美黄色片欧美黄色片| 亚洲激情在线av| 美女大奶头视频| 国产精品综合久久久久久久免费| 成人特级黄色片久久久久久久| 变态另类成人亚洲欧美熟女| 亚洲人成伊人成综合网2020| 欧美xxxx性猛交bbbb| 好男人在线观看高清免费视频| 亚洲欧美日韩卡通动漫| 日韩欧美免费精品| 我要搜黄色片| 最近视频中文字幕2019在线8| 欧美性猛交黑人性爽| 国产精品三级大全| 日韩欧美免费精品| 床上黄色一级片| 亚洲av美国av| 成人特级av手机在线观看| 欧美性猛交╳xxx乱大交人| 国产成人影院久久av| 久久天躁狠狠躁夜夜2o2o| 黄色视频,在线免费观看| 亚洲最大成人av| 亚洲av日韩精品久久久久久密| 97超视频在线观看视频| 日本精品一区二区三区蜜桃| 亚洲av电影在线进入| 国产欧美日韩精品亚洲av| 日本精品一区二区三区蜜桃| 丰满的人妻完整版| 久久这里只有精品中国| 一区福利在线观看| 一a级毛片在线观看| 国产高清激情床上av| 久久亚洲精品不卡| 一进一出好大好爽视频| 他把我摸到了高潮在线观看| 国产精品久久久久久亚洲av鲁大| 成年版毛片免费区| 一个人看视频在线观看www免费| 国内精品久久久久精免费| 欧美黄色淫秽网站| 男女视频在线观看网站免费| 欧美在线黄色| 色噜噜av男人的天堂激情| 国产蜜桃级精品一区二区三区| 国产精品,欧美在线| www.999成人在线观看| 日日干狠狠操夜夜爽| 欧美成人a在线观看| 国产三级黄色录像| 亚洲av美国av| 欧美一区二区亚洲| 51午夜福利影视在线观看| 在线观看美女被高潮喷水网站 | 国产免费一级a男人的天堂| 三级男女做爰猛烈吃奶摸视频| 久久精品国产亚洲av涩爱 | 九九久久精品国产亚洲av麻豆| 亚洲国产精品999在线| 在线观看一区二区三区| www.999成人在线观看| 日韩欧美免费精品| АⅤ资源中文在线天堂| 免费高清视频大片| av在线天堂中文字幕| 日韩欧美免费精品| 亚洲精华国产精华精| 夜夜爽天天搞| 最后的刺客免费高清国语| av女优亚洲男人天堂| 亚洲一区高清亚洲精品| 99热这里只有是精品在线观看 | 国产欧美日韩精品一区二区| 亚洲国产精品合色在线| 亚洲国产日韩欧美精品在线观看| 国产精品亚洲美女久久久| 久久久久国产精品人妻aⅴ院| av天堂中文字幕网| 不卡一级毛片| 动漫黄色视频在线观看| 国产高清视频在线播放一区| 欧美在线黄色| 97超视频在线观看视频| 精品国产三级普通话版| 国产高清激情床上av| 日韩中字成人| 亚洲综合色惰| 国产成年人精品一区二区| a级毛片a级免费在线| 桃红色精品国产亚洲av| 男人狂女人下面高潮的视频| 国产久久久一区二区三区| 亚洲av成人精品一区久久| a级一级毛片免费在线观看| 免费av观看视频| 成人高潮视频无遮挡免费网站| 国语自产精品视频在线第100页| 亚洲成人免费电影在线观看| 精品一区二区三区人妻视频| 男人舔奶头视频| 精品日产1卡2卡| 人妻久久中文字幕网| 国产 一区 欧美 日韩| 中文字幕久久专区| netflix在线观看网站| 人人妻人人看人人澡| 欧美日本视频| 日韩精品中文字幕看吧| 亚洲中文字幕日韩| 此物有八面人人有两片| 午夜福利18| 国产精品av视频在线免费观看| 亚州av有码| 国产精品爽爽va在线观看网站| 久久精品综合一区二区三区| 一区二区三区免费毛片| 国产真实乱freesex| 亚洲国产色片| 国内久久婷婷六月综合欲色啪| 午夜视频国产福利| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品在线美女| 在线十欧美十亚洲十日本专区| 国产精品久久久久久久电影| 久久久久久九九精品二区国产| 午夜a级毛片| 俄罗斯特黄特色一大片| 高清毛片免费观看视频网站| 久久午夜亚洲精品久久| 色精品久久人妻99蜜桃| 一本一本综合久久| 色综合欧美亚洲国产小说| 久久国产精品人妻蜜桃| 亚洲国产日韩欧美精品在线观看| 亚洲经典国产精华液单 | 久久精品国产清高在天天线| 久久精品综合一区二区三区| 一个人观看的视频www高清免费观看| 桃红色精品国产亚洲av| 好男人在线观看高清免费视频| 国产精品久久久久久久久免 | 中文字幕精品亚洲无线码一区| 欧美另类亚洲清纯唯美| 国产激情偷乱视频一区二区| 成人国产综合亚洲| ponron亚洲| 波多野结衣高清无吗| 国内揄拍国产精品人妻在线| 亚洲成人中文字幕在线播放| 国产一区二区三区在线臀色熟女| 精品欧美国产一区二区三| 毛片女人毛片| 国产69精品久久久久777片| 高清在线国产一区| 天堂av国产一区二区熟女人妻| 亚洲色图av天堂| 色噜噜av男人的天堂激情| 亚洲黑人精品在线| 国产久久久一区二区三区| 窝窝影院91人妻| 麻豆成人av在线观看| 色在线成人网| 美女黄网站色视频| 99国产精品一区二区蜜桃av| 国产视频内射| 国产av麻豆久久久久久久| 精品不卡国产一区二区三区| 91字幕亚洲| 深夜精品福利| 国产亚洲精品久久久com| 亚洲一区高清亚洲精品| 99热这里只有是精品50| 精品一区二区三区av网在线观看| 国产探花在线观看一区二区| 亚洲精品一卡2卡三卡4卡5卡| 国产三级中文精品| 亚洲综合色惰| 日韩av在线大香蕉| 性欧美人与动物交配| 午夜免费男女啪啪视频观看 | 精品久久久久久久久久久久久| 国产欧美日韩精品一区二区| 亚洲中文字幕日韩| 国产91精品成人一区二区三区| 五月玫瑰六月丁香| 夜夜爽天天搞| 久久草成人影院| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 夜夜看夜夜爽夜夜摸| 色噜噜av男人的天堂激情| 欧美性猛交黑人性爽| 国产精品久久久久久人妻精品电影| 国产成人aa在线观看| 国产v大片淫在线免费观看| 午夜视频国产福利| av在线蜜桃| 亚洲av免费在线观看| 亚洲精品影视一区二区三区av| 亚洲国产色片| 国产av一区在线观看免费| 欧美激情久久久久久爽电影| 午夜免费成人在线视频| 少妇的逼好多水| 俺也久久电影网| 99热这里只有精品一区| 中文字幕久久专区| 999久久久精品免费观看国产| 日韩成人在线观看一区二区三区| 性欧美人与动物交配| 亚洲av电影不卡..在线观看| 午夜精品一区二区三区免费看| 亚洲第一区二区三区不卡| 一进一出好大好爽视频| 午夜免费激情av| 精品久久久久久久久久免费视频| 老熟妇仑乱视频hdxx| 听说在线观看完整版免费高清| 一卡2卡三卡四卡精品乱码亚洲| 国产老妇女一区| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 欧美一区二区精品小视频在线| 美女被艹到高潮喷水动态| 不卡一级毛片| 欧美黑人欧美精品刺激| 桃红色精品国产亚洲av| 国产69精品久久久久777片| netflix在线观看网站| 男人舔奶头视频| 久久久久久国产a免费观看| 99在线视频只有这里精品首页| 亚洲av成人精品一区久久| 老司机深夜福利视频在线观看| 午夜两性在线视频| 一区福利在线观看| 久久6这里有精品| 成年免费大片在线观看| 成人欧美大片| 亚洲 欧美 日韩 在线 免费| 免费av毛片视频| 午夜亚洲福利在线播放| 亚洲成av人片免费观看| 欧美色欧美亚洲另类二区| 免费一级毛片在线播放高清视频| 亚洲 欧美 日韩 在线 免费| 国产精品av视频在线免费观看| 精品免费久久久久久久清纯| 色视频www国产| 在线观看舔阴道视频| 国产午夜精品久久久久久一区二区三区 | 丰满人妻一区二区三区视频av| 久久久久国内视频| 校园春色视频在线观看| 99视频精品全部免费 在线| 国产91精品成人一区二区三区| 在现免费观看毛片| 国产精品人妻久久久久久| 久久久久亚洲av毛片大全| 美女xxoo啪啪120秒动态图 | 蜜桃亚洲精品一区二区三区| 午夜日韩欧美国产| 国产精品久久久久久久电影| 国产精品自产拍在线观看55亚洲| 国产日本99.免费观看| 免费观看精品视频网站| 亚洲中文字幕一区二区三区有码在线看| 亚洲av成人av| 精品无人区乱码1区二区| 欧美绝顶高潮抽搐喷水| 久久精品国产亚洲av香蕉五月| 午夜福利视频1000在线观看| 免费无遮挡裸体视频| 淫秽高清视频在线观看| 国内精品久久久久精免费| 一a级毛片在线观看| 国产美女午夜福利| 国产亚洲精品av在线| 一级a爱片免费观看的视频| 精品久久久久久久末码| 亚洲av成人精品一区久久| 国内少妇人妻偷人精品xxx网站| 一个人免费在线观看的高清视频| 可以在线观看的亚洲视频| 18+在线观看网站| 久久国产乱子伦精品免费另类| 欧美zozozo另类| 免费观看的影片在线观看| 欧美一级a爱片免费观看看| 亚洲 国产 在线| 色综合站精品国产| 国产精品久久久久久久久免 | 亚洲不卡免费看| 精品国产亚洲在线| 国内毛片毛片毛片毛片毛片| 亚洲久久久久久中文字幕| 国产午夜福利久久久久久| 久久精品综合一区二区三区| 亚洲美女搞黄在线观看 | 99久久久亚洲精品蜜臀av| 亚洲av电影在线进入| 首页视频小说图片口味搜索| 女同久久另类99精品国产91| 精品久久国产蜜桃| 首页视频小说图片口味搜索| 午夜视频国产福利| avwww免费| 给我免费播放毛片高清在线观看| 亚洲av中文字字幕乱码综合| 精品久久久久久,| 国产精品爽爽va在线观看网站| 亚洲欧美激情综合另类| 又黄又爽又免费观看的视频| 久久这里只有精品中国| 青草久久国产| 人人妻人人看人人澡| 国产午夜精品久久久久久一区二区三区 | 国产一区二区激情短视频| 亚洲人成伊人成综合网2020| 久久精品国产亚洲av涩爱 | 精品一区二区免费观看| 最后的刺客免费高清国语| 久久久久久久亚洲中文字幕 | 日韩欧美精品v在线| 性色av乱码一区二区三区2| 高清在线国产一区| 婷婷精品国产亚洲av| 亚洲熟妇熟女久久| 女人被狂操c到高潮| 在线观看av片永久免费下载| 一a级毛片在线观看| 少妇被粗大猛烈的视频| 欧美乱妇无乱码| 全区人妻精品视频| 日韩高清综合在线| av福利片在线观看| 天天躁日日操中文字幕| 国产精品国产高清国产av| 亚洲性夜色夜夜综合| 日日摸夜夜添夜夜添小说| 男女床上黄色一级片免费看| 真实男女啪啪啪动态图| 欧美午夜高清在线| 亚洲第一电影网av| av欧美777| 欧美另类亚洲清纯唯美| 亚洲国产精品sss在线观看| 99久久精品国产亚洲精品| 久久精品人妻少妇| 欧美日韩福利视频一区二区| 国产一级毛片七仙女欲春2| 高清毛片免费观看视频网站| av福利片在线观看| 毛片一级片免费看久久久久 | 美女cb高潮喷水在线观看| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 欧美精品啪啪一区二区三区| 一个人免费在线观看电影| 欧美日韩国产亚洲二区| 乱人视频在线观看| 9191精品国产免费久久| 噜噜噜噜噜久久久久久91| 在线国产一区二区在线| 久久久久久久午夜电影| 国产av麻豆久久久久久久| 国产一区二区在线av高清观看| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 日韩亚洲欧美综合| 12—13女人毛片做爰片一| 国产亚洲精品久久久com| 永久网站在线| 日本 欧美在线| 午夜日韩欧美国产| 天堂动漫精品| 夜夜躁狠狠躁天天躁| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人成电影免费在线| 国产一级毛片七仙女欲春2| 国产精品嫩草影院av在线观看 | 亚洲久久久久久中文字幕| 色综合亚洲欧美另类图片| 中文字幕人妻熟人妻熟丝袜美| 91麻豆精品激情在线观看国产| 天堂影院成人在线观看| 日韩欧美三级三区| 国产色婷婷99|