Strong coupling between height of gaps and thickness of thermal boundary layer in partitioned convection system

A direct numerical simulation(DNS) method is used to calculate the partitioned convection system with Ra number ranging from 10~7 to 2×10~9.Using the boundary layer thickness to normalize the height of gaps d, we find a strong consistency between the variation of the TD number(the average value of the temperature in each heat transfer channel is averaged after taking the absolute values) with the change of the height of gaps and the variation of the TD number with the change of Ra number in partitioned convection.For a given thickness of partition, heights of gaps are approximately equal to 0.5 or 1 time of the thermal boundary layer thickness λ_θ at different Ra numbers.TD number representing temperature characteristics is almost the constant value, which means that TD number is a function of d/λ_θ only.Analysis of local temperature field of area in gaps shows that the temperature distribution in the gaps are basically the same when d/λ_θ is certain.The heat transfer Nu number of the system at d/λ_θ≈ 0.5 is larger than that of d/λ_θ≈ 1, both of them have the same scaling law with Ra number and N_u～Ra~(0.25).     

Analytical study of Cattaneo–Christov heat flux model for a boundary layer flow of Oldroyd-B fluid

We investigate the Cattaneo–Christov heat flux model for a two-dimensional laminar boundary layer flow of an incompressible Oldroyd-B fluid over a linearly stretching sheet. Mathematical formulation of the boundary layer problems is given. The nonlinear partial differential equations are converted into the ordinary differential equations using similarity transformations. The dimensionless velocity and temperature profiles are obtained through optimal homotopy analysis method(OHAM). The influences of the physical parameters on the velocity and the temperature are pointed out. The results show that the temperature and the thermal boundary layer thickness are smaller in the Cattaneo–Christov heat flux model than those in the Fourier's law of heat conduction.     

Comprehensive wind correction for a Rayleigh Doppler lidar from atmospheric temperature and pressure influences and Mie contamination

A correction considering the effects of atmospheric temperature, pressure, and Mie contamination must be performed for wind retrieval from a Rayleigh Doppler lidar(RDL), since the so-called Rayleigh response is directly related to the convolution of the optical transmission of the frequency discriminator and the Rayleigh–Brillouin spectrum of the molecular backscattering. Thus, real-time and on-site profiles of atmospheric pressure, temperature, and aerosols should be provided as inputs to the wind retrieval. Firstly, temperature profiles under 35 km and above the altitude are retrieved, respectively,from a high spectral resolution lidar(HSRL) and a Rayleigh integration lidar(RIL) incorporating to the RDL. Secondly,the pressure profile is taken from the European Center for Medium range Weather Forecast(ECMWF) analysis, while radiosonde data are not available. Thirdly, the Klett–Fernald algorithms are adopted to estimate the Mie and Rayleigh components in the atmospheric backscattering. After that, the backscattering ratio is finally determined in a nonlinear fitting of the transmission of the atmospheric backscattering through the Fabry–Perot interferometer(FPI) to a proposed model.In the validation experiments, wind profiles from the lidar show good agreement with the radiosonde in the overlapping altitude. Finally, a continuous wind observation shows the stability of the correction scheme.     

Numerical study of heat-transfer in two-and quasi-two-dimensional Rayleigh–Bénard convection

A detailed comparative numerical study between the two-dimensional(2 D) and quasi-two-dimensional(quasi-2 D)turbulent Rayleigh–B′enard(RB) convection on flow state, heat transfer, and thermal dissipation rate(TDR) is made. The Rayleigh number(Ra) in our simulations ranges up to 5×10~(10) and Prandtl number(Pr) is fixed to be 0.7. Our simulations are conducted on the Tianhe-2 supercomputer. We use an in-house code with high parallelization efficiency, based on the extended PDM–DNS scheme. The comparison shows that after a certain Ra, plumes with round shape, which is called the "temperature islands", develop and gradually dominate the flow field in the 2 D case. On the other hand, in quasi-2 D cases, plumes remain mushroom-like. This difference in morphology becomes more significant as Ra increases, as with the motion of plumes near the top and bottom plates. The exponents of the power-law relation between the Nusselt number(Nu) and Ra are 0.3 for both two cases, and the fitting pre-factors are 0.099 and 0.133 for 2 D and quasi-2 D respectively,indicating a clear difference in magnitude of the heat transfer rate between two cases. To understand this difference in the magnitude of Nu, we compare the vertical profile of the horizontally averaged TDR for both two cases. It is found that the profiles of both cases are nearly the same in the bulk, but they vary near boundaries. Comparing the bifurcation height zb with the thermal boundary layer thickness dq, it shows that zb δ_θ(3 D) δ_θ(2 D) and all three heights obey a universal power-law relation z ~Ra~(-0.30). In order to quantify the difference further, we separate the domain by zb, i.e., define the area between two zb(near top and bottom plates respectively) as the "mid region" and the rest as the "side region", and integrate TDR in corresponding regions. By comparing the integral it is found that most of the difference in TDR between two cases, which is connected to the heat transfer rate, occurs within the thermal boundary layers. We also compare the ratio of contributions to total heat transfer in BL–bulk separation and side–mid separation.     

Electron temperature fluctuation in the HT-7 Tokamak plasma observed by electron cyclotron emission imaging

The fluctuation of the electron temperature has been measured by using the electron cyclotron emission imaging in the Hefei Tokamak-7 (HT-7) plasma. The electron temperature fluctuation with a broadband spectrum shows that it propagates in the electron diamagnetic drift direction, and the mean poloidal wave-number $\bar{k}_{\theta}$is calculated to be about 1.58 cm^-1, or \bar{k}_θ ρ_s≈ 0.34. It indicates that the fluctuation should come from the electron drift wave turbulence. The linear global scaling of the electron temperature fluctuation with the gradient of electron temperature is consistent with the mixing length scale qualitatively. Evolution of spectrum of the fluctuation during the sawtooth oscillation phases is investigated, and the fluctuation is found to increase with the gradient of electron temperature increasing during most phases of the sawtooth oscillation. The results indicate that the electron temperature gradient is probably the driver of the fluctuation enhancement. The steady heat flux driven by electron temperature fluctuation is estimated and compared with the results from power balance estimation.     

Heat transfer analysis in the flow of Waiters＇ B fluid with a convective boundary condition

Radiative heat transfer in the steady two-dimensional flow of Walters＇ B fluid with a non-uniform heat source/sink is investigated. An incompressible fluid is bounded by a stretching porous surface. The convective boundary condition is used for the thermal boundary layer problem. The relevant equations are first simplified under usual boundary layer assumptions and then transformed into a similar form by suitable transformations. Explicit series solutions of velocity and temperature are derived by the homotopy analysis method （HAM）. The dimensionless velocity and temperature gradients at the wall are calculated and discussed.     

Short-Time Resistively-Shunted Junction Dynamic Study on Two-Dimensional Fully Frustrated XY Model

Short-time dynamic behaviour of two-dimensional fully frustrated XY model near the Ising-like phase-transition temperature is studied with resistively shunted junction dynamics.The transition temperature Tc and the dynamic and static critical exponents z,2β/ν and ν are estimated under both a periodic boundary condition (PBC) and a fluctuating twist boundary condition (FTBC) using the short-time dynamic scaling analysis.The results show that Tc,2β/ν and ν are almost independent of the boundary conditions,whereas z depends on the boundary condition.The value of z=2.18(5) of the PBC is roughly equal to that of the FFXY model with Monte Carlo dynamics under the PBC,while z=2.05(4) of the FTBC is in agreement with that of the Coulomb gas model.     

Effects of Prandtl number in two-dimensional turbulent convection

We report a numerical study of the Prandtl-number(Pr) effects in two-dimensional turbulent Rayleigh-Benard convection.The simulations were conducted in a square box over the Pr range from 0.25 to 100 and over the Rayleigh number(Ra) range from 107 to 1010.We find that both the strength and the stability of the large-scale flow decrease with the increasing of Pr,and the flow pattern becomes plume-dominated at high Pr.The evolution in flow pattern is quantified by the Reynolds number(Re),with the Ra and the Pr scaling exponents varying from 0.54 to 0.67 and-0.87 to-0.93,respectively.It is further found that the non-dimensional heat flux at small Ra diverges strongly for different Pr,but their difference becomes marginal as Ra increases.For the thermal boundary layer,the spatially averaged thicknesses for all the Pr numbers can be described by δ_θ～Ra~(-0.30) approximately,but the local values vary a lot for different Pr,which become more uniform with Pr increasing.     

He and H sequential implantation induced surface blistering and the exfoliation of Si covered with an oxide layer

Si wafers with a 220 nm top oxide layer were sequentially implanted at room temperature with 40 keV He and 35 keV H ions at a fluence of 5×10¹⁶ /cm² and 1×10¹⁶ /cm² , respectively. Techniques of scanning electron microscopy, atomic force microscopy and cross-sectional transmission electron microscopy （XTEM） were used to characterize the thermal evolution of surface damage as well as defect microstructures. Surface blisters as well as the localized exfoliation （～0.42 μm in depth） have been observed for samples annealed at temperatures of 500℃ and above. XTEM observations reveal a variety of defect microstructures, including cavities, platelets, nanometer or micrometer sized cracks and dislocations. The platelets and cracks are mainly distributed at a depth of about 0.42 μm parallel to the sample surface, which are responsible for the occurrence of the observed surface features. The relationship between surface damage and defect microstructures is described in detail.     

Analysis of Heat Transfer Phenomenon in Magnetohydrodynamic Casson Fluid Flow Through Cattaneo–Christov Heat Diffusion Theory

Heat transport phenomenon of two-dimensional magnetohydrodynamic Casson fluid flow by employing Cattaneo–Christov heat diffusion theory is described in this work. The term of heat absorption/generation is incorporated in the mathematical modeling of present flow problem. The governing mathematical expressions are solved for velocity and temperature profiles using RKF 45 method along with shooting technique. The importance of arising nonlinear quantities namely velocity, temperature, skin-friction and temperature gradient are elaborated via plots. It is explored that the Casson parameter retarded the liquid velocity while it enhances the fluid temperature. Further, we noted that temperature and thickness of temperature boundary layer are weaker in case of Cattaneo–Christov heat diffusion model when matched with the profiles obtained for Fourier's theory of heat flux.     

Geometric quantities of lower doubly excited bound states of helium

Expectation values of single electron and interelectronic geometric quantities such as $\langle r\rangle$, $\langle r_{12}\rangle$, $\langle r_<\rangle$, $\langle r_>\rangle$, $\langle \cos\theta_{12}\rangle$ and $\langle \theta_{12}\rangle$ are calculated for doubly excited $2{\rm p}n{\rm p}\,{}^1P^{\,\rm e}\,(3\leq n\leq5),\, 2{\rm p}n{\rm p}\,{}^3\!P^{\,\rm e}\,(2\leq n\leq5)$ and $2{\rm p}n{\rm d}\,{}^{1,3}D^{\,\rm o}\,(3\leq n\leq5)$ states of helium using Hylleraas-$B$-spline basis set. The energy levels converge to at least 10 significant digits in our calculations. The extrapolated values of geometric quantities except for $\langle \theta_{12}\rangle$ reach 10 significant digits as well; $\langle \theta_{12}\rangle$ reaches at least 7 significant digits using a multipole expansion approach. Our results provide a precise reference for future research.     

激光冷却SH~–阴离子的理论研究

本文采用多组态相互作用及Davidson修正方法和全电子基组计算了SH~-阴离子的X~1∑~+,a~3∏和A~1∏态的势能曲线、电偶极矩和跃迁偶极矩.计算的光谱常数与实验值及已有的理论值符合得很好.在计算中考虑了自旋-轨道耦合效应.计算得到a~3∏_1(v'=0)?X~1∑_(0+)~+(v"=0)和A~1∏_1(v'=0)?X~1Σ_(0+)~+(v"=0)跃迁具有高对角分布的弗兰克-康登因子,分别为0.9990和0.9999;计算得到a~3∏_1和A~1∏_1态的自发辐射寿命分别为1.472和0.188 ms.A~1∏_1?X~1∑_(0+)~+跃迁存在中间态a~3∏_(0+)和a~3∏_1,但中间态对激光冷却SH~-阴离子的影响可以忽略.分别利用a~3∏_1(v'=0)? X~1∑_(0+)~+(v"=0)和A~1∏_1(v'=0)? X~1∑_(0+)~+(v"=0)跃迁构建了准闭合的能级系统,冷却所需的激光波长分别为492.27和478.57 nm.最后预测了激光冷却SH~-阴离子能达到的多普勒温度和反冲温度.这些结果为进一步实验提供了理论参数.     

Covalent hadronic molecules induced by shared light quarks

After examining Feynman diagrams corresponding to the ${\bar{D}}^{(* )}{{\rm{\Sigma }}}_{c}^{(* )}$, ${\bar{D}}^{(* )}{{\rm{\Lambda }}}_{c}$, ${D}^{(* )}{\bar{K}}^{* }$, and ${D}^{(* )}{\bar{D}}^{(* )}$ hadronic molecular states, we propose a possible binding mechanism induced by shared light quarks. This mechanism is similar to the covalent bond in chemical molecules induced by shared electrons. We use the method of QCD sum rules to calculate its corresponding light-quark-exchange diagrams, and the obtained results indicate a model-independent hypothesis: the light-quark-exchange interaction is attractive when the shared light quarks are totally antisymmetric so they obey the Pauli principle. We build a toy model with four parameters to formulize this picture and estimate binding energies of some possibly-existing covalent hadronic molecules. A unique feature of this picture is that the binding energies of the (I)J^P = (0)1⁺ $D{\bar{B}}^{* }/{D}^{* }\bar{B}$ hadronic molecules are much larger than those of the (I)J^P = (0)1⁺ ${{DD}}^{* }/\bar{B}{\bar{B}}^{* }$ ones, while the (I)J^P = (1/2)1/2⁺ $\bar{D}{{\rm{\Sigma }}}_{c}/\bar{D}{{\rm{\Sigma }}}_{b}/B{{\rm{\Sigma }}}_{c}/B{{\rm{\Sigma }}}_{b}$ hadronic molecules have similar binding energies.     

Electron tunnelling phase time and dwell time through an associated delta potential barrier

The electron tunnelling phase time τP and dwell time τD through an associated delta potential barrier U(x) = ξδ(x) are calculated and both are in the order of 10^-17～10^-16s. The results show that the dependence of the phase time on the delta barrier parameter ξ can be described by the characteristic length lc = h^2/meξ and the characteristic energy Ec=meξ^2/h^2 of the delta barrier, where me is the electron mass, lc and Ec are assumed to be the effective width and height of the delta barrier with lcEc=ξ, respectively. It is found that TD reaches its maximum and τD = τp as the energy of the tunnelling electron is equal to Ec/2, i.e. as lc =λDB, λDB is de Broglie wave length of the electron.     

Searching for isotropic stochastic gravitational-wave background in the international pulsar timing array second data release

We search for isotropic stochastic gravitational-wave background (SGWB) in the International Pulsar Timing Array second data release. By modeling the SGWB as a power-law, we find very strong Bayesian evidence for a common-spectrum process, and further this process has scalar transverse (ST) correlations allowed in general metric theory of gravity as the Bayes factor in favor of the ST-correlated process versus the spatially uncorrelated common-spectrum process is 30 ± 2. The median and the 90% equal-tail amplitudes of ST mode are ${{ \mathcal A }}_{\mathrm{ST}}={1.29}_{-0.44}^{+0.51}\times {10}^{-15}$, or equivalently the energy density parameter per logarithm frequency is ${{\rm{\Omega }}}_{\mathrm{GW}}^{\mathrm{ST}}={2.31}_{-1.30}^{+2.19}\times {10}^{-9}$, at frequency of 1 year⁻¹. However, we do not find any statistically significant evidence for the tensor transverse (TT) mode and then place the 95% upper limits as ${{ \mathcal A }}_{\mathrm{TT}}\lt 3.95\times {10}^{-15}$, or equivalently ${{\rm{\Omega }}}_{\mathrm{GW}}^{\mathrm{TT}}\lt 2.16\times {10}^{-9}$, at frequency of 1 year⁻¹.     

Parity-violating neutron spin rotation in hydrogen and deuterium

We calculate the (parity-violating) spin-rotation angle of a polarized neutron beam through hydrogen and deuterium targets, using pionless effective field theory up to next-to-leading order. Our result is part of a program to obtain the five leading independent low-energy parameters that characterize hadronic parity violation from few-body observables in one systematic and consistent framework. The two spin-rotation angles provide independent constraints on these parameters. Our result for np spin rotation is $\frac{1} {\rho }\frac{{d\varphi _{PV}^{np} }} {{dl}} = \left[ {4.5 \pm 0.5} \right] rad MeV^{ - \frac{1} {2}} \left( {2g^{\left( {^3 S_1 - ^3 P_1 } \right)} + g^{\left( {^3 S_1 - ^3 P_1 } \right)} } \right) - \left[ {18.5 \pm 1.9} \right] rad MeV^{ - \frac{1} {2}} \left( {g_{\left( {\Delta {\rm I} = 0} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} - 2g_{\left( {\Delta {\rm I} = 2} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} } \right)$\frac{1} {\rho }\frac{{d\varphi _{PV}^{np} }} {{dl}} = \left[ {4.5 \pm 0.5} \right] rad MeV^{ - \frac{1} {2}} \left( {2g^{\left( {^3 S_1 - ^3 P_1 } \right)} + g^{\left( {^3 S_1 - ^3 P_1 } \right)} } \right) - \left[ {18.5 \pm 1.9} \right] rad MeV^{ - \frac{1} {2}} \left( {g_{\left( {\Delta {\rm I} = 0} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} - 2g_{\left( {\Delta {\rm I} = 2} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} } \right), while for nd spin rotation we obtain $\frac{1} {\rho }\frac{{d\varphi _{PV}^{nd} }} {{dl}} = \left[ {8.0 \pm 0.8} \right] rad MeV^{ - \frac{1} {2}} g^{\left( {^3 S_1 - ^1 P_1 } \right)} + \left[ {17.0 \pm 1.7} \right] rad MeV^{ - \frac{1} {2}} g^{\left( {^3 S_1 - ^3 P_1 } \right)} + \left[ {2.3 \pm 0.5} \right] rad MeV^{ - \frac{1} {2}} \left( {3g_{\left( {\Delta {\rm I} = 0} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} - 2g_{\left( {\Delta {\rm I} = 1} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} } \right)$\frac{1} {\rho }\frac{{d\varphi _{PV}^{nd} }} {{dl}} = \left[ {8.0 \pm 0.8} \right] rad MeV^{ - \frac{1} {2}} g^{\left( {^3 S_1 - ^1 P_1 } \right)} + \left[ {17.0 \pm 1.7} \right] rad MeV^{ - \frac{1} {2}} g^{\left( {^3 S_1 - ^3 P_1 } \right)} + \left[ {2.3 \pm 0.5} \right] rad MeV^{ - \frac{1} {2}} \left( {3g_{\left( {\Delta {\rm I} = 0} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} - 2g_{\left( {\Delta {\rm I} = 1} \right)}^{\left( {^1 S_0 - ^3 P_0 } \right)} } \right), where the g ^(X-Y), in units of $MeV^{ - \frac{3} {2}}$MeV^{ - \frac{3} {2}}, are the presently unknown parameters in the leading-order parity-violating Lagrangian. Using naıve dimensional analysis to estimate the typical size of the couplings, we expect the signal for standard target densities to be $\left| {\frac{{d\varphi _{PV} }} {{dl}}} \right| \approx \left[ {10^{ - 7} \ldots 10^{ - 6} } \right]\frac{{rad}} {m}$\left| {\frac{{d\varphi _{PV} }} {{dl}}} \right| \approx \left[ {10^{ - 7} \ldots 10^{ - 6} } \right]\frac{{rad}} {m} for both hydrogen and deuterium targets. We find no indication that the nd observable is enhanced compared to the np one. All results are properly renormalized. An estimate of the numerical and systematic uncertainties of our calculations indicates excellent convergence. An appendix contains the relevant partial-wave projectors of the three-nucleon system.     

The splitting of low-lying states for hydroxyl molecule under spin--orbit coupling

The splitting of potential energy curves for the states $X^{2}\Pi _{3/2}$, $^{2}\Pi _{1/2}$ and $A^{2}\Sigma ^{ +}$ of hydroxyl OH under spin--orbit coupling (SOC) has been calculated by using the SO multi-configuration quasi-degenerate perturbation theory (SO-MCQDPT). Their Murrell--Sorbie (M--S) potential functions have been derived, then, the spectroscopic constants for $X^{2}\Pi _{3/2}$,$^{ 2}\Pi _{1/2}$ and $A^{2}\Sigma ^{ + }$ have been derived from the M--S function. The calculated dissociation energies for the three states are $D_{0}$[OH($X^{2}\Pi _{3/2})$]=34966.632cm$^{-1}$, $D_{0}$[OH($^{2}\Pi _{1/2})$]=34922.802cm$^{-1}$, and $D_{0}$[OH($A^{2}\Sigma ^{ + })$]=17469.794cm$^{-1}$, respectively. The vertical excitation energy $\nu [ {{ }^2\Pi _{1/2} ( {\nu = 0} ) \to {X}{ }^2\Pi _{3/2} ( {\nu = 0} )} ] = 139.6{\rm cm}^{-{\rm 1}}$. All the spectroscopic data for the $X^{2}\Pi _{3/2}$ and $^{2}\Pi _{1/2 }$ are given for the first time except the dissociation energy of $X^{2}\Pi _{3/2}$.     

Magnetization and magnetic phase diagrams of a spin-1/2 ferrimagnetic diamond chain at low temperature

We used the Jordan–Wigner transform and the invariant eigenoperator method to study the magnetic phase diagram and the magnetization curve of the spin-1/2 alternating ferrimagnetic diamond chain in an external magnetic field at finite temperature. The magnetization versus external magnetic field curve exhibits a 1/3 magnetization plateau at absolute zero and finite temperatures, and the width of the 1/3 magnetization plateau was modulated by tuning the temperature and the exchange interactions. Three critical magnetic field intensities H_(CB), H_(CE) and H_(CS) were obtained, in which the H_(CB) and H_(CE) correspond to the appearance and disappearance of the 1/3 magnetization plateau, respectively, and the higher H_(CS) correspond to the appearance of fully polarized magnetization plateau of the system. The energies of elementary excitation ωσ,k(σ = 1, 2, 3) present the extrema of zero at the three critical magnetic fields at 0 K, i.e., [hω_(3,k)(HCB)]_(min)= 0, [hω_(2,k)(H_(CE))]_(max)= 0 and [hω _(2,k)(H_(CS))]_(min)= 0, and the magnetic phase diagram of magnetic field versus different exchange interactions at 0 K was established by the above relationships. According to the relationships between the system's magnetization curve at finite temperatures and the critical magnetic field intensities, the magnetic field-temperature phase diagram was drawn. It was observed that if the magnetic phase diagram shows a three-phase critical point, which is intersected by the ferrimagnetic phase, the ferrimagnetic plateau phase, and the Luttinger liquid phase, the disappearance of the1/3 magnetization plateau would inevitably occur. However, the 1/3 magnetization plateau would not disappear without the three-phase critical point. The appearance of the 1/3 magnetization plateau in the low temperature region is the macroscopic manifestations of quantum effect.     

Constraints on $H_0$ from WMAP and BAO Measurements *

We report the constraints of $H_0$ obtained from Wilkinson Microwave Anisotropy Probe (WMAP) 9-year data combined with the latest baryonic acoustic oscillations (BAO) measurements. We use the BAO measurements from 6dF Galaxy Survey (6dFGS), the SDSS DR7 main galaxies sample (MGS), the BOSS DR12 galaxies, and the eBOSS DR14 quasars. Adding the recent BAO measurements to the cosmic microwave background (CMB) data from WMAP, we constrain cosmological parameters $\Omega_m=0.298\pm0.005$, $H_0=68.36^{+0.53}_{-0.52} {\rm km}\cdot {\rm s}^{-1}\cdot {\rm Mpc}^{-1}$, $\sigma_8=0.8170^{+0.0159}_{-0.0175}$ in a spatially flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, and $\Omega_m=0.302\pm0.008$, $H_0=67.63\pm1.30 {\rm km}\cdot{\rm s}^{-1}\cdot {\rm Mpc}^{-1}$, $\sigma_8=0.7988^{+0.0345}_{-0.0338}$ in a spatially flat $w$CDM model, respectively. Our measured $H_0$ results prefer a value lower than 70 ${\rm km}\cdot {\rm s}^{-1}\cdot{\rm Mpc}^{-1}$, consistent with the recent data on CMB constraints from Planck (2018), but in $3.1$ and $3.5\sigma$ tension with local measurements of SH0ES (2018) in $\Lambda$CDM and $w$CDM framework, respectively. Our results indicate that there is a systematic tension on the Hubble constant between SH0ES and the combination of CMB and BAO datasets.