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

    Closed corner divertor with B × ?B away from the divertor: a promising divertor scenario for tokamak power exhaust

    2023-11-16 05:37:58徐國盛,賈國章,丁根凡
    Plasma Science and Technology 2023年10期

    Abstract

    Keywords: divertor concept, closed divertor corner, E × B drifts, SOLPS-ITER simulation,EAST tokamak

    1.Introduction

    Controlling the heat flux onto the divertor target plates and the resulting material erosion is one of the most critical challenges that hinders the further development of fusion energy.Therefore, considerable effort has been made to develop innovative divertor solutions in the tokamak fusion community.There are generally two approaches for advanced divertor solutions [1-6].One approach is to use advanced poloidal magnetic configurations to widen the flux expansion and increase the connection length[1,2].The other approach is to develop optimized geometric structures to improve the divertor closure and trap neutrals and impurities to enhance dissipation and facilitate divertor detachment[3-5].Detached or partially detached divertor operations are desirable for a steady-state fusion reactor to avoid excessive heat deposition and material erosion at the divertor target plate.The dynamics of deuterium neutrals(molecules and atoms)and their spatial distribution in the divertor chamber of tokamaks are critical to achieve divertor detachment.A high divertor neutral pressure,which can be generally realized by increasing the divertor closure,is of great importance for particle and power exhausts[5].Several innovative divertors with closed geometric structures have been proposed,including the V-shape divertor for SlimCS [3], the SAS divertor developed in DIII-D [4, 5]and the recently installed right-angle closed divertor in EAST[6].The closed-divertor strategy aims at achieving divertor regimes with high dissipation and stable detachment over the entire divertor target at low levels of fueling, low upstream densities and impurity influx to meet the requirement of compatibility with high-performance steady-state core energy confinement scenarios in fusion reactors.

    Figure 1.Schematic diagram of drifts in SOL and divertor regions for(a)B×?B towards the lower divertor,i.e.favorable Bt direction,and(b)B×?B away from the lower divertor,i.e.unfavorable Bt direction,with green arrows for Eθ×B drifts,blue arrows for Er×B drifts,purple arrows for Pfirsch-Schlüter flows and red arrows for parallel transport from the upstream SOL.The Pfirsch-Schlüter flow directs to lower divertor for upwards B × ?B, and changes its direction for downwards B × ?B.The Pfirsch-Schlüter flow is always in the same direction of the plasma current, thus it reverses poloidally with the change of Bt direction rather than Ip direction.

    Influences of vertical and horizontal divertor target plates on the recycling neutrals were investigated earlier in ASDEX and JET [7] where corner pumps were applied.A general feature of vertical divertors is the accessibility to the partial detachment state: the plasma detaches mainly near the strike point,whereas in the far scrape-off layer(SOL)the plasma is still attached with rather high electron temperature Tet.For horizontal target plates with a pump duct at the SOL side,however,heat fluxes near the strike point are typically higher than those for vertical targets, since the recycled neutrals cannot be concentrated effectively to achieve cooling plasma conditions.

    Upgrading of the previous lower graphite divertor to tungsten with a steady-state heat flux handling capacity of~10 MW m-2is a critical step to transition to a tokamak with full-metal plasma-facing components for EAST.One of the primary principles for designing the new tungsten divertor concept on EAST was devoted to develop and validate divertor schemes that have high dissipations and are fully detached at low upstream densities with low levels of impurity seeding, which are compatible with advanced longpulse H-mode core plasma scenarios [6].A new tungsten divertor was developed and installed recently in EAST,which features a right-angled divertor corner on the low field side that joins the horizontal and vertical divertor plates to enhance the retention of recycled neutrals, as shown in figure 1.The new EAST divertor has a more closed structure than the JET divertor where the pumping entrance is located in between the vertical and inclined slope target plates [8].Directions of plasma drifts are illustrated in figure 1 for cases with B×?B directed away from or towards the new lower divertor.The original idea of designing the new divertor is the corner effect,i.e.as the strike point locates at the horizontal plate, the recycling particles are reflected largely to the SOL, and then restricted by the vertical plate and thus trapped near the corner,so that Tein the SOL can be reduced significantly due to the build-up of high neutral pressure in the corner area.The principal contribution of this new manuscript is to study the differences between the favorable and unfavorable Btdirections and the impact of different strike point locations through SOLPS-ITER simulations, while only cases with a favorable Btdirection were investigated in [6].

    The rest of this paper is organized as follows.In section 2, the principle of the closed corner divertor is introduced.In section 3, preliminary results from EAST discharges demonstrating the corner effect are presented.In section 4, the simulation results carried out with SOLPSITER simulations are shown.Summaries and discussions are presented in section 5.

    2.Principle of the closed corner divertor

    The ‘corner effect’ facilitates the local buildup of high neutral pressure near the corner, and thus strongly dissipative divertor regimes can be acquired so as to maintain stable detachment along the entire divertor target at relatively lower separatrix densities and impurity seeding rates,as required for compatibility with advanced core confinements.Specifically, when the outer strike point is located at the horizontal target, particularly near the corner, a large fraction of recycled neutrals is reflected towards the SOL.

    Figure 2.(a) Divertor Langmuir probes implemented in the lower diverter target at the low field side.The black, red and blue curves demonstrate the magnetic separatrices, magnetic reconstructed with EFIT, at three time slices (3.98 s, 5.90 s and 7.34 s), respectively.The blue points indicate the locations of divertor Langmuir probes.(b)Electron temperature and(c)ion saturation current measured with divertor Langmuir probes during the scanning of strike point from the horizontal towards the vertical target plate.

    The vertical target plate will act as a baffle.These neutrals then become ionized and move along field lines back to the horizontal plate or are reflected by the vertical plate,increasing the effective dwelling time of neutrals and causing them to be trapped in the corner area.As a result, this trapping of recycling neutrals dramatically enhances divertor closure and reduces the number of recycled particles that can penetrate into the pedestal from the lower divertor.The leakage of neutrals into the PFR will be less compared to a vertical divertor geometry where the recycled neutrals will fly directly towards the PFR.In contrast, as the strike point locates on vertical target, the plate directs recycled neutrals towards the strike point and the private flux region (PFR)where the electron temperature and density are much lower than the SOL,so that they cannot be fully re-ionized and can fuel the pedestal from the PFR through the X-point area.Furthermore, the recycled neutral particles from the outer divertor target area are significantly more than that with the strike point on the horizontal target as there are more particles in the far SOL hitting the vertical plate.These neutral particles readily move across the separatrix, increasing the pedestal fueling [9].Therefore, the operation with strike point being on the horizontal target exhibits a much stronger capability of confining the recycling particles inside the divertor and reducing the pedestal fueling from the new divertor.The new divertor based on closed divertor concept in EAST provides a test bed for achieving steady-state detachment operations with high heating powers.The ‘corner effect’ is confirmed in preliminary experimental results and reasonable agreement with SOLPS-ITER simulations including drift effects is obtained.

    3.Experimental measurements

    The recently installed lower tungsten divertor in EAST was operated successfully during the recent experimental campaign (May-July, 2021) including discharges with injected heating power as high as 7.5 MW.A series of dedicated divertor shots have been carried out, such as long-pulse detachment feedback control discharges and high-power discharges to validate the design idea of this tungsten divertor and investigate the divertor plasma physics involved.

    For instance, during the L-mode discharge #98338 [6],the strike point position of outer divertor was scanned slowly to the vertical target from the horizontal target plate as indicated in figure 2(a),where three typical positions of the outer strike point at three different time slices are shown.The corresponding equilibria are calculated by magnetic reconstruction with EFIT.At t = 3.98, 5.90 and 7.34 s, the outer strike point was situated far from and close to the corner at the horizontal target plate,and close to the bottom of the vertical target plate, respectively.The main operation parameters in this discharge are Ip≈ 500.0 kA, Bt≈ 2.5 T, nel≈2.1×1019m-3, WMHD≈115.0 kJ, and q95≈6.3 under double null (DN) magnetic configuration with upwards ion B×?B drift,where nel,WMHDand q95denote line-averaged plasma density, stored energy and edge safety factor,respectively.The corresponding triangularity reduces gradually from 0.63 to 0.5 as the outer strike point scans from horizontal to vertical target plates.The heating power is 2.3 MW in total, which includes 0.5 MW LHW at 2.45 GHz,1.4 MW LHW at 4.6 GHz and 0.4 MW ECW.No extrinsic impurity was injected in this shot.

    Figure 3.Computational meshes in modeling(quadrangles in sky blue for plasmas,extended triangles in lavender for neutrals)for magnetic configurations at(a)3.98 s,(b)5.90 s,and(c)7.34 s of the EAST shot#98338.D+flows are estimated at interfaces numbered with S1-S5.

    The divertor closure can be influenced by the distance between the corner and the strike point.Figure 2(b)illustrates Teevolution at the divertor target measured with divertor Langmuir probes 10-15.The peak Teis situated in the SOL region near the strike point.The maximum Tedecreases as the strike point approaches the corner along the horizontal target plate.The lowest maximum Teis obtained as the strike point locates in the vicinity of the corner on the horizontal plate(measured by probe 13).The highest maximum Teis measured as the strike point locates at the vertical plate.Therefore, the corner effect is clearly evidenced by this experimental result.

    In the following section, SOLPS-ITER simulations including drift effects are performed using the three different divertor plasma configurations corresponding to t = 3.98,5.90 and 7.34 s.Both upward and downward B × ?B drift directions are modeled to gain more insight into the geometric effects and drift effects.

    4.Modeling results by using SOLPS-ITER

    4.1.Basic setups in simulations

    The meshes used in SOLPS-ITER simulation corresponding to t=3.98,5.90 and 7.34 s are illustrated in figures 3(a)-(c)and are termed as horizontal, near-corner and vertical cases,respectively.In the modeling, only deuterium speciesis included.Both upward and downward ion B × ?B drifts are considered for each configuration, and thus 6 cases in total are modelled.Fixed ionic deuterium density at core-edge interface (CEI) is applied and its value is adjusted to acquire similar electron density at the upstream separatrix,neu≈1.0×1019m-3,for upward and downward B × ?B directions to make reasonable comparisons with respect to the divertor plasma states.D2puffing is switched off and at CEI boundary the D+flux is self-adjusted to obtain the preset density.The radiated power loss by intrinsic impurities is about 0.8 MW in the core region for this discharge,and thus the total power that enters CEI is set as 1.5 MW.The power is partitioned equally between ions and electrons.

    Cross-field drifts,especially E×B drift,would generally influence the divertor asymmetry through generating significant convective particle flux [10].However, it is worth noting that incorporating full drifts and currents into SOLPSITER modeling requires the use of a smaller time step and can be challenging, especially for cases with ion B × ?B drift pointed away from the X-point, to converge due to the associated numerical instabilities.Here,drifts and currents are fully activated except for the fraction accounting for the E×B drift which is set as 0.6.Recycling coefficients at first wall and divertor target plates are assumed as 1.0.The ratios for released thermal D2molecules and reflected fast D atoms due to the wall recycling are calculated with the plasmasurface interaction model included in the SOLPS-ITER package [11].Leakage boundary conditions are applied at PFR and SOL boundaries.Sheath conditions for ion and electron temperatures,as well as for ion densities are set at the target plate boundary.The radial heat conductivity and particle diffusivity are set as constant,χi=χe=1.2 m2s-1,Dn=0.5 m2s-1.The resultant near-SOL λqis 6.0-8.2 mm with these transport coefficients.The experimental λqvalue varies in the range of 3-8 mm when scanning the strike point positions,which is slightly smaller than the simulated λq.Due to incomplete activation of drifts in the simulations, we did not precisely match the simulated profiles with the measured profiles in the experiments.However, this small discrepancy does not affect the main physical conclusions of the work.

    4.2.Simulation results

    Figure 4.The simulated 2D patterns of the neutral pressure contributed by deuterium atoms and molecules in ((a), (d))horizontal, ((b), (e))near-corner and((c),(f))vertical cases.The subfigures in the first row correspond to simulations with upwards B×?B while the results with downwards B × ?B are shown in the second row.

    Table 1.D+flux that passes across the five interfaces indicated in figure 3.For S1,S3 and S5,the flux is driven by poloidal E×B drift.For S2 and S4, the flux is resulted from radial E×B drift.The unit is 1020 s-1.

    Convective flow due to classic drifts is dominated by E×B drift in the divertor region.Therefore, we firstly calculate the D+flux generated by E×B drift at five interfaces, S1-S5as shown in figure 3, to evaluate its contribution to particle exchange between inner and outer divertor regions.S1and S5are selected at the inner and outer divertor entrances,respectively.S2and S4are along the divertor legs.S3separates the inner PFR and outer PFR.According to SOLPSITER conventions, positive directions are poloidally clockwise towards the outer divertor from the inner divertor, and radially outwards to SOL edge from the core.The fluxes through the interfaces, which are contributed by E×B drift,are summarized in table 1.For example, in the Horizontal case with downwards B×?B,the poloidal D+flux entering the outer divertor through S5is 1.26×1021s-1and that leaving the inner divertor through S1is 1.11×1021s-1.The poloidal D+flux through S3, -1.83×1021s-1,indicates the transferring of ionic deuterium to the inner divertor from the outer divertor, acting as a particle sink at the outer divertor and a particle source at the inner divertor.Table 1 demonstrates that this situation applies also to the near-corner and vertical cases with downwards B×?B.As a result,the inner divertor becomes denser and colder while the outer divertor has a reduced density and increased temperature in configurations with downwards B × ?B.As the direction of magnetic field is reversed(upwards B×?B),the drift fluxes at the interfaces change sign which makes the inner divertor a particle source and the outer divertor a particle sink,facilitating outer divertor detachment.The contribution from the Pfirsch-Schlüter flow is relatively smaller compared to the E×B drift (not shown in the table).

    Figure 5.The simulated 2D patterns of the ionization source for D+in((a),(d))horizontal,((b),(e))near-corner and((c),(f))vertical cases.The subfigures in the first row correspond to simulations with upwards B×?B while the results with downwards B×?B are shown in the second row.

    In tokamaks, neutrals are produced mostly from plasma ions through recycling at plasma-facing components.These neutral particles re-ionize and charge-exchange with plasma ions,which disperse heat and momentum.The plasma state in divertor relates strongly to the neutral density and pressure near the target plate [12-14].The combination of the geometric effects of the new EAST divertor and the drift effects with B×?B pointed away from the divertor generates a promising divertor scenario to facilitate detachment and heat flux mitigation by making use of the recycling neutrals.Patterns of the neutral pressure distribution(including deuterium molecules and atoms) at the same upstream electron densityneu≈1.0×1019m-3for the six cases are demonstrated in figure 4.For horizontal target plate cases in figures 4(a) and(d), the divertor corner locates on the SOL side and a neutral cushion is formed in the area between the corner and the outer strike point.As a result,the electron temperature Teis reduced as will be shown later.The maximum neutral pressure is about 9.3 Pa for upwards B× ?B and 7.2 Pa for downwards B × ?B.When the strike point moves closer to the corner,the maximum neutral pressure increases to as high as 11.5 Pa as indicated in figure 4(b) with upwards B × ?B.By changing the Btdirection, although the closure of corner is weakened with downwards B×?B,the peak value still exceeds 10.0 Pa,as seen in figure 4(e),due to the geometric effects.As the outer strike point moves further onto the vertical divertor plate, the recycled neutrals are directed towards the strike point and the PFR where the corner is situated.As some of the particles are lost to PFR, the resulting neutral pressure as shown in figure 4(c)is maximized close to the strike point but with a maximum of only~5.0 Pa, much smaller than the Horizontal or the near-corner case.

    The distribution of ionization particle source for the 6 cases is illustrated in figure 5.When the strike point is located on the horizontal plate,recycled neutrals are reflected towards the SOL region, creating an ionization region in front of the outer target plate between the closed corner and the strike point.This ionization region provides access to plasma conditions with high density and low temperature, which can facilitate approaching detachment.As the strike point approaches to the corner, the ionization intensity increases,especially with upwards B×?B in figure 5(b).As the strike point is situated on the vertical target,the ionization intensity is much lower in the vicinity of the strike point compared with the strike point on the horizontal target plate.

    As the inner strike point moves towards the slant reflector plate from the vertical plate, influences of the corner on neutral concentrations were also indicated in the simulations of ITER divertor [15] leading to an enhanced detachment degree and mitigations of the heat load.

    The primary philosophy of the conventional vertical divertor target is to re-ionize the recycled neutrals in the vicinity of the strike point to efficiently reduce the electron temperature Teand heat flux qtthere, such as the ITER-like divertor.However, the impinging particles are largely reflected as neutrals towards the PFR by the vertical target,and most of the particles cannot return to the SOL.As a result, it is difficult to concentrate the recycled neutrals towards the SOL which typically leads to relatively lower neutral density and pressure and higher Tein the far-SOL.While for the horizontal divertor especially with a closed corner, incoming particles are largely directed towards the SOL region where most of them are ionized or reflected back again by the vertical plate.Therefore, particles are trapped within the SOL region near the corner, facilitating achievement of high neutral pressure with enhanced dissipation and reduced Te.If the closed corner is replaced with a pumping duct, such as the outer divertor of JET [16], the divertor closure will be greatly relaxed.For example, simulation results during the conceptual design stage of the new EAST divertor demonstrate that high neutral pressures will almost disappear as the pumping duct is positioned at the lower end of the vertical plate [17].

    Figure 6.The simulated profiles ofalong the outer target plates for horizontal(black square),near-corner(red dot),and vertical(blue star)cases with B×?B directed upwards(left)and downwards(right).The position of the outer strike point is indicated by the dashed black line.

    In figure 6,the spatial profiles of ne,Te,j‖,outand qtalong the outer divertor target are plotted.For upwards B×?B,the maximum electron density in the near-corner case is about 1.5 times larger than the horizontal case and 3 times higher than the vertical case.Correspondingly, the near-corner case has the lowest electron temperature and the largest particle flux,while the vertical case has the highest electron temperature and the lowest particle flux.The simulated Teandj‖,outare in reasonable agreement with the measurements as shown in figure 2(b).The discrepancy between the simulations and experiments could be ascribed to the incomplete activation of drift effects and the uncertainties in radial transport coefficients.The heat flux in the vertical case is lower due to several factors, including a relatively low particle flux, a misalignment between the maximum particle flux and the maximum electron temperature,and a more grazing incidence angle of field lines in the vertical case.For downwards B×?B,the outer diverter acts as a particle source while the inner divertor as a particle sink, resulting in the reduction of plasma density and the increase of Tein the outer SOL as shown in figures 6(e)-(h).

    5.Summary and discussions

    For the divertor in future steady-state tokamak fusion reactors, not only exhausts of particle and power but also controlling of tungsten impurity source is critical to obtain stationary high-performance operations [18].To keep the tungsten sputtering erosion and tungsten concentration in the core plasma at an acceptable level, it is required to achieve a divertor target temperature of Tet< 5 eV and density net>1021m-3[19].Although the typical vertical-target divertor,for instance the ITER divertor, has the advantages of relatively smaller field-line incident angle, lower heat flux in the vicinity of strike point and better particle exhaust capability,however it has the disadvantage of too high target electron temperature Tetin the far SOL, typically much larger than 10 eV, which is unacceptable for tungsten impurity sputtering.The resulting erosion rate is not serious for today’s devices, but will be unacceptable for future fusion reactors[19, 20].The high Tetis mainly because the vertical target is unfavorable for the concentration of particles in the SOL.To solve this problem,we proposed the concept of closed corner divertor [6] and the divertor operational scenario which combines the ‘corner effect’ and the effect of E×B drifts with B × ?B away from the divertor.As indicated in the SOLPS-ITER simulation, as the outer strike point locates on the horizontal target, the target plate will reflect a significant number of recycled neutrals into the SOL where they are re-ionized and move along field lines back to and recycle again at the horizontal plate or are reflected again by the vertical plate.Thus, these particles are trapped within the corner area, dramatically enhancing the divertor closure and reducing the number of recycled particles that can penetrate into the pedestal from the divertor.As a result of the concentration of recycling particles in the SOL, the temperature at the outer target Tetis effectively reduced, as demonstrated by the EAST experiments in figure 2.The present work does not cover the conditions of detachment either in modeling or in experiment aspects.However, radiated power loss due to impurities can be used to achieve detachment.The required levels of impurity seeding for different cases (such as those with varying strike point locations and Btdirections) may differ and depend on the characteristics of the background divertor plasma.To elaborate further, the effectiveness of impurity seeding for achieving detachment may depend on the properties of the divertor plasma.In particular, a divertor plasma with high density/neutral pressure and low temperature (near-corner case with an unfavorable Btdirection) can facilitate high power exhaust,which in turn might necessitate lower levels of impurity seeding.Therefore, the near-corner case with an unfavorable Btdirection is likely to satisfy the compatibility requirement between a detached divertor and a high-performance core.

    Furthermore, we notice that this divertor scenario has a significantly larger particle flux along the entire horizontal plate compared with the conventional vertical target scenario,as indicated in figure 6(c), mainly due to the much higher density, indicating that this divertor scenario may have a stronger particle exhaust capability, which is extremely important for fusion reactors to remove the helium ash.

    The E×B drifts play critical roles in enhancing the outer divertor density and particle flux.With B×?B away from the divertor, the E×B drifts in the divertor move particles from the inner divertor to the outer divertor and the Pfirsch-Schlüter flows move particles from the upstream towards the outer divertor(the former dominates over the latter), as indicated in figure 1, which significantly increase the particle flux and density in the outer divertor.This effect has been demonstrated by the EAST experiments.A relatively lower Tetand higher jsare always observed for B×?B away from the divertor under the same parameter conditions.As a result, the detachment density threshold is much lower and the detachment is more easily accessible with B × ?B away from the divertor in comparison with B × ?B pointed to the divertor [21].

    The SAS divertor in DIII-D [4] is a highly closed divertor, and exhibits advantages in improving neutral trapping and facilitations of divertor plasma detachment.As compared with the SAS divertor in DIII-D,the EAST‘closed corner’divertor design is more conductive for actively watercooled W/Cu plasma-facing components, provides great engineering reliability, and does not require precise control over the strike point location, which can be challenging to achieve in EAST.The divertor scenario proposed here may potentially provide a promising divertor solution for steadystate operations of future tokamak fusion reactors with much higher power fluxes.However, the power threshold for L-H transition is usually higher with B × ?B away from the divertor, which is a major unfavorable factor.Partly for this reason,ITER is designed to operate with downwards B×?B direction.In addition, a relatively lower H98factor is usually obtained with B × ?B directed away from the divertor compared with B × ?B towards the divertor [22], but impurity exhaust capability is usually stronger,which is more important for steady-state high-performance operations in a tokamak with full-metal plasma-facing components.As the controlling of heat flux and impurity concentration are more serious challenges for future fusion reactors,they should have higher priorities than other factors,thus this divertor scenario may be a solution that deserves serious considerations.

    Acknowledgments

    This work was funded by the National Magnetic Confinement Fusion Program of China (Nos.2019YFE03030000,2019YFE03080500 and 2022YFE03060004), National Natural Science Foundation of China (No.U19A20113).

    ORCID iDs

    Lingyi MENG (孟令義)https://orcid.org/0000-0002-8633-5383

    亚洲成人中文字幕在线播放| 三级男女做爰猛烈吃奶摸视频| 国产高清三级在线| 精品一区二区三区视频在线| 偷拍熟女少妇极品色| 激情五月婷婷亚洲| 国产精品人妻久久久影院| 日韩制服骚丝袜av| 国产精品av视频在线免费观看| 日韩伦理黄色片| 久久精品久久久久久噜噜老黄| 亚洲成人精品中文字幕电影| 欧美日韩精品成人综合77777| 高清午夜精品一区二区三区| 国产有黄有色有爽视频| 中文字幕久久专区| 夜夜爽夜夜爽视频| 国产欧美日韩精品一区二区| 成人欧美大片| 高清视频免费观看一区二区| 亚洲色图av天堂| 日本免费在线观看一区| 欧美最新免费一区二区三区| 97超碰精品成人国产| a级一级毛片免费在线观看| 久久久精品免费免费高清| 人人妻人人看人人澡| 日韩在线高清观看一区二区三区| 中文资源天堂在线| 色5月婷婷丁香| 欧美+日韩+精品| 国产成人freesex在线| 最近中文字幕2019免费版| 天堂中文最新版在线下载 | 精品人妻偷拍中文字幕| 夜夜看夜夜爽夜夜摸| 欧美xxⅹ黑人| 交换朋友夫妻互换小说| 只有这里有精品99| 哪个播放器可以免费观看大片| 国产真实伦视频高清在线观看| 婷婷色麻豆天堂久久| 男女啪啪激烈高潮av片| 91aial.com中文字幕在线观看| 少妇人妻一区二区三区视频| 久久国产乱子免费精品| 国产真实伦视频高清在线观看| 精品国产三级普通话版| 中国三级夫妇交换| 综合色av麻豆| 免费观看性生交大片5| 免费观看性生交大片5| av国产久精品久网站免费入址| 2021少妇久久久久久久久久久| 春色校园在线视频观看| 亚洲av免费高清在线观看| 日韩 亚洲 欧美在线| 韩国高清视频一区二区三区| 亚洲成人一二三区av| 日本三级黄在线观看| 老女人水多毛片| 国产色婷婷99| 免费黄频网站在线观看国产| 欧美3d第一页| 婷婷色综合www| 少妇人妻精品综合一区二区| 啦啦啦啦在线视频资源| 水蜜桃什么品种好| 51国产日韩欧美| 美女视频免费永久观看网站| 久久久久久久大尺度免费视频| 最近最新中文字幕大全电影3| 日韩亚洲欧美综合| 少妇 在线观看| 欧美成人午夜免费资源| 亚洲av成人精品一区久久| 国产伦精品一区二区三区四那| 国产探花极品一区二区| 免费av不卡在线播放| 免费观看av网站的网址| 国产永久视频网站| 日韩中字成人| xxx大片免费视频| 成人漫画全彩无遮挡| 免费观看无遮挡的男女| 国产亚洲av嫩草精品影院| 在线免费十八禁| 99热这里只有是精品在线观看| 99热这里只有精品一区| 一级爰片在线观看| 国产又色又爽无遮挡免| 人人妻人人澡人人爽人人夜夜| 欧美日韩一区二区视频在线观看视频在线 | 国产人妻一区二区三区在| 啦啦啦啦在线视频资源| 老女人水多毛片| 国产毛片在线视频| 成人亚洲欧美一区二区av| 成人毛片a级毛片在线播放| 国产午夜精品久久久久久一区二区三区| 国产成人免费观看mmmm| 18禁裸乳无遮挡动漫免费视频 | 国产成人a区在线观看| 欧美高清成人免费视频www| 一个人看视频在线观看www免费| 九九在线视频观看精品| 亚洲精品视频女| av国产免费在线观看| av免费在线看不卡| 舔av片在线| 亚洲婷婷狠狠爱综合网| 黄色视频在线播放观看不卡| 中文字幕制服av| 国产成人a区在线观看| 搞女人的毛片| 永久网站在线| 欧美日韩视频精品一区| 少妇人妻久久综合中文| 国产精品爽爽va在线观看网站| 少妇猛男粗大的猛烈进出视频 | 亚洲第一区二区三区不卡| 免费看a级黄色片| 美女国产视频在线观看| 精品一区在线观看国产| 亚洲国产精品成人综合色| 亚洲国产色片| 久久久久久久大尺度免费视频| 可以在线观看毛片的网站| 免费大片黄手机在线观看| 国产成人精品久久久久久| 精华霜和精华液先用哪个| 日韩av免费高清视频| 99九九线精品视频在线观看视频| 男插女下体视频免费在线播放| 最近的中文字幕免费完整| 精品亚洲乱码少妇综合久久| 18+在线观看网站| 国产高清不卡午夜福利| 一区二区三区精品91| 国产精品国产三级国产av玫瑰| 黄色怎么调成土黄色| 国产人妻一区二区三区在| 99久久人妻综合| 人人妻人人爽人人添夜夜欢视频 | 成人国产av品久久久| 国产免费视频播放在线视频| 国内精品宾馆在线| 在线免费观看不下载黄p国产| 中文字幕人妻熟人妻熟丝袜美| 高清视频免费观看一区二区| 人体艺术视频欧美日本| av国产免费在线观看| 国产成人福利小说| 一二三四中文在线观看免费高清| 国产黄色视频一区二区在线观看| 国产乱人偷精品视频| 男人添女人高潮全过程视频| 亚洲国产日韩一区二区| 亚洲精品aⅴ在线观看| 少妇人妻 视频| 欧美一区二区亚洲| 18禁在线播放成人免费| 亚洲国产最新在线播放| 久久久久网色| 尤物成人国产欧美一区二区三区| 尾随美女入室| 国产女主播在线喷水免费视频网站| 黄片wwwwww| 国产亚洲午夜精品一区二区久久 | 国产探花极品一区二区| 精品国产露脸久久av麻豆| 亚洲欧美日韩另类电影网站 | 一区二区三区免费毛片| 99热全是精品| av专区在线播放| 成年女人看的毛片在线观看| 国产精品秋霞免费鲁丝片| 成人漫画全彩无遮挡| 日本一二三区视频观看| 欧美丝袜亚洲另类| 亚洲av成人精品一区久久| 国产白丝娇喘喷水9色精品| 亚洲一区二区三区欧美精品 | 中文字幕人妻熟人妻熟丝袜美| 亚洲,欧美,日韩| 中文字幕av成人在线电影| 亚洲av免费在线观看| 天美传媒精品一区二区| 三级国产精品片| 亚洲三级黄色毛片| 亚洲最大成人中文| 亚洲激情五月婷婷啪啪| 97超视频在线观看视频| 视频中文字幕在线观看| 国产精品嫩草影院av在线观看| 高清毛片免费看| 国产欧美亚洲国产| 婷婷色麻豆天堂久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品久久午夜乱码| 熟女人妻精品中文字幕| 精品人妻一区二区三区麻豆| 欧美+日韩+精品| 2021天堂中文幕一二区在线观| 国产一区二区亚洲精品在线观看| 欧美日韩国产mv在线观看视频 | 老司机影院成人| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 亚洲精品aⅴ在线观看| 国产国拍精品亚洲av在线观看| 国产老妇伦熟女老妇高清| 亚洲av成人精品一区久久| 色网站视频免费| 在线观看一区二区三区| 成人国产麻豆网| 久久久国产一区二区| 精品久久久久久久久av| 久久综合国产亚洲精品| 狂野欧美白嫩少妇大欣赏| 青春草视频在线免费观看| 99视频精品全部免费 在线| 五月伊人婷婷丁香| 69人妻影院| 国产精品久久久久久久久免| 国产成人freesex在线| 国产精品国产三级专区第一集| 免费不卡的大黄色大毛片视频在线观看| 一级毛片我不卡| 男女那种视频在线观看| 99视频精品全部免费 在线| 久久99热这里只有精品18| 女人十人毛片免费观看3o分钟| 乱系列少妇在线播放| 可以在线观看毛片的网站| 五月开心婷婷网| 黄片无遮挡物在线观看| 久久久久久国产a免费观看| 日本欧美国产在线视频| 久久精品国产亚洲av天美| 亚洲av二区三区四区| 亚洲精品日韩av片在线观看| 欧美一区二区亚洲| 好男人在线观看高清免费视频| 亚洲人成网站高清观看| 三级经典国产精品| 在线观看一区二区三区| 日韩成人av中文字幕在线观看| 91精品国产九色| 内射极品少妇av片p| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频 | 国产 一区精品| 综合色丁香网| 久热久热在线精品观看| 在线 av 中文字幕| 国产伦精品一区二区三区四那| 在线观看av片永久免费下载| 国产黄a三级三级三级人| 亚洲精品第二区| 欧美高清成人免费视频www| 国产极品天堂在线| 久久ye,这里只有精品| 亚洲欧美日韩另类电影网站 | 熟女人妻精品中文字幕| 久热这里只有精品99| 天天躁夜夜躁狠狠久久av| 欧美少妇被猛烈插入视频| 免费黄色在线免费观看| 日韩电影二区| 三级国产精品欧美在线观看| 国产成人精品婷婷| 五月天丁香电影| 久热这里只有精品99| 只有这里有精品99| 69人妻影院| 欧美极品一区二区三区四区| 精品酒店卫生间| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 夜夜爽夜夜爽视频| 1000部很黄的大片| 中国三级夫妇交换| 毛片一级片免费看久久久久| 99热这里只有精品一区| 亚洲经典国产精华液单| 精品一区二区免费观看| 国产在视频线精品| 国产综合懂色| 国产av码专区亚洲av| 国产精品国产三级国产av玫瑰| 在线 av 中文字幕| 国产欧美日韩精品一区二区| 亚洲人成网站在线观看播放| 中文天堂在线官网| 亚洲精品国产成人久久av| 免费看av在线观看网站| 中国美白少妇内射xxxbb| 婷婷色综合大香蕉| 卡戴珊不雅视频在线播放| 成人亚洲精品av一区二区| h日本视频在线播放| 亚洲精品乱久久久久久| 国内精品宾馆在线| 国产欧美日韩精品一区二区| 久久99精品国语久久久| 69人妻影院| 97在线人人人人妻| 国产精品久久久久久久电影| 九色成人免费人妻av| 久久久久性生活片| 搞女人的毛片| 亚州av有码| 亚洲成人久久爱视频| 不卡视频在线观看欧美| 在线观看一区二区三区| 高清欧美精品videossex| 成年女人在线观看亚洲视频 | 日韩欧美精品免费久久| 老女人水多毛片| 婷婷色麻豆天堂久久| 免费观看性生交大片5| 涩涩av久久男人的天堂| 大话2 男鬼变身卡| 欧美日韩视频高清一区二区三区二| 18禁动态无遮挡网站| 国产高清不卡午夜福利| 亚洲国产欧美在线一区| 人妻少妇偷人精品九色| 内射极品少妇av片p| 日韩不卡一区二区三区视频在线| 一本一本综合久久| 亚洲精品一二三| 亚洲熟女精品中文字幕| 丰满少妇做爰视频| 精品久久久久久久人妻蜜臀av| 成人黄色视频免费在线看| 春色校园在线视频观看| 欧美97在线视频| 午夜福利在线在线| 91久久精品国产一区二区成人| 日日啪夜夜爽| 交换朋友夫妻互换小说| 观看美女的网站| 老师上课跳d突然被开到最大视频| 亚洲av日韩在线播放| 国产亚洲最大av| 极品教师在线视频| 国产老妇女一区| 久久久久久久午夜电影| 成年免费大片在线观看| 成人高潮视频无遮挡免费网站| 亚州av有码| 日韩一区二区三区影片| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| 一本一本综合久久| 国产精品人妻久久久久久| 亚洲第一区二区三区不卡| 九九久久精品国产亚洲av麻豆| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 菩萨蛮人人尽说江南好唐韦庄| 国产精品女同一区二区软件| 狂野欧美白嫩少妇大欣赏| 老女人水多毛片| 国内精品美女久久久久久| 国国产精品蜜臀av免费| 免费黄网站久久成人精品| 国产永久视频网站| 毛片一级片免费看久久久久| 久久这里有精品视频免费| 精品人妻视频免费看| 国产免费一区二区三区四区乱码| 成人欧美大片| 26uuu在线亚洲综合色| 少妇人妻精品综合一区二区| 夫妻性生交免费视频一级片| 久久ye,这里只有精品| 亚洲不卡免费看| 插逼视频在线观看| 激情五月婷婷亚洲| 久久精品国产鲁丝片午夜精品| 三级国产精品欧美在线观看| 欧美日韩视频高清一区二区三区二| 内地一区二区视频在线| 你懂的网址亚洲精品在线观看| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 亚洲欧美一区二区三区黑人 | 午夜福利视频1000在线观看| av在线亚洲专区| av播播在线观看一区| 欧美成人一区二区免费高清观看| 国产色婷婷99| 日韩中字成人| 国产老妇伦熟女老妇高清| 国产伦理片在线播放av一区| 亚洲av成人精品一二三区| 午夜爱爱视频在线播放| 国产毛片a区久久久久| 国产黄a三级三级三级人| 国产精品国产三级国产专区5o| 久久久久久久久久人人人人人人| 男女边吃奶边做爰视频| 丝瓜视频免费看黄片| 精品人妻偷拍中文字幕| 欧美日韩视频高清一区二区三区二| 国产日韩欧美在线精品| 亚洲欧美日韩无卡精品| 毛片女人毛片| 精品少妇黑人巨大在线播放| xxx大片免费视频| 国产欧美日韩一区二区三区在线 | 国产一区二区在线观看日韩| 国产亚洲av嫩草精品影院| 免费看不卡的av| 免费不卡的大黄色大毛片视频在线观看| 色视频www国产| 身体一侧抽搐| 欧美xxxx黑人xx丫x性爽| 久热这里只有精品99| 国内揄拍国产精品人妻在线| 在现免费观看毛片| 特级一级黄色大片| 久久久久久久久大av| 亚洲精品中文字幕在线视频 | 久久ye,这里只有精品| 最近最新中文字幕大全电影3| 免费观看a级毛片全部| 久久久久九九精品影院| 精品酒店卫生间| 在线天堂最新版资源| 国产黄色免费在线视频| 国产伦精品一区二区三区视频9| 亚洲美女视频黄频| 18禁动态无遮挡网站| 国产91av在线免费观看| av网站免费在线观看视频| 18禁在线无遮挡免费观看视频| 国产午夜精品一二区理论片| 七月丁香在线播放| 99热这里只有精品一区| av在线观看视频网站免费| 久久久色成人| 亚洲av二区三区四区| 99九九线精品视频在线观看视频| 真实男女啪啪啪动态图| 91aial.com中文字幕在线观看| 亚洲欧美中文字幕日韩二区| 亚洲无线观看免费| 又粗又硬又长又爽又黄的视频| 中文字幕免费在线视频6| 天天躁夜夜躁狠狠久久av| 日韩欧美精品免费久久| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 久久精品久久精品一区二区三区| a级一级毛片免费在线观看| 中国三级夫妇交换| 熟妇人妻不卡中文字幕| 只有这里有精品99| 丝袜脚勾引网站| 亚洲国产高清在线一区二区三| 少妇猛男粗大的猛烈进出视频 | 国产男女内射视频| 亚洲天堂国产精品一区在线| 日本色播在线视频| 国产亚洲一区二区精品| av播播在线观看一区| 少妇 在线观看| 欧美精品一区二区大全| av国产久精品久网站免费入址| 伊人久久精品亚洲午夜| 久久ye,这里只有精品| 午夜免费观看性视频| 啦啦啦在线观看免费高清www| av线在线观看网站| 久久精品国产自在天天线| 啦啦啦啦在线视频资源| 日韩av不卡免费在线播放| 熟女人妻精品中文字幕| 纵有疾风起免费观看全集完整版| 欧美成人一区二区免费高清观看| 国产精品一区二区三区四区免费观看| 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 国产一区二区亚洲精品在线观看| 国产免费又黄又爽又色| 国产精品久久久久久精品电影| 国产一级毛片在线| 国产午夜精品久久久久久一区二区三区| 午夜福利视频1000在线观看| 欧美区成人在线视频| 六月丁香七月| 女人久久www免费人成看片| 国产精品女同一区二区软件| 免费av观看视频| 一级片'在线观看视频| 国产色婷婷99| 免费av毛片视频| 日本三级黄在线观看| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 日日啪夜夜撸| 亚洲欧美一区二区三区国产| av在线蜜桃| 在线观看人妻少妇| 99热国产这里只有精品6| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕| 亚洲欧美中文字幕日韩二区| 国产精品一区二区性色av| 日韩欧美精品免费久久| 成年版毛片免费区| 久久精品久久精品一区二区三区| 99re6热这里在线精品视频| 国产综合精华液| 日韩在线高清观看一区二区三区| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| 亚洲内射少妇av| 国产免费又黄又爽又色| 大码成人一级视频| 国产成人a区在线观看| 亚洲自偷自拍三级| 丝袜喷水一区| 国产精品国产三级国产av玫瑰| 别揉我奶头 嗯啊视频| 看非洲黑人一级黄片| 国产欧美日韩一区二区三区在线 | 国产男女内射视频| 丝袜喷水一区| 99久久九九国产精品国产免费| 在线观看免费高清a一片| 国产精品一区二区三区四区免费观看| 国产成人a∨麻豆精品| 大香蕉久久网| 91久久精品国产一区二区成人| 国产欧美亚洲国产| 少妇猛男粗大的猛烈进出视频 | 国模一区二区三区四区视频| 国产高清不卡午夜福利| 国产在视频线精品| 亚洲国产高清在线一区二区三| 亚洲最大成人手机在线| 欧美xxⅹ黑人| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 久久久久九九精品影院| 大话2 男鬼变身卡| 男人爽女人下面视频在线观看| 在线看a的网站| 五月天丁香电影| 91久久精品电影网| 五月天丁香电影| 91久久精品电影网| 最近的中文字幕免费完整| 国产黄色免费在线视频| 建设人人有责人人尽责人人享有的 | 一本一本综合久久| 我的老师免费观看完整版| 99久国产av精品国产电影| 制服丝袜香蕉在线| 欧美老熟妇乱子伦牲交| 99九九线精品视频在线观看视频| 亚洲在久久综合| 色播亚洲综合网| 成人亚洲欧美一区二区av| 午夜福利视频1000在线观看| 国产高潮美女av| 97精品久久久久久久久久精品| 亚洲真实伦在线观看| 麻豆精品久久久久久蜜桃| 91aial.com中文字幕在线观看| 晚上一个人看的免费电影| 久久久成人免费电影| 国产精品一二三区在线看| 三级国产精品欧美在线观看| 黄色视频在线播放观看不卡| 久久精品人妻少妇| 热re99久久精品国产66热6| 黄片wwwwww| 97在线人人人人妻| 2018国产大陆天天弄谢| 午夜免费观看性视频| 久久久精品欧美日韩精品| 欧美人与善性xxx| 国产成人a∨麻豆精品| 大陆偷拍与自拍| 男插女下体视频免费在线播放| 欧美xxⅹ黑人| 亚洲成人久久爱视频| 汤姆久久久久久久影院中文字幕| 免费看日本二区| 成人国产麻豆网| 亚洲av成人精品一区久久| 97在线视频观看| 欧美人与善性xxx| 18禁在线无遮挡免费观看视频| 制服丝袜香蕉在线| 少妇 在线观看| 精华霜和精华液先用哪个| 久久精品国产鲁丝片午夜精品| 精品国产乱码久久久久久小说| 又爽又黄a免费视频| 国产成人福利小说| 欧美高清成人免费视频www| 嫩草影院新地址| 亚洲色图av天堂| 国产亚洲精品久久久com| 日韩av不卡免费在线播放| 蜜桃亚洲精品一区二区三区|