• <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í)
    精品国产三级普通话版| 精品国内亚洲2022精品成人| 国产精品国产高清国产av| 精品人妻偷拍中文字幕| 我的女老师完整版在线观看| 岛国在线免费视频观看| 国产精品美女特级片免费视频播放器| 十八禁国产超污无遮挡网站| 亚洲色图av天堂| 亚洲精品自拍成人| 波多野结衣高清作品| 91久久精品国产一区二区成人| 一进一出抽搐动态| 日韩视频在线欧美| 亚洲中文字幕日韩| 国产真实伦视频高清在线观看| 国产精品三级大全| 午夜福利视频1000在线观看| 在线播放国产精品三级| 99久久九九国产精品国产免费| 欧美日韩在线观看h| 成年女人永久免费观看视频| 国产精品嫩草影院av在线观看| 又爽又黄a免费视频| 一个人看的www免费观看视频| 精品久久久久久久久亚洲| 青青草视频在线视频观看| 午夜激情欧美在线| 国产av麻豆久久久久久久| 色哟哟哟哟哟哟| 爱豆传媒免费全集在线观看| 日韩欧美精品免费久久| 亚洲国产精品久久男人天堂| 国产精品.久久久| 少妇人妻一区二区三区视频| 听说在线观看完整版免费高清| 国产美女午夜福利| 国产欧美日韩精品一区二区| 12—13女人毛片做爰片一| 亚洲国产欧美在线一区| 99久久精品国产国产毛片| 一本精品99久久精品77| 小说图片视频综合网站| 日韩欧美一区二区三区在线观看| 色播亚洲综合网| 晚上一个人看的免费电影| 精品少妇黑人巨大在线播放 | eeuss影院久久| 国产精品女同一区二区软件| 久久人人爽人人片av| 亚洲av中文字字幕乱码综合| 又粗又爽又猛毛片免费看| 变态另类丝袜制服| 99在线视频只有这里精品首页| 国产成人精品一,二区 | 韩国av在线不卡| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 国产av不卡久久| 好男人在线观看高清免费视频| 天堂网av新在线| 国产极品精品免费视频能看的| 在线免费观看的www视频| 亚洲va在线va天堂va国产| 中文字幕人妻熟人妻熟丝袜美| 日本成人三级电影网站| 国产精品久久久久久精品电影| 欧美成人a在线观看| 日本在线视频免费播放| 大香蕉久久网| 男女下面进入的视频免费午夜| 欧美激情在线99| 成年av动漫网址| 三级经典国产精品| 亚洲欧洲国产日韩| videossex国产| 国产亚洲av片在线观看秒播厂 | 观看免费一级毛片| 亚洲第一区二区三区不卡| 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 欧美最黄视频在线播放免费| 国内少妇人妻偷人精品xxx网站| 国产精品三级大全| 久久久久性生活片| 中文字幕av在线有码专区| 久久久午夜欧美精品| 狂野欧美激情性xxxx在线观看| 韩国av在线不卡| 色5月婷婷丁香| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 蜜桃久久精品国产亚洲av| 亚洲av二区三区四区| 久久久久久伊人网av| 午夜免费激情av| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 国产一区二区激情短视频| 午夜激情福利司机影院| 亚洲精华国产精华液的使用体验 | 直男gayav资源| 欧美日韩综合久久久久久| 我的女老师完整版在线观看| 精品一区二区三区视频在线| 国产av麻豆久久久久久久| 老司机影院成人| 欧美潮喷喷水| 日产精品乱码卡一卡2卡三| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 国产高清不卡午夜福利| 爱豆传媒免费全集在线观看| 国产人妻一区二区三区在| 少妇的逼水好多| 日本与韩国留学比较| 69人妻影院| 午夜福利成人在线免费观看| 在线a可以看的网站| 我的女老师完整版在线观看| 精品一区二区三区视频在线| 97超碰精品成人国产| 久久99蜜桃精品久久| 狠狠狠狠99中文字幕| 久久久久久久久久久丰满| 欧美日韩一区二区视频在线观看视频在线 | 国产成人a∨麻豆精品| 亚洲天堂国产精品一区在线| 精品一区二区三区人妻视频| 久久久久久大精品| 久久草成人影院| 婷婷精品国产亚洲av| 99久久中文字幕三级久久日本| 精品欧美国产一区二区三| 免费人成在线观看视频色| 日韩人妻高清精品专区| 观看美女的网站| 午夜精品国产一区二区电影 | 日日撸夜夜添| 亚洲国产精品成人综合色| 大型黄色视频在线免费观看| 亚洲乱码一区二区免费版| 亚洲综合色惰| 亚洲av一区综合| 人妻少妇偷人精品九色| 亚洲第一电影网av| 精品少妇黑人巨大在线播放 | 男人舔奶头视频| 欧美精品国产亚洲| 夜夜夜夜夜久久久久| 一级av片app| 在线播放国产精品三级| 亚洲无线观看免费| 亚洲激情五月婷婷啪啪| 五月伊人婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 国产精品嫩草影院av在线观看| 精品99又大又爽又粗少妇毛片| av在线亚洲专区| 我的老师免费观看完整版| 日韩欧美一区二区三区在线观看| 搡老妇女老女人老熟妇| 精品人妻熟女av久视频| 日韩欧美国产在线观看| 国内精品一区二区在线观看| 久久99蜜桃精品久久| 成人特级av手机在线观看| 久久午夜亚洲精品久久| 日韩欧美国产在线观看| 黄色日韩在线| 男女做爰动态图高潮gif福利片| 岛国在线免费视频观看| 精品久久久久久久久久免费视频| 日日啪夜夜撸| 99国产精品一区二区蜜桃av| 午夜精品一区二区三区免费看| 国产av在哪里看| 99久久无色码亚洲精品果冻| 欧美一级a爱片免费观看看| 亚洲乱码一区二区免费版| av卡一久久| 黄色欧美视频在线观看| 村上凉子中文字幕在线| 国产成人影院久久av| 赤兔流量卡办理| 男女下面进入的视频免费午夜| 男人和女人高潮做爰伦理| 亚洲欧美成人精品一区二区| 久久精品国产亚洲av香蕉五月| 亚洲四区av| 中国美白少妇内射xxxbb| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 亚洲,欧美,日韩| 久久久久久久久大av| 国产精品美女特级片免费视频播放器| 国产精品国产高清国产av| 国产国拍精品亚洲av在线观看| 日本黄色片子视频| 亚洲成人av在线免费| 男女边吃奶边做爰视频| 中文字幕免费在线视频6| 国内精品美女久久久久久| 少妇裸体淫交视频免费看高清| 久久久久久久久久久丰满| 精品人妻偷拍中文字幕| av卡一久久| 看片在线看免费视频| 久久精品国产自在天天线| 午夜精品在线福利| 精品少妇黑人巨大在线播放 | 国产精品久久久久久精品电影小说 | 精品不卡国产一区二区三区| 国产视频首页在线观看| 久久精品国产亚洲av天美| 成人毛片a级毛片在线播放| 搡老妇女老女人老熟妇| 久久久精品大字幕| 久久国产乱子免费精品| 在现免费观看毛片| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 成人二区视频| 国产免费一级a男人的天堂| 国产黄色小视频在线观看| 老女人水多毛片| 男人舔奶头视频| 男女视频在线观看网站免费| 欧美最黄视频在线播放免费| eeuss影院久久| 在现免费观看毛片| 国产精品一区二区三区四区免费观看| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 国产高清视频在线观看网站| 青春草国产在线视频 | а√天堂www在线а√下载| 久久人妻av系列| 久久99精品国语久久久| 亚洲欧美清纯卡通| 亚洲欧美日韩东京热| 国产成年人精品一区二区| 国产精品一二三区在线看| 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 亚洲美女搞黄在线观看| 在线观看av片永久免费下载| 性插视频无遮挡在线免费观看| 精品久久久久久久人妻蜜臀av| 少妇丰满av| 久99久视频精品免费| 欧美不卡视频在线免费观看| 九色成人免费人妻av| 欧美高清性xxxxhd video| 欧美一区二区精品小视频在线| 嫩草影院入口| 欧美日本视频| av又黄又爽大尺度在线免费看 | 精品久久久噜噜| 国产一区二区在线观看日韩| 久久精品国产亚洲网站| 99视频精品全部免费 在线| 观看美女的网站| 国产精品野战在线观看| 国产精品久久久久久精品电影| 国产亚洲欧美98| 国产亚洲精品久久久久久毛片| 校园人妻丝袜中文字幕| 国产在线男女| 在线免费观看的www视频| 特大巨黑吊av在线直播| 91午夜精品亚洲一区二区三区| videossex国产| 特级一级黄色大片| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 精华霜和精华液先用哪个| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 一个人观看的视频www高清免费观看| 日韩欧美在线乱码| 欧美精品国产亚洲| 国产av一区在线观看免费| 欧美日本亚洲视频在线播放| 国产精品美女特级片免费视频播放器| 在线天堂最新版资源| 国产黄片视频在线免费观看| 国产午夜精品一二区理论片| 深夜精品福利| 国产久久久一区二区三区| 久久人人精品亚洲av| 国产成人aa在线观看| 亚洲国产高清在线一区二区三| 精品免费久久久久久久清纯| 国产在线男女| 久久国产乱子免费精品| 色播亚洲综合网| 国产av一区在线观看免费| 在线观看66精品国产| 丰满的人妻完整版| 国产精品久久久久久av不卡| 级片在线观看| 伦理电影大哥的女人| 国产精品,欧美在线| 国产真实乱freesex| 麻豆成人av视频| 久久99热这里只有精品18| 一本一本综合久久| 欧美高清成人免费视频www| 少妇人妻一区二区三区视频| 亚洲高清免费不卡视频| 国产欧美日韩精品一区二区| 深夜a级毛片| 男插女下体视频免费在线播放| 国产欧美日韩精品一区二区| 我的女老师完整版在线观看| 亚洲国产精品久久男人天堂| 男人舔女人下体高潮全视频| 美女高潮的动态| 中文字幕av成人在线电影| 男人的好看免费观看在线视频| 免费无遮挡裸体视频| 久久综合国产亚洲精品| 国产午夜精品一二区理论片| 床上黄色一级片| 国产精品麻豆人妻色哟哟久久 | 蜜臀久久99精品久久宅男| 在线免费十八禁| 免费观看在线日韩| 欧美激情在线99| 精品久久久噜噜| 亚洲成人精品中文字幕电影| 一边摸一边抽搐一进一小说| 欧美最黄视频在线播放免费| 此物有八面人人有两片| 一个人免费在线观看电影| 亚洲不卡免费看| 国产成人精品久久久久久| 乱人视频在线观看| 九九爱精品视频在线观看| 人妻夜夜爽99麻豆av| 美女xxoo啪啪120秒动态图| 日本三级黄在线观看| 男的添女的下面高潮视频| 亚洲av第一区精品v没综合| 亚洲在线自拍视频| 免费观看a级毛片全部| 国产成人影院久久av| 91精品一卡2卡3卡4卡| 久久久久久久久大av| a级毛色黄片| 欧美区成人在线视频| 最近视频中文字幕2019在线8| 精品久久久噜噜| 国产一区亚洲一区在线观看| 国产综合懂色| 99热只有精品国产| 中文字幕人妻熟人妻熟丝袜美| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲丝袜综合中文字幕| 又粗又硬又长又爽又黄的视频 | 日韩欧美一区二区三区在线观看| 别揉我奶头 嗯啊视频| 日韩精品有码人妻一区| 午夜精品国产一区二区电影 | 成年免费大片在线观看| 国产极品天堂在线| 国产视频内射| 亚洲真实伦在线观看| 免费在线观看成人毛片| 2022亚洲国产成人精品| 日本撒尿小便嘘嘘汇集6| 亚洲精品自拍成人| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久亚洲av鲁大| 美女 人体艺术 gogo| 在线观看av片永久免费下载| 国产麻豆成人av免费视频| 日韩人妻高清精品专区| 久久精品国产亚洲网站| 人人妻人人澡欧美一区二区| 亚洲一区高清亚洲精品| 深夜a级毛片| 91精品国产九色| 亚洲成人中文字幕在线播放| 特大巨黑吊av在线直播| 国内揄拍国产精品人妻在线| 校园春色视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲不卡免费看| 九九久久精品国产亚洲av麻豆| 久久久久久久久久久丰满| 久久久成人免费电影| 观看免费一级毛片| 国产男人的电影天堂91| 免费人成在线观看视频色| 精品久久久久久久久av| 国产午夜精品一二区理论片| av在线天堂中文字幕| 尾随美女入室| 国产黄色小视频在线观看| 天天一区二区日本电影三级| 99久久无色码亚洲精品果冻| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 波野结衣二区三区在线| 亚洲精品国产av成人精品| 成人特级黄色片久久久久久久| 久久久久九九精品影院| 我的老师免费观看完整版| 午夜爱爱视频在线播放| 亚洲精品久久久久久婷婷小说 | 极品教师在线视频| 日韩精品青青久久久久久| 春色校园在线视频观看| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 亚洲精品乱码久久久v下载方式| 中文字幕av成人在线电影| 好男人视频免费观看在线| 亚洲精华国产精华液的使用体验 | 午夜久久久久精精品| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站| 亚洲av成人精品一区久久| 哪个播放器可以免费观看大片| 可以在线观看的亚洲视频| 在线观看免费视频日本深夜| 91在线精品国自产拍蜜月| 看片在线看免费视频| 丰满人妻一区二区三区视频av| 在线免费观看不下载黄p国产| 欧美区成人在线视频| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 亚洲欧美清纯卡通| 男女啪啪激烈高潮av片| 国产单亲对白刺激| 日日撸夜夜添| 国产伦精品一区二区三区四那| 国产精品爽爽va在线观看网站| 久久久久久久久久成人| 亚洲av熟女| 一级黄色大片毛片| 一本一本综合久久| 免费看a级黄色片| 精品日产1卡2卡| 国产精品久久久久久精品电影| 精品久久久噜噜| 成人特级av手机在线观看| 久久久久久久久久成人| 久久午夜福利片| 国内精品美女久久久久久| 99热6这里只有精品| 成人三级黄色视频| 中文字幕人妻熟人妻熟丝袜美| 又粗又硬又长又爽又黄的视频 | 久久久久久久久久黄片| 少妇的逼水好多| 69人妻影院| 精品少妇黑人巨大在线播放 | 国产精品爽爽va在线观看网站| 一本精品99久久精品77| 99久久精品热视频| 午夜激情欧美在线| 日本免费一区二区三区高清不卡| 国产亚洲精品久久久com| 国产精品一区二区在线观看99 | 日韩大尺度精品在线看网址| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| av在线观看视频网站免费| 狂野欧美白嫩少妇大欣赏| 悠悠久久av| 老司机影院成人| 你懂的网址亚洲精品在线观看 | 亚洲国产精品成人久久小说 | 亚洲欧美日韩高清专用| 在线a可以看的网站| 免费搜索国产男女视频| 欧美激情国产日韩精品一区| 99九九线精品视频在线观看视频| 成人午夜精彩视频在线观看| АⅤ资源中文在线天堂| 联通29元200g的流量卡| 黄色视频,在线免费观看| h日本视频在线播放| 国产亚洲av嫩草精品影院| 日日啪夜夜撸| 99久久人妻综合| 久久精品影院6| 91狼人影院| 亚洲五月天丁香| 悠悠久久av| 久久久国产成人精品二区| 国产蜜桃级精品一区二区三区| 国产探花在线观看一区二区| 最近中文字幕高清免费大全6| 一个人看的www免费观看视频| 好男人在线观看高清免费视频| 久久人人精品亚洲av| 在线播放国产精品三级| 91久久精品国产一区二区三区| 人妻制服诱惑在线中文字幕| 国产亚洲av嫩草精品影院| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲av天美| 九色成人免费人妻av| 国产精品乱码一区二三区的特点| 亚洲最大成人手机在线| 亚洲欧美日韩无卡精品| 男女那种视频在线观看| 青青草视频在线视频观看| 美女脱内裤让男人舔精品视频 | 亚洲中文字幕一区二区三区有码在线看| 久久精品夜色国产| 热99在线观看视频| 免费一级毛片在线播放高清视频| 狠狠狠狠99中文字幕| 国产极品精品免费视频能看的| 成年免费大片在线观看| 色吧在线观看| 在线播放无遮挡| 一边摸一边抽搐一进一小说| 久久欧美精品欧美久久欧美| 久久久久久久久久久丰满| 国产精品日韩av在线免费观看| 老熟妇乱子伦视频在线观看| 国产人妻一区二区三区在| 少妇被粗大猛烈的视频| 国产日韩欧美在线精品| а√天堂www在线а√下载| 午夜福利在线观看免费完整高清在 | 全区人妻精品视频| 麻豆精品久久久久久蜜桃| www.色视频.com| kizo精华| videossex国产| av天堂在线播放| 在线播放国产精品三级| 黄色视频,在线免费观看| 人体艺术视频欧美日本| 亚洲一区高清亚洲精品| 国产一区亚洲一区在线观看| 天天一区二区日本电影三级| 欧美色视频一区免费| 国产一区二区在线观看日韩| 美女被艹到高潮喷水动态| 丝袜美腿在线中文| 天天躁夜夜躁狠狠久久av| 国产久久久一区二区三区| 一级毛片电影观看 | 午夜福利高清视频| 国产精品久久电影中文字幕| 亚洲精品国产成人久久av| 变态另类成人亚洲欧美熟女| 99热这里只有是精品在线观看| 久久中文看片网| 精品少妇黑人巨大在线播放 | 欧美最黄视频在线播放免费| 69av精品久久久久久| 99久久久亚洲精品蜜臀av| 亚洲久久久久久中文字幕| 久久久久久大精品| 九九在线视频观看精品| 亚洲精品456在线播放app| 国产亚洲91精品色在线| 成年版毛片免费区| 国产美女午夜福利| 亚洲av男天堂| 欧美激情国产日韩精品一区| 麻豆成人午夜福利视频| 国产高清激情床上av| 亚洲欧美成人精品一区二区| 免费观看的影片在线观看| 少妇人妻一区二区三区视频| 成人一区二区视频在线观看| 91久久精品国产一区二区成人| 欧美+亚洲+日韩+国产| 级片在线观看| 日本av手机在线免费观看| 一区二区三区高清视频在线| 亚洲久久久久久中文字幕| a级毛片免费高清观看在线播放| 三级男女做爰猛烈吃奶摸视频| 亚洲成a人片在线一区二区| 深夜精品福利| 免费无遮挡裸体视频| 久久人人爽人人片av| 黑人高潮一二区| 国内精品一区二区在线观看| 少妇人妻精品综合一区二区 | 日韩中字成人| 一级毛片aaaaaa免费看小| videossex国产| 精品人妻视频免费看| 啦啦啦观看免费观看视频高清| 男女视频在线观看网站免费| 蜜臀久久99精品久久宅男| 日韩视频在线欧美| 毛片一级片免费看久久久久| 人人妻人人澡欧美一区二区| 欧美一级a爱片免费观看看| 国产成人a∨麻豆精品| 亚洲乱码一区二区免费版| 免费观看精品视频网站| 欧美成人一区二区免费高清观看| 亚洲成人精品中文字幕电影| 日韩一区二区视频免费看| 国模一区二区三区四区视频| 97人妻精品一区二区三区麻豆| 国产在线男女| 禁无遮挡网站| 男人舔女人下体高潮全视频| 亚洲欧美成人精品一区二区|