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

    Effect of Pore Distribution on Melting Behavior of Paraffin in Fractal Metal Foam?

    2018-11-24 07:36:02JiaYiZheng鄭佳宜TaoWang王濤TaoTaoChen陳濤濤andYanShunYu余延順
    Communications in Theoretical Physics 2018年10期
    關(guān)鍵詞:王濤

    Jia-Yi Zheng(鄭佳宜), Tao Wang(王濤),Tao-Tao Chen(陳濤濤),and Yan-Shun Yu(余延順)

    1MIIT Key Laboratory of Thermal Control of Electronic Equipment,School of Energy and Power Engineering,Nanjing University of Science&Technology,Nanjing 210094,China

    2Shanghai Institute of Satellite Engineering,Shanghai 200240,China

    AbstractAn unsteady model of melting heat transfer in metal foam filled with paraffin is developed and numerically analyzed.In the model,the pore distribution of metal foam is described by Sierpinski fractal.By this fractal description, six types of metal foams with the identical porosity,different pore distributions are reconstructed.The effect of pores distribution on the heat transfer performance of the paraffin/copper foam composite is emphatically investigated.The solid-liquid distribution,the evolution of the melting front,the dynamic temperature response and the total melting time in these six fractal structures are all examined and analyzed.The results indicate that,the pores distribution has a significant impact on the melting heat transfer of the paraffin in the metal foam.When the first-level fractal pore is located near the heat source boundary,the thermal hysteresis effect of the paraffin on the solid matrix heat transfer is more significant.If the cross-sectional length of the pores along the normal direction of heat transfer is larger,the heat transfer barrier of the solid matrix is greater.It is favorable for the paraffin heat transfer when the pores specific surface area is larger.

    Key words:fractal,foam,pores distribution,melting,heat transfer

    1 Introduction

    Nowadays the energy storage technology becomes a hot research topic on the background of the worldwide shortages of energy supplies.Energy storage technology allows excess thermal energy to be collected for later use,and the stored energy to be released when necessary,which compensates for the intermittent and instability of renewable energy(e.g.,solar energy).Phase change material(PCM),a kind of promising energy storage material,has widespread applications in the fields of building energy conservation,battery thermal management,[1]heating and air-conditioning,because of its good chemical stability and high energy storage density.[2]However,the thermal conductivity of PCM is generally low,which limits its further development in the field of energy storage technology.Metal foam can effectively enhance the thermal performance of PCM,i.e.,it not only accelerates the phase change rate of PCM,but also changes the evolution of the melting or solidification front.The pore structure of the metal foam is closely related to its heat transfer performance.Especially,the effect of the pore distribution on the phase change behavior of PCM in metal foam still remains less understood.

    There are a number of approaches that are adopted in order to enhance the thermalconductivity of PCM by experimental observations and numerical simulations,including dispersing expanded graphitein PCM,[3?4]compositing PCM/carbon foam,[5]or compositing PCM/metal foam,etc.Among them,metal foam is a common and effective additive for its large thermal conductivity to improve the average thermal performance of PCM and metal foam.A quantitative experiment was carried out to show that the addition of metal foam can shorten the heat storage time(74.4%of pure paraffin)and improve the temperature uniformity of the PCM.[6]Zhao[7]also made a quantitative experiment to prove that the presence of metal foam also can increase the overall heat transfer rate by 3–10 times compared with the pure paraffin wax,and the reason for the difference in heat transfer rate may attribute to the structures and properties of metal foam.The pore structure of metal foam(e.g.,porosity,pore size,PPI(number of pores per inch),pore distribution)is irregular and complex,which aroused scientific interest in recent years and became a hot research topic.Jin[8]carried out a visualized study to explore the pore size effects on the melting heat transfer of paraffin wax saturated in copper foam.Sriharsha[9]also investigated the effect of pore size and porosity on the performance of PCM infiltrated metal foam,especially when the pore size reduces to less than 100μm.Deng[10]studied the role of porosity,fractal dimension and pore distribution on the melting behaviors of PCM in porous metal foam.The above researchers all considered the pore structure of metal foam,which plays an important role in improving the thermal conductivity of PCM.

    For the complex porous structure of metal foam,the structural parameters of metal foam can be obtained by geometry reconstruction,CT scan,[5]etc.Although the CT scan can accurately construct two-dimensional or three-dimensional pore structure of metal foam,it has the drawbacks of high cost,low efficiency and complicated correction function.Hence,Feng[11]and Giorgio[12]constructed the open-cell metal foam using sphere-centered tetrakaidecahedron and carried out the heat transfer simulation between the metal foam and PCM in porescale.In addition,a great number of scientific researches proved that the microstructures of naturally formed and partially artificially made porous media are fractal.[13]They all have the features of multi-scales and self-similar(strict[14?15]or statistical[16?18]),thus,they can be characterized by fractal geometry.Fractals has plenty of important applications in fields of single phase flow,multiphase flow,heat transfer,reaction,et al.[19?22]Zhang[23]introduced the Sierpinski fractal to construct a porous metal foam and numerically analyzed the solidification heat transfer processes of PCM in fractals,the results show that the presence of fractal metal foam provides a fast heat flow channel both horizontally and vertically,which induces the enhancement of heat transfer between PCM and metal foam.Deng[10]characterized the pore distribution of metal foam by fractal Brownian motion and investigated the melting behaviors of PCM in metal foam with a particular focus on the role of pore distribution.Therefore,the fractal geometry is a useful and effective method to characterize the pore structure of metal foam and investigate the role of pore structure on melting behaviors of PCM embedded in metal foam in detail and depth.

    In the previous researches,the porosity and pore size were generally regarded as the key parameters to characterize the pore structure of porous media.[5,8?9]However,the recent researches indicate that the fluid flow and heat transfer in porous structures are also affected by the pores distribution even if the porosity is identical.[10,16]It is noteworthy that the pores distribution is a crucial parameter that should not be neglected when the melting heat transfer performance is studied.For the random pores distribution of metal foam,the pores inevitably contact each other,which thereby forms the open-cell or closed-cell metal foam.Therefore,six types of metal foams with diverse pores distributions(including open-cell and closed-cell metal foams)are reconstructed using Sierpinski fractal.Based on the fractal description,an unsteady melting heat transfer model of paraffin in copper foam is developed and numerically analyzed.We mainly focus on the role of pores distribution on melting heat transfer performance of the paraffin/copper foam composite.In addition,the solid-liquid distribution,the evolution of the melting front,the dynamic temperature response,and the total melting time of these six fractal structures are compared and analyzed.

    2 Porous Structure Reconstruction

    The fractal geometry with a feature of multiscale and self-similar is commonly introduced to characterize the pore structures of porous media. A typical twodimensional Sierpinski carpet is used to reconstruct the geometrical structures of porous metal foam.As shown in Fig.1,the construction process of the Sierpinski carpet can be described as follows:firstly,the reconstruction begins with a square that is divided into 9 congruent subsquares in a 3-by-3 grid,and then the central sub-square is removed where PCM can be filled in.After that,the same procedure is applied recursively at the remaining 8 sub-squares.In this way,different levels of Sierpinski carpets are reconstructed.The above recursive algorithm is theoretically infinite.Considering the computational expense and the exhibition effect,the third level fractal is chosen to study the influence of pores distribution on the thermal properties of metal foam filled with PCM.It is easy to obtain the fractal dimension of the Sierpinski carpet D=ln8/ln3=1.893.It is worth noting that the black and white regions represent the solid matrix and pores,respectively.In the simulation,the solid matrix is made of copper,while the pores are filled with paraffin.

    Fig.1 Sierpinski fractals.

    Fig.2 Six types of Sierpinski fractals.

    In order to investigate the role of pores distribution on heat transfer processes of copper foams filled with paraffin,six types of pores configurations of copper foams are reconstructed(see Fig.2).It is obviously shown that the porosity and the fractal dimension of the six Sierpinski fractals are identical,while the pores and solid matrix distributions of which are totally different.That is,each pore is distributed independently in the solid matrix,which forms the closed-cell structure(see Fig.2(a)),while some pores are in touch with each other,and even some pores are in touch with the boundary of the solid matrix,which forms the open-cell structures.Therefore,the specific surfaces of the six Sierpinski fractals are not exactly the same.

    3 Melting Performances of Paraffin in Copper Foams

    3.1 Theoretical Model

    For purpose of theoretically studying the melting performances of paraffin in the copper foams,a twodimensional unsteady heat transfer model accompanied with melting of paraffin embedded in pores of copper foams is established.As shown in Fig.3,except for the constant wall temperature boundary(Th)at the bottom,the other boundaries are adiabatic.Because if the left boundary is chosen as the heat source boundary,the left part of pores filled with paraffin melt first.The left boundary heat can produce local convection of the liquid paraffi n,that is,the upper left corner of the solid paraffin melts faster than that of the lower left corner.But if the thermal boundary is below,the effect of the liquid paraffin convection can be eliminated.The length of the square computational domain H is 300 mm.The coefficient of thermal conductivity(λ),the specific heat(cp)and density(ρ)of the copper foam and paraffin,as well as the phase transition temperature(Tm)and latent heat(Lat)of paraffin are listed in Table 1,respectively.

    Fig.3 Schematic of heat transfer in metal foam filled with PCM.

    Table 1 Physical properties of copper foam and paraffin.

    In order to simplify the melting process of composite paraffin/copper foam,some appropriate assumptions are made as follows:

    (i)The physical properties of copper foam is constant,the phase transformation occurs in a temperature interval between Tm?T/2 and Tm+?T/2,where?T=2 K.

    (ii)The paraffin is fully filled in pores of the copper foams,i.e.,there is no gap between the paraffin and copper foam,or the contact thermal resistance between the paraffin and copper foam can be ignored.

    (iii)The melting process is divided into three stages,i.e.,the solid zone,mushy zone and liquid zone.The physical property parameters do not change with the temperature in the solid zone or liquid zone,and vary linearly with the temperature in the mushy zone.

    Heat conduction with phase change is more complicated than general variable coefficient heat conduction,because there is a time-varying melting interface in the mushy region.Plenty of methods are applied to study the phase change problem.[24?26]The apparent heat capacity method is an effective method that treats the melting heat of paraffin as a large heat capacity within a small temperature range.The magnitude of the apparent heat capacity is determined by the latent heat of melting and the range of the melting temperature.This method can transform the governing equations and conservation con-ditions of each phase into a single phase nonlinear heat conduction equation applied to the whole region.Based on the above assumptions,the model of the unsteady heat transfer accompanied with the melting of paraffin in copper foam is established.The governing equations are as follows:

    where ρ is the density,T is the temperature,λ is the thermal conductivity.The physical property parameters of ρ,λ and cpin the mushy zone are acquired via the apparent heat capacity method,which are shown as follows:

    where cpis the equivalent volume specific heat capacity,when the temperature is less than Tm??T/2,cprepresents the specific heat of solid paraffin.When the temperature is more than Tm+?T/2,cprepresents the specific heat of liquid paraffin.When the temperature is between Tm??T/2 and Tm+?T/2,the latent heat of the paraffi n is regarded as an apparent heat capacity between solid paraffin and liquid paraffin.

    where θ is a smoothed function,representing the fraction of solid paraffin.That is,θ is equal to 1 when the paraffin is solid,while it is 0 when the paraffin is liquid.αmis the mass fraction.The sub-script s and l represent solid state and liquid state of paraffin,respectively.

    3.2 Initial and Boundary Conditions

    At the initial stage,there is no temperature difference between the copper foam and paraffin,they are in a thermal equilibrium with the ambient temperature,T0=293.15 K.As shown in Fig.3,when the melting process begins,the bottom boundary is set as a constant high temperature Th,while the upper,left and right boundaries are assumed to be adiabatic.If it is not specified,the initial temperature difference(Th?T0)is set as 30 K.

    Initial condition:

    Boundary conditions:

    There is a thermal balance between the heat conduction of the copper foam and the melting heat of paraffin in the pores,and they are coupled via the solid boundary of the pores,i.e.,the continuities of temperature and heatflux occur at the solid boundary of the pores.

    where ? represents the interface between the copper foam and paraffin,n is the normal direction of heat flux.

    3.3 Numerical Simulation

    The finite element method is applied to solve the heat conduction equation and apparent heat capacity equations.According to the pore characteristics of the Sierpinski fractals,two-dimensional structured meshes(quadrilaterals elements)are built to all the computational region.Considering the existence of a large temperature gradient in the pores filled with PCM during the heat transfer process,it is necessary to refine meshes in the pores region.Meanwhile,a mesh independence test is conducted with a series of mesh sizes,because the mesh size greatly affects the accuracy of the numerical simulation and the speed of calculation.

    3.4 Model Verification

    In order to verify the above theoretical model,the predicted results are compared with Ref.[23],which investigated the solidification heat transfer process of PCM in the standard Sierpinski fractal metal foam(see Fig.4(a)).As shown in Fig.4(b),the liquid fraction of PCM in our model agrees well with that in Ref.[23],which verifies the rationality of the theoretical model in the paper.

    Fig.4 Comparison of liquid fraction between our model and Ref.[23].

    4 Results and Discussions

    Based on the above theoretical model and the construction of different pore structures,the heat transfer characteristics of composite paraffin/copper foam are analyzed and discussed,especially the role of pores distribution on heat transfer process is investigated.Then,we will discuss from the two aspects,the melting process of the paraffin in copper foam and its influence factors analysis.

    4.1 Melting Behaviors

    Fig.5 (Color online)Solid-liquid distribution of paraffin in pores of copper foam during melting process.

    Copper foam can significantly enhance the heat transfer of paraffin,because the thermal conductivity of copper foam is much greater than that of paraffin.Figure 5 shows the phase interface evolution of the different fractals.The six types of constructed fractals have the same porosity,fractal dimension and different pores distributions.Obviously,the melting characteristics of the paraffin in pores of the copper foams are quite different.The melting rate of paraffin can reflect the heat-conducting property of fractal porous structures,which can be analyzed by the position of the maximum pore and its specific surface area.According to the distance between the maximum pore and the bottom constant temperature boundary,it is divided into low(Figs.5(b)and 5(e)),middle(Fig.5(d))and high(Figs.5(c)and 5(f))positions.Similarly,according to the specific surface area of the maximum pore,it is divided into big(Figs.5(e)and 5(f)),medium(Figs.5(b),5(c),and 5(d))and small(Fig.5(a))pores.When t=10000 s,the solid fraction of Fig.5(f)is the highest,because the big pores locate at the high position,that is,the largest area of pores is filled with lower thermal conductivity of paraffin at the position where there is a relatively small temperature gradient.The solid fraction of Fig.5(e)is followed by the solid fraction of Fig.5(f),because although the specific surface area of Fig.5(e)is equal to that of Fig.5(f),the maximum pores of Fig.5(e)is located on the thermal boundary where the temperature gradient is the largest,thus the paraffin can perform the phase transition at the maximum efficiency.Then the solid fractions of Figs.5(c)and 5(d)are approximately the same and followed by the solid fraction of Fig.5(e),for the maximum pore area of Figs.5(c)and 5(d)is less than that of Figs.5(e)and 5(f),the contact area of paraffi n and copper foam is larger than that of Figs.5(c)and 5(d),so the efficiency of heat exchange between the copper foam and paraffin is higher.As for which heat transfer effect of Figs.5(c)and 5(d)is better,further data analysis is illustrated later.Finally the solid fractions of Figs.5(a)and 5(b)are comparable and the lowest,for although the surface of the maximum pore of Fig.5(b)is larger than that of Fig.5(a),the maximum pore of Fig.5(b)located at the boundary with the maximum temperature gradient,so the solid fractions of Figs.5(a)and 5(b)are basically the same.Based on the above analysis,we can conclude that the melting rate of paraffin is lower when the specific surface area of the maximum pore is larger and the distance away from the constant temperature boundary is farther.

    As shown in Fig.6,the paraffin melting front in each pore moves inward during the melting process.The paraffi n in the third level pores from the bottom to top boundary melts first,then the paraffin in the second level pores from the bottom to top boundary melts,the paraffin in the first level pores melts lastly.Noteworthily,the pores distribution plays an important role in the paraffin melting process.When t=2500 s,the blue regions are almost hardly been seen in the second level pores in Figs.6(a)and 6(b),from which we can infer that the pores structures of Figs.6(a)and 6(b)is very beneficial for heat conduction.When t=10000 s,it is illustrated that the size of the low temperature zone is the biggest in Fig.6(f),which indicates that the thermal conductivity of pores structures in Fig.6(f)is the worst.According to this rule,the thermal conductivity of pores structures in Fig.6(e)is better than that of Fig.6(f),the thermal conductivity of pores structures in Figs.6(c)and 6(d)are almost the same,which are better than that of Fig.6(e).

    Fig.6 (Color online)Temperature evolution of the six fractals during melting process.

    4.2 Influence Factors Analysis

    In order to get more detailed and deeper understanding of the paraffin melting behaviors in copper foams with different pores structure,more data is collected for comparative analysis.Figure 7 shows the temperature evolution of the representative points(1~13 points)in the physical model(see Fig.3).Why these thirteen points?Firstly,along the direction of temperature gradient,three typical regions(see Figs.3(a),3(b),and 3(c))are selected,which contain the second-level and third-level fractals.Then,similarly along the direction of temperature gradient,twelve typical positions(see No.2–No.13 in Fig.3)are selected in the above mentioned three typical regions.Each typical area contains two points in the center of the second-level and third-level fractals,and the other two points are in the copper foam.In addition to the twelve positions,the last point(see No.1 in Fig.3)is located in the center of the first-level fractal.

    As shown in Fig.7,according to the characteristics of these curves,the thirteen curves can be divided into four groups for analysis and discussion.Curve 1 is much different from other curves,because it has the largest pore area where the maximum amount of paraffin is filled,so the heat transfer is the slowest and the temperature lags behind that of in the metal skeleton and the other second or third-level pores.It is seen that when t=5000 s,the metal skeleton and all other second or third-level pores almost reach thermal equilibrium,only point 1 has not yet completed phase change process,there is also a temperature difference of about 23?C.However,it can be deduced that if the time is long enough,the temperature of all points will reach thermal equilibrium.Curves 5,9,13 represent the center temperature of the second-level fractal pores.It is seen that along the temperature gradient,the time of phase change heat storage of curves 5,9,13 becomes shorter,which are about 2250 s,1750 s,and 1250 s,respectively.This is due to the heat capacity of the paraffin is much higher than that of the copper foam,thus the hysteresis of heat transfer of the paraffin is significant.Therefore,the three curves are quite different in the figure.

    Fig.7(Color online)Temperature evolution of the representative points in physical model.

    Another group curves 3,7,11 represent the center temperature of the third-level fractal pores.Similar to the second-level fractal pores,along the temperature gradient,the time of phase change heat storage of curves 3,7,11 becomes shorter,which are about 600 s,550 s,and 250 s,respectively.This is because the pore area of the thirdlevel fractal pores are smaller than that of the second-level fractal pores,therefore,the corresponding paraffin content and heat capacity are also smaller,and the time of phase change heat storage of the third-level fractal pores is shorter than that of the second-level fractal pores.The last group curves 4,2/6,8/10,12 are all selected in the copper foam.Along the direction of temperature gradient,points 4 and 2 belong to the region(a),points 6 and 8 belong to the region(b),points 10 and 12 belong to the region(c)(see Fig.3).Owing to the thermal conductivity of the copper foam is as high as 398 W·(m·K)?1,so the curves 4 and 2,curves 6 and 8,curves 10 and 12 almost coincide,respectively.Furthermore,with the end of phase change process in the second-level fractal pores,the temperature difference between the paraffin in the secondlevel fractal pores and the corresponding metal zone will be reduced,and finally reach the thermal equilibrium(at about 4500 s).

    The thermal conductivities of the copper foams with different pores structures are not easy to distinguish by the previous cloud images,such as Figs.5(a)and 5(b),Figs.5(c)and 5(d),Figs.6(a)and 6(b),Figs.6(c)and 6(d).But it can be clearly compared with the line graph that is shown in Fig.8.Figure 8 shows the average temperature evolution of the paraffin and copper foam in the six types of fractal structures(see Fig.2).Due to the thermal conductivity of copper is far greater than that of paraffin,so the twelve curves in Fig.8 are mainly divided into two parts,the above hollow symbols are for the copper foam,the lower solid symbols are for the paraffin.It is illustrated that the heat transfer rates of the copper foams are large at the initial stage(0~1200 s),and the time of thermal balance of the copper foams(about 4000 s)is much shorter than that of the paraffin(>10000 s).

    Fig.8 (Color online)Effect of pores distribution on average temperature of paraffin and copper foam.

    For heat transfer of solid matrix,the heat transfer rate and temperature of curves Figs.(8c)and(8f)are quite similar and relatively large in the initial stage(0~4000 s).The reason is that the first-level fractal pore of Figs.2(c)and 2(f)that block the heat transfer are located far away from the constant temperature boundary,so the heat passes rapidly through the copper foam.On the contrary,the first-level fractal pore of Figs.2(b)and 2(e)that block the heat transfer are distribute on the constant temperature boundary,which is the big obstacle to the heat transfer,therefore,the heat transfer rate and temperature of curves Figs.(8b)and(8e)are smaller than other porous structures in the initial stage(0~4000 s).However,the cross section length of the first-level fractal pore of Fig.2(e)in the normal direction of heat flux is greater than that of Fig.2(b),so the heat transfer rate of Fig.2(e)is less than that of Fig.2(b).Next,Figs.2(a)and 2(d)with the first-level fractal pore locating in the middle are analyzed.Due to the feature of porous structure of Fig.2(d)is each level of fractal pores are connected,while the feature of porous structure of Fig.2(a)is that pores at each level are independent of each other.Therefore,the cross section length of the fractal pores of Fig.2(d)in the normal direction of heat transfer is greater than that of Fig.2(a),the barrier effect of pore structure of Fig.2(d)on heat transfer is greater than that of Fig.2(a),that is,the heat transfer efficiency of Fig.2(a)is greater than Fig.2(b).It is worth noting that after 5000 s,the six types of solid matrix reach thermal equilibrium.In fact,great heating rate of the solid matrix demand for the first-level fractal pore is far away from the heat source boundary.The transversal length of pores in the normal direction of heat flux increases,which will result in the greater barriers to heat transfer of solid matrix.

    For heat transfer of paraffin,there is minor differences between Figs.2(a)and 2(b),2(c)and 2(d)on the heat transfer performance of paraffin in the cloud chart.It is clearly seen that the heat transfer performance of paraffin of curve(8a)is slightly better than that of(8b).Similarly,the heat transfer performance of paraffin of curve(8c)is slightly better than that of(8d).It is because Fig.2(b)has the connection of each level of fractal pores,resulting in the specific surface area of Fig.2(b)is less than that of Fig.2(a).Although the specific surface areas of Figs.2(c)and 2(d)are identical,the contact length of Fig.2(c)with the heat flux direction is longer than that of Fig.2(d).If the specific surface area of pores is large,the heat is more likely to be transferred from the surface of the pores to the paraffin,which is conducive to the storage of heat in the paraffin.The average paraffin temperature of Fig.2(f)is the lowest,because the paraffin in the first-level fractal pore is at the top of porous media where the temperature gradient is relatively small,and the area of the connected pores is the largest.Although the connected pore area of Fig.2(e)is equal to that of Fig.2(f),Fig.2(e)is at the constant temperature boundary with the maximum temperature gradient.It is interesting that before 2500 s,the heat transfer performance of the paraffin in Fig.2(f)is better than that in Fig.2(e),but after 2500 s,the heat transfer performance of the paraffin in Fig.2(e)is better than that in Fig.2(f).The reason is that the solid matrix of Fig.2(f)has the best thermal conductivity,which will bring heat quickly to the paraffin.While along the heatflux upward,the specific surface area of the lower part of Fig.2(f)is greater than that of Fig.2(e).Therefore,the heat will quickly heat up the paraffin in the second and third-level fractal pores.After 2500 s,when all the paraffin in the second and third-level fractal pores is melt,the first-level fractal pore of Fig.2(e)is at the constant temperature boundary with the maximum temperature gradient,so the paraffin heat transfer of curve Fig.8(e)catches up that of Fig.8(f).Therefore,we can conclude that the first-level fractal pore is at the boundary with the maximum temperature gradient,and the larger the specific surface area of pores,the more favorable for the heat transfer of paraffin.

    Fig.9(Color online)Temperature evolution of the firstlevel fractal central point in the six fractal structures.

    The heat storage effect of the six fractal structures can be determined through temperature variation of the paraffin at the center of the first-level fractal pore.Figure 9 illustrates the temperature evolution of the first-level fractal central point in the six fractal structures.When the constant temperature heat source is loaded,the paraffi n temperature at the first-level center of each fractal increases.The trend of the increase and the average temperature variation of the paraffin in all pores(the solid symbol curve in Fig.8)are basically the same,but the degree of the variation is different.Before t=12750 s,curves Figs.9(a)and 9(b)almost coincide. Whereas,curves Figs.9(c),9(d)and 9(e)almost coincide in the whole computational time.This is because although thefirst-level pores of Figs.2(c),2(d)are different from that of Fig.2(e)in shape,the size of the unmelt region is quite the same,thus the temperature variations are also quite the same.It is worth noting that curves Figs.9(e)and 9(f)do not intersect as solid symbol curves Figs.(8e)and(8f),because the temperatures of the paraffin in the firstlevel fractal pore linearly increase when they are observed individually.In addition,there is the most unmelt paraffi n in the first-level fractal pore of Fig.2(f),which results in the biggest temperature lag.Therefore,curve Fig.9(f)is the lowest.After t=12750 s,curves Figs.9(a)and 9(b)steep rise,because the paraffin in Figs.2(a)and 2(b)completely finish the phase transition at that time,that is,the paraffin melts from solid to liquid.Due to the thermal conductivity of the liquid paraffin is much lower than that of the copper foam,and the liquid paraffin of Fig.2(b)is in the thermal boundary with the largest temperature gradient,the temperature variation of curve Fig.9(b)is not larger than that of curve Fig.9(a).

    Fig.10(Color online)Comparison of total melting time in the six fractal structures.

    Under the same constant temperature boundary condition,if more time is required for the paraffin to complete the melting process,it indicates the thermal conductivity of porous fractal structure is poor and the energy storage efficiency of the paraffin/copper foam composite is low.Figure 10 shows the comparison of the total melting time in the six paraffin/copper foam composites.In fact,the histogram(Fig.10)corresponds to the previous cloud chart(Figs.5 and 6)and line graph(Figs.8 and 9),which can accurately and quantitatively compare the heat transfer performance of the six fractal structures.It is illustrated that Fig.2(f)has the largest pores connectivity area,that is,the pores specific surface area of Fig.2(f)is the smallest,and the first-level fractal pore of Fig.2(f)is at the top region far from the constant temperature boundary,so the paraffin has the longest melting time(t=48204.3 s),which is nearly 5 times as much as Fig.2(a)and more than 2 times as much as Fig.2(e).The pores connectivity area of Fig.2(e)is the same as that of Fig.2(f),but the first-level fractal pore of Fig.2(e)is at the constant temperature boundary with the maximum temperature gradient,therefore the melting time of the paraffin in Fig.2(e)(t=20080.4 s)is much smaller than that of Fig.2(f).The pores connectivity area of Fig.2(c)is the same as that of Fig.2(d),and less than that of Figs.2(e)and 2(f),but because of the difference positions of the first-level fractal pore,that is,the first-level pore of Fig.2(c)is at the top,while that of Fig.2(d)is in the middle,so the melting time of the paraffin in Fig.2(c)is more than that of Fig.2(d),and less than that of Figs.2(e)and 2(f).The pores connectivity area of Figs.2(a)is less than Fig.2(b),and the pores connectivity area of Fig.2(b)is equal to that of Figs.2(c)and 2(d),which indicates the pore specific surface area of Fig.2(a)is the largest,and the heat transfer performance of Fig.2(a)is the best,but due to the first-level pore of Fig.2(a)is in the middle,while that of Fig.2(b)is at the bottom with the maximum temperature gradient,therefore the melting time of the paraffin in Fig.2(b)(t=10897.2 s)is slightly longer than that of Fig.2(a),but less than that of others.

    5 Conclusion

    An unsteady melting heat transfer model of paraffin in copper foam is developed and numerically analyzed.Six types of copper foam structures are reconstructed by Sierpinski fractal geometry.The effect of pores distribution on the heat transfer performance of the paraffin/copper foam composite is emphatically studied.The solid-liquid distribution,the evolution of the melting front,the dynamic temperature response,and the total melting time in these six fractal structures are all compared and investigated.The main conclusions are as follows:

    (i)Although the porosity and the fractal dimension of the six Sierpinski fractals are identical,their pores and solid matrix distributions,and the paraffin melting behaviors are totally different.The influence of the pores distribution on the melting heat transfer is significant.

    (ii) If the first-level fractal pore is located far from the heat source boundary,the effect of the paraffin on the heat transfer of solid matrix is weaker,and the temperature rising rate of the solid matrix is greater.If the cross-sectional length of the pores along the normal direction of heat transfer is larger,the heat transfer barrier of the solid matrix is greater.

    (iii)It is favorable for the paraffin heat transfer when the first-level fractal pore is closer to the heat source boundary with the largest temperature gradient,and the pores specific surface area is larger,that is,the pores connectivity area is smaller.The melting time of paraffin in Fig.2(a)is the shortest,while that of Fig.2(f)is the longest.

    猜你喜歡
    王濤
    綿師學(xué)人
    ——王濤
    Review of a direct epitaxial approach to achieving micro-LEDs
    “雞兔同籠”問題解法探析及思考
    Nonlinear excitation of a geodesic acoustic mode by reversed shear Alfvén eignemodes
    Transition to chaos in lid–driven square cavity flow?
    王濤油畫作品
    大眾文藝(2020年23期)2021-01-04 08:48:40
    王濤 李佳星作品
    大眾文藝(2020年22期)2020-12-13 11:37:16
    ONE-DIMENSIONAL VISCOUS RADIATIVE GAS WITH TEMPERATURE DEPENDENT VISCOSITY?
    Study on the Reduced Traffic Congestion Method Based on Dynamic Guidance Information?
    王濤作品
    一边摸一边抽搐一进一出视频| av女优亚洲男人天堂| 久久97久久精品| 美女扒开内裤让男人捅视频| 国产麻豆69| 成人黄色视频免费在线看| 亚洲精品自拍成人| 99精国产麻豆久久婷婷| 人人妻人人爽人人添夜夜欢视频| 一级爰片在线观看| www.自偷自拍.com| 少妇 在线观看| 在现免费观看毛片| 91精品伊人久久大香线蕉| 9色porny在线观看| 多毛熟女@视频| 国产亚洲欧美精品永久| 免费不卡黄色视频| 精品一区二区三区四区五区乱码 | 天天躁日日躁夜夜躁夜夜| 色综合欧美亚洲国产小说| 中文字幕色久视频| 色婷婷av一区二区三区视频| 免费在线观看黄色视频的| 操出白浆在线播放| 国产欧美日韩综合在线一区二区| 免费高清在线观看视频在线观看| 精品福利永久在线观看| 大陆偷拍与自拍| 日韩,欧美,国产一区二区三区| 1024香蕉在线观看| 国产精品久久久人人做人人爽| 国产又色又爽无遮挡免| 久久狼人影院| 欧美97在线视频| 国产又爽黄色视频| 精品亚洲成国产av| 97在线人人人人妻| 新久久久久国产一级毛片| 久久人人97超碰香蕉20202| 国产精品蜜桃在线观看| 国产欧美日韩一区二区三区在线| 999精品在线视频| 在线观看www视频免费| 韩国高清视频一区二区三区| 国产精品国产三级国产专区5o| 国产老妇伦熟女老妇高清| 亚洲欧美精品综合一区二区三区| 亚洲国产精品成人久久小说| 久久精品亚洲熟妇少妇任你| 下体分泌物呈黄色| 97精品久久久久久久久久精品| 一二三四在线观看免费中文在| 一区二区三区乱码不卡18| 国产在视频线精品| www日本在线高清视频| 色精品久久人妻99蜜桃| 亚洲一码二码三码区别大吗| 亚洲精品日韩在线中文字幕| 中文欧美无线码| 亚洲色图 男人天堂 中文字幕| 国产精品一国产av| 日韩一卡2卡3卡4卡2021年| 国产成人欧美在线观看 | 最近手机中文字幕大全| 国产精品一区二区在线不卡| 国产成人欧美在线观看 | 国产av一区二区精品久久| 亚洲人成网站在线观看播放| 久久久欧美国产精品| 蜜桃在线观看..| av视频免费观看在线观看| 亚洲精华国产精华液的使用体验| 欧美黑人欧美精品刺激| 黄色视频不卡| 看十八女毛片水多多多| 美女脱内裤让男人舔精品视频| 色94色欧美一区二区| 男女之事视频高清在线观看 | 丝袜喷水一区| 亚洲情色 制服丝袜| 在现免费观看毛片| 久久久久国产一级毛片高清牌| 操美女的视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 777米奇影视久久| 国产精品 欧美亚洲| 成人18禁高潮啪啪吃奶动态图| 国产成人免费无遮挡视频| 五月天丁香电影| 久久久精品区二区三区| videos熟女内射| 一区二区日韩欧美中文字幕| 19禁男女啪啪无遮挡网站| 制服丝袜香蕉在线| 天天躁夜夜躁狠狠久久av| 如日韩欧美国产精品一区二区三区| 十分钟在线观看高清视频www| 美女扒开内裤让男人捅视频| 看免费av毛片| 国产精品熟女久久久久浪| 激情视频va一区二区三区| 蜜桃国产av成人99| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲,欧美精品.| 美女高潮到喷水免费观看| 天天操日日干夜夜撸| 美女主播在线视频| 一边亲一边摸免费视频| 亚洲精品久久午夜乱码| 日韩免费高清中文字幕av| 国产 精品1| 午夜精品国产一区二区电影| 精品一区二区三卡| 久久精品国产亚洲av涩爱| 狠狠婷婷综合久久久久久88av| 1024视频免费在线观看| 日本vs欧美在线观看视频| 男女高潮啪啪啪动态图| 日日爽夜夜爽网站| 午夜激情av网站| 乱人伦中国视频| 日韩一区二区视频免费看| 国产一区二区激情短视频 | 欧美精品亚洲一区二区| 美女视频免费永久观看网站| 精品一区二区三卡| 亚洲精品美女久久av网站| av网站免费在线观看视频| 又大又黄又爽视频免费| 日韩中文字幕视频在线看片| 久久青草综合色| 亚洲精品第二区| 久久热在线av| 欧美日韩视频高清一区二区三区二| 日韩免费高清中文字幕av| 国产精品久久久久久久久免| 日韩一区二区视频免费看| 人人妻人人澡人人爽人人夜夜| 老鸭窝网址在线观看| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 女性被躁到高潮视频| 亚洲成人一二三区av| 日日啪夜夜爽| tube8黄色片| 午夜激情久久久久久久| 桃花免费在线播放| 午夜久久久在线观看| 王馨瑶露胸无遮挡在线观看| 国产亚洲一区二区精品| 国产亚洲av片在线观看秒播厂| avwww免费| 大陆偷拍与自拍| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 久久久久精品国产欧美久久久 | 免费少妇av软件| 久久精品人人爽人人爽视色| a级毛片在线看网站| 久久久精品国产亚洲av高清涩受| 成人三级做爰电影| 日日爽夜夜爽网站| 老司机深夜福利视频在线观看 | 精品久久久久久电影网| 中文字幕最新亚洲高清| 9191精品国产免费久久| 下体分泌物呈黄色| 国产精品久久久久久精品古装| 国产在线一区二区三区精| 秋霞伦理黄片| 天天添夜夜摸| 国产精品成人在线| 国产成人免费观看mmmm| 美女扒开内裤让男人捅视频| 精品少妇久久久久久888优播| 一区二区三区乱码不卡18| 国产精品国产三级国产专区5o| 日本91视频免费播放| 久久亚洲国产成人精品v| 国产在线免费精品| 高清黄色对白视频在线免费看| 国产 一区精品| 精品午夜福利在线看| 精品国产国语对白av| 欧美日韩亚洲高清精品| 新久久久久国产一级毛片| 一本一本久久a久久精品综合妖精| 一级毛片黄色毛片免费观看视频| 国产伦人伦偷精品视频| 亚洲欧美一区二区三区黑人| 亚洲熟女精品中文字幕| 国产激情久久老熟女| 国产成人av激情在线播放| 超碰成人久久| 飞空精品影院首页| 天天影视国产精品| 啦啦啦 在线观看视频| 色精品久久人妻99蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 欧美另类一区| 女人爽到高潮嗷嗷叫在线视频| 国产精品99久久99久久久不卡 | 亚洲国产欧美一区二区综合| 少妇人妻 视频| 制服人妻中文乱码| 欧美日韩精品网址| 欧美在线黄色| 黑人巨大精品欧美一区二区蜜桃| 97在线人人人人妻| av在线老鸭窝| 1024香蕉在线观看| 精品福利永久在线观看| 狂野欧美激情性xxxx| 亚洲国产最新在线播放| 狠狠精品人妻久久久久久综合| 国产毛片在线视频| 国产成人精品久久二区二区91 | 好男人视频免费观看在线| 999精品在线视频| 国产一卡二卡三卡精品 | 国产成人啪精品午夜网站| 亚洲自偷自拍图片 自拍| 国产深夜福利视频在线观看| 午夜免费男女啪啪视频观看| 国产av精品麻豆| 99热国产这里只有精品6| 国产一区二区 视频在线| 亚洲精品成人av观看孕妇| 国产免费又黄又爽又色| 老汉色av国产亚洲站长工具| 交换朋友夫妻互换小说| 女人精品久久久久毛片| 人妻人人澡人人爽人人| 久久久久久久大尺度免费视频| 熟女av电影| 国产福利在线免费观看视频| 99精品久久久久人妻精品| 欧美久久黑人一区二区| 一级a爱视频在线免费观看| 国产乱人偷精品视频| 看免费成人av毛片| 只有这里有精品99| 国产日韩一区二区三区精品不卡| 国产精品99久久99久久久不卡 | 成人国产av品久久久| 久久 成人 亚洲| 久久精品久久久久久久性| 久热这里只有精品99| av网站免费在线观看视频| 一级片'在线观看视频| 超碰成人久久| 日本91视频免费播放| 香蕉丝袜av| 在线观看免费日韩欧美大片| 精品久久蜜臀av无| 亚洲四区av| avwww免费| 亚洲精品,欧美精品| 亚洲精品久久午夜乱码| 青春草国产在线视频| 高清视频免费观看一区二区| 最黄视频免费看| 欧美人与性动交α欧美软件| 久久精品aⅴ一区二区三区四区| 在线观看免费日韩欧美大片| 嫩草影院入口| 精品亚洲成国产av| 精品国产一区二区三区久久久樱花| 波多野结衣一区麻豆| 丰满乱子伦码专区| 水蜜桃什么品种好| 国产不卡av网站在线观看| 免费高清在线观看视频在线观看| 9色porny在线观看| 999久久久国产精品视频| 性高湖久久久久久久久免费观看| 亚洲美女视频黄频| 高清欧美精品videossex| 搡老岳熟女国产| 色网站视频免费| 日本vs欧美在线观看视频| 美女福利国产在线| 国产深夜福利视频在线观看| 黄色一级大片看看| 人人妻人人澡人人看| 又黄又粗又硬又大视频| 少妇被粗大的猛进出69影院| 欧美日韩亚洲综合一区二区三区_| 色吧在线观看| 色婷婷久久久亚洲欧美| 一个人免费看片子| 欧美精品一区二区大全| 黑人猛操日本美女一级片| 亚洲三区欧美一区| 亚洲第一青青草原| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 久久人人97超碰香蕉20202| 亚洲成色77777| 18禁裸乳无遮挡动漫免费视频| 老司机亚洲免费影院| 国产又爽黄色视频| 国产成人精品久久久久久| 2018国产大陆天天弄谢| 建设人人有责人人尽责人人享有的| 80岁老熟妇乱子伦牲交| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| www日本在线高清视频| 午夜免费观看性视频| 亚洲av日韩在线播放| 国产亚洲欧美精品永久| 在线免费观看不下载黄p国产| 18禁国产床啪视频网站| 久久99精品国语久久久| 久久久久久久精品精品| 三上悠亚av全集在线观看| 亚洲av电影在线进入| 成人国语在线视频| 亚洲自偷自拍图片 自拍| 一级爰片在线观看| 天天躁夜夜躁狠狠久久av| 欧美另类一区| 国产精品免费视频内射| 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 午夜日本视频在线| 最新在线观看一区二区三区 | 中文字幕人妻熟女乱码| 热re99久久国产66热| 熟妇人妻不卡中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱久久久久久| 国产 精品1| 免费在线观看完整版高清| 波多野结衣一区麻豆| 色播在线永久视频| 九九爱精品视频在线观看| 90打野战视频偷拍视频| 最新在线观看一区二区三区 | 精品久久久久久电影网| 热re99久久国产66热| av在线老鸭窝| 亚洲一区二区三区欧美精品| xxxhd国产人妻xxx| 在线观看人妻少妇| 黑丝袜美女国产一区| 可以免费在线观看a视频的电影网站 | 亚洲精品成人av观看孕妇| 久久鲁丝午夜福利片| 一级毛片 在线播放| 久久人人爽人人片av| videosex国产| 午夜久久久在线观看| 超碰97精品在线观看| 欧美最新免费一区二区三区| 一边亲一边摸免费视频| 久久久精品区二区三区| 成人亚洲欧美一区二区av| 高清av免费在线| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 午夜激情久久久久久久| 婷婷成人精品国产| 色94色欧美一区二区| 国产精品人妻久久久影院| 亚洲欧美中文字幕日韩二区| 亚洲成人国产一区在线观看 | 国产精品秋霞免费鲁丝片| 亚洲一区中文字幕在线| 尾随美女入室| 女人精品久久久久毛片| 97精品久久久久久久久久精品| 19禁男女啪啪无遮挡网站| 亚洲一级一片aⅴ在线观看| 侵犯人妻中文字幕一二三四区| av有码第一页| 大码成人一级视频| 久热爱精品视频在线9| 国产在线免费精品| 亚洲av福利一区| 亚洲精品中文字幕在线视频| 亚洲成人av在线免费| 亚洲国产最新在线播放| 日韩电影二区| 91aial.com中文字幕在线观看| 我要看黄色一级片免费的| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 秋霞在线观看毛片| 女性生殖器流出的白浆| 国产激情久久老熟女| 亚洲成人一二三区av| 狠狠婷婷综合久久久久久88av| 久久久久久久久免费视频了| 国产日韩欧美亚洲二区| 1024视频免费在线观看| av卡一久久| 人妻 亚洲 视频| 人人妻人人澡人人看| 亚洲欧美日韩另类电影网站| 国产精品一区二区在线观看99| svipshipincom国产片| 午夜精品国产一区二区电影| 男女下面插进去视频免费观看| 黄片播放在线免费| 国产成人系列免费观看| 亚洲人成电影观看| 欧美人与性动交α欧美软件| 1024香蕉在线观看| 最近最新中文字幕免费大全7| 一级黄片播放器| 99国产综合亚洲精品| 国产亚洲最大av| 999久久久国产精品视频| 男的添女的下面高潮视频| 亚洲精品av麻豆狂野| 国产精品香港三级国产av潘金莲 | 久久久久精品久久久久真实原创| 丝袜脚勾引网站| 免费高清在线观看日韩| 青春草亚洲视频在线观看| 日本91视频免费播放| 黄片无遮挡物在线观看| 老司机影院成人| 久久久久精品人妻al黑| 啦啦啦在线观看免费高清www| 免费黄网站久久成人精品| 日韩大码丰满熟妇| av不卡在线播放| 精品人妻熟女毛片av久久网站| 多毛熟女@视频| 精品亚洲乱码少妇综合久久| a级毛片在线看网站| 99久久99久久久精品蜜桃| 高清视频免费观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产成人一精品久久久| 亚洲一区二区三区欧美精品| 一本久久精品| 捣出白浆h1v1| 啦啦啦在线观看免费高清www| 中文字幕亚洲精品专区| 久久免费观看电影| 国产男女内射视频| 日韩制服丝袜自拍偷拍| 日韩一本色道免费dvd| 最近手机中文字幕大全| 亚洲中文av在线| 制服人妻中文乱码| 欧美97在线视频| 日韩大片免费观看网站| 亚洲专区中文字幕在线 | 久久久久精品性色| 大香蕉久久网| 久久热在线av| 亚洲国产看品久久| 亚洲av成人精品一二三区| 国产精品久久久久久久久免| 一级片免费观看大全| 欧美 亚洲 国产 日韩一| 国产男女超爽视频在线观看| 精品一区二区三区av网在线观看 | 久久精品aⅴ一区二区三区四区| 成人国语在线视频| 韩国精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇人妻精品综合一区二区| 国产成人精品久久久久久| 美女福利国产在线| 亚洲在久久综合| 国产日韩一区二区三区精品不卡| 国产精品香港三级国产av潘金莲 | 成人三级做爰电影| av在线app专区| 日韩伦理黄色片| 天堂俺去俺来也www色官网| 亚洲精品视频女| 精品国产国语对白av| 久久鲁丝午夜福利片| 国产视频首页在线观看| 免费少妇av软件| 性少妇av在线| 九九爱精品视频在线观看| 免费观看a级毛片全部| 青春草国产在线视频| 男女边摸边吃奶| 水蜜桃什么品种好| 在线看a的网站| 中文字幕高清在线视频| 91精品三级在线观看| 国产亚洲最大av| 免费黄频网站在线观看国产| 日本wwww免费看| 晚上一个人看的免费电影| 大香蕉久久成人网| 久久人人爽人人片av| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 男的添女的下面高潮视频| 99久国产av精品国产电影| 男女免费视频国产| 一本大道久久a久久精品| 亚洲成人av在线免费| 国产1区2区3区精品| 精品少妇黑人巨大在线播放| 9色porny在线观看| 久久亚洲国产成人精品v| 久久久久精品国产欧美久久久 | 母亲3免费完整高清在线观看| 亚洲四区av| 国产成人精品久久久久久| 久久精品亚洲熟妇少妇任你| 成人亚洲欧美一区二区av| 99九九在线精品视频| 午夜福利在线免费观看网站| 在线观看免费高清a一片| 国产精品.久久久| 精品少妇内射三级| 人体艺术视频欧美日本| 欧美人与性动交α欧美软件| 免费观看av网站的网址| 色婷婷久久久亚洲欧美| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 伊人亚洲综合成人网| 一二三四在线观看免费中文在| 日韩一本色道免费dvd| 欧美日韩精品网址| 伊人亚洲综合成人网| 丝袜脚勾引网站| 国产精品 国内视频| 97人妻天天添夜夜摸| 在线看a的网站| 国产老妇伦熟女老妇高清| 免费女性裸体啪啪无遮挡网站| 高清欧美精品videossex| 成人国语在线视频| 自线自在国产av| 亚洲国产av新网站| 午夜精品国产一区二区电影| 国产极品粉嫩免费观看在线| av线在线观看网站| 波多野结衣一区麻豆| 女人被躁到高潮嗷嗷叫费观| 美女主播在线视频| 国产乱人偷精品视频| 色吧在线观看| 精品国产一区二区三区四区第35| 可以免费在线观看a视频的电影网站 | 老汉色av国产亚洲站长工具| av在线播放精品| 亚洲中文av在线| 亚洲国产欧美网| 国产精品av久久久久免费| 老司机亚洲免费影院| 中文天堂在线官网| 国产 精品1| 欧美亚洲 丝袜 人妻 在线| 悠悠久久av| 久久久精品国产亚洲av高清涩受| 三上悠亚av全集在线观看| 99热国产这里只有精品6| 国产深夜福利视频在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲欧美色中文字幕在线| 国产精品欧美亚洲77777| 国产午夜精品一二区理论片| 成人影院久久| 久久av网站| 少妇猛男粗大的猛烈进出视频| 亚洲av电影在线进入| 一本大道久久a久久精品| 国产成人精品无人区| 最近的中文字幕免费完整| 在线观看免费日韩欧美大片| 在线亚洲精品国产二区图片欧美| 国产成人免费观看mmmm| 久久久久久久国产电影| 成人国产av品久久久| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 黄色视频不卡| 国精品久久久久久国模美| 人体艺术视频欧美日本| 老司机影院成人| 亚洲一区二区三区欧美精品| 久久久久久免费高清国产稀缺| 99香蕉大伊视频| 中文字幕制服av| 日韩中文字幕视频在线看片| 亚洲精品久久成人aⅴ小说| 在线观看国产h片| 黄色一级大片看看| 91精品伊人久久大香线蕉| 黄频高清免费视频| 一本久久精品| 老鸭窝网址在线观看| 亚洲欧美日韩另类电影网站| 婷婷色综合大香蕉| 高清在线视频一区二区三区| 男女边吃奶边做爰视频| 亚洲成av片中文字幕在线观看| 大片免费播放器 马上看| 人体艺术视频欧美日本| 老司机在亚洲福利影院| 亚洲一区二区三区欧美精品| 晚上一个人看的免费电影| 午夜老司机福利片| 高清在线视频一区二区三区| 久久久久视频综合| 赤兔流量卡办理| 亚洲,一卡二卡三卡| 夫妻性生交免费视频一级片|