• <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一区| 黄片小视频在线播放| 欧美老熟妇乱子伦牲交| 99精国产麻豆久久婷婷| 一区二区三区精品91| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美一区二区综合| 久久99精品国语久久久| 少妇的丰满在线观看| 2018国产大陆天天弄谢| 久久精品久久精品一区二区三区| 免费日韩欧美在线观看| 纯流量卡能插随身wifi吗| 在线观看人妻少妇| 国产色婷婷99| 欧美老熟妇乱子伦牲交| 国产免费福利视频在线观看| 男的添女的下面高潮视频| 黄片无遮挡物在线观看| 大码成人一级视频| 久久国产亚洲av麻豆专区| 国产精品嫩草影院av在线观看| 亚洲熟女毛片儿| 777久久人妻少妇嫩草av网站| 91精品伊人久久大香线蕉| 久久国产亚洲av麻豆专区| 日日爽夜夜爽网站| 9色porny在线观看| 国产成人91sexporn| 人成视频在线观看免费观看| 午夜av观看不卡| 成人国产av品久久久| 色综合欧美亚洲国产小说| 久久久久精品国产欧美久久久 | 岛国毛片在线播放| 国产av码专区亚洲av| 亚洲激情五月婷婷啪啪| 日韩中文字幕视频在线看片| 午夜福利视频精品| 少妇 在线观看| 自线自在国产av| 久久久国产欧美日韩av| 80岁老熟妇乱子伦牲交| 香蕉丝袜av| 亚洲色图 男人天堂 中文字幕| xxx大片免费视频| 成人18禁高潮啪啪吃奶动态图| 丝袜美足系列| 国产精品久久久久久精品古装| 美女午夜性视频免费| 91精品三级在线观看| 男女免费视频国产| 桃花免费在线播放| 美女福利国产在线| 亚洲av综合色区一区| www.av在线官网国产| 国产免费现黄频在线看| 黄片无遮挡物在线观看| www.熟女人妻精品国产| 一本大道久久a久久精品| 免费日韩欧美在线观看| 国产野战对白在线观看| 黄色一级大片看看| 只有这里有精品99| 大片电影免费在线观看免费| 少妇人妻精品综合一区二区| 午夜免费观看性视频| 永久免费av网站大全| 看非洲黑人一级黄片| av网站在线播放免费| 免费日韩欧美在线观看| 极品少妇高潮喷水抽搐| 老司机影院毛片| 亚洲美女视频黄频| 久久国产亚洲av麻豆专区| 国产日韩欧美视频二区| 午夜福利,免费看| 午夜福利网站1000一区二区三区| 搡老乐熟女国产| 我要看黄色一级片免费的| 亚洲成人免费av在线播放| 国产精品二区激情视频| 这个男人来自地球电影免费观看 | 亚洲精品国产av蜜桃| avwww免费| 国产免费现黄频在线看| 日韩精品有码人妻一区| 色视频在线一区二区三区| 丝袜脚勾引网站| 男女午夜视频在线观看| 久久女婷五月综合色啪小说| 国产激情久久老熟女| 亚洲欧美一区二区三区久久| 韩国精品一区二区三区| 国产一区二区激情短视频 | 波多野结衣av一区二区av| 亚洲精品美女久久久久99蜜臀 | 一级毛片 在线播放| 精品一区二区三区av网在线观看 | 在线观看www视频免费| 欧美激情 高清一区二区三区| 久久毛片免费看一区二区三区| 国产一区二区激情短视频 | 国产精品免费视频内射| 1024香蕉在线观看| 这个男人来自地球电影免费观看 | 狂野欧美激情性bbbbbb| 男女午夜视频在线观看| 国产精品av久久久久免费| 99香蕉大伊视频| 亚洲国产精品一区三区| 亚洲av综合色区一区| 十八禁网站网址无遮挡| 一区在线观看完整版| 色视频在线一区二区三区| 精品少妇一区二区三区视频日本电影 | 久久久久网色| 久热这里只有精品99| 成年美女黄网站色视频大全免费| 日本wwww免费看| 黄色怎么调成土黄色| 国产精品国产三级国产专区5o| 在线观看免费视频网站a站| 波野结衣二区三区在线| 成人三级做爰电影| 9热在线视频观看99| 免费看不卡的av| 伊人亚洲综合成人网| 两个人看的免费小视频| 如何舔出高潮| 免费观看av网站的网址| 国产男人的电影天堂91| 亚洲精品在线美女| 九草在线视频观看| 国产一区二区在线观看av| 少妇猛男粗大的猛烈进出视频| 欧美成人精品欧美一级黄| 国产不卡av网站在线观看| 成人黄色视频免费在线看| 国产成人av激情在线播放| 精品国产超薄肉色丝袜足j| a级片在线免费高清观看视频| 大香蕉久久成人网| 天美传媒精品一区二区| 免费看av在线观看网站| 欧美激情极品国产一区二区三区| 操出白浆在线播放| 色播在线永久视频| 精品亚洲成a人片在线观看| 中文字幕人妻丝袜一区二区 | 天天躁日日躁夜夜躁夜夜| 久久精品国产a三级三级三级| 亚洲av日韩精品久久久久久密 | 久久毛片免费看一区二区三区| 制服丝袜香蕉在线| 久久久久国产一级毛片高清牌| 亚洲欧美清纯卡通| 日韩视频在线欧美| 亚洲av日韩在线播放| 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说| 美女中出高潮动态图| 成人亚洲欧美一区二区av| 51午夜福利影视在线观看| 91国产中文字幕| 精品一区二区三区四区五区乱码 | 亚洲七黄色美女视频| 色网站视频免费| 亚洲av电影在线进入| 制服丝袜香蕉在线| 成人毛片60女人毛片免费| 亚洲精品久久久久久婷婷小说| bbb黄色大片| 青青草视频在线视频观看| 热99国产精品久久久久久7| 麻豆精品久久久久久蜜桃| 国产亚洲av片在线观看秒播厂| 国产成人午夜福利电影在线观看| 欧美激情 高清一区二区三区| 精品人妻熟女毛片av久久网站| 日本黄色日本黄色录像| 欧美日韩亚洲综合一区二区三区_| 免费观看av网站的网址| 成人黄色视频免费在线看| 中国三级夫妇交换| 精品一区在线观看国产| 亚洲国产欧美网| 亚洲欧美中文字幕日韩二区| 另类精品久久| 日本欧美视频一区| 亚洲成人免费av在线播放| 99re6热这里在线精品视频| 久久精品国产亚洲av高清一级| 在线观看免费午夜福利视频| 亚洲成国产人片在线观看| 黄色怎么调成土黄色| 好男人视频免费观看在线| 亚洲精品日本国产第一区| 日韩 亚洲 欧美在线| 免费看av在线观看网站| 人人妻,人人澡人人爽秒播 | 99九九在线精品视频| 久久久久精品性色| 天堂俺去俺来也www色官网| 成人午夜精彩视频在线观看| 亚洲中文av在线| 欧美精品人与动牲交sv欧美| av免费观看日本| 各种免费的搞黄视频| 色视频在线一区二区三区| 亚洲精品第二区| 国产色婷婷99| 老司机影院毛片| 国产一区二区激情短视频 | 亚洲色图综合在线观看| 国产老妇伦熟女老妇高清| 亚洲av成人不卡在线观看播放网 | 一区在线观看完整版| 在线观看三级黄色| 伦理电影免费视频| 97精品久久久久久久久久精品| 国产欧美日韩一区二区三区在线| 欧美97在线视频| 国产精品.久久久| 精品国产一区二区三区四区第35| 少妇人妻精品综合一区二区| 午夜精品国产一区二区电影| 极品少妇高潮喷水抽搐| 在线精品无人区一区二区三| 国产成人免费无遮挡视频| 一二三四在线观看免费中文在| 久久亚洲国产成人精品v| tube8黄色片| 欧美最新免费一区二区三区| 日韩av在线免费看完整版不卡| 日韩制服骚丝袜av| 9色porny在线观看| 国产无遮挡羞羞视频在线观看| 久久久精品免费免费高清| 成人亚洲欧美一区二区av| 亚洲熟女精品中文字幕| 国产精品欧美亚洲77777| 一边摸一边抽搐一进一出视频| 国产日韩欧美亚洲二区| 久久精品熟女亚洲av麻豆精品| 久久国产精品大桥未久av| 波多野结衣一区麻豆| 日本av手机在线免费观看| 国产日韩欧美视频二区| 美女午夜性视频免费| 黄片小视频在线播放| 中文天堂在线官网| 成人国产av品久久久| 亚洲五月色婷婷综合| 久久久国产欧美日韩av| 啦啦啦 在线观看视频| 成人毛片60女人毛片免费| 成人手机av| 久久ye,这里只有精品| 91精品伊人久久大香线蕉| 免费黄色在线免费观看| 午夜福利视频在线观看免费| 亚洲在久久综合| 欧美精品一区二区免费开放| 日韩av不卡免费在线播放| 日本一区二区免费在线视频| 少妇人妻 视频| 午夜av观看不卡| 国产激情久久老熟女| 在线观看免费高清a一片| 一级a爱视频在线免费观看| 免费黄频网站在线观看国产| 色视频在线一区二区三区| 下体分泌物呈黄色| 欧美精品高潮呻吟av久久| 精品人妻熟女毛片av久久网站| 亚洲欧美精品自产自拍| 亚洲欧洲国产日韩| 中文字幕av电影在线播放| 青春草国产在线视频| 一本一本久久a久久精品综合妖精| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 叶爱在线成人免费视频播放| 91老司机精品| 国产精品香港三级国产av潘金莲 | av网站在线播放免费| 久久久精品区二区三区| 免费观看人在逋| 色视频在线一区二区三区| 国产野战对白在线观看| 少妇人妻久久综合中文| 热re99久久精品国产66热6| 人人妻人人澡人人爽人人夜夜| 国产高清国产精品国产三级| 久久国产精品男人的天堂亚洲| 欧美精品一区二区免费开放| 婷婷色综合www| 在线精品无人区一区二区三| 一级毛片电影观看| 韩国精品一区二区三区| 欧美日韩视频精品一区| 亚洲欧美一区二区三区久久| 国产成人精品福利久久| 亚洲欧洲精品一区二区精品久久久 | 超碰97精品在线观看| 久久国产精品大桥未久av| 男女国产视频网站| 欧美av亚洲av综合av国产av | 综合色丁香网| 一级爰片在线观看| 视频区图区小说| 久热爱精品视频在线9| 老司机影院毛片| 亚洲av电影在线观看一区二区三区| 午夜精品国产一区二区电影| 这个男人来自地球电影免费观看 | 国产探花极品一区二区| 久久久久视频综合| 婷婷成人精品国产| 成人毛片60女人毛片免费| 色视频在线一区二区三区| 日韩视频在线欧美| 少妇人妻久久综合中文| 久久久久久免费高清国产稀缺| 一级毛片黄色毛片免费观看视频| 国精品久久久久久国模美| 久久婷婷青草| 各种免费的搞黄视频| 亚洲精品国产av蜜桃| 黄片无遮挡物在线观看| 国产精品99久久99久久久不卡 | 美国免费a级毛片| 欧美黑人欧美精品刺激| 国产片内射在线| 男人添女人高潮全过程视频| 色婷婷av一区二区三区视频| 久久久久久久大尺度免费视频| 国产一区二区三区av在线| 99久国产av精品国产电影| 久久国产精品大桥未久av| 黑人欧美特级aaaaaa片| 亚洲精品美女久久久久99蜜臀 | 久久久久国产一级毛片高清牌| 婷婷色av中文字幕| 成人黄色视频免费在线看| 日本av手机在线免费观看| 久久久久久久国产电影| 亚洲国产精品999| 成人亚洲精品一区在线观看| av免费观看日本| 丰满少妇做爰视频| 国产探花极品一区二区| 免费av中文字幕在线| 亚洲欧美精品综合一区二区三区| 在线天堂最新版资源| 国产激情久久老熟女| 亚洲精品一区蜜桃| 欧美日韩视频精品一区| 777久久人妻少妇嫩草av网站| 久久精品人人爽人人爽视色| 国产精品女同一区二区软件| 一区二区日韩欧美中文字幕| 狂野欧美激情性xxxx| 中文字幕亚洲精品专区| 两性夫妻黄色片| 国产精品一区二区精品视频观看| 亚洲欧美精品自产自拍| 精品久久久精品久久久| 欧美老熟妇乱子伦牲交| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 多毛熟女@视频| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人看| 激情视频va一区二区三区| 久久99一区二区三区| 热99国产精品久久久久久7| 国产高清不卡午夜福利| 精品卡一卡二卡四卡免费| 国产一区二区三区综合在线观看| √禁漫天堂资源中文www| 国产野战对白在线观看| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 亚洲色图 男人天堂 中文字幕| 美女视频免费永久观看网站| 欧美精品一区二区大全| 国产免费一区二区三区四区乱码| 国产一区有黄有色的免费视频| 在线观看免费午夜福利视频| 晚上一个人看的免费电影| 婷婷色麻豆天堂久久| 亚洲一区中文字幕在线| 波野结衣二区三区在线| 久久青草综合色| 欧美97在线视频| 国产精品成人在线| 亚洲精品中文字幕在线视频| 最黄视频免费看| 亚洲五月色婷婷综合| 精品国产乱码久久久久久男人| 夜夜骑夜夜射夜夜干| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 婷婷色麻豆天堂久久| 亚洲精品在线美女| 亚洲欧美成人综合另类久久久| 高清在线视频一区二区三区| 一级片'在线观看视频| 一区福利在线观看| 黄片播放在线免费| 久久久国产精品麻豆| 国语对白做爰xxxⅹ性视频网站| 免费女性裸体啪啪无遮挡网站| 美女福利国产在线| 亚洲成色77777| 黄色怎么调成土黄色| 久热这里只有精品99| 一本—道久久a久久精品蜜桃钙片| 如何舔出高潮| 亚洲天堂av无毛| 欧美成人午夜精品| av在线观看视频网站免费| 成年人午夜在线观看视频| 国产熟女欧美一区二区| 亚洲婷婷狠狠爱综合网| 欧美xxⅹ黑人| 不卡av一区二区三区| 国产成人精品久久久久久| 一本大道久久a久久精品| 国产成人一区二区在线| 9色porny在线观看| 亚洲国产精品999| 观看美女的网站| 肉色欧美久久久久久久蜜桃| 国产男人的电影天堂91| 国产精品蜜桃在线观看| 国产精品偷伦视频观看了| 午夜激情av网站| 天堂中文最新版在线下载| 99九九在线精品视频| 午夜福利免费观看在线| 9色porny在线观看| 亚洲精华国产精华液的使用体验| 各种免费的搞黄视频| 国产深夜福利视频在线观看| 亚洲精品,欧美精品| 亚洲一区中文字幕在线| 男女下面插进去视频免费观看| 亚洲成人手机| 国产精品一区二区在线观看99| 久久韩国三级中文字幕| 国产在线视频一区二区| 国产成人av激情在线播放| 久久亚洲国产成人精品v| 国产激情久久老熟女| 国产日韩欧美在线精品| 大片电影免费在线观看免费| 秋霞伦理黄片| 美女扒开内裤让男人捅视频| 国产又爽黄色视频| 在线观看国产h片| 免费看不卡的av| svipshipincom国产片| 欧美老熟妇乱子伦牲交| 我的亚洲天堂| 国产av精品麻豆| 亚洲欧美一区二区三区久久| 成人亚洲精品一区在线观看| 国产精品熟女久久久久浪| 在线看a的网站| 色婷婷久久久亚洲欧美| 亚洲精品,欧美精品| 精品国产超薄肉色丝袜足j| 在线观看免费日韩欧美大片| 亚洲av成人不卡在线观看播放网 | 一级爰片在线观看| 久久热在线av| 99国产精品免费福利视频| 亚洲国产最新在线播放| 制服丝袜香蕉在线| 老司机深夜福利视频在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 高清av免费在线| 欧美黑人精品巨大| 久久久久精品国产欧美久久久 | 大码成人一级视频| 午夜影院在线不卡| www.精华液| 亚洲人成网站在线观看播放| 国产精品一区二区在线不卡| 久久av网站| 黄色怎么调成土黄色| 不卡视频在线观看欧美| 我的亚洲天堂| 国产成人精品福利久久| 黄色毛片三级朝国网站| 亚洲第一区二区三区不卡| 热99国产精品久久久久久7| 国产女主播在线喷水免费视频网站| 欧美精品亚洲一区二区| 日韩av在线免费看完整版不卡| 精品少妇一区二区三区视频日本电影 | 亚洲第一区二区三区不卡| 国产精品国产三级国产专区5o| 亚洲色图 男人天堂 中文字幕| 老司机靠b影院| 成人手机av| 国产亚洲一区二区精品| 免费在线观看黄色视频的| 亚洲国产精品国产精品| 国产成人欧美在线观看 | 亚洲综合精品二区| 精品福利永久在线观看| 麻豆精品久久久久久蜜桃| 一区二区av电影网| 一级爰片在线观看| 国产精品三级大全| av有码第一页| 久久精品国产综合久久久| 人人妻人人澡人人爽人人夜夜| 欧美人与性动交α欧美软件| 亚洲七黄色美女视频| 色婷婷久久久亚洲欧美| 美女主播在线视频| 熟妇人妻不卡中文字幕| 成人手机av| 国产欧美日韩一区二区三区在线| 七月丁香在线播放| 一级爰片在线观看| 日本av手机在线免费观看| 90打野战视频偷拍视频| 婷婷色综合www| 中文欧美无线码| 纯流量卡能插随身wifi吗| 精品一区二区三区四区五区乱码 | 黄色视频在线播放观看不卡| 欧美精品一区二区免费开放| 狠狠精品人妻久久久久久综合| 女人精品久久久久毛片| 亚洲精品一区蜜桃| 久久午夜综合久久蜜桃| 美女高潮到喷水免费观看| 91精品伊人久久大香线蕉| 制服丝袜香蕉在线| 国产老妇伦熟女老妇高清| 波野结衣二区三区在线| 中国国产av一级| 考比视频在线观看| 精品一区二区免费观看| 少妇 在线观看| 免费在线观看完整版高清| 欧美中文综合在线视频| 无遮挡黄片免费观看| 亚洲av成人精品一二三区| 日韩 亚洲 欧美在线| 国产成人精品久久久久久| 欧美日韩国产mv在线观看视频| 国产在线免费精品| 成年女人毛片免费观看观看9 | av网站免费在线观看视频| 久久久久精品性色| 视频在线观看一区二区三区| 久久精品久久久久久久性| 亚洲国产精品一区二区三区在线| 欧美日韩av久久| 国产女主播在线喷水免费视频网站| 国产在视频线精品| 国产av精品麻豆| 美女福利国产在线| 精品人妻熟女毛片av久久网站| 国产精品三级大全| av在线app专区| av福利片在线| 如何舔出高潮| 这个男人来自地球电影免费观看 | 亚洲,欧美,日韩| av天堂久久9| 久久久精品免费免费高清| 2018国产大陆天天弄谢| 黄色 视频免费看| 男的添女的下面高潮视频| 亚洲自偷自拍图片 自拍| 深夜精品福利| 久久久精品94久久精品| 男人添女人高潮全过程视频| 国产97色在线日韩免费| 亚洲av在线观看美女高潮| 人人妻人人澡人人看| 午夜精品国产一区二区电影| 超色免费av| 丰满饥渴人妻一区二区三| 汤姆久久久久久久影院中文字幕| 日韩 欧美 亚洲 中文字幕| 大码成人一级视频| 亚洲第一av免费看| 亚洲精品乱久久久久久| 国产精品久久久久久久久免| 两性夫妻黄色片| 国产乱人偷精品视频| 另类精品久久| 日本av免费视频播放| 91成人精品电影| 国产亚洲最大av| 看十八女毛片水多多多| 成年人午夜在线观看视频| 两性夫妻黄色片| av在线app专区| 激情五月婷婷亚洲| 久久久久久久久久久免费av| 曰老女人黄片| 欧美日韩一级在线毛片| 久久精品熟女亚洲av麻豆精品|