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

    Toroidal component of velocity for geodesic acoustic modes in the edge plasmas of the J-TEXT tokamak

    2021-10-31 08:15:08YaoyuXIE謝耀禹KaijunZHAO趙開君ZhipengCHEN陳志鵬JiaqiDONG董家齊KimitakaITOHZhongyongCHEN陳忠勇YuejiangSHI石躍江YonghuaDING丁永華JunCHENG程鈞LongwenYAN嚴(yán)龍文HaiLIU劉海ZhifengCHENG程芝峰ZhoujunYANG楊州軍NengchaoWANG王能超LuWANG王璐JianqiangXU許健強(qiáng)Yunfe
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:王璐龍文永華

    Yaoyu XIE (謝耀禹), Kaijun ZHAO (趙開君), Zhipeng CHEN (陳志鵬),Jiaqi DONG (董家齊), Kimitaka ITOH, Zhongyong CHEN (陳忠勇),Yuejiang SHI (石躍江), Yonghua DING (丁永華), Jun CHENG (程鈞),Longwen YAN (嚴(yán)龍文), Hai LIU (劉海), Zhifeng CHENG (程芝峰) ,Zhoujun YANG (楊州軍), Nengchao WANG (王能超), Lu WANG (王璐) ,Jianqiang XU (許健強(qiáng)), Yunfeng LIANG (梁云峰) and J-TEXT Team,9

    1 School of Nuclear Science and Engineering, East China University of Technology, Nanchang 330013,People’s Republic of China

    2 International Joint Research Laboratory of Magnetic Confinement Fusion and Plasma Physics,State Key Laboratory of Advanced Electromagnetic Engineering and Technology,School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, People’s Republic of China

    3 Southwestern Institute of Physics, Chengdu 610041, People’s Republic of China

    4 Institute of Science and Technology Research, Chubu University, Aichi 487-8501, Japan

    5 Research Center for Plasma Turbulence, Kyushu University, Kasuga 816-8580, Japan

    6 National Institute for Fusion Science, Toki 509-5292, Japan

    7 ENN Science and Technology Development Co., Ltd, and Hebei Key Laboratory of Compact Fusion Langfang 065001, People’s Republic of China

    8 Institute of Fusion Science, School of Physical Science and Technology, Southwest Jiaotong University,Chengdu 610031, People’s Republic of China

    Abstract The toroidal component of the velocity for geodesic acoustic modes (GAMs) is first demonstrated.Multiple Langmuir probe arrays set up near the top tokamak of the J-TEXT were utilized for this study.A significant peak at the GAM frequency is observed in Mach number fluctuations.The toroidal velocity for the GAMs is estimated as~10–100 m s?1 and increases with the poloidal velocity.The ratio of toroidal component to the poloidal one of the velocity is mainly located in the interval between 0.3 and 1.0.With higher safety factors q,the ratio almost does not change with decreasing the safety factor, whereas it goes up sharply at low q.The coherencies between poloidal electric fields and Mach number fluctuations in turbulence frequency bands are also evaluated, and are higher than those between radial electric fields and Mach number fluctuations.

    Keywords: zonal flow, geodesic acoustic mode, tokamak, return flow, Mach number

    1.Introduction

    Zonal flows (ZFs) self-generated by turbulence widely exist in nature and laboratory plasmas.Typical examples include space plasmas and magnetically confined laboratories,as well as atmospheres of stars.This study in the fusion field aims at understanding the nonlinear mechanism responsible for the coherent structure formation and for anomalous cross-field transport induced by turbulent fluctuations.

    It has been recognized in the last decade that the E×B sheared flows can reduce or even suppress the microscale turbulence and thus transport due to turbulence.Generally,the sheared flows exist in the forms of the mean flows due to ion pressure gradient and the zonal flows driven by turbulence[1–4].The sheared flow formation due to turbulence nonlinear interaction via Reynolds stress has been well documented[3].This suggests that a feedback loop exists between turbulence and sheared flows.The turbulence-flow system can organize and regulate itself [5, 6].

    Two kinds of ZFs have been detected in toroidal plasmas,called the low-frequency zonal flows(LFZFs)[4,7]with zero frequency and the geodesic acoustic modes (GAMs) [8, 9]having finite frequency,whereγiis the adiabatic index,TiandTeare ion and electron temperature, respectively, andMiandRare ion mass and major radius of the torus, respectively.The ratio of the toroidal component to the poloidal one of the velocity is expected to be ?2qcosθ for the LFZFs (whereqis the safety factor and θ is the poloidal angle), and that for the GAMs isq?1cosθ [3].

    In the last two decades,the LFZFs and GAMs have been studied extensively [10–24].The common features of GAMs and LFZFs have been obtained in a great number of devices.They all show the poloidal and toroidal symmetries, and the finite radial wave vectors [10–15].The density fluctuation of the GAMs in the top plasmas of the tokamak is higher than that in the middle plane [16].The bicoherence analysis suggests that the nonlinear three-wave coupling between the ZF and ambient turbulence(AT)is a possible mechanism for the ZF formation [10, 14].The direct measurements of nonlinear energy transfer among LFZFs, GAMs and turbulence show that the kinetic energy is transferred from turbulence into LFZFs and GAMs [17].The comparison of turbulence stress force and dynamic force of GAMs suggests that GAMs are driven by turbulence stresses dominantly[18].The role of the LFZFs during low to high confinement mode(L–H)transition and transport barrier formation is examined [19–21].Nevertheless, the toroidal components of the velocity for the ZFs have not been studied experimentally.

    Here, the toroidal component of the velocity for the GAMs is presented.The toroidal velocity is ~10–100 m s?1and increases with the poloidal velocity.The ratio of the toroidal velocity to the poloidal one is ~0.3–1.0.In the region of 3.5

    The structure of the remaining work is as below.Section 2 gives the experiment arrangement.Section 3 provides the experimental results, involving the properties of GAMs, toroidal component of the velocity for the GAMs,phase shifts between radial electric field and Mach number fluctuations in turbulence frequency band, statistical analysis of the velocities for the GAMs, nonlinear coupling and turbulence modulation etc.Sections 4 and 5 present the discussion and conclusion, respectively.

    2.Experimental setup

    The experiments were performed in Ohmic plasmas with limiter configuration on the J-TEXT tokamaks.The major radius isR=1.05 m and the minor radius isa=0.255 m for the J-TEXT tokamak [22].The toroidal magnetic field isBt=1.5–2.0 T and the plasma current isIp=150–180 kA.The line averaged electron density and the edge safety factor areNe=(1–3)×1019m?3andqa=3.2–4.3, respectively.The charging gas is hydrogen and the plasma shape is circular.In this experiment, a fast reciprocating probe array is used to measure the floating potentials, density, electron temperature and Mach number,etc[23].This probe array has two steps and 12 tips as shown in figure 1.It was set up on the top of the J-TEXT tokamak.Because of Shafranov displacement, the probe array is located at the position with a poloidal angle of θ~85°.On the other hand,the ratio of the toroidal component to the poloidal one of the velocity is expected as ~q?1cosθ for the GAM.These made it possible to measure the toroidal flows of the GAMs.To avoid the shadow,pins 2 and 4,and 6 and 8,should not be too close to their connecting lines with respect to magnetic field line;there is a small deviation of the angle.Assuming that they are shadowed, the measured density will be underestimated.For this experiment, the shadow is examined through changing the angle, and it does not occur.The diameter and length of the tips are 2 mm and 3 mm, respectively.The sampling rate for the digitizer is 2 MHz.

    Figures 2(a)–(d)describe the movement of the Langmuir probe array, floating potentials and ion saturated current measured by the Mach probe at the up and down sides,respectively.The Δrmeans the distance between the last closed flux surface (LCFS) and the measurement point.The negative sign of the Δrindicates inside the LCFS.The data were analyzed when the probes stayed inside the plasmas stably.

    Figure 1.Langmuir probe array,Vf (1, 3, 5, 7), ?V (2, 6), +V (4, 8), Idown (9, 11), Iup (10, 12).

    Figure 2.(a) Movement of the Langmuir probe arrays, (b) floating potentials, ion saturated currents at (c) up and (d) down sides.

    Figure 3.(a)The auto-power spectra of the potential fluctuations.(b)Phase shifts between potentials at two tips with poloidal spans of 6 mm(blue solid line)and with radial spans of 3 mm(red dashed line).

    3.Experimental results

    3.1.Properties of GAMs

    The auto-power spectrum of the potential fluctuations is described in figure 3(a).The resolution of the frequency is 0.5 kHz for this analysis.A low-frequency fluctuation is located in the frequency band of 0–4 kHz and the GAM is at the frequency off~15 kHz.The half width of the spectra is~5–8 kHz and thus the calculated lifetime is~120–200 μs for the GAMs.The phase shifts between potentials at two positions with spans ofθd=6 mm in the poloidal direction and ofdr=3 mm in the radial direction are given in figure 3(b).At the GAM frequency band,the phase shift is near zero,while at the turbulence frequency band, the linear dispersion relation is clearly demonstrated.The poloidal and radial wave vectors for the GAM are estimated askθ= Δφ1,3/dθ~0.3±0.1 cm?1andkr= Δφ3,5/dr~1.4±0.1 cm?1,respectively,where the phase shiftφΔ1,3is estimated fromδVf1andδVf3whileφΔ3,5fromδVf3andδVf5.TheδVfi(i=1, 2, …) refers to the folating potential fluctuations at thei-tip.In this paper,the error-bars are given by evaluating standard deviation.Considering the interference of turbulence,the poloidal wave vector is expected to be small.Although a poloidal component is involved in the measuredkrdue to the probe geometry, the estimatedkris still reliable for the GAMs due to its poloidal symmetry.For the turbulence, an accuratekrvalue can also be obtained with the formula k=kr+kθ, where k is the wave vector measured by the probe.This observation is in agreement with the theory prediction of GAMs,i.e,thekθis~0 and thekris finite for the GAMs.The peak in the frequency band of~30–100 kHz is ambient turbulence.

    The finite frequency is one of the particular features of GAMs.Figure 4 shows the relationship between the GAM frequency and the square root of the electron temperature measured by Langmuir probe array shot by shot.The measured GAM frequencyfGAMis proportional to the square root of the electron temperature.The electron temperature is detected by standard four tips.The GAM frequency is also calculated with the formulaassuming that the electron and ion temperature are identical in the edge plasmas.The scalingfGAM~1.3fthholds for the results, indicating that the measured frequencies are in good agreement with the prediction of the GAM theory.

    Figure 4.The measured GAM frequency (red square) and the calculated GAM frequency (green cycle).

    Figure 5.(a) The auto-power spectra of the density fluctuations,(b) the coherency spectra between the density and floating potential fluctuations, and (c) their phase shifts.

    The auto-power spectra of the density fluctuations are given in figure 5(a).The small peak is clearly shown at the frequency of ~15 kHz.The amplitudes of the density fluctuations at the GAM frequency and the average density are

    δn=5×1015m?3andn0=1×1018m?3, respectively.Theδn n0/is~0.5%.For the GAMs, the formula[24] should be satisfied.Here,qrandρsare the radial wave vector and the ion Larmor radius,respectively.In the top plasmas,θ=sin1,the value of the right-hand side is estimated as 0.7% and close toδn n0/ ~0.5%.

    The phase relationship between density fluctuations and potentials is also evaluated.Figures 5(b) and (c) present the coherency spectra and the phase differences between density fluctuations and potentials, respectively.The coherencies are significant above the noise level at the GAM frequency,indicating that the density fluctuations are well correlated with the potential fluctuations for the GAMs.The phase shifts at the GAM frequency are about 0.6π.

    Figure 6.(a)The auto-power spectra of(a)Er,(b)Mach number and(c) Eθ fluctuations, (d) the coherency spectra between Eθ and Mach number fluctuations (blue solid line), and between Er and Mach number fluctuations (red dashed line) (qa=4.1).

    3.2.Toroidal component of the velocity for the GAMs

    The toroidal component of the velocity for the GAMs deserves to be measured due to its relevance to the return flows owning to the plasma incompressibility.The Mach probe is parallel to the magnetic field line for the Mach number measurement.Due to the small angle between the parallel and toroidal directions,which is less than 5°, the toroidal component of the parallel velocity is approximately equivalent to the parallel velocity.Therefore,it is reasonable to estimate the toroidal velocity via the parallel Mach number.Figures 6(a)–(c) provide the auto-power spectra of the radial electric field fluctuations, Mach numbers,and poloidal electric field fluctuations atqa=4.1, respectively.The Mach numbers can be estimated as=Mmach0.4 ln (Iup/Idown)[25].From the measured Mach number and electron temperature, the toroidal velocity is evaluated asVφ≈MmachCs.The ion sound speed is written aswhereTe(mi) is the electron temperature (ion mass).A sharp peak at the GAM frequency occurs in the spectra of the radial electric feild fulctuationsδEr=(δVf5?δVf3)/dr.Small peaks at the GAM frequency are also observed in the Mach numberδM//and poloidal electric field fluctuationsδEθ=(δVf1?δVf3)/dθas shown in fgiures 6(b) and (c).The peak inδ θEis smaller than that inδErfor about an order of magnitude and may come from the fact that the two tips to measure the poloidal electric field are not placed on the same magnetic surface.Removing the turbulence interference, the poloidal velocity at the GAM frequency is calculated asVθGAM= ?ErGAM/Bt~60 m s?1.The folw velocity at the GAM frequency from the Mach number spectra is evaluated as

    Figure 7.The auto-power spectra of(a)Er,(b)Mach number and(c)Eθ fluctuations, (d) the coherency spectra between Eθ and Mach number fluctuations (blue solid line), and between Er and Mach number fluctuations (red dashed line) (qa=3.4).

    φVGAM~35 m s?1.The ratio of the toroidal component of velocities to poloidal one is~0.6.Note that for the GAMs,δEθ=0,thereby theδErdoes not contain a poloidal component of electric feild caused by the probe geometry.However, for turbulence,theδErcould be overestimated because theδEθis not equal to zero.

    The correlation between turbulence components in different directions is a mystery.Figure 6(d) shows the coherency spectra between the poloidal electric field fluctuations and Mach numbers.The high coherencies betweenδ θEandδM//in the turbulence frequency region of~30–100 kHz indicate thatδ θEandδM//have strong correlation for turbulence.In contrast, the coherencies between the radial electric field fluctuationsδErand Mach numberδM//in the turbulence frequency band are much lower than those ofδEθandδM//.This indicates that the correlation betweenδErandδM//is weak for turbulence.The strong correlation betweenδEθandδM//in the turbulence frequency band is required for the generation of the parallel stresses〈δV rδM//〉 .However, the coherence betweenδErandδM//is higher than that betweenδEθandδM//in the GAM frequency band.The possible reason is that the poloidal component of the velocity is dominant for the GAMs.

    The toroidal component of the velocity for the GAMs has also been detected atqa=3.4.Figures 7(a)–(c) show the auto-power spectra and the coherency spectra between fluctuations, respectively.Similarly, the sharp peak at the GAM frequency is observed in the spectra of theδErwhile small peaks appear in the spectra of theδE θandδM//.TheVθGAMis evaluated as~100 m s?1, respectively.The estimated toroidal velocityVφGAMat the GAM frequency is~80 m s?1from the Mach number spectra.The ratios of these velocities are calculated asVφGAM/VθGAM~0.8.The coherencies betweenδErandδM//are also lower (higher) than those ofδ θEandδM//in the turbulence (GAM) frequency band atqa=3.4.

    Figure 8.(a)The toroidal component of velocity as a function of poloidal velocity,and(b)the ratio of the toroidal velocity to the poloidal one for the GAMs at various edge safety factors; the dashed line denotes the theoretical prediction.

    Figure 9.Phase shifts between poloidal electric field and Mach number fluctuations at qa=3.4 (solid line) and 4.1 (dashed line).

    Note that the electron temperature fluctuations are generally considered to be small and can be neglected in the lowtemperature edge plasmas.In addition, for the measurement of the electron temperature fluctuations,the standard four tips are often used, while for the floating potential fluctuations,just one tip is needed.This indicates that the measurements for the electron temperature and floating potential fluctuations are not at the same position.Therefore, commonly, only floating potentials are utilized for the estimation of smallerscale electric fluctuations.

    3.3.Statistic analysis of the velocities for GAMs

    The poloidal and toroidal velocities for the GAMs are analyzed statistically near 20 shots.The toroidal component of the velocity as a function of the poloidal velocity for the GAMs is shown in figure 8(a).One can see that the toroidal velocity is less than the poloidal velocity in most cases and it tends to increase with the poloidal velocity.Here, except for the edge safety factors, other parameters are kept constant.The ratio of the toroidal component to the poloidal one of velocities for the GAMs as a function of the edge safety factor is given in figure 8(b).The ratio is mainly localized in the interval between 0.3 and 0.8.With a lower edge safety factor,the ratio tends to increase.The observed ratio ofVφGAM/VθGAM=0.3–1.0 is larger than the theoretical expectation ofθ qcos / =0.02–0.03.The possible reason is that the nonlinear effect is not taken into account for the equationVφGAM/VθGAM~q?1cosθ.The nonlinear coupling between the toroidal flow and turbulence will be discussed further.Note that the ratio rises rapidly when theqais reduced to 3.2, while in the interval betweenq=3.5 and 4.5, theVφGAM/VθGAMchanges slightly.The larger ratio at lowqmight be explained by the appearance of magnetic islands in the edge or the higher Ohmic heating power.Both islands and heating power could enhance turbulence and consequently drive flows.

    The observation is not consistent with the theoretical expectation [24], i.e.,VφGAM/VθGAM~qa?1cosθ.The approximate evaluation ofθwill be given later in the discussion section.The toroidal flows come from the return flow because of plasma incompressibility.In toroidal plasma, the GAMs directed mainly in the poloidal direction can induce the compression of plasma.However, to maintain incompressibility,a return flow along the magnetic field line will be generated to compensate for this compression.

    3.4.Cross phase between poloidal electric field and Mach number fluctuations

    The toroidal flow can be driven by turbulent nonlinear interaction via parallel Reynolds stress.The value of Reynolds stress〈δVrδM//〉 is influenced by the amplitudes and the cross correlation (or cross phase) of theδVrandδM//.The radial velocity fluctuation is written asδVr=δE θ/Bt.The stronger correlation between theδ θEandδM//is observed as above.Here,the cross phase of theδ θEandδM//is examined.Figure 9 gives the phase shifts betweenδ θEand Mach number atqa=3.4 and 4.1.One can expect that the cross phase is almost constant and about 0.6π in the turbulence frequency bands for two cases.The invariant phase between theδ θEandδM//is beneficial to the generation of the〈δVrδM//〉 .In the calculation for the cross phase,δ θEandδM//are selected at the same frequency.Therefore, the discussion corresponds to the toroidal flow, which has a zero frequency, but not for GAM because of its finite frequency.Generally, the momentum and energy conservations should be satisfied for the nonlinear coupling, i.e.,k=k1±k2andf=f1±f2,whereki(i=1, 2) refers to the wave vector.Assumingf1=f2,thenf=0.However, the mechanism is similar to that of the GAM generation.The nonlinear coupling between GAMs and turbulence will be given in next section.

    Figure 10.The auto-bicoherence spectrum of Mach numbers.

    3.5.Nonlinear coupling

    Nonlinear coupling is considered to be an important mechanism for the GAM generation and needs to be exploited further.Figure 10 gives the auto-bicoherence spectra of the Mach numbersδM//.The auto-bicoherences are calculated aswhere=B f3( )The values ofaboutf3=f1±|f2|~15 kHz,f2~±15 kHz andf1~15 kHz are higher.The cross-bicoherence spectra of the Mach numbersδM//and the turbulent velocityδVrare shown in figure 11(a),here,andTwolines aboutf3=f1±|f2|~0–10 kHz and ~15 kHz are apparently demonstrated.The linef3<10 kHz may come from the LFZFs.The summed bicoherence is given in figure 11(b)and calculated aswhereNis the total number ofThe noise level is expressed as 1/N0; hereN0is the realization number.This coupling analysis suggests that it is possible that the components at the GAMs frequency may come from the nonlinear interaction between turbulence and GAMs via three-wave coupling.It also has the implication of those flows being driven by parallel stresses.

    Figure 11.(a) The cross-bicoherence spectrum ofδVr andδM/ / ,and (b) summed bicoherence (the red dashed line is the noise level).

    Figure 12.(a)The auto-power spectrum of the turbulence envelopes in the frequency bands of ~50–200 kHz, (b) the coherency spectra between potential fluctuations and turbulence envelopes,and (c)the corresponding phase shifts.

    3.6.Turbulence modulation

    Figures 12(a)–(c) present the auto-power spectra of the turbulence envelope from the floating potential fluctuations, the coherency spectra between turbulence envelope and floating potentials, and their phase shifts, respectively.The envelope can be calculated with Hilbert transform from the higher frequency turbulence.The turbulence in the frequency band of~50–300 kHz is selected for this evaluation.A peak is clearly shown at the GAM frequency of~18 kHz in the power spectra.The high coherency at the GAM frequency band indicates that the GAM is well correlated with turbulence.This observation apparently shows that the turbulence is modulated by GAMs.The phase shifts between envelopes and potentials are~π at the GAM frequency.In addition,the coherencies in the frequency band of 0–4 kHz are also high and the phase shifts are~π.This may come from the fact that the LFZF modulates turbulence.The observation is consistent with the previous measurements in the middle plane of the HL-2A tokamak, where the turbulence envelope is well correlated with GAMs and the corresponding phase shifts are near π [14].The invariance of the phase shifts between turbulence envelope and GAMs indicates that the GAM effects on turbulence are the same in the middle plane and near the top magnetic surface of tokamaks.

    4.Discussion

    The toroidal component of velocity for the GAMs is estimated asVφGAM=10 ?100m s?1in the top of the J-TEXT tokamaks.This is not consistent with the theoretical expectation, where theVφGAM=0in the top plasmas of tokamaks [3].

    The most possible error resource is the movement of the magnetic axis for the present measurement.This can mean that the Mach probes and two tips to measure the poloidal electric field are not placed on the same magnetic surface[24,26].The error introduced by the misalignment is described by the ratiosinδsinθ,whereMint(MGAM)denotes the misalignment interference to Mach number (the Mach number for the GAM), and thelandδare the distance between two Mach probes and the mismatch angle, respectively [24].Considering Shafranov displacement, the plasma movement is ~2 cm and the poloidal angleθis estimated as~85°.This indicates that the probe is localized near the top of the magnetic surface.In this case, theφVGAMis expected asThe mismatch angle between the Mach probe array and magnetic surface due to the installment is less than 5°.Here, the observed toroidal velocity increases with the poloidal velocity.Meanwhile,the ratio0.3?1.0 is dominant and larger than the theoretical expectationcosθ~0.02 ?0.03.In the region of 3.5

    Note that here theTe=Tiis assumed.In most cases,theTishould be smaller than theTe.Therefore, the assumption may lead to the overestimation of the toroidal velocity.In addition, the calculationVφGAM=M//Cscan also enlarge the error bar.However, these two error sources are difficult to remove.

    5.Conclusion

    The toroidal component of the velocity for the GAMs is first demonstrated.Multiple Langmuir probe arrays set up near the top tokamak of the J-TEXT were used for this investigation.A significant peak at the GAM frequency is observed in Mach number fluctuations and well correlated with the GAMs.Similar to the previous measurements in a great number of devices,the GAM frequency is proportional to the square root of the electron temperature.The phase shift between densities and potentials at the GAM frequency is close to half π.Turbulence envelope analysis suggests that the GAM effects on turbulence are the same in the middle plane and near the top magnetic surface of tokamaks.The poloidal and toroidal velocities for the GAMs are estimated as~20–130 m s?1and~10–100 m s?1, respectively.The ratio of toroidal component to the poloidal one of the velocity is mainly 0.3–1.0.The toroidal velocity is also observed to increase with the poloidal velocity.The estimatedVφGAM/VθGAMalmost does not change in the interval between 3.5 and 4.5,while it rises rapidly at lowq.The coherencies between poloidal electric fields and Mach number fluctuations in turbulence frequency band are also evaluated and higher than those betweenδErandδM//.The high coherency betweenδE θandδM//is favorable for the generation of the turbulent Reynolds stress.

    Acknowledgments

    This work is supported by National Natural Science Foundation of China (Nos.12075057, 11775069, 11320101005,51821005 and 11875020);Jiangxi Provincial Natural Science Foundation (No.20202ACBL201002) and Doctoral Foundation (Nos.DHBK2017134 and DHBK 2018059); and Grant-in-Aid for Scientific Research of JSPS (Nos.15H02155, 15H02335, 16H02442).

    ORCID iDs

    猜你喜歡
    王璐龍文永華
    Relativistic Hartree–Fock model and its recent progress on the description of nuclear structure*
    Bandgap evolution of Mg3N2 under pressure:Experimental and theoretical studies
    Free-boundary plasma equilibria with toroidal plasma flows
    Non-Hermitian quasicrystal in dimerized lattices?
    How To Get Along With Your Friends Better
    Improved Fibroblast Adhesion and Proliferation by Controlling Multi-level Structure of Polycaprolactone Microfiber
    Club Recruitment
    A m,p-Laplacian Parabolic Equation with Nonlinear Absorption and Boundary Flux
    直徑不超過2的無爪圖的2—因子
    脾踩踏板有利于學(xué)習(xí)
    日本与韩国留学比较| av在线播放精品| 亚洲最大成人av| 久久韩国三级中文字幕| 午夜爱爱视频在线播放| videos熟女内射| 人人妻人人看人人澡| 成年女人在线观看亚洲视频 | 国产精品人妻久久久影院| 亚洲精品成人久久久久久| 中文字幕av成人在线电影| 日韩人妻高清精品专区| 三级国产精品欧美在线观看| 91久久精品电影网| 国产高潮美女av| 日韩一本色道免费dvd| 亚洲精品自拍成人| 2021天堂中文幕一二区在线观| 男女边摸边吃奶| 永久免费av网站大全| 22中文网久久字幕| 国产精品无大码| 国产成年人精品一区二区| 男女国产视频网站| 天美传媒精品一区二区| 欧美丝袜亚洲另类| 国产在线男女| 亚洲精品乱码久久久久久按摩| 美女黄网站色视频| 禁无遮挡网站| 国产精品麻豆人妻色哟哟久久 | 免费观看的影片在线观看| 久久久久久久大尺度免费视频| 看黄色毛片网站| 免费观看在线日韩| 色综合亚洲欧美另类图片| 亚洲精品久久午夜乱码| 美女高潮的动态| 啦啦啦韩国在线观看视频| 欧美人与善性xxx| 欧美激情久久久久久爽电影| 97人妻精品一区二区三区麻豆| 男女边吃奶边做爰视频| 99久国产av精品国产电影| 亚洲美女搞黄在线观看| 少妇高潮的动态图| 国产精品福利在线免费观看| 最近视频中文字幕2019在线8| 成人性生交大片免费视频hd| 丝袜喷水一区| 亚洲欧美一区二区三区黑人 | 尤物成人国产欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 日韩不卡一区二区三区视频在线| 日本免费a在线| 熟女人妻精品中文字幕| 成人欧美大片| 肉色欧美久久久久久久蜜桃 | 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久成人| 日韩欧美精品v在线| 激情五月婷婷亚洲| 非洲黑人性xxxx精品又粗又长| 99热这里只有是精品50| 成人毛片60女人毛片免费| 精品一区二区免费观看| 成人午夜高清在线视频| 国产91av在线免费观看| 亚洲成色77777| 国产精品av视频在线免费观看| 国产黄片视频在线免费观看| 亚洲av不卡在线观看| 麻豆成人av视频| 男女啪啪激烈高潮av片| 色哟哟·www| 97热精品久久久久久| 中文字幕制服av| 高清视频免费观看一区二区 | 九九爱精品视频在线观看| 伦精品一区二区三区| 欧美精品一区二区大全| 一级毛片电影观看| 麻豆乱淫一区二区| 国产探花极品一区二区| 深夜a级毛片| 亚洲av一区综合| 国产精品伦人一区二区| 国产成人一区二区在线| 欧美日韩在线观看h| 久久综合国产亚洲精品| 国产一级毛片在线| 精品酒店卫生间| 日韩精品有码人妻一区| 久久午夜福利片| 日本三级黄在线观看| 少妇高潮的动态图| 亚洲精品乱久久久久久| 亚洲欧美中文字幕日韩二区| 欧美日本视频| 欧美bdsm另类| 国产免费视频播放在线视频 | 草草在线视频免费看| 性插视频无遮挡在线免费观看| 亚洲美女搞黄在线观看| 国产av码专区亚洲av| 高清日韩中文字幕在线| 亚洲精品日本国产第一区| 国产探花在线观看一区二区| 国产黄色免费在线视频| 舔av片在线| 亚洲av一区综合| av福利片在线观看| 精品不卡国产一区二区三区| 国产老妇女一区| 边亲边吃奶的免费视频| 99九九线精品视频在线观看视频| 一区二区三区四区激情视频| 能在线免费观看的黄片| 欧美丝袜亚洲另类| 直男gayav资源| 波多野结衣巨乳人妻| 成人亚洲精品一区在线观看 | 亚洲欧洲日产国产| 深爱激情五月婷婷| 午夜福利视频1000在线观看| 国产成人aa在线观看| av女优亚洲男人天堂| 国产一区二区在线观看日韩| 亚洲av国产av综合av卡| 我的老师免费观看完整版| 亚洲精品自拍成人| 国产精品av视频在线免费观看| 麻豆国产97在线/欧美| or卡值多少钱| 嫩草影院入口| 51国产日韩欧美| 极品少妇高潮喷水抽搐| 亚洲国产日韩欧美精品在线观看| 熟女电影av网| 97人妻精品一区二区三区麻豆| 色综合站精品国产| 日本一二三区视频观看| 亚洲三级黄色毛片| 精品人妻视频免费看| 日日摸夜夜添夜夜爱| 亚洲精品乱码久久久v下载方式| 国产男女超爽视频在线观看| 久久久精品94久久精品| 久久久久久久国产电影| 我的老师免费观看完整版| a级毛片免费高清观看在线播放| 午夜激情福利司机影院| 国产男女超爽视频在线观看| 久久久精品94久久精品| 丝瓜视频免费看黄片| 九九在线视频观看精品| 亚洲国产欧美在线一区| 国产视频内射| 国产一区二区三区综合在线观看 | av线在线观看网站| 亚洲成色77777| 自拍偷自拍亚洲精品老妇| 免费大片18禁| 精品欧美国产一区二区三| 成人美女网站在线观看视频| 自拍偷自拍亚洲精品老妇| 丝袜美腿在线中文| 精品国内亚洲2022精品成人| 神马国产精品三级电影在线观看| 久久精品久久精品一区二区三区| 精品一区二区三卡| 国产一区二区在线观看日韩| 久久久久精品性色| 天天一区二区日本电影三级| 五月玫瑰六月丁香| 欧美 日韩 精品 国产| 黄色日韩在线| 少妇的逼好多水| 麻豆av噜噜一区二区三区| 九草在线视频观看| 夫妻午夜视频| 18禁动态无遮挡网站| 全区人妻精品视频| 舔av片在线| 晚上一个人看的免费电影| 国产精品1区2区在线观看.| 亚洲av.av天堂| 亚洲欧美清纯卡通| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 午夜福利视频1000在线观看| 精品少妇黑人巨大在线播放| 免费av毛片视频| 熟女人妻精品中文字幕| 国产老妇女一区| 精品一区在线观看国产| 中文在线观看免费www的网站| 国产黄色免费在线视频| 亚洲最大成人手机在线| 九草在线视频观看| 嫩草影院精品99| 国产久久久一区二区三区| 国产精品一区二区性色av| 免费大片18禁| 国产伦精品一区二区三区四那| 看非洲黑人一级黄片| 在现免费观看毛片| 黄色日韩在线| 嫩草影院新地址| 天堂网av新在线| 国产午夜精品久久久久久一区二区三区| 一区二区三区乱码不卡18| 国产在线男女| 国产精品一及| 欧美成人精品欧美一级黄| 精品熟女少妇av免费看| 一区二区三区四区激情视频| 免费看不卡的av| 亚洲伊人久久精品综合| 午夜激情欧美在线| 寂寞人妻少妇视频99o| 亚洲av一区综合| 国产色爽女视频免费观看| 亚洲自拍偷在线| 午夜福利视频1000在线观看| 网址你懂的国产日韩在线| 一级毛片电影观看| 国产免费又黄又爽又色| 婷婷六月久久综合丁香| 亚洲成人一二三区av| 肉色欧美久久久久久久蜜桃 | 免费看av在线观看网站| 两个人的视频大全免费| 日本av手机在线免费观看| 99热这里只有是精品50| 婷婷六月久久综合丁香| 国产精品久久久久久久久免| 亚洲精品aⅴ在线观看| 午夜福利视频精品| 不卡视频在线观看欧美| 淫秽高清视频在线观看| 免费少妇av软件| 亚洲欧美日韩无卡精品| 精品欧美国产一区二区三| 日韩视频在线欧美| 22中文网久久字幕| 国产伦精品一区二区三区四那| 精品国内亚洲2022精品成人| 亚洲最大成人中文| 国产有黄有色有爽视频| 久久久精品94久久精品| 国产免费视频播放在线视频 | 久久午夜福利片| 一本一本综合久久| 中文在线观看免费www的网站| 国产精品不卡视频一区二区| 伦理电影大哥的女人| 国产成人精品久久久久久| 亚洲欧美精品专区久久| 黑人高潮一二区| 一边亲一边摸免费视频| 国产综合精华液| 神马国产精品三级电影在线观看| 午夜精品国产一区二区电影 | 欧美日韩视频高清一区二区三区二| 一区二区三区免费毛片| 少妇熟女欧美另类| 日本三级黄在线观看| kizo精华| 在线a可以看的网站| 婷婷色麻豆天堂久久| 菩萨蛮人人尽说江南好唐韦庄| 日本猛色少妇xxxxx猛交久久| 18禁在线无遮挡免费观看视频| 免费在线观看成人毛片| .国产精品久久| 精品一区二区免费观看| 一区二区三区四区激情视频| 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 日本-黄色视频高清免费观看| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 天堂中文最新版在线下载 | 国产精品一区二区在线观看99 | 国产淫片久久久久久久久| 精品久久久噜噜| av免费观看日本| 国产黄片视频在线免费观看| 少妇的逼水好多| 哪个播放器可以免费观看大片| 在现免费观看毛片| 免费看光身美女| 色综合站精品国产| 伊人久久国产一区二区| 夜夜看夜夜爽夜夜摸| 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 人人妻人人澡欧美一区二区| 国产黄色小视频在线观看| 精品久久久久久成人av| 我的老师免费观看完整版| 一级毛片久久久久久久久女| 日韩成人伦理影院| 精品久久久久久电影网| 婷婷色av中文字幕| 国产免费福利视频在线观看| 卡戴珊不雅视频在线播放| 国产亚洲av片在线观看秒播厂 | 一级毛片 在线播放| 精品一区二区三区视频在线| av.在线天堂| 国产精品日韩av在线免费观看| 午夜激情福利司机影院| 亚洲精品日韩在线中文字幕| 18禁在线无遮挡免费观看视频| 伊人久久精品亚洲午夜| 免费不卡的大黄色大毛片视频在线观看 | 成年人午夜在线观看视频 | 亚洲成色77777| 国产成年人精品一区二区| 欧美xxxx性猛交bbbb| 国产午夜精品一二区理论片| 一边亲一边摸免费视频| 中文字幕制服av| 美女脱内裤让男人舔精品视频| 啦啦啦中文免费视频观看日本| 国产高清三级在线| 精品人妻视频免费看| a级毛片免费高清观看在线播放| 免费黄色在线免费观看| 国产精品一区二区在线观看99 | 日韩欧美一区视频在线观看 | 麻豆成人午夜福利视频| 国产在线男女| 亚洲最大成人av| 男人和女人高潮做爰伦理| 精品少妇黑人巨大在线播放| 日本免费a在线| 最近最新中文字幕大全电影3| 日韩制服骚丝袜av| 99热6这里只有精品| 1000部很黄的大片| 成人午夜高清在线视频| 黄色配什么色好看| 麻豆国产97在线/欧美| 欧美97在线视频| 自拍偷自拍亚洲精品老妇| 成人性生交大片免费视频hd| 欧美日韩综合久久久久久| 91在线精品国自产拍蜜月| 亚洲在线自拍视频| 亚洲精品,欧美精品| 一级黄片播放器| 国产精品.久久久| 免费观看无遮挡的男女| av专区在线播放| 2018国产大陆天天弄谢| 成年av动漫网址| 亚洲电影在线观看av| 一区二区三区乱码不卡18| 少妇的逼水好多| 亚洲伊人久久精品综合| 观看免费一级毛片| 蜜臀久久99精品久久宅男| 亚洲欧美日韩卡通动漫| 中文乱码字字幕精品一区二区三区 | 欧美一级a爱片免费观看看| 看免费成人av毛片| av网站免费在线观看视频 | 亚洲精品国产成人久久av| av女优亚洲男人天堂| 老女人水多毛片| 人体艺术视频欧美日本| 久久精品夜夜夜夜夜久久蜜豆| 一本久久精品| 久久精品夜色国产| 免费人成在线观看视频色| 少妇丰满av| 成人一区二区视频在线观看| 国产av码专区亚洲av| 高清日韩中文字幕在线| 免费观看无遮挡的男女| 国产精品久久久久久久久免| 黄色一级大片看看| av免费观看日本| 日韩成人av中文字幕在线观看| 日韩一区二区三区影片| 蜜桃亚洲精品一区二区三区| 久久久久九九精品影院| 秋霞伦理黄片| 亚洲精品亚洲一区二区| 国产视频首页在线观看| 国产黄色视频一区二区在线观看| 秋霞伦理黄片| 国产精品国产三级专区第一集| 免费观看的影片在线观看| 少妇的逼水好多| 免费大片黄手机在线观看| 黄片wwwwww| 在线免费观看不下载黄p国产| 亚洲无线观看免费| 日韩一本色道免费dvd| 日韩精品有码人妻一区| 大片免费播放器 马上看| 最近手机中文字幕大全| 国产色爽女视频免费观看| 成人亚洲精品一区在线观看 | 亚洲av在线观看美女高潮| 久久久久精品性色| 久久久久久久亚洲中文字幕| 熟女电影av网| 成人国产麻豆网| 亚洲国产最新在线播放| 成人高潮视频无遮挡免费网站| 国产真实伦视频高清在线观看| 亚洲国产欧美人成| 欧美 日韩 精品 国产| 亚洲丝袜综合中文字幕| 国产女主播在线喷水免费视频网站 | 又大又黄又爽视频免费| 我要看日韩黄色一级片| 晚上一个人看的免费电影| 亚洲激情五月婷婷啪啪| 色吧在线观看| 非洲黑人性xxxx精品又粗又长| 身体一侧抽搐| 国产成人精品婷婷| videossex国产| 成年女人在线观看亚洲视频 | 中文字幕免费在线视频6| 内地一区二区视频在线| 九九在线视频观看精品| 国产精品一区二区性色av| 国产精品av视频在线免费观看| 亚洲精品成人av观看孕妇| 久久亚洲国产成人精品v| 日本色播在线视频| 一级二级三级毛片免费看| 丝袜喷水一区| 午夜福利在线在线| 丰满乱子伦码专区| 久久精品久久久久久久性| 国产探花极品一区二区| 欧美日韩在线观看h| 天堂√8在线中文| 大香蕉97超碰在线| 国产v大片淫在线免费观看| 伦理电影大哥的女人| 亚洲欧美精品专区久久| 91久久精品国产一区二区成人| 成人二区视频| 欧美最新免费一区二区三区| 看非洲黑人一级黄片| 嫩草影院新地址| 亚洲av电影不卡..在线观看| 在线播放无遮挡| 人妻夜夜爽99麻豆av| 中文精品一卡2卡3卡4更新| av女优亚洲男人天堂| 久久这里有精品视频免费| 91久久精品电影网| 亚洲精品乱久久久久久| 亚洲国产精品成人久久小说| 久久久久久久亚洲中文字幕| 汤姆久久久久久久影院中文字幕 | 欧美另类一区| 久久久久性生活片| 国产熟女欧美一区二区| 午夜激情久久久久久久| 日韩一区二区三区影片| 天天一区二区日本电影三级| 成人国产麻豆网| 欧美日韩精品成人综合77777| 免费播放大片免费观看视频在线观看| 一边亲一边摸免费视频| 男人爽女人下面视频在线观看| 大又大粗又爽又黄少妇毛片口| av在线蜜桃| 搡老乐熟女国产| 国产黄色视频一区二区在线观看| 国内精品美女久久久久久| 我的女老师完整版在线观看| 国产免费视频播放在线视频 | 免费观看精品视频网站| 国产一级毛片七仙女欲春2| 日韩成人av中文字幕在线观看| 啦啦啦啦在线视频资源| 人妻制服诱惑在线中文字幕| 2021天堂中文幕一二区在线观| 少妇的逼好多水| 午夜激情福利司机影院| 网址你懂的国产日韩在线| 美女xxoo啪啪120秒动态图| 亚洲久久久久久中文字幕| 国产熟女欧美一区二区| 一个人观看的视频www高清免费观看| 国产在视频线在精品| 22中文网久久字幕| 一区二区三区乱码不卡18| 成人一区二区视频在线观看| 熟女电影av网| 久久久久久久午夜电影| 中文在线观看免费www的网站| 我的老师免费观看完整版| 国产精品日韩av在线免费观看| 一本一本综合久久| 看免费成人av毛片| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频 | 日日摸夜夜添夜夜添av毛片| 高清日韩中文字幕在线| 狂野欧美白嫩少妇大欣赏| h日本视频在线播放| 国产视频首页在线观看| kizo精华| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 搡女人真爽免费视频火全软件| 成人欧美大片| 97在线视频观看| 免费观看精品视频网站| 久久6这里有精品| 伊人久久精品亚洲午夜| 好男人在线观看高清免费视频| 三级国产精品欧美在线观看| 高清欧美精品videossex| 欧美日韩一区二区视频在线观看视频在线 | 久久精品夜夜夜夜夜久久蜜豆| 久久精品熟女亚洲av麻豆精品 | 偷拍熟女少妇极品色| 三级男女做爰猛烈吃奶摸视频| 一区二区三区免费毛片| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 中文天堂在线官网| 精品久久久久久成人av| 最新中文字幕久久久久| 赤兔流量卡办理| 97在线视频观看| 午夜免费男女啪啪视频观看| videossex国产| 精品不卡国产一区二区三区| a级毛色黄片| 久久99精品国语久久久| 97精品久久久久久久久久精品| 日本wwww免费看| 日日啪夜夜爽| 欧美性感艳星| 伦精品一区二区三区| 2021少妇久久久久久久久久久| 欧美xxⅹ黑人| 99久国产av精品国产电影| av免费观看日本| 狠狠精品人妻久久久久久综合| 国产一区二区亚洲精品在线观看| 久久这里有精品视频免费| 国产av不卡久久| 看免费成人av毛片| 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 成人av在线播放网站| 91午夜精品亚洲一区二区三区| 国产乱人偷精品视频| 七月丁香在线播放| 精品亚洲乱码少妇综合久久| 欧美性感艳星| 成人高潮视频无遮挡免费网站| 简卡轻食公司| 久久热精品热| av黄色大香蕉| 免费不卡的大黄色大毛片视频在线观看 | 亚洲人成网站在线观看播放| 亚洲高清免费不卡视频| 欧美日韩综合久久久久久| 一二三四中文在线观看免费高清| 国内精品一区二区在线观看| 亚洲av免费在线观看| 亚洲国产精品成人久久小说| 午夜福利高清视频| 三级毛片av免费| www.色视频.com| 国产成人精品一,二区| 啦啦啦中文免费视频观看日本| 国产毛片a区久久久久| 18禁在线播放成人免费| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 免费看光身美女| 日日干狠狠操夜夜爽| 男女边摸边吃奶| 大香蕉97超碰在线| 国内揄拍国产精品人妻在线| 亚洲久久久久久中文字幕| 搡老乐熟女国产| 人妻制服诱惑在线中文字幕| 午夜激情福利司机影院| 亚洲乱码一区二区免费版| 啦啦啦韩国在线观看视频| 搡老乐熟女国产| 高清视频免费观看一区二区 | 国产午夜精品一二区理论片| 99九九线精品视频在线观看视频| 亚洲第一区二区三区不卡| 亚洲av免费高清在线观看| 边亲边吃奶的免费视频| 又爽又黄a免费视频| 国产精品人妻久久久影院| 三级经典国产精品| 精品久久久精品久久久| 中文字幕av在线有码专区| 日韩三级伦理在线观看| 国产在线男女|