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

    Effect of single-legged coil on 3D plasma boundary corrugation in EAST

    2021-10-31 08:15:08QibinLUAN欒其斌WendaZHANG張聞達YoujieDENG鄧有杰YanfeiWANG王燕飛LiLI李莉YueqiangLIU劉鉞強XiaojiangHUANG黃曉江andFangchuanZHONG鐘方川
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:李莉

    Qibin LUAN (欒其斌), Wenda ZHANG (張聞達), Youjie DENG (鄧有杰),Yanfei WANG (王燕飛), Li LI (李莉), Yueqiang LIU (劉鉞強),Xiaojiang HUANG (黃曉江) and Fangchuan ZHONG (鐘方川)

    1 Faculty of Electronic information and Electrical Engineering, Dalian University of Technology, Dalian 116024, People’s Republic of China

    2 School of Physics, Dalian University of Technology, Dalian 116024, People’s Republic of China

    3 College of Science, Donghua University, Shanghai 201620, People’s Republic of China

    4 General Atomics, San Diego, CA 92186-5608, United States of America

    Abstract A single-legged coil behind the lower divertor and covering a 120°toroidal angle is utilized in a recent EAST discharge, for the purpose of increasing the wetted area of the divertor surface by locally modifying the magnetic field near the X-point.The plasma response, in particular, the plasma boundary surface corrugation due to the single-legged coil current, is modeled by the updated MARS-F code, by computing the plasma displacement for all important toroidal harmonics (n=1, 2, 4 and 5) associated with the partial toroidal coverage by the coil.The plasma response produced by the single-legged coil is found to be non-local and is of the kinkpeeling type.For a reference EAST plasma with a lower single-null magnetic configuration,the plasma boundary corrugation near the X-point, produced by the upper single-legged coil, is about twice as large as that produced by the lower single-legged coil,despite the proximity of the latter to the X-point.

    Keywords: single-legged coil, plasma response, EAST

    1.Introduction

    High heat flux on the plasma-facing components, due to various plasma perturbations, is a challenging issue for modern tokamaks, especially for future fusion reactors [1].Increasing the wetted area of the divertor is a plausible way to enhance the ability of the divertor to endure high heat flux from the plasma region.

    A new type of divertor has recently been developed in EAST, which is referred to as the Fishtail divertor [2].This divertor concept consists of a Fishtail-coil and the conventional divertor in EAST, where the Fishtail-coil is a singlelegged coil installed behind the lower divertor.In recent EAST discharge, this Fishtail-coil (single-legged coil) only covers a 120° toroidal angle.The other part of the (closed)coil current circuit is located outside the device far away from the plasma.Since the Fishtail-coil is close to the divertor, it can produce a magnetic field, with moderate coil current comparable to the equilibrium magnetic field near the divertor region.In other words,the magnetic footprint on the divertor surface can be modified by the Fishtail-coil.In particular,it is expected that the wetted area on the divertor surface, due to heat flux, can be widened.This is the main purpose of installing the Fishtail-coil in EAST experiments.

    However, besides changes to the vacuum magnetic field near the divertor, the plasma boundary can also be modified by the single-legged coil current, manifested as the boundary surface corrugation due to the plasma displacement in the edge region.The latter is due to the plasma response to the magnetic field perturbation produced by the coil current.As shown in many previous experimental and theoretical studies[3–6], the plasma response plays an important role in controlling the edge-localized modes (ELMs) by resonant magnetic perturbations(RMPs).It is therefore of high interest to investigate how the EAST plasma will respond to the magnetic field produced by the single-legged coil current,assuming the maximal current level as designed for experiments.

    Previous work [5, 7–9] has shown that 3D plasma boundary corrugation is an important measure of the plasma response to the externally applied magnetic field perturbation.In this work, we compute the plasma boundary displacement due to the single-legged coil current, utilizing the MARS-F code[10].Due to partial coverage of the toroidal angle by the single-legged coil, field perturbations with multiple toroidal harmonics are generated, and the plasma response to the dominant toroidal harmonics needs to be modeled.Because of the special coil current arrangement for the single-legged coil(compared to the conventional ELM control coils), the coil model within MARS-F has to be updated first,which we will also describe in this work.

    The paper is organized as follows.Section 2 describes the plasma and coil models assumed in this study, including the adopted EAST equilibrium.Section 3 reports the modeling results of the plasma boundary corrugation due to the singlelegged coil as designed in EAST.Section 4 reports the modeling results, hypothetically assuming that the singlelegged coil were installed at the upper divertor.Section 5 summarizes the work.

    2.Models and formulation

    In this section, we introduce the single-legged coil design in EAST and explain the corresponding numerical model implemented in MARS-F.We start by describing the plasma equilibrium adopted in this study.

    2.1.Plasma equilibrium

    The equilibrium is reconstructed from the EAST discharge 52340 at 3450 ms, representing a reference plasma from the ELM control experiments [11].Note that no Fishtail divertor had been installed when this discharge was carried out.Radial profiles of the plasma safety factor and the equilibrium pressure are shown in figure 1, whereψPis the normalized equilibrium poloidal magnetic flux.The safety factor values at the magnetic axis and at the 95% poloidal flux are 0.933 and 5.059, respectively.

    2.2.Single-legged coil configuration

    As mentioned before,the single-legged coil in EAST includes a single coil located at (1.87 ?1.09 m) on the (R,Z) poloidal plane and just behind the lower divertor[2].The coil location is indicted in figure 2(a),where for reference,we also plot the upper and lower rows of the ELM control coils in EAST.

    The ELM control coils in EAST are composed of eight window-frame coils in each row,with each coil covering 37°along the toroidal angle.This single-legged coil is very different—it in fact consists of only a single coil (with single turn) in the poloidal cross-section (figure 2(a)).It extends to 120° along the toroidal angle.The circuit is closed by wires that are installed farther away from the plasma.Therefore,these additional wire segments do not provide a significant contribution to the vacuum field perturbation near the plasma.

    The peculiar geometry of the single-legged coil generates multiple toroidal harmonics for the given coil current, as shown in figure 2(b).Besides then=0 component, which contributes to the equilibrium field, the main sidebands are the toroidal mode harmonics ofn=1, 2, 4 and 5.The toroidal spectrum of the coil current is calculated using the following formula:

    whereH(φ) is the Heaviside function equal to 1 for 0≤φ≤2π/ 3,and 0 elsewhere.The main non-axisymmetric sidebands are also documented in table 1.The actual current flowing in the single-legged coil is fixed at the maximal level of 5 kAt(single turn).Note that then=3k(k=1,2, 3, …) components of the coil current vanish due to the designed toroidal coverage(120°)of the single-legged coil,as is also evident from equation (1).Note also that the toroidal Fourier components of the coil current,I,ncoilare generally complex numbers.The phase of the complex number cannot be ignored when computing the total plasma response.

    Table 1.Dominant non-axisymmetric toroidal spectrum of thesingle-legged coil current.Shown are the coil current amplitude for each nth toroidal component, current amplitude ratio in percentage and current ratio including the phase.The actual current flowing in the single-legged coil is fixed at 5 kAt.

    In this work, we focus on modeling the plasma response and resulting boundary corrugation to the non-axisymmetric vacuum magnetic field components produced by the singlelegged coil.A DC current is assumed in the modeling,although AC current has been applied in experiments.In terms of the plasma response, the main AC effect is the reduction of the actual field reaching the plasma, due to shielding by various conducting structures surrounding the coil.Some of the conducting structures (e.g.the passive stabilizing plate) are of 3D geometry and cannot be properly modeled by MARS-F.However, this is not an issue with the DC single-legged coil current.

    We model the linear plasma response to the single-legged coil current, where the linear superposition principle applies.This means that we will compute the plasma response to each of the toroidal Fourier harmonics of the coil current with MARS-F, then sum all the response fields together to obtain the total field.

    2.3.MARS-F formulation and single-legged coil model

    In this work,we utilize the MARS-F code[10]to compute the linear plasma response to the single-legged coil.In the plasma region, MARS-F solves the linearized single fluid, resistive and full magnetohydrodynamic (MHD) equations in general toroidal geometry.For the response problem, the main equations read as follows:

    where the variablesξ,v, j,b andpare the plasma displacement, perturbed velocity, current, magnetic field and pressure,respectively.ρ,B,J,Pand Ω denote the equilibrium density, magnetic field, current, pressure and angular frequency of the toroidal rotation, respectively.Ris the plasma major radius,φ?the unit vector along the toroidal angle andnthe toroidal harmonic number.

    The single-legged coil current, as a source term, enters into the MHD equation via Ampere’s law ?× b =jRMP.The previous MARS-F implementation always assumes that the coil current flows in a closed circuit.In other words,there are always two parallel legs of the coil along the toroidal/poloidal angle, where currents flow in opposite directions.However, the single-legged coil has only one leg along the toroidal angle for the current-carrying segment that we model.Therefore, we need to update the code in order to accommodate the one-leg geometry.

    To model the conventional window-frame ELM control coils, MARS-F assumes the following distribution for the toroidal component of the coil current density along the poloidal angleχ,

    To model the single-legged coil,we only assume one leg with the coil current density represented by.The MARS-F code is updated to include this single-leg option.Obviously,the new implementation can also be used to recover the conventional window-frame coil, by assuming two sets of coils with opposite signs of current.This consistency check has been carried out using the upper set of ELM control coils shown in figure 2(a) as a test example.First, the perturbed magnetic field produced by the upper row is computed with standard MARS-F implementation.Next,two perturbed fields are computed with the new implementation, assuming the single-leg geometry and opposite coil currents.These two fields are then superposed and compared with that from the first step.The results are identical (figure 3), with or without the inclusion of the plasma response.Note that the amplitude of the coil current is assumed to be 1 kAt in this test,and the toroidal spectrum is chosen to ben=1–3.

    For example, figure 3(a) compares radial profiles of the amplitude of the poloidal Fourier harmonicsm=1–8 of then=1 radial field perturbation, between the two MARS-F implementations.Both the vacuum field(solid)and total field perturbation are compared together with the resistive plasma response(dashed).The corresponding curves overlap between the two implementations.The amplitude of the resonant harmonics is extracted from figure 3(a) and plotted in

    figure 3(b)forn=1 RMP.The results again overlap between the two implementations.Similar tests have also been performed forn=2 and 3 RMP, with good agreement, as reported in figures 3(c) and (d), respectively.

    3.Single-legged coil-induced plasma boundary corrugation

    We compute the plasma response to the single-legged coil,with the coil current model described in the previous section.We focus on the plasma boundary corrugation,i.e.the normal displacement of the plasma boundary surface as one of the key indicators of the plasma response.This choice is motivated by the following considerations.

    First, according to previous studies [5–7, 12, 13], the normal displacement of the plasma boundary near the X-point(for divertor plasmas),ξxpt,together with the resonant radial magnetic field perturbation near the plasma edge,,serves as important figures of merit for predicting ELM control by RMPs in tokamak experiments.These two figures of merit,which often predict similar optimal ELM control coil configurations, have been successfully applied to a range of ELM control (either mitigation or suppression) experiments including MAST [12], ASDEX Upgrade [6, 7], EAST [14],DIII-D[15,16]and HL-2A[17].The same criteria have also recently been applied to ELM control optimization in DEMO[18] and ITER [8, 9].

    Second, partial toroidal coverage of the single-legged coil implies the co-existence of several toroidal harmonics of comparable amplitude (figure 2(b)).The radial locations of resonant radial magnetic perturbations(near the plasma edge)are obviously different for differentn, making it difficult to employ thecriterion.The X-point displacement criterionξxpt,on the other hand, is defined for the plasma boundary surface (which is independent of the toroidal mode number)and can thus be easily superposed for RMP fields with multiple toroidal harmonics.This superposition is important in order to investigate the plasma response to the single-legged coil with partial toroidal coverage.

    We start by reporting the MARS-F computed plasma response for each individual-nRMP field produced by the single-legged coil current (at 5 kAt).Figure 4 shows results for then=1 component.For comparison, we consider here both figures of merit mentioned before.Figures 4(a) and (b)plot the radial distribution of the poloidal Fourier harmonics of the resonant radial magnetic field, for both then= 1 vacuum field and the totaln=1 field including the plasma response.Compared to the vacuum field, a strong screening of the resonant field perturbation occurs, in both the plasma core and edge regions.Figure 4(c) plots the normal displacement amplitude of the plasma surface along the geometric poloidal angle.The X-point is located at ?103°.We are mostly interested in the local maximum of the plasma surface displacement near the X-point, which is about 1.07 mm in our case.The overall structure of the plasma displacement, including that inside the plasma, is shown in figure 4(d).An interesting observation is that large plasma displacement occurs near both the top and bottom of the plasma(as also confirmed by figure 4(c)),despite the fact that the single-legged coil is located only near the lower divertor.This means that the plasma response is not a local phenomenon.The computed displacement is strongest near the plasma edge, indicating a peeling-like response.

    Figure 1.Radial profiles of (a) the safety factor and (b) the equilibrium pressure normalized by /μ 0,for the EAST discharge 52340 at 3450 ms.Here,ψp is the normalized equilibrium poloidal magnetic flux.

    Figure 2.(a)Poloidal cross-section of the plasma boundary(black curve),the ELM control coils(blue squares)and the single-legged coil(red circle).(b) Toroidal spectrum of the single-legged coil current, normalized by the actual current flowing in the coil.

    Figure 3.Comparison of the resonant radial magnetic field perturbation between the standard and new implementations of MARS-F for the coil current.Both the vacuum field and total field perturbation are compared together with the plasma response: (a) radial profiles of the m=1–8 poloidal Fourier harmonics of the n=1 radial field perturbation, (b)–(d) amplitude of the pitch resonant radial magnetic field components at the corresponding rational surface locations for the n=1,2,3 perturbation,respectively.The upper row of the ELM control coils in EAST is considered, with a coil current of 1 kAt.Vertical dashed lines in (a) indicate the radial location of the n=1 rational surfaces.

    Figure 4.MARS-F computed plasma response to the n=1 single-legged coil current.(a) Radial distribution of the poloidal Fourier harmonics of the resonant radial magnetic field perturbation,with the vacuum field(in blue)being compared to the total field perturbation(in red)including the plasma response,(b)amplitude of the resonant radial magnetic field components at the corresponding rational surfaces,for the vacuum field(in blue)and the total field(in red),respectively,(c)amplitude of the normal displacement of the plasma boundary surface along the geometric poloidal angle, (d) amplitude of the normal displacement at the poloidal cross-section of the plasma.Vertical dashed l ines in(a)indicate the radial location of the n=1 rational surfaces.X-point corresponds to the poloidal angle of ?103°in(c).Single-leggedcoil current amplitude for n=1 harmonic? is assumed to be 5 kAt.

    Figure 5.MARS-F computed plasma response to the (a)–(b) n=2, (c)–(d) n=4 and (e)–(f) n=5 components of the single-legged coil current.Amplitude of the normal displacement of the plasma boundary surface along the geometric poloidal angle (left panels) and normal displacement at the poloidal cross-section(right panels).X-point corresponds to the poloidal angle of ?103°in the left panels.Single-legged coil current amplitude for each harmonic is assumed to be 5 kAt.

    Figure 6.Comparison of the overall plasma boundary corrugation amplitude along both the poloidal and toroidal angles,produced by(a)the single-legged coil and(b)the lower row of ELM control coils.Here,the amplitude of the coil current is fixed at 5 kAt in both(a)and(b),and the same set of weighting factors, Cn from equation (10), are assumed for the toroidal sidebands of the single-legged coil and the ELM control coils.

    Figure 7.MARS-F computed plasma response to the(a)–(b)n=1,(c)–(d)n=2,(e)–(f)n=4 and(g)–(h)n=5 components of the upper single-legged coil current.Amplitude of the normal displacement of the plasma boundary surface along the geometric poloidal angle (left panels) and the normal displacement at the poloidal cross-section (right panels) are plotted.X-point corresponds to the poloidal angle of?103° in the left panels.Single-legged coil current amplitude for each harmonic is assumed to be 5 kAt.

    Figure 8.Computed plasma boundary corrugation along both the poloidal and toroidal angles, assuming the upper single-legged coil with 5 kAt current.The n=1, 2, 4, 5 toroidal harmonics of the plasma response are superposed.

    The normal displacement amplitude forn=2, 4 and 5 single-legged coil currents is summarized in figure 5.The corresponding coil current amplitude is again fixed at 5 kAt,in order to compare it with the results for then=1 coil current.The MARS-F computed local maxima of the plasma surface displacement near the X-point are 0.757, 0.099 and 0.022 mm, with then=2, 3 and 5 single-legged coil currents,respectively.It is evident that,with the same amplitude of the coil current for each toroidal harmonic, the plasma response to then=1 (5) coil configuration produces the largest (smallest) X-point displacement.

    Several other interesting observations can also be made from figure 5.First, a relatively large plasma displacement is generally computed near the top and bottom of the torus.This is despite the fact that the external field perturbation is produced by the single-legged coil current located only near the bottom (behind the lower divertor).This again confirms the non-local nature of the plasma response,for all the low-ncoil current components.

    Next, although the plasma surface displacement near the X-point,∣ξxpt∣,is weakest for then=5 coil current, the ratio ofξmid∣∣toξxpt∣∣is largest.Here,ξmid∣∣denotes the plasma surface displacement at the outboard mid-plane.This is a clear indication of the kink-like plasma response [12] to then=5 single-legged coil current.

    As mentioned before, the plasma response to each individual toroidal component of the single-legged coil current can be superposed, in order to construct the overall response to the ‘physical’ coil current that contains a range of toroidal spectra.In particular,the overall plasma boundary corrugation amplitude,which now varies along both poloidal and toroidal angles, can be obtained utilizing the coil current weighting factors documented in table 1

    whereCnis the corresponding (complex) weighting coefficient for eachn-component of the single-legged coil current.

    In order to better understand the efficiency of the singlelegged coil geometry (in terms of the generated plasma corrugation magnitude),we compare it with that produced by the lower row of the ELM control coils in EAST, but assuming the same coil current and same sideband distribution as that of the single-legged coil.The plasma displacements from different toroidal harmonics are superposed following equation (10), and compared between the single-legged and RMP coils in figure 6.Two interesting observations can be made here.First, the maximum displacement along the toroidal angle and near the X-point (at ?103° poloidal angle),∣ξxpt∣ ,is 0.35 mm due to the single-legged coil current(figure 6(a)).This value is about four times smaller than that(1.27 mm) produced by the lower row of the ELM control coils with the same current configuration.This difference is due to two factors: (i) the lower row of ELM control coils is closer to the plasma than the single-legged coil; (ii)the ELM control coils geometrically cover a finite poloidal angle due to their window-frame structure.

    Next, we find it interesting that despite being a singlelegged coil (with a single turn) located behind the lower divertor, the single-legged coil is capable of generating a plasma response that is global, i.e.not just located near the X-point.Large displacement is also triggered near the top of the plasma, which is far from the single-legged coil location.This global response effect is similar to that produced by the ELM control coils.In fact, maximum displacement occurs near the top of the plasma at the poloidal angle of 106°,with both types of coil geometry.The corresponding values are 0.42 mm with the single-legged coil and 1.5 mm with the ELM control coils.

    4.Effect of upper single-legged coil on plasma boundary corrugation

    Due to engineering constraints, the single-legged coil is applied for an upper single-null magnetic configuration in EAST experiments,in contrast to the lower single-null plasma that we adopt in this study.Therefore, to better simulate the experimental situation, we assume that the single-legged coil is located near the upper top,at(R=1.87 m,Z=1.09 m)in the MARS-F model,while keeping the same lower single-null plasma equilibrium.In what follows,we investigate the effect of this new single-legged coil configuration on the plasma response.

    We again start by comparing the plasma response level among various toroidal sidebands of the single-legged coil current, assuming the same current of 5 kAt for eachn-harmonic.The computed normal displacements for then=1,2,4, 5 components, plotted in figure 7, show that the local maxima of the plasma surface corrugation near the X-point are 2.022, 0.831, 0.11 and 0.024 mm, respectively, i.e.sharply decreasing with increasingn-number.

    Compared to the results reported in section 3, we find that the lower and upper single-legged coils produce a similar plasma response to each toroidal harmonic.The strongest plasma surface distortion generally occurs near the top and bottom of the torus.The plasma displacement near the(lower single-null) X-point, produced by the upper single-legged coil,is slightly larger than that produced by the lower singlelegged coil.This again indicates the non-local nature of the plasma response.Despite a generally small displacement, the displacement ratio of the outboard mid-plane to the X-point,∣ξmid∣ /∣ξxpt∣,is largest for then=5 toroidal component.

    Figure 8 plots the computed plasma boundary corrugation along the poloidal and toroidal angles, produced by the upper single-legged coil, after superposition of then=1, 2,4, 5 contributions according to equation (10).Compared to the lower single-legged coil (figure 6(a)), the upper singlelegged coil produces about double (0.651 mm versus 0.344 mm) the plasma displacement near the X-point,assuming the same (lower single-null) plasma and same coil current of 5 kAt.Note that this is despite the much larger distance between the upper single-legged coil and the X-point.However,the X-point displacement with both singlelegged coils is still lower than that produced by the lower row of ELM control coils (figure 6(b)) with the same coil current configuration.

    5.Conclusion and discussion

    The Fishtail divertor, composed of a single-legged coil and the lower divertor, is designed in EAST for the purpose of modifying the local magnetic field structure near the X-point.However, our modeling shows that the single-legged coil is capable of triggering a more global plasma response.

    The single-legged coil is modeled as a single-current filament in the updated MARS-F code, covering a 120° toroidal angle.This finite toroidal coverage is realized by considering a range of important toroidal sidebandsn=1, 2, 4,5.Modeling result shows that the plasma response amplitude,in terms of the boundary corrugation, drops sharply with increasing toroidal mode number, assuming the same coil current.On the other hand, based on the computed response for each individual toroidal harmonic, the total plasma response to the single-legged coil current can be reconstructed via linear superposition with proper weighting factors, producing plasma surface corrugation that varies along both poloidal and toroidal angles.

    The total response shows that the single-legged coil triggers large plasma displacement near the top and bottom of the torus.This holds, independent of whether the single-legged coil is placed near the top or bottom of the plasma,for the same lower single-null equilibrium.Thus, the modeling results show that (i) the plasma response to the single-legged coil current is non-local and (ii) the total response is predominantly of the kink-peeling type.

    Another interesting result is that we find that for a discharge with the lower single-null magnetic configuration, the upper single-legged coil produces about double the X-point displacement compared to the lower single-legged coil,despite the fact that the latter is much closer to the X-point.There seems to be a certain correlation between the X-point displacement and the magnetic footprint on the divertor surface[19,20].The last result here suggests that a single-legged coil, placed opposite to the X-point on the poloidal plane (as in the EAST experiment),may be more effective in producing a wider magnetic footprint.Direct modeling of the magnetic footprint due to the single-legged coil,based on the MARS-F computed plasma response field, is the next step in our investigation.

    Acknowledgments

    This work is funded by the Fundamental Research Funds for the Central Universities (No.2232021G-10), National Natural Science Foundation of China (Nos.12075053, 11505021 and 11975068) and the National Key R&D Program of China(Nos.2017YFE030052, 2017YFE0301100, 2017YFE0301104,2017YFE0300500 and 2017YFE0300501).The work was also supported by the U.S.DoE Office of Science (Nos.DEFG02–95ER54309 and DEFC02–04ER54698).This paper was prepared as an account of work sponsored by an agency of the United States Government.Neither the United States Government nor any agency thereof,nor any of their employees,makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy,completeness,or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights.Reference herein to any specific commercial product,process,or service by trade name, trademark, manufacturer, or otherwise,does not necessarily constitute or imply its endorsement,recommendation or favouring by the United States Government or any agency thereof.The views and opinions of the authors expressed herein do not necessarily represent or reflect those of the United States Government or any agency thereof.The views and opinions expressed herein do not necessarily reflect those of the European Commission.

    猜你喜歡
    李莉
    Impact renaming non-alcoholic fatty liver disease to metabolic associated fatty liver disease in prevalence, characteristics and risk factors
    THE REGULARITY CRITERIA OF WEAK SOLUTIONS TO 3D AXISYMMETRIC INCOMPRESSIBLE BOUSSINESQ EQUATIONS?
    李莉作品(一)
    大眾文藝(2021年17期)2021-09-29 03:03:20
    李莉作品(二)
    大眾文藝(2021年17期)2021-09-29 03:03:20
    Dynamic and inner-dressing control of four-wave mixing in periodically-driven atomic system?
    故鄉(xiāng)一把土
    當代音樂(2019年4期)2019-06-11 14:54:47
    裝錯芯片的機器人
    Non-Markovianity Measure Based on Brukner–Zeilinger Invariant Information for Unital Quantum Dynamical Maps?
    樹葉上的優(yōu)點
    Numerical simulation of viscous flow past an oscillating square cylinder using a CIP-based model*
    一进一出抽搐gif免费好疼| 欧美色欧美亚洲另类二区 | 久久人妻av系列| 国产私拍福利视频在线观看| 免费在线观看黄色视频的| 日本 av在线| 亚洲成人免费电影在线观看| 波多野结衣巨乳人妻| 18禁裸乳无遮挡免费网站照片 | 午夜福利18| 伊人久久大香线蕉亚洲五| 黄色 视频免费看| 亚洲欧美激情在线| 欧美一级毛片孕妇| 欧美激情高清一区二区三区| 亚洲久久久国产精品| 成人国产一区最新在线观看| 亚洲中文日韩欧美视频| 9色porny在线观看| 在线av久久热| 男女下面插进去视频免费观看| 午夜a级毛片| 12—13女人毛片做爰片一| 日韩国内少妇激情av| 91大片在线观看| 非洲黑人性xxxx精品又粗又长| 嫁个100分男人电影在线观看| 国产精品日韩av在线免费观看 | 久久中文字幕人妻熟女| 欧美日本亚洲视频在线播放| 999精品在线视频| 亚洲在线自拍视频| 午夜久久久在线观看| 九色国产91popny在线| 亚洲在线自拍视频| 看黄色毛片网站| 国产成人精品久久二区二区免费| 少妇的丰满在线观看| 欧美色欧美亚洲另类二区 | 亚洲成av人片免费观看| 久久九九热精品免费| 男女下面插进去视频免费观看| 国产成人精品久久二区二区91| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲三区欧美一区| 操出白浆在线播放| 亚洲人成电影免费在线| 欧美成狂野欧美在线观看| 91av网站免费观看| av在线天堂中文字幕| 亚洲第一青青草原| 国产欧美日韩精品亚洲av| 色哟哟哟哟哟哟| 亚洲欧美激情在线| 国产亚洲精品av在线| 亚洲在线自拍视频| 国产区一区二久久| 国产av又大| 久久久国产成人精品二区| 99香蕉大伊视频| 好男人在线观看高清免费视频 | 757午夜福利合集在线观看| 中文字幕最新亚洲高清| 多毛熟女@视频| 无人区码免费观看不卡| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 欧洲精品卡2卡3卡4卡5卡区| 日本vs欧美在线观看视频| 老司机午夜十八禁免费视频| 精品国产美女av久久久久小说| 久久精品aⅴ一区二区三区四区| 欧美日本亚洲视频在线播放| 日韩免费av在线播放| 999精品在线视频| 丝袜在线中文字幕| 女人被躁到高潮嗷嗷叫费观| 午夜精品国产一区二区电影| 精品久久久精品久久久| svipshipincom国产片| 伦理电影免费视频| www.精华液| 精品欧美一区二区三区在线| 日本一区二区免费在线视频| 成年女人毛片免费观看观看9| 国产精品 欧美亚洲| 一级毛片精品| 欧美激情久久久久久爽电影 | 亚洲欧美日韩另类电影网站| 999精品在线视频| 亚洲国产毛片av蜜桃av| 精品一品国产午夜福利视频| 天天躁夜夜躁狠狠躁躁| www日本在线高清视频| 麻豆一二三区av精品| 精品一区二区三区四区五区乱码| 免费搜索国产男女视频| 国产成人精品无人区| 久久精品人人爽人人爽视色| 变态另类丝袜制服| 久久香蕉国产精品| 成在线人永久免费视频| 亚洲精品国产精品久久久不卡| 日韩欧美三级三区| 亚洲激情在线av| 国产精品一区二区三区四区久久 | 国产精品久久电影中文字幕| 欧美大码av| 欧美乱色亚洲激情| 国产亚洲av高清不卡| 亚洲中文字幕日韩| 成人三级黄色视频| 午夜福利18| 日本欧美视频一区| 亚洲第一av免费看| 最好的美女福利视频网| 最好的美女福利视频网| 巨乳人妻的诱惑在线观看| av视频免费观看在线观看| 久久国产精品男人的天堂亚洲| 国产一区在线观看成人免费| 亚洲激情在线av| 19禁男女啪啪无遮挡网站| 成在线人永久免费视频| 亚洲精品中文字幕一二三四区| 免费人成视频x8x8入口观看| 久久午夜亚洲精品久久| 性少妇av在线| 成年人黄色毛片网站| 久久久久久久久久久久大奶| 精品国内亚洲2022精品成人| 别揉我奶头~嗯~啊~动态视频| 久久久水蜜桃国产精品网| 97碰自拍视频| 精品欧美一区二区三区在线| 天堂√8在线中文| 波多野结衣巨乳人妻| 欧美成人性av电影在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲熟女毛片儿| 正在播放国产对白刺激| 女人精品久久久久毛片| aaaaa片日本免费| 视频在线观看一区二区三区| 香蕉丝袜av| 亚洲成av片中文字幕在线观看| 乱人伦中国视频| 亚洲国产高清在线一区二区三 | 99re在线观看精品视频| 国产精品免费视频内射| 免费久久久久久久精品成人欧美视频| 黄色片一级片一级黄色片| 1024视频免费在线观看| 久久久久久久精品吃奶| 美女高潮到喷水免费观看| 日韩欧美一区视频在线观看| 在线观看www视频免费| 又黄又爽又免费观看的视频| 日日摸夜夜添夜夜添小说| 色尼玛亚洲综合影院| 高清在线国产一区| 中文字幕色久视频| 国产精品秋霞免费鲁丝片| 如日韩欧美国产精品一区二区三区| 久久久久久久久久久久大奶| 中文字幕人妻熟女乱码| 成人免费观看视频高清| av欧美777| 自线自在国产av| 好男人电影高清在线观看| 一区二区三区高清视频在线| 变态另类成人亚洲欧美熟女 | 亚洲av熟女| 久久伊人香网站| 777久久人妻少妇嫩草av网站| 久久久久亚洲av毛片大全| 免费高清视频大片| 国产三级黄色录像| 99riav亚洲国产免费| 国产单亲对白刺激| 麻豆一二三区av精品| 国产精品1区2区在线观看.| 国产野战对白在线观看| 免费高清视频大片| 一边摸一边抽搐一进一小说| 亚洲欧美精品综合久久99| 一级毛片精品| 亚洲欧美激情在线| 丰满的人妻完整版| 18禁美女被吸乳视频| 欧美色欧美亚洲另类二区 | 国产av精品麻豆| 国产精品永久免费网站| 男女午夜视频在线观看| a在线观看视频网站| 在线国产一区二区在线| 岛国视频午夜一区免费看| 国产亚洲欧美98| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 日日爽夜夜爽网站| 国产精品,欧美在线| 日韩一卡2卡3卡4卡2021年| 曰老女人黄片| 99国产综合亚洲精品| 国产精品永久免费网站| av天堂在线播放| 久久精品成人免费网站| 久久精品亚洲精品国产色婷小说| 久久中文看片网| 国产97色在线日韩免费| 日本黄色视频三级网站网址| 久久人妻熟女aⅴ| 夜夜躁狠狠躁天天躁| 亚洲国产毛片av蜜桃av| 热99re8久久精品国产| 在线观看66精品国产| 国产一区二区三区在线臀色熟女| 久久热在线av| 日韩 欧美 亚洲 中文字幕| 欧美日韩一级在线毛片| av免费在线观看网站| 国产精品日韩av在线免费观看 | 女生性感内裤真人,穿戴方法视频| 精品日产1卡2卡| 久久婷婷成人综合色麻豆| av视频在线观看入口| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 国产成人av教育| 成年女人毛片免费观看观看9| 久久精品亚洲精品国产色婷小说| 乱人伦中国视频| 天堂动漫精品| 久久久精品欧美日韩精品| 男人操女人黄网站| 亚洲国产精品久久男人天堂| 九色亚洲精品在线播放| 日本 av在线| 久久精品成人免费网站| 亚洲天堂国产精品一区在线| 999久久久精品免费观看国产| 亚洲激情在线av| 国产精华一区二区三区| 人人妻人人爽人人添夜夜欢视频| 午夜福利视频1000在线观看 | 最好的美女福利视频网| 99在线人妻在线中文字幕| 国产精品 国内视频| 亚洲av成人一区二区三| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 欧美日韩福利视频一区二区| 亚洲电影在线观看av| 亚洲一码二码三码区别大吗| 波多野结衣巨乳人妻| 91麻豆精品激情在线观看国产| 亚洲成国产人片在线观看| 国产一级毛片七仙女欲春2 | svipshipincom国产片| 亚洲少妇的诱惑av| 精品熟女少妇八av免费久了| 动漫黄色视频在线观看| 在线观看午夜福利视频| 免费无遮挡裸体视频| 俄罗斯特黄特色一大片| 亚洲av片天天在线观看| 久久香蕉国产精品| 啦啦啦观看免费观看视频高清 | 亚洲色图综合在线观看| 深夜精品福利| 久久人妻av系列| 国产不卡一卡二| 国产成人一区二区三区免费视频网站| 亚洲精品国产区一区二| 69av精品久久久久久| 男男h啪啪无遮挡| 亚洲成人国产一区在线观看| 亚洲精品国产一区二区精华液| 国产精品影院久久| 国产成人一区二区三区免费视频网站| 老熟妇仑乱视频hdxx| 麻豆av在线久日| 免费在线观看视频国产中文字幕亚洲| 国产精品日韩av在线免费观看 | 亚洲无线在线观看| 国产精品电影一区二区三区| 午夜福利,免费看| 亚洲精品国产一区二区精华液| 亚洲少妇的诱惑av| 极品人妻少妇av视频| 亚洲精品在线美女| 在线观看免费日韩欧美大片| 亚洲欧美日韩另类电影网站| 1024视频免费在线观看| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 久久精品亚洲熟妇少妇任你| 搡老妇女老女人老熟妇| 搞女人的毛片| 黄色视频,在线免费观看| 一级毛片精品| 女生性感内裤真人,穿戴方法视频| 九色国产91popny在线| 欧美黑人精品巨大| 欧洲精品卡2卡3卡4卡5卡区| 久久亚洲精品不卡| 亚洲第一欧美日韩一区二区三区| 亚洲少妇的诱惑av| 女警被强在线播放| 黄色 视频免费看| 在线观看午夜福利视频| 一二三四在线观看免费中文在| 嫩草影院精品99| 国产精品野战在线观看| 美女高潮到喷水免费观看| 91老司机精品| 欧美黑人欧美精品刺激| 两个人看的免费小视频| 在线永久观看黄色视频| 日韩欧美国产在线观看| 欧美老熟妇乱子伦牲交| 午夜久久久久精精品| 精品人妻在线不人妻| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| 性少妇av在线| 久久久国产欧美日韩av| 中文字幕色久视频| 国产xxxxx性猛交| 女人被狂操c到高潮| 国产男靠女视频免费网站| 欧美日韩一级在线毛片| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 日韩中文字幕欧美一区二区| 在线观看免费日韩欧美大片| 国产99白浆流出| 18美女黄网站色大片免费观看| 国产精品乱码一区二三区的特点 | 日韩精品中文字幕看吧| 日韩欧美免费精品| 久久国产亚洲av麻豆专区| 岛国视频午夜一区免费看| 国产av一区在线观看免费| 波多野结衣av一区二区av| 三级毛片av免费| 男人舔女人的私密视频| 精品无人区乱码1区二区| 亚洲 国产 在线| 欧美色视频一区免费| 动漫黄色视频在线观看| 99在线视频只有这里精品首页| 久久久久久免费高清国产稀缺| 亚洲国产精品成人综合色| 色婷婷久久久亚洲欧美| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 性色av乱码一区二区三区2| 国产精品野战在线观看| 亚洲av电影在线进入| 国产1区2区3区精品| 美女高潮到喷水免费观看| 国产精品野战在线观看| 精品久久久久久久毛片微露脸| 丝袜美足系列| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 俄罗斯特黄特色一大片| 久久精品亚洲熟妇少妇任你| av网站免费在线观看视频| 久久久久国产精品人妻aⅴ院| 18禁美女被吸乳视频| 国产精品乱码一区二三区的特点 | 国产一区二区三区视频了| 国产精华一区二区三区| 一区二区三区高清视频在线| 久久这里只有精品19| 男人舔女人的私密视频| 国产亚洲欧美98| 亚洲成人精品中文字幕电影| 色av中文字幕| 久久天堂一区二区三区四区| 夜夜躁狠狠躁天天躁| 国产欧美日韩精品亚洲av| 免费在线观看日本一区| 90打野战视频偷拍视频| 国产精品亚洲美女久久久| 国产亚洲精品久久久久5区| 日韩国内少妇激情av| 满18在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 电影成人av| 女同久久另类99精品国产91| 久久性视频一级片| 亚洲成人久久性| 国产熟女午夜一区二区三区| 波多野结衣高清无吗| e午夜精品久久久久久久| 久久人人精品亚洲av| 国产97色在线日韩免费| 国产精品香港三级国产av潘金莲| 亚洲天堂国产精品一区在线| 天堂√8在线中文| 91精品国产国语对白视频| 琪琪午夜伦伦电影理论片6080| 一区二区三区国产精品乱码| 俄罗斯特黄特色一大片| 手机成人av网站| 麻豆av在线久日| 精品国产亚洲在线| 黄色片一级片一级黄色片| 亚洲少妇的诱惑av| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 我的亚洲天堂| 久久久国产精品麻豆| 国产熟女xx| 精品免费久久久久久久清纯| 久久国产精品影院| 欧美成狂野欧美在线观看| 纯流量卡能插随身wifi吗| 国产成人av教育| 久久久国产欧美日韩av| 欧美国产日韩亚洲一区| 中文字幕人妻丝袜一区二区| 午夜日韩欧美国产| 大型av网站在线播放| 麻豆成人av在线观看| 国产亚洲欧美精品永久| av视频免费观看在线观看| 久久热在线av| 91成年电影在线观看| 黑人操中国人逼视频| 免费女性裸体啪啪无遮挡网站| 成年版毛片免费区| 成人手机av| 天天躁夜夜躁狠狠躁躁| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 精品国内亚洲2022精品成人| 免费看美女性在线毛片视频| 亚洲色图av天堂| www.www免费av| 女人被躁到高潮嗷嗷叫费观| 两人在一起打扑克的视频| 午夜影院日韩av| 国产精品久久电影中文字幕| 天堂动漫精品| 香蕉久久夜色| 亚洲五月天丁香| 欧美成人午夜精品| 欧美日韩福利视频一区二区| 美女高潮喷水抽搐中文字幕| 人人澡人人妻人| 91成年电影在线观看| 欧美国产精品va在线观看不卡| 久久久久久久久久久久大奶| 国产成人欧美| 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片 | 亚洲激情在线av| 国产激情欧美一区二区| 黄频高清免费视频| 亚洲中文字幕日韩| 免费在线观看亚洲国产| 国产真人三级小视频在线观看| 成人亚洲精品av一区二区| 久久精品国产清高在天天线| 波多野结衣巨乳人妻| 51午夜福利影视在线观看| 精品不卡国产一区二区三区| 人人妻人人澡人人看| 亚洲第一av免费看| 精品一区二区三区av网在线观看| 日韩视频一区二区在线观看| 欧美不卡视频在线免费观看 | 老鸭窝网址在线观看| 亚洲一区二区三区色噜噜| 免费在线观看亚洲国产| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 极品教师在线免费播放| 又紧又爽又黄一区二区| 欧美日韩一级在线毛片| 国产真人三级小视频在线观看| 国产日韩一区二区三区精品不卡| 亚洲狠狠婷婷综合久久图片| 欧美精品亚洲一区二区| 午夜精品久久久久久毛片777| 国产三级在线视频| 国产精品一区二区在线不卡| 一个人免费在线观看的高清视频| 激情在线观看视频在线高清| 别揉我奶头~嗯~啊~动态视频| 波多野结衣av一区二区av| 国产精品一区二区在线不卡| 国产一区二区三区在线臀色熟女| 一二三四在线观看免费中文在| av福利片在线| avwww免费| 脱女人内裤的视频| 亚洲一码二码三码区别大吗| 欧美亚洲日本最大视频资源| 老司机午夜福利在线观看视频| 亚洲,欧美精品.| 黄色毛片三级朝国网站| 一个人免费在线观看的高清视频| 黑人操中国人逼视频| 美女午夜性视频免费| 在线天堂中文资源库| 两个人看的免费小视频| 在线免费观看的www视频| 日韩欧美国产在线观看| 亚洲欧美日韩无卡精品| 免费不卡黄色视频| 日韩精品青青久久久久久| 黄色片一级片一级黄色片| 久久 成人 亚洲| 日韩 欧美 亚洲 中文字幕| 色综合欧美亚洲国产小说| av天堂在线播放| 操美女的视频在线观看| 波多野结衣高清无吗| 国产麻豆69| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 此物有八面人人有两片| 99久久国产精品久久久| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 久久久久久免费高清国产稀缺| 午夜福利视频1000在线观看 | 国产精品爽爽va在线观看网站 | 日本撒尿小便嘘嘘汇集6| 免费少妇av软件| 日韩大尺度精品在线看网址 | 国产一区二区三区视频了| 国产亚洲欧美98| 亚洲第一青青草原| 搞女人的毛片| 欧美色视频一区免费| 国产真人三级小视频在线观看| 久久青草综合色| 99国产精品一区二区蜜桃av| 天天躁夜夜躁狠狠躁躁| 亚洲第一电影网av| 亚洲性夜色夜夜综合| 一区二区三区激情视频| 国产成人欧美在线观看| 长腿黑丝高跟| 国产乱人伦免费视频| 国产麻豆69| 国产精品,欧美在线| 精品卡一卡二卡四卡免费| 99riav亚洲国产免费| 色在线成人网| 一级毛片女人18水好多| 国产精品久久久久久人妻精品电影| 亚洲专区中文字幕在线| 亚洲精品国产色婷婷电影| 黑丝袜美女国产一区| 在线观看免费视频日本深夜| 国产97色在线日韩免费| 99热只有精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 人人澡人人妻人| 精品免费久久久久久久清纯| 国产精品综合久久久久久久免费 | 午夜影院日韩av| 亚洲欧美日韩高清在线视频| 亚洲国产精品成人综合色| 亚洲av成人不卡在线观看播放网| 91大片在线观看| 搡老妇女老女人老熟妇| 纯流量卡能插随身wifi吗| www.熟女人妻精品国产| 欧美激情高清一区二区三区| 亚洲成人久久性| 国产精品 欧美亚洲| 99国产综合亚洲精品| 999精品在线视频| www.精华液| av在线天堂中文字幕| 国产99白浆流出| 午夜久久久在线观看| 成在线人永久免费视频| 身体一侧抽搐| svipshipincom国产片| 久久中文字幕人妻熟女| 国产成人欧美| 高清在线国产一区| 国产伦人伦偷精品视频| 亚洲全国av大片| 国产人伦9x9x在线观看| 精品免费久久久久久久清纯| 999久久久国产精品视频| 亚洲天堂国产精品一区在线| 亚洲精品中文字幕一二三四区| 久久人人97超碰香蕉20202| 久久人妻av系列| cao死你这个sao货| 亚洲精品一区av在线观看| 淫妇啪啪啪对白视频| 欧美在线一区亚洲| 亚洲av成人不卡在线观看播放网| 久久久久久久久免费视频了| 久久亚洲真实| 91成人精品电影| 午夜免费鲁丝| 中亚洲国语对白在线视频| 亚洲国产精品999在线|