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

    Numerical investigation of velocity ratio effect in combined wall and offset jet flows *

    2019-01-05 08:08:50NidhalHnaienSalwaMarzoukHabibBenAissiaJacquesJay

    Nidhal Hnaien, Salwa Marzouk , Habib Ben Aissia, Jacques Jay

    1. Unité de Métrologie et des systèmes énergétiques, Ecole National d'ingénieurs de Monastir, Tunisie

    2. Centre de Thermique de Lyon (CETHIL–UMR CNRS 5008), INSA de Lyon, France

    Abstract: In the present work, we propose to numerically study a combined turbulent wall and offset jet flow (this combination will be denoted WOJ). Several turbulence models were tested in this study: the standard -ωk, SST -ωk, standard -kε, RNG, and realizable -kε model. A parametric study was performed to determine the effect of offset ratio H and the velocity ratio r on the longitudinal and transverse positions of the merge point (MP), the combined point (CP), the upper vortex center (UVC) and the lower vortex center (LVC). Correlations that predict the position of these characteristic points of the WOJ flow as a function of H and r have been provided. Results show that any increase in the velocity ratio causes a displacement of the MP, CP, UVC and LVC further upstream along the longitudinal direction. Concerning the transverse positions, the increase of velocity ratio results in a deviation of the merge point (MP) and the lower vortex center (LVC) toward the strong jet (LWJ) whereas the transverse position of combined point(CP) and the upper vortex center (UVC) is almost independent of the velocity ratio.

    Key words: Finite volume, combined point, merge point, offset jet, velocity ratio

    Introduction

    The multi-jet flows play an important role in various industrial applications such as boiler, burners,fuel injection systems, vertical takeoff and landing of air planes, ejectors, heating and air conditioning systems ...etc. The combined wall and offset jet flow(wall offset jet (WOJ)) (Fig. 1) is frequently encountered in many industrial processes especially in the heat exchanger and the wastewater evacuation process which makes interesting the study of the WOJ flow dynamic characteristics. The first reported work on this kind of flow was that of Wang and Tan[1], who have experimentally studied the WOJ flow characteristics using particle image velocimetry (PIV). They considered water as the working fluid. The Reynolds number and the offset ratio were respectively 10 000 and 2. The different flow characteristics such as the average velocity, the shear stress and the jet halfwidth at different shear layers are also measured.Results show that the zone in the vicinity of nozzleplate is characterized by a periodic large-scale Karman-like vortex shedding. This phenomenon results in a periodic interaction between the wall jet(LWJ) and the offset jet (UOJ). Vishnuvardhanarao and Das[2]numerically studied the heat transfer mechanism in the WOJ flow using the standard -kε turbulence model. The flow is assumed steady and incompressible. The Reynolds numbers considered were taken between 10 000 and 40 000 and the Prandtl number was =0 7Pr .. These authors considered two different thermal boundary conditions:constant wall heat flux and constant wall temperature.Vishnuvardhanarao and Das[2]have found that far downstream of the nozzle plate, the local Nusselt number Nu increases with increasing the mass flow ejected by elevating either the wall jet (LWJ) or the offset jet (UOJ) initial velocity. They also showed that the local Nusselt number and consequently the heat transfer exchanged between the wall and the whole flow is higher when the wall is under constant heat flux boundary condition. Vishnuvardhanarao and Das[2]also noticed that for an UOJ initial velocity of 0.25 m/s and that of LWJ equal to 1 m/s which corresponds to a velocity ratio of 4, the average Nusselt number is maximum and it is higher by 5%compared with the reference value.

    The numerical study of the turbulent WOJ flow is also achieved by Kumar and Das[3]using the high Reynolds number standard -kε turbulence model.The Reynolds number, the jet spacing and the turbulence intensity adopted were respectively 20 000,9 and 5%. The Kumar and Das[3]numerical code was validated on the experimental results of Pelfrey and Liburdy[4]and Vishnuvardhanarao and Das[5]for a single offset jet (SOJ) flow. The Kumar and Das[3]numerical results were compared with various available data in the literature. A discussion was also performed on the pressure and the turbulent kinetic energy fields. The maximum velocity decay and the half width evolution are also presented. Li et al.[6]have numerically studied the WOJ flow using the large eddy simulation (LES) method. Governing equations describing this flow are discretized using the finite volume method with an orthogonal mesh of 6.17×106nodes. The longitudinal velocity and turbulence intensity evolution for different longitudinal positions are presented and their distribution is in good agreement with experimental results of Wang and Tan[1]. Li et al.[6]showed that the interaction between the LWJ and UOJ is greater in the zone near the nozzle plate (the converging zone) and that these two jets begin the merging process at the merge point.In the second study, Li et al.[7]have numerically investigated the interaction between the LWJ and the UOJ in a WOJ flow for different offset ratio H and velocity ratio r using several turbulence models.The comparison between different turbulence models shows that the realizable -kε model is the most effective for the prediction of this type of flow. Li et al.[7]have noticed that for a fixed offset ratio, the flow pattern change from a single offset jet to that of a single wall jet flow when increasing the velocity ratio.They have also showed that when increasing the velocity ratio while maintaining constant the UOJ equal to 1m/s, the WOJ flow maximum velocity decay rate increase.

    Mondal and Das[8]numerically studied using two-dimensional unsteady RANS equations the nozzle spacing effect on the zone in the vicinity of the nozzle plate (the converging zone) of a WOJ flow. The numerical simulation was performed for different nozzle spacing and nozzle width. These authors noted that when the non-dimensional nozzle spacing was taken between 0.7 and 2.1, the converging region is characterized by the presence of large scale Karmanlike vortices. However, outside this range, two steady vortices with opposite direction are formed in the converging zone. The experimental setup of Pelfrey and Liburdy[4]on a SOJ flow with a Reynolds number of 15 000 and an offset ratio of 7 was recently numerically simulated by Kumar[9]. The nozzle width was d =12. 5 mm and the nozzle outlet velocity was 17.4 m/s. After validating the geometric configuration of Pelfrey and Liburdy[4], Kumar[9]added in parallel to the offset jet a second wall jet to form a WOJ flow.The purpose of his study was to investigate the effect of the wall jet (LWJ) addition on the whole WOJ flow structure. The numerical analysis performed by Kumar[9]compares the SOJ flow characteristic with those of WOJ flow for an offset ratio H ranging from 3 to 15 at an interval of 2. Correlations that predict the longitudinal and transverse positions of the merge point, combined points and the recirculation vortex centers are also provided with respect to a single parameter which is the offset ratio H.

    Through all the work that we have just mentioned, it appears that the majority of studies on the WOJ flow are concentrated mainly in the determination of the different zones of the flow field and on the offset ratio H effect on the MP, CP,UVC and LVC positions seen their important roles in the flow development, lunching and interruption of the merging process. In most WOJ flow studies,correlations that predict the positions of these various characteristic points are limited to a single parameter which is the offset ratio H. However, to our knowledge, there is no significant attempt to study the velocity ratio effect on the positions of these characteristic points for heights offset ratios H. This study aims to fill these gaps. It is important to note that this is the first time that correlations that predict the longitudinal and transverse positions of the merge point, the combined points and the vortices centers are provided as a function of two parameters together which are the offset ratio H and the velocity ratio r. Thus, the main objective of the present study is to investigate the simultaneous effect of the velocity ratio and the offset ratio on the WOJ flow characteristics. This study allowed us to provide correlations that predict the longitudinal and transverse positions of the merge point, the combined point and the vortex centers. Indeed, these correlations may be useful in the quick estimation of these different characteristics points that manage this type of flow in various industrial applications.

    1. Mathematical formulation

    1.1 . Geometric configuration

    The geometric configuration of Kumar[9]is used to study the influence of the velocity ratio in the WOJ flow (Fig. 1).

    Both nozzles are identical, rectangular, each ofaligned in the direction of the transverse flow (along y axis). The offset ratio H is set equal to 9, the turbulence intensity and the Reynolds number at each nozzle outlet were respectively I = 5% and Re = 15 000 which relate to an initial velocity value u0=17.4 m/s. The computed region dimensions are chosen so as not to affect the flow spreading. Several configurations were tested for finally adopt the following dimensions 100d and 50d respectively along the longitudinal and the transverse directions.Figure 2 shows the domain dimensions as well as the location of the Cartesian coordinates system adopted.

    Fig. 1 Schematic diagram of combined wall and offset jet flow(WOJ)

    Fig. 2 Grid size and boundary conditions

    1.2 Hypothesis

    The system of equations governing the flow is written in the Cartesian coordinate system whose origin σ is located on the horizontal wall axis (Fig.1). The following hypothesis will be also adopted:

    (1) The flow is two-dimensional and isothermal.

    (2) The work fluid is air assumed incompressible and the thermo-physical properties are constant.

    (3) The flow is supposed in a steady state.

    (4) The jet is emitted in the longitudinal direction.

    (5) The flow is assumed to be turbulent and fully developed.

    1.3 Governing equations

    The Reynolds-average Navier-Stocks (RANS)equations can be written in Cartesian tensor form as follows:

    Additional terms now appear which represent the effect of turbulence. These Reynolds stresses must be modeled in order to close Eq. (2). One common method employs the Boussinesq approximation to relate the Reynolds stress to the mean velocity gradient as follows

    The turbulent kinetic energy k and its specific dissipation rate “ω” noted w are obtained from the following equations:

    For k-ω model with high Reynolds number,α = 1. The dissipation of the turbulent kinetic energyand that of its specific dissipation rateare given as follows:

    1.4 Grid distribution and boundary conditions

    A non uniform grid is adopted along the longitudinal and the transverse directions. Indeed, a fine grid is used near the nozzle plate and a little looser further. As shown in Fig. 2, an uniform grid is used for the [A], [B] and [C] segments with a spacing s = 0.05. The [D] segment has a non uniform grid with spacing of 0.5 and an expansion ratio e =1.0450 which result in a nodes number of 280 on these segments ([A], [B], [C] and [D]). Note that the grid spacing s is dimensionless based on the nozzle width d and the expansion ratio e is a non-dimensional coefficient. For [E] segment, the grid is non uniform with s = 0.50 and e = 1.0185 which give the node number of 200. It is noted that all opposite segments have the same grid size. Consequently we obtained a total number of quadratic cells over the entire area as follows

    This particular choice of nodes number will be discussed later in the grid size sensitivity section. To complete the problem, besides the equations mentioned above, it is necessary to take into account the boundary conditions which are summarized in Table 1.

    Table 1 Boundary and emission conditions

    1.5 Numerical method

    In this work, the transport equation associated with the boundary and emission conditions are solved numerically by the finite volume method and by using the computational fluid dynamics software “Fluent”.The computed region is divided into finite number of sub-regions called “control volume”. The resolution method is to integrate on each control volume the transport equations such as the momentum conservation, mass conservation, the turbulence kinetic energy k and the specific dissipation rate of the turbulence kinetic energy w. These equations are discretized using the second order upwind. The coupling velocitypressure is based on the SIMPLEC algorithm. The disposal method of Gauss-Seidel associated with a relaxation technique is used to solve the resulting tri-diagonal matrix. A non uniform grid is adopted in the longitudinal and transverse directions, a fine grid is used near the nozzles and a little looser further. The convergence of the global solution is obtained when the normalized residuals fall below 10-4. We have verified that the increase in this accuracy had practically no influence on the results. The normalized residual is expressed as followsis the variable such as u, v, k, q or ω and,are the coefficients of the discretized equations andis the source term.

    2. Results and discussion

    2.1 Velocity field validation and grid sensitivity study

    Figure 3(a) shows the prediction of the longitudinal mean velocity transverse distribution on the plane X =7 using the standard k-omega turbulence model for H = 9 and R e = 15 000.

    Fig. 3 Transverse distributions of the longitudinal velocities at X=7

    Fig. 4 Longitudinal velocity U, static pressure P and turbulent kinetic energy K contours for =9H and =Re 15000

    As shown in Fig. 3(a), the numerical results are given using three different grid sizes such as(180/125), (280/200) and (330/333). These grids contain respectively 22 500, 56 000 and 109 800 quadratic cells. It seems clear that the velocity profile predicted by the grid (280/200) and (330/333) are almost identical while a remarkable difference was found between the prediction by the grids (180/125)and (280/200). It seems important to note that the grid(180/125) has a maximum discrepancy of 7.3%compared with the results of Kumar[9]while the grid(280/200) and (330/333) have almost the same maximum discrepancies approximately of 1.3%compared with Kumar[9]results. Thus, the grid(280/200) is sufficient to obtain a numerical solution independent of the grid size and validated by the numerical results of Kumar[9]. We presented in Fig.3(b), the transverse distribution of the longitudinal velocity on the plane =7X with (280/200) grid size using five different turbulence models: standard

    k -ω , SST k -ω , standard k -ε , RNG k -ε and realizable k -ε . The longitudinal velocity profile predicted by the various turbulence models are in agreement with the numerical results of Kumar[9],except the k -ε realizable turbulence models witch over-predicts the velocity profile for 2.5 ≤Y ≤ 5 .Whereas for Y > 5 , the standard k -ω model provides the best agreement with the Kumar[9] results.Thus, the standard k -ω turbulence model with(280/200) grid size will be adopted in this numerical study.

    2.2 Flow characteristics

    We present in Fig. 4 the longitudinal velocity U contours (Fig. 4(a)), the static pressure P contours(Fig. 4(b)) and the turbulent kinetic energy K contours (Fig. 4(c)) in the three WOJ flow zones namely the convergence zone forthe merging zone for (and the combined zonefor an offset ratio =H 9, a Reynolds number R e = 15000 and a velocity ratio =1r .

    The convergence zone is characterized by a negative longitudinal velocity U which highlights a recirculation flow (Fig. 4(a)) and by a negative static pressure P defining a depression zone (Fig. 4(b)).This last zone is characterized by an increase in the turbulent kinetic energy due to the strong interaction between the LWJ and the UOJ that attract and mix at the merge point defined at X =10.48. At this merge point, the static pressure P (Fig. 4(b)) and the turbulent kinetic energy K (Fig. 4(c)) have high values while the longitudinal velocity is zero (Fig.4(a)). Downstream of the merge point is located the merging zone. In this zone occurs the two jets (LWJ and UOJ) merging process. This region is characterized by a gradual increase of the longitudinal velocity (Fig. 4(a)), by a decrease in the static pressure(Fig. 4((b)) and by a decrease in the turbulent kinetic energy (Fig. 4(c)). Beyond the merging zone, we found the combined zone defined for X≥ 1 8.32, in this zone, the two jets combine and their resulting flow behaves like a single wall jet. In this latter zone(combined zone), the longitudinal velocity, the static pressure and the kinetic turbulent energy decrease continuously according to the longitudinal distances(Fig. 4).

    2.3 Offset ratio H effect (=1)r

    To highlight the offset ratio H effect on the MP, CP, UVC and LVC positions, we present in Fig.5 the longitudinal velocity contours U =0, the transverse velocity contours V =0 and the dU / d Y =0 contours for different offset ratio H and for a velocity ratio =1r .

    Fig. 5 U = 0, V =0 and d U / d Y =0 contours for r = 1,Re = 15000 and for different offset ratios

    Fig. 6 Static pressure (1000)P contours for =1r and for two different offset ratios

    The =0U and =0V contours intersection defines three characteristic points of the WOJ flow in the recirculation zone namely the merge point (MP),the upper vortex center (UVC) and the lower vortex center (LVC)[9]while the combined point is identified0 contour extremity (Fig. 5). It is clear from this figure that the MP, CP, UVC and LVC longitudinal positions increase when increasing the offset ratio H. In addition, these characteristic points move transversely away from the wall (higher Y values) when the offset ratio H increases.

    This effect of the offset ratio H on the various characteristic points longitudinal positions can be explained from Fig. 6. This figure in which are plotted the static pressure contours for two different offset ratio H = 5 (Fig. 6(a)) and H = 15 (Fig. 6(b))shows that increasing H results in a decrease in the recirculation zone depression compared with the atmospheric pressure (Pminis -110×10-3and -64×10-3respectively for H = 5 and H =15). This results in lower deviation of the UOJ toward the LWJ and consequently towards the horizontal wall for H = 15 compared to H =5. Therefore, the various characteristic points (MP, CP, UVC and LVC) are less attracted toward the nozzle plate (defined at X =0)and less deflected toward the horizontal wall (defined at Y =0).

    Figure 7 shows the variation of the merge point(Fig. 7(a)) and combined point (Fig. 7(b)) longitudinal position as a function of the offset ratio H for a velocity ratio =1r . The experimental and numerical results for the WOJ flow[1,3,9]and for two parallel jets flow[10-14]are used for comparison. As already noticed in Fig. 5, the Fig. 7 also shows that increasing H results in an increase inpandaccording to linear functions defined for =1r and for 5≤≤H 15 by the following equations:

    Equation (17) shown in Fig. 7(a) and which predictmpX is consistent with the results of Wang and Tan[1], Kumar and Das[3], Kumar[9], Nasr and Lai[10], Durve et al.[12], Hnaien et al.[13]and Wang et al.[14]while certain difference is noticed when comparing with the study of Anderson and Spall[11]forpredicted by Eq. (18) (Fig. 7(b)),our numerical results are in good agreement with those of Anderson and Spall[11], Kumar and Das[3],Kumar[9], Durve et al.[12], Hnaien et al.[13]while they differ from other studies[10,1,14]. The gap observed between the numerical simulation and experimental results are probably due to two reasons: (1) In addition to the experimental techniques variety for each study,the initial and emissions conditions adopted are not completely identical. (2) The geometric configurations of the WOJ flow and that of two parallel jets flow free parallel jet (FPJ) are not identical. Indeed, in the WOJ configuration, the LWJ is a wall jet that develops along a horizontal wall located at =0y ,while in the FPJ configuration, the wall jet becomes a free jet (absence of the horizontal wall).

    The evolution of the merge and combined points transverse positions with respect to the offset ratio H is shown in Fig. 8 for a velocity ratio =1r . It is noted thatandincrease linearly when increasing H following respectively these equations:

    Fig. 8 Transverse positions of merge and combined points with respect to H for =1r and for FPJ and WOJ flows

    Figure 9 shows the U = 0, V =0 andcontours for an offset ratio H = 9 and for different velocity ratios.

    The analysis of these figures shows that when increasing the velocity ratio, the MP, CP, UVC and LVC approach the nozzle plate defined at =0X while deviating toward the strong jet (LWJ) and therefore toward the wall defined at =0Y . This remark is in good agreement with our previous work[13], with that of Elbanna and Sabbaght[15], Durve et al.[12]in the case of two parallel jets flow (FPJ).These authors noted that when increasing the velocity ratio, the lower jet (with lower outlet velocity) is attracted toward the strong jet (with higher outlet velocity). The flow pattern presented by the stream functions is also showed in Fig. 9 for the same offset ratio =9H and for a velocity ratio =1r . These figures clearly show the two counter rotated vortices

    2.4 Velocity ratio effect

    Fig. 9 U = 0, V =0 and dU / dY =0 contours and streamline pattern for H = 9 and for different velocity ratios

    (UVC and LVC) for the different considered velocity ratios.

    The deviation phenomenon of the flow characteristic points (MP, CP, UVC and LVC) can be explained from the static pressure contours for =9H and for two different velocity ratios (see Fig. 10). The pressure reduction rate produced by the LWJ or the UOJ depends on the amount of the entrained fluid which in turn depends on the outlet velocity. For high velocity ratios, the jet with lower outlet velocity (UOJ)produces low pressure reduction rate compared with that produced by the jet with high outlet velocity(LWJ).

    Fig. 10 Static pressure (1000 )P contours for =9H and for two different velocity ratios

    Therefore, as shown in Fig. 10(b), for a high velocity ratio, the pressure at the UOJ peripheryis greater than that at the LWJ peripherywhereas when r=1 (Fig. 10(a)), the two jets have almost the same pressures at their peripheries (69- ×Therefore, when the velocity ratio increases, the lower jet (UOJ) with higher pressure at its periphery is deflected toward the strong jet (LWJ) with lower pressure at its periphery.

    Figure 9 also shows that increasing the velocity ratio causes a reduction of the merging zone length(this length is defined as the difference between the merge and combined points longitudinal positions).This is obvious since by increasing r, the UOJ outlet velocity decreases while keeping constant that of the LWJ, therefore the UOJ mass flow decreases and rapidly mixed with that of the LWJ. It also seems important to note that the UVC and LVC longitudinal positions decrease and approach the same location when elevating the velocity ratio. Indeed, for =2r ,the UVC and the LVC have almost the same longitudinal position

    2.4.1 Merge and combined points

    Figure 11 shows the merge point (Fig. 11(a)) and the combined point (Fig. 11(b)) longitudinal positions as a function of the offset ratio H for different velocity ratio r.

    Fig. 11 Longitudinal positions of the merge and combined points

    This figure shows that increasing H leads to an increase inandfollowing polynomial functions for different velocity ratios. For fixed value of H,anddecrease for higher value of velocity ratio (Fig. 11). Indeed, for >1r , the MP and CP move further upstream along the longitudinal direction indicating an acceleration in the merging process launch. This figure shows also that the slope ofandevolution as a function of H decreases as the velocity ratio increases. This means that the H effect on the MP and CP longitudinal positions is less noticed for high velocity ratios. This can be explained by the fact that the transverse velocity amplitude in >1r case is greater than that in the =1r case. On the other hand, when r is increased from 1 to 2,decreases by 23%, 32.1%and 41.5% respectively for =5H , =7H and(Fig. 11(a)), likewisedecreases by 49.8%, 50.5% and 59.1% (Fig. 11(b)). We can then conclude that the velocity ratio effect on the MP and CP longitudinal positions is more accentuated for higher offset ratios H. This can be explained by the fact that the MP and CP upstream displacement in the longitudinal direction when increasing the velocity ratio are damped by the shearing effect of the horizontal wall. Or from Fig. 5, we note that increasing the offset ratio H result that the MP and CP transversally move away from the horizontal wall which weakened the wall shearing effect on these points. Indeed, when the offset ratio H increases from 5 to 15, the MP and CP transversally deviate from the wall respectively by 41.6% and 64.6%. Thus for elevated values of H and r, the merge and combined points displacement becomes less damped by the horizontal wall shearing effect.

    Fig. 12 Scatter plots relative to longitudinal positions of the merge and combined points

    Knowing the merge and combined points positions is essential to have a complete idea about the lunching and interruption of the mixing process in the WOJ flow. These positions are strongly affected by the offset ratio H and the velocity ratio r.Therefore, it is desirable to determine an estimating method of the merge and combined points positions based on these important parameters (Hand r).

    In this regard, we proposed empirical correlations that predict the longitudinal positions of these points.Fig. 12 present the scatter plots relating tocurves as a function of Ln()Hfor an offset ratio H between 5 and 15 and a velocity ratio r between 1 and 2, the scatter plots may be represented by the following linear functions:

    Equations (21), (22) allow us to deduce correlations that attach the MP and CP longitudinal positions to H and r. These correlations are as follows:

    Figures 14(a), 14(b) respectively show the merge and combined point’s transverse evolution for 515H≤≤ and 12r≤≤. These figures reveal thatandincrease when increasing H for the various considered values of velocity ratios. It is clear from Fig. 14(a) that thecurve slope as a function of H slightly decreases when r increases, which means that the offset ratio H effect on theevolution is slightly less observed for high velocity ratios. For fix value of the offset ratio H, the increase in the velocity ratio results in adecrease,therefore, the merge point is laterally deflected towards the strong jet (LWJ) and subsequently towards the horizontal wall. When the velocity ratio increases from 1 to 2, the merge point undergoes a lateral deviation of 8.6% and 12.2% respectively for H =5 and H = 15. Therefore the effect of the velocity ratio inmpY is sparsely influenced by the offset ratio H. For an offset ratio =9H , when the velocity ratio r varied from 1 to 2, the merge point is deviated laterally by 13.1%. On the other hand, in the numerical study of Durve et al.[12]on two parallel jets flow (FPJ) and respecting the same Cartesian coordinates origin adopted in this present study, the merge point is deviated by 43% for the same velocity ratio r interval (from 1 to 2). Then we can conclude that the velocity ratio effect is more pronounced in FPJ flow compared to WOJ flow. This is obvious since the horizontal wall in the WOJ flow partially damp the lateral deviation of the merge point toward the strong jet (LWJ).

    Fig. 13 longitudinal positions of the merge and combined points compared to correlations (23) and (24)

    For 1.5r≤ , the combined point transverse positionincreases linearly as a function of H with an increase rate which is independent of the velocity ratio r (Fig. 14(b)). For fixed offset ratio H, when the velocity ratio increases, the combined point transverse position slightly decreases. Indeed, for r≤ 1 .5, the Ycpcurves for different r may be considerate similar and represented by a single linear equation

    For >1.5r , the combined point transverse positiondecreases as r increases for <11H ,which mean that the merge point slightly deviate toward the strong jet (LWJ). While for 11H= , the transverse position of the CP is the same for >1.5r ,while for >11H ,increases when increasing the velocity ratio. Thus, the combined point transversely deviates this time toward the lower jet (UOJ).

    Fig. 14 Transverse positions of the merge and combined points

    From Eq. (26), we can deduce the following correlation which attaches the merge point transverse position to H and r

    Fig. 15 Scatter plots relative to transverse position of the merge point and comparison with correlation (27)

    2.4.2 Vortices centers

    In WOJ flow, two vortices with opposite directions are observed in the recirculation zone: an upper vortex (UVC) in the side of the UOJ and another lower (LVC) in the side of the LWJ, these vortices(Fig. 9) are responsible for the UOJ deviation toward the LWJ and consequently toward the horizontal wall.

    In Fig. 16, we plot the curves which represent the UVC (Fig. 16(a)) and LVC (Fig. 16(b)) longitudinal positions as a function of H for different r.

    Fig. 16 Longitudinal positions of the upper and lower vortex centers

    The figure analysis shows that increasing the offset ratio H result in longitudinal displacement of the UVC and LVC further downstream. In addition,the H effect on UVC and LVC longitudinal positions is less observed for high velocity ratios. This results a reduction in the(Fig. 16(a)) and(Fig. 16(b)) increase rate for higher value of velocity ratio r. For =5H , we note that the UVC and LVC longitudinal positions are almost equal and independent of the velocity ratio. This observation is in good agreement with Kumar[9]who found thatandare almost the same for low offset ratio H. For H >5, the increase in velocity ratio causes longitudinal displacement of the UVC and LVC further upstream (anddecrease). The same figure shows that an increase in the velocity ratio r from 1 to 2 lead to a longitudinal displacement of the UVC further upstream by 40.6% and 49.9%respectively for H=9 and H=15. On the other hand,the LVC positions shift further upstream by 20.1%and 33.3% in the same interval of r. These percentages allows us to conclude that the velocity ratio effect onandeis more pronounced for high offset ratio H and thatis more influenced by the velocity ratio variation compared to. This can be explained by the fact that the UVC and LVC movement is damped by two type of resistance; one is applied by the horizontal wall on the lower vortex and the other applied by the surrounding fluid on the upper vortex. The resistance exerted by the wall on the lower vortex is greater than that exerted by the surrounding fluid on the upper vortex. Therefore, the UVC is more influenced by the velocity ratio compared to the LVC.

    Figures 17(a), 17(b) show respectively the UVC and LVC transverse positions with respect to H for different r.

    Fig. 17 Transverse positions of the upper and lower vortex centers

    It is clear from Fig. 17(a) for all velocity ratios considered that the UVC and LVC transverse positions increase linearly when increasing the offset ratio H. The same figure shows that the H effect on the UVC transverse location is independent of r,this is clear from the constant value of theuvcY curve slope regardless of the velocity ratio. On the other hand, for fixed offset ratio H, the UVC transverse position slightly decreases when increasing r.Thereafter, we can considerate that theuvcY curves as a function of H are the same for all velocity ratios considered and can be described by the flowing equation

    The Fig. 17(b) analysis shows an increase in the transverse position of the LVC when H increases for all velocity ratios. The decrease in thecurve slope for high value of r indicates a weaker H effect on the LVC transverse displacement for higher values of r. For fixed offset ratio H, increasing the velocity ratio results in a reduction of the LVC transverse position which reflect a LVC deviation toward the strong jet (LWJ) and subsequently toward the horizontal wall.

    3. Conclusions

    The present work has numerically investigated turbulent flow in which a wall and an offset jet flow(WOJ) combined. This study is based on the finite volume method using the computational fluid dynamics software “Fluent”. The objective of this study is to examine the effect of the velocity ratio r and the offset ratio H on the dynamic characteristics of the WOJ flow. The discussion focuses mainly on the -ωk model's validity to predict the characteristics points of WOJ flow such as the merge point(MP), the combined point (CP), the upper vortex center (UVC) and the lower vortex center (LVC). The present work can be used for a quick estimate of these points' position that enormously affect the WOJ flow in different industrial applications. The results derived from this study may be summarized as follows:

    (1) For a fixed velocity ratio, increasing the offset ratio H results in a longitudinal displacement of the MP, CP, UVC and LVC further downstream(,,andincrease). The H effect on the MP, CP, UVC and LVC longitudinal position is less observed for high velocity ratios.Correlation that predict the MP and CP longitudinal positions are respectively given by Eqs. (23), (24).

    (2) For a fixed velocity ratio, when increasing H, the MP, CP, UVC and LVC move further away from the wall in the transverse direction (andincrease). Correlation that predicts the MP transverse position is given by Eq. (27).

    (3) The H effect on the MP and LVC transverse positions is less observed when increasing the velocity ratio while the increase in the velocity ratio r has virtually no effect on theandevolution respectively for 12r≤≤ and 1 5r.≤which are given as a function of the offset ratio respectively by Eqs. (25), (28). For a velocity ratio beyond =1 5r ., the H effect onevolution increase for an offset ratio value above H = 11.

    (4) For a fixed offset ratio H, increasing the velocity ratio results in a longitudinal displacement of the MP, CP, UVC and LVC further upstream in the longitudinal direction (,,anddecrease. This decrease is more observed for higher offset ratio H.

    (5) For a fixed offset ratio H, increasing the velocity ratio produces a transverse deviation of the MP and LVC toward the strong jet (LWJ) and subsequently toward the horizontal wall. This deviation is more intense for higher offset ratio H.On the other hand, the UVC transverse position is almost independent of the velocity ratio (slightly decreases as r increases). For the CP, its transverse position is also independent of the velocity ratio only when 1 5r.≤ (a slight decrease is noted as r increases). Beyond =1 5r . and for an offset ratio more than =11H , increasing the velocity ratio causes a deviation of the CP this time towards the lower jet (UOJ).

    (6) The velocity ratio effect on the MP transverse position is more observed in the FPJ flow compared to the WOJ flow since the horizontal wall shearing effect in the WOJ flow partially damp the transverse displacement of the MP when changing the velocity ratio.

    (7) The UVC longitudinal displacement is more affected by the velocity ratio variation compared with that of LVC because of the higher resistance applied by the horizontal wall on the lower vortex compared with that applied by the surrounding fluid on the upper vortex.

    Acknowledgement

    The helpful comments and suggestions of the reviewers are gratefully acknowledged by the authors.

    中文字幕亚洲精品专区| 午夜福利影视在线免费观看| 国产片特级美女逼逼视频| 亚洲精品中文字幕在线视频| 亚洲,欧美精品.| 最近最新中文字幕大全免费视频 | 成人国产av品久久久| 黄片小视频在线播放| 在线看a的网站| 国产97色在线日韩免费| 99九九在线精品视频| 亚洲国产成人一精品久久久| 国产成人一区二区在线| 国产精品无大码| 搡老乐熟女国产| 一级毛片我不卡| 久久99热这里只频精品6学生| 国产片内射在线| 人妻人人澡人人爽人人| 男人操女人黄网站| 97精品久久久久久久久久精品| 人妻一区二区av| 亚洲成人国产一区在线观看 | 国产女主播在线喷水免费视频网站| 91aial.com中文字幕在线观看| 亚洲精品乱久久久久久| av电影中文网址| 欧美日韩亚洲高清精品| 国产免费又黄又爽又色| 国产成人精品无人区| 国产熟女午夜一区二区三区| 天堂8中文在线网| 丝袜脚勾引网站| 亚洲美女搞黄在线观看| 大香蕉久久网| 少妇被粗大猛烈的视频| 美女福利国产在线| 免费黄网站久久成人精品| 欧美人与善性xxx| 母亲3免费完整高清在线观看| 操出白浆在线播放| av国产精品久久久久影院| 少妇人妻精品综合一区二区| 亚洲人成网站在线观看播放| 国产精品国产三级专区第一集| 亚洲成国产人片在线观看| 午夜福利视频精品| www.自偷自拍.com| 看免费成人av毛片| 99精品久久久久人妻精品| 亚洲av电影在线观看一区二区三区| 女性被躁到高潮视频| 女人高潮潮喷娇喘18禁视频| 夫妻性生交免费视频一级片| 色婷婷av一区二区三区视频| 女人爽到高潮嗷嗷叫在线视频| 久久久久久人人人人人| 啦啦啦 在线观看视频| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产一区二区精华液| 王馨瑶露胸无遮挡在线观看| 午夜福利视频在线观看免费| 亚洲欧美一区二区三区黑人| 中文精品一卡2卡3卡4更新| 精品卡一卡二卡四卡免费| 欧美 日韩 精品 国产| 久久久欧美国产精品| 国产女主播在线喷水免费视频网站| 中国三级夫妇交换| 亚洲国产最新在线播放| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 久久久国产一区二区| 一区二区三区乱码不卡18| 亚洲美女搞黄在线观看| 在线天堂最新版资源| 亚洲国产av影院在线观看| 啦啦啦在线免费观看视频4| 国产亚洲av高清不卡| 精品一品国产午夜福利视频| 色精品久久人妻99蜜桃| 中国国产av一级| 狂野欧美激情性bbbbbb| 欧美老熟妇乱子伦牲交| 男人操女人黄网站| 日日摸夜夜添夜夜爱| 老司机靠b影院| 亚洲精品美女久久久久99蜜臀 | svipshipincom国产片| 亚洲四区av| 涩涩av久久男人的天堂| 国产成人精品久久久久久| 成人国语在线视频| 亚洲精品国产av成人精品| 日日啪夜夜爽| 午夜福利网站1000一区二区三区| svipshipincom国产片| 国产精品99久久99久久久不卡 | 午夜免费鲁丝| 如何舔出高潮| 亚洲人成77777在线视频| 人人妻人人澡人人看| 丁香六月欧美| 天天躁狠狠躁夜夜躁狠狠躁| 久久99一区二区三区| 久久久久视频综合| 午夜激情av网站| 两个人免费观看高清视频| 亚洲av欧美aⅴ国产| 亚洲av福利一区| 色婷婷久久久亚洲欧美| 免费看不卡的av| 麻豆av在线久日| 久久久久精品国产欧美久久久 | 国产精品 欧美亚洲| 国产 精品1| 精品一品国产午夜福利视频| 国产精品无大码| 色视频在线一区二区三区| 999精品在线视频| 看免费成人av毛片| 色吧在线观看| 国产老妇伦熟女老妇高清| 精品国产一区二区三区四区第35| 亚洲一区二区三区欧美精品| 成人漫画全彩无遮挡| 少妇 在线观看| 国产午夜精品一二区理论片| 观看av在线不卡| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 在线天堂中文资源库| 亚洲欧洲国产日韩| 亚洲四区av| 久久婷婷青草| 一区二区三区精品91| 9色porny在线观看| av线在线观看网站| 午夜福利网站1000一区二区三区| 国产精品一区二区在线不卡| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91 | 国产在线一区二区三区精| 午夜日本视频在线| 18在线观看网站| 国产精品熟女久久久久浪| 国产国语露脸激情在线看| 久久精品国产亚洲av涩爱| 亚洲av在线观看美女高潮| 2021少妇久久久久久久久久久| 亚洲国产精品成人久久小说| 曰老女人黄片| 777久久人妻少妇嫩草av网站| 国产精品久久久久久精品电影小说| 麻豆乱淫一区二区| 一边摸一边做爽爽视频免费| 丝瓜视频免费看黄片| 国产午夜精品一二区理论片| 街头女战士在线观看网站| 国产成人a∨麻豆精品| 精品亚洲成国产av| 欧美亚洲 丝袜 人妻 在线| 91精品国产国语对白视频| 国产日韩欧美在线精品| 免费观看a级毛片全部| 亚洲欧美一区二区三区久久| 欧美老熟妇乱子伦牲交| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 国产高清不卡午夜福利| 中文字幕最新亚洲高清| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 国产日韩一区二区三区精品不卡| 亚洲国产精品国产精品| 亚洲人成电影观看| 久久精品国产综合久久久| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到| av国产久精品久网站免费入址| 欧美在线黄色| 精品国产一区二区久久| 少妇人妻精品综合一区二区| 精品国产乱码久久久久久男人| 国产精品三级大全| 精品亚洲乱码少妇综合久久| 免费人妻精品一区二区三区视频| 最近手机中文字幕大全| 99re6热这里在线精品视频| 亚洲国产av影院在线观看| 老司机在亚洲福利影院| 国产成人精品在线电影| 极品人妻少妇av视频| 国产精品一区二区精品视频观看| 亚洲在久久综合| 免费不卡黄色视频| 国产无遮挡羞羞视频在线观看| 国产精品av久久久久免费| √禁漫天堂资源中文www| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| 精品少妇内射三级| 国产精品国产av在线观看| 日本一区二区免费在线视频| 日韩中文字幕欧美一区二区 | 久久毛片免费看一区二区三区| 亚洲精品日韩在线中文字幕| 大片免费播放器 马上看| 成年动漫av网址| 中文天堂在线官网| 看免费av毛片| 亚洲视频免费观看视频| 成年女人毛片免费观看观看9 | 极品少妇高潮喷水抽搐| 如日韩欧美国产精品一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久人人做人人爽| 日本一区二区免费在线视频| 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 国产成人啪精品午夜网站| 久久久久久久大尺度免费视频| 亚洲精品美女久久av网站| 啦啦啦在线观看免费高清www| 飞空精品影院首页| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 精品人妻熟女毛片av久久网站| 一本色道久久久久久精品综合| 国产日韩欧美视频二区| 成人毛片60女人毛片免费| 久久精品久久精品一区二区三区| 久热这里只有精品99| 最黄视频免费看| 老司机亚洲免费影院| 日本黄色日本黄色录像| 伦理电影免费视频| 精品视频人人做人人爽| 日韩大片免费观看网站| 飞空精品影院首页| 午夜日本视频在线| a级毛片黄视频| 成人亚洲精品一区在线观看| 国产亚洲最大av| 国产乱人偷精品视频| 亚洲av中文av极速乱| av电影中文网址| 99re6热这里在线精品视频| av网站在线播放免费| 老汉色∧v一级毛片| 欧美激情高清一区二区三区 | 丝袜美腿诱惑在线| 黄色怎么调成土黄色| 一区福利在线观看| 18禁裸乳无遮挡动漫免费视频| 啦啦啦啦在线视频资源| 伊人亚洲综合成人网| 亚洲成人手机| 国产欧美亚洲国产| 两个人看的免费小视频| 国产又色又爽无遮挡免| 国产男女内射视频| 波多野结衣av一区二区av| 亚洲美女黄色视频免费看| 老司机亚洲免费影院| 伦理电影大哥的女人| 在线观看国产h片| 国产免费又黄又爽又色| 18禁裸乳无遮挡动漫免费视频| 国产精品蜜桃在线观看| 免费人妻精品一区二区三区视频| 久久久久国产精品人妻一区二区| 最近最新中文字幕大全免费视频 | 久久久久久人妻| 久久久亚洲精品成人影院| 久久精品国产综合久久久| 亚洲男人天堂网一区| 日韩av免费高清视频| 伊人久久国产一区二区| av福利片在线| 欧美成人精品欧美一级黄| 国产精品 国内视频| 亚洲欧美清纯卡通| 五月开心婷婷网| 久久人人爽av亚洲精品天堂| 国产午夜精品一二区理论片| 国产成人啪精品午夜网站| 国产野战对白在线观看| 51午夜福利影视在线观看| 久久综合国产亚洲精品| 1024香蕉在线观看| 男女国产视频网站| 午夜福利乱码中文字幕| 国产亚洲欧美精品永久| 国产xxxxx性猛交| 精品久久久久久电影网| 超碰97精品在线观看| 麻豆av在线久日| 国产精品嫩草影院av在线观看| 色94色欧美一区二区| 9色porny在线观看| 亚洲人成网站在线观看播放| 一级毛片电影观看| 考比视频在线观看| 欧美少妇被猛烈插入视频| 一个人免费看片子| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| 少妇人妻精品综合一区二区| 亚洲精品视频女| 美女国产高潮福利片在线看| 亚洲第一区二区三区不卡| 一二三四中文在线观看免费高清| 国产在线视频一区二区| 精品第一国产精品| 无遮挡黄片免费观看| 日本wwww免费看| 亚洲av综合色区一区| 亚洲精品久久久久久婷婷小说| a 毛片基地| 亚洲七黄色美女视频| 啦啦啦视频在线资源免费观看| 亚洲第一青青草原| 国产精品女同一区二区软件| 一本一本久久a久久精品综合妖精| 啦啦啦在线观看免费高清www| 91成人精品电影| 欧美人与性动交α欧美精品济南到| 国产亚洲一区二区精品| 91国产中文字幕| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 女性被躁到高潮视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲综合精品二区| 精品免费久久久久久久清纯 | videos熟女内射| 操出白浆在线播放| www.自偷自拍.com| 99热网站在线观看| 涩涩av久久男人的天堂| av片东京热男人的天堂| 免费黄频网站在线观看国产| 18禁裸乳无遮挡动漫免费视频| 十八禁高潮呻吟视频| 国产精品秋霞免费鲁丝片| 黄色 视频免费看| 如日韩欧美国产精品一区二区三区| 高清欧美精品videossex| 国产毛片在线视频| 国产黄频视频在线观看| 欧美少妇被猛烈插入视频| 国产 一区精品| av女优亚洲男人天堂| 精品国产露脸久久av麻豆| 丰满迷人的少妇在线观看| 精品国产超薄肉色丝袜足j| 久久人妻熟女aⅴ| 久久青草综合色| 国产免费又黄又爽又色| 高清黄色对白视频在线免费看| 18在线观看网站| 美女主播在线视频| 一二三四在线观看免费中文在| 19禁男女啪啪无遮挡网站| 午夜激情久久久久久久| 韩国av在线不卡| h视频一区二区三区| 视频在线观看一区二区三区| xxxhd国产人妻xxx| 97精品久久久久久久久久精品| 欧美日本中文国产一区发布| 亚洲天堂av无毛| 欧美乱码精品一区二区三区| 黄色 视频免费看| 99国产综合亚洲精品| 亚洲欧洲日产国产| 妹子高潮喷水视频| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线观看免费高清www| 丰满饥渴人妻一区二区三| 久久av网站| 亚洲免费av在线视频| 老司机影院成人| 欧美日本中文国产一区发布| 18禁裸乳无遮挡动漫免费视频| 熟妇人妻不卡中文字幕| 久久免费观看电影| 久久精品久久精品一区二区三区| 亚洲在久久综合| 啦啦啦中文免费视频观看日本| 不卡视频在线观看欧美| 一二三四中文在线观看免费高清| 99精国产麻豆久久婷婷| 女性被躁到高潮视频| 欧美精品人与动牲交sv欧美| 青春草国产在线视频| 欧美精品一区二区免费开放| 免费黄网站久久成人精品| 菩萨蛮人人尽说江南好唐韦庄| 97在线人人人人妻| 爱豆传媒免费全集在线观看| 欧美黑人精品巨大| 亚洲情色 制服丝袜| 亚洲精品美女久久久久99蜜臀 | 无限看片的www在线观看| 肉色欧美久久久久久久蜜桃| av电影中文网址| 久久精品国产亚洲av高清一级| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 一二三四在线观看免费中文在| 超色免费av| 亚洲精品在线美女| 最黄视频免费看| 亚洲伊人色综图| 七月丁香在线播放| 免费少妇av软件| 欧美日韩亚洲国产一区二区在线观看 | 男女高潮啪啪啪动态图| 人人妻人人爽人人添夜夜欢视频| 校园人妻丝袜中文字幕| 夫妻性生交免费视频一级片| 亚洲精品日本国产第一区| 日韩av在线免费看完整版不卡| 韩国高清视频一区二区三区| 少妇人妻 视频| 电影成人av| 老司机亚洲免费影院| a级毛片在线看网站| 国产日韩欧美视频二区| av免费观看日本| 999久久久国产精品视频| 亚洲国产av新网站| 午夜影院在线不卡| 一级毛片电影观看| 成年人免费黄色播放视频| 精品国产国语对白av| 久久人妻熟女aⅴ| 国产黄色视频一区二区在线观看| 中文欧美无线码| 久久久久久人人人人人| 亚洲成国产人片在线观看| 日韩一区二区视频免费看| 中文字幕最新亚洲高清| 亚洲国产毛片av蜜桃av| 色婷婷av一区二区三区视频| 精品一区二区三卡| 中文字幕最新亚洲高清| 国产精品.久久久| 夜夜骑夜夜射夜夜干| av网站免费在线观看视频| 国产精品一二三区在线看| 操出白浆在线播放| 免费少妇av软件| 国产亚洲最大av| 成年人免费黄色播放视频| 亚洲欧美一区二区三区久久| 国产精品 国内视频| 超碰97精品在线观看| 亚洲国产av影院在线观看| 欧美在线黄色| 青青草视频在线视频观看| 90打野战视频偷拍视频| 亚洲欧美中文字幕日韩二区| 免费高清在线观看视频在线观看| www.自偷自拍.com| 国产男女内射视频| 日本黄色日本黄色录像| 免费在线观看视频国产中文字幕亚洲 | 亚洲少妇的诱惑av| 2018国产大陆天天弄谢| 中文字幕人妻丝袜一区二区 | 欧美精品一区二区大全| 少妇 在线观看| 中文字幕人妻熟女乱码| 无遮挡黄片免费观看| www.熟女人妻精品国产| 性高湖久久久久久久久免费观看| 欧美人与性动交α欧美软件| 亚洲人成电影观看| 亚洲精品一区蜜桃| 久久久精品94久久精品| 亚洲中文av在线| 国产精品熟女久久久久浪| 免费高清在线观看视频在线观看| 乱人伦中国视频| 久久久久国产一级毛片高清牌| 国产精品久久久人人做人人爽| av在线播放精品| 亚洲第一区二区三区不卡| 韩国av在线不卡| 午夜老司机福利片| 日韩电影二区| 午夜免费鲁丝| 中文字幕人妻丝袜制服| 国产亚洲一区二区精品| 女人精品久久久久毛片| 18禁国产床啪视频网站| 精品少妇内射三级| 久久免费观看电影| 久久人人97超碰香蕉20202| 中文字幕色久视频| 男女床上黄色一级片免费看| 中文字幕精品免费在线观看视频| 午夜日本视频在线| 老熟女久久久| 男女国产视频网站| 午夜福利在线免费观看网站| 久久久久国产精品人妻一区二区| 亚洲美女视频黄频| 考比视频在线观看| 午夜日韩欧美国产| 久久天躁狠狠躁夜夜2o2o | 国产在线免费精品| 国产精品久久久久久精品电影小说| 亚洲色图 男人天堂 中文字幕| 无遮挡黄片免费观看| 精品第一国产精品| 你懂的网址亚洲精品在线观看| 亚洲欧美激情在线| 免费不卡黄色视频| 亚洲,欧美精品.| 亚洲国产欧美网| 亚洲精品自拍成人| 国产在线一区二区三区精| 久久久精品94久久精品| 亚洲人成77777在线视频| 久久久久国产精品人妻一区二区| 日韩制服骚丝袜av| 考比视频在线观看| av国产久精品久网站免费入址| 亚洲精品自拍成人| 亚洲少妇的诱惑av| 成人漫画全彩无遮挡| 丝袜美腿诱惑在线| 另类精品久久| 亚洲成国产人片在线观看| 国产成人一区二区在线| 2021少妇久久久久久久久久久| 制服诱惑二区| 男人添女人高潮全过程视频| 丁香六月欧美| 一级,二级,三级黄色视频| av线在线观看网站| 成人18禁高潮啪啪吃奶动态图| 国产av码专区亚洲av| 麻豆精品久久久久久蜜桃| 国产成人91sexporn| 国产精品国产av在线观看| 99九九在线精品视频| 一区二区三区精品91| 天天躁狠狠躁夜夜躁狠狠躁| xxxhd国产人妻xxx| 欧美精品一区二区大全| www.自偷自拍.com| 制服人妻中文乱码| 人体艺术视频欧美日本| 另类精品久久| 久久久国产一区二区| 伦理电影大哥的女人| 97精品久久久久久久久久精品| 日韩一卡2卡3卡4卡2021年| a级片在线免费高清观看视频| 亚洲综合精品二区| 哪个播放器可以免费观看大片| 美女中出高潮动态图| 国产成人a∨麻豆精品| 日韩大码丰满熟妇| av在线老鸭窝| 校园人妻丝袜中文字幕| av网站免费在线观看视频| 乱人伦中国视频| 精品第一国产精品| 亚洲av日韩在线播放| 女人被躁到高潮嗷嗷叫费观| 久久久国产一区二区| 久久精品人人爽人人爽视色| 侵犯人妻中文字幕一二三四区| 国产97色在线日韩免费| 精品人妻熟女毛片av久久网站| 欧美日韩精品网址| 操美女的视频在线观看| 精品少妇黑人巨大在线播放| 色婷婷久久久亚洲欧美| 精品国产一区二区久久| 午夜福利,免费看| 亚洲一区中文字幕在线| 日韩制服骚丝袜av| 99国产精品免费福利视频| 五月天丁香电影| 日韩不卡一区二区三区视频在线| 在线观看免费视频网站a站| av国产久精品久网站免费入址| 欧美精品一区二区免费开放| 成年动漫av网址| 丰满少妇做爰视频| 一区在线观看完整版| 精品人妻在线不人妻| 亚洲国产欧美网| 亚洲色图 男人天堂 中文字幕| 秋霞伦理黄片| 国产免费一区二区三区四区乱码| 如日韩欧美国产精品一区二区三区| 亚洲精品久久成人aⅴ小说| 精品久久久精品久久久| 一区二区三区精品91| 国产淫语在线视频| 啦啦啦在线免费观看视频4| 亚洲精品久久久久久婷婷小说| 国产xxxxx性猛交| 色播在线永久视频| 亚洲熟女精品中文字幕| 赤兔流量卡办理|