概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generali...概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation,GDEE).对于GDEE的求解通常采用有限差分法,然而,由于GDEE是初始条件间断的变系数一阶双曲偏微分方程,通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象.文章从全局逼近的角度出发,基于Chebyshev拟谱法为GDEE构造了全局插值格式,解决了数值色散、数值耗散以及网格敏感性问题.考虑GDEE的系数在每个时间步长均为常数,推导了GDEE在每一个时间步长内时域上的序列矩阵指数解.由于序列矩阵指数解形式上是解析的,从而很好地克服了数值稳定性问题.两个数值算例表明,通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合,且数值耗散和数值色散现象几乎可以忽略.此外,拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制,因此该方法具有良好的数值稳定性和计算效率.展开更多
In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell＇s equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not req...In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell＇s equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not require permittivity averaging so that better accuracy can be achieved. Our numerical results on optical force computation show that the standard FDTD method converges linearly, while the proposed method achieves quadratic convergence and better accuracy.展开更多
The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing co...The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.展开更多
Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principl...Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principle. Secondly, a study for the requirements demonstration on the next-generation GOCE Follow-On satellite gravity gradiometry system is developed using different satellite orbital altitudes and measurement accuracies of satellite gravity gradiometer by the new analytical error model of SGG. The research results show that it is preferable to design satellite orbital altitudes of 300 km–400km and choose the measurement accuracies of 10-13/s2 –10-15/s2 from satellite gravity gradiometer. Finally, the complementarity of the four-stage satellite gravity missions, including past CHAMP, current GRACE, and GOCE, and next-generation GOCE Follow-On, is contrastively demonstrated for precisely recovering the Earth’s full-frequency gravitational field with high spatial resolution.展开更多
In this article, we discuss nonsymmetric solutions of the colored Yang-Baxter equation dependent on spectral as well as colored parameters and give all seven-vertex solutions by Wu＇s method. It is also proved that th...In this article, we discuss nonsymmetric solutions of the colored Yang-Baxter equation dependent on spectral as well as colored parameters and give all seven-vertex solutions by Wu＇s method. It is also proved that the solutions are composed of six groups of basic solutions up to five solution transformations. Moreover, al l solutions can be classified into two categories called Baxter type and free-fermion type.展开更多
This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton）...This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.展开更多
The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic mod...The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic models in this field, which are widely accepted and proven to be valid in many experiments and engineering. Although the contact models have evolved considerably in recent years, the verifications of the models are most based on the indirect methods such as electrical conductivity and contact stiffness, because of the lack of effective methods to directly measure the variation of contact surface. In this paper, the total reflection(TR) method is introduced into the verification of contact models.An experiment system based on TR method is constructed to measure the real contact area of two PMMA specimens. The comparison analysis between the results of experiment and models suggests that the experiment result has the same trend with simulation, the MB model has better agreement with the experimental result because this method can take into account the variation of radius and the merging of asperities, while the GW model has a huge deviation because of the dependence on resolution and the lack of considering the variation of radius and asperity's merging process. Taking the interaction of asperities into account could give a better result that is closer to the experiment. Our results and analysis prove that the experimental methods in this paper could be used as a more direct and valid method to quantitatively measure the real contact area and to verify the contact models.展开更多
In this work,we investigate a classical pseudomonotone and Lipschitz continuous variational inequality in the setting of Hilbert space,and present a projection-type approximation method for solving this problem.Our me...In this work,we investigate a classical pseudomonotone and Lipschitz continuous variational inequality in the setting of Hilbert space,and present a projection-type approximation method for solving this problem.Our method requires only to compute one projection onto the feasible set per iteration and without any linesearch procedure or additional projections as well as does not need to the prior knowledge of the Lipschitz constant and the sequentially weakly continuity of the variational inequality mapping.A strong convergence is established for the proposed method to a solution of a variational inequality problem under certain mild assumptions.Finally,we give some numerical experiments illustrating the performance of the proposed method for variational inequality problems.展开更多
A computational method of constraint stabilization and correction is introduced. The method is based on the Baumgart＇s one-step method. Constraint conditions are addressed to stabilize and correct the solution. Two e...A computational method of constraint stabilization and correction is introduced. The method is based on the Baumgart＇s one-step method. Constraint conditions are addressed to stabilize and correct the solution. Two examples are given to illustrate the results of the method.展开更多
文摘概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation,GDEE).对于GDEE的求解通常采用有限差分法,然而,由于GDEE是初始条件间断的变系数一阶双曲偏微分方程,通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象.文章从全局逼近的角度出发,基于Chebyshev拟谱法为GDEE构造了全局插值格式,解决了数值色散、数值耗散以及网格敏感性问题.考虑GDEE的系数在每个时间步长均为常数,推导了GDEE在每一个时间步长内时域上的序列矩阵指数解.由于序列矩阵指数解形式上是解析的,从而很好地克服了数值稳定性问题.两个数值算例表明,通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合,且数值耗散和数值色散现象几乎可以忽略.此外,拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制,因此该方法具有良好的数值稳定性和计算效率.
基金supported by the Air Force Office of Scientific Research (AFOSR) under Grant numbers FA9550-04-1-0213 and FA9550-07-1-0010
文摘In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell＇s equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not require permittivity averaging so that better accuracy can be achieved. Our numerical results on optical force computation show that the standard FDTD method converges linearly, while the proposed method achieves quadratic convergence and better accuracy.
基金the National Natural Science Foundation of China(No.U2032141)the Open Project of Guangxi Key Laboratory of Nuclear Physics and Nuclear Technology(No.NLK2022-02)+4 种基金the Central Government Guidance Funds for Local Scientific and Technological Development,China(Guike ZY22096024)the Natural Science Foundation of Henan Province(No.202300410479)the Guizhou Provincial Science and Technology Projects(No.ZK[2022]203)the Foundation of Fundamental Research for Young Teachers of Zhengzhou University(No.JC202041041)the Physics Research and Development Program of Zhengzhou University(No.32410217).
文摘The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.
基金Project supported by the Main Direction Program of Knowledge Innovation of Chinese Academy of Sciences for Distinguished Young Scholar (Grant No. KZCX2-EW-QN114)the National Natural Science Foundation of China for Young Scholar (Grant Nos. 41004006, 41131067, 11173049, and 41202094)+5 种基金the Merit-based Scientific Research Foundation of the State Ministry of Human Resources and Social Security of China for Returned Overseas Chinese Scholars(Grant No. 2011)the Open Research Fund Program of the Key Laboratory of Computational Geodynamics of Chinese Academy of Sciences (Grant No. 2011-04)the Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, China (Grant No. 11-01-02)the Open Research Fund Program of the Key Laboratory of Geo-Informatics of National Administration of Surveying, Mapping and Geoinformation of China(Grant No. 201322)the Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Grant No. PLN1113)the Foundation of State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing (Grant No. PRP/open-1206)
文摘Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principle. Secondly, a study for the requirements demonstration on the next-generation GOCE Follow-On satellite gravity gradiometry system is developed using different satellite orbital altitudes and measurement accuracies of satellite gravity gradiometer by the new analytical error model of SGG. The research results show that it is preferable to design satellite orbital altitudes of 300 km–400km and choose the measurement accuracies of 10-13/s2 –10-15/s2 from satellite gravity gradiometer. Finally, the complementarity of the four-stage satellite gravity missions, including past CHAMP, current GRACE, and GOCE, and next-generation GOCE Follow-On, is contrastively demonstrated for precisely recovering the Earth’s full-frequency gravitational field with high spatial resolution.
基金supported by NKBRPC(2004CB31800, 2006CB805905)Knowledge Innovation Funds of CAS (KJCX3-SYW-S03)
文摘In this article, we discuss nonsymmetric solutions of the colored Yang-Baxter equation dependent on spectral as well as colored parameters and give all seven-vertex solutions by Wu＇s method. It is also proved that the solutions are composed of six groups of basic solutions up to five solution transformations. Moreover, al l solutions can be classified into two categories called Baxter type and free-fermion type.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10871117 and 10571110)
文摘This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.
基金Project supported by the National Natural Science Foundation of China(Grant No.11872033)the Beijing Natural Science Foundation,China(Grant No.3172017)
文摘The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic models in this field, which are widely accepted and proven to be valid in many experiments and engineering. Although the contact models have evolved considerably in recent years, the verifications of the models are most based on the indirect methods such as electrical conductivity and contact stiffness, because of the lack of effective methods to directly measure the variation of contact surface. In this paper, the total reflection(TR) method is introduced into the verification of contact models.An experiment system based on TR method is constructed to measure the real contact area of two PMMA specimens. The comparison analysis between the results of experiment and models suggests that the experiment result has the same trend with simulation, the MB model has better agreement with the experimental result because this method can take into account the variation of radius and the merging of asperities, while the GW model has a huge deviation because of the dependence on resolution and the lack of considering the variation of radius and asperity's merging process. Taking the interaction of asperities into account could give a better result that is closer to the experiment. Our results and analysis prove that the experimental methods in this paper could be used as a more direct and valid method to quantitatively measure the real contact area and to verify the contact models.
基金funded by National University ofCivil Engineering(NUCE)under grant number 15-2020/KHXD-TD。
文摘In this work,we investigate a classical pseudomonotone and Lipschitz continuous variational inequality in the setting of Hilbert space,and present a projection-type approximation method for solving this problem.Our method requires only to compute one projection onto the feasible set per iteration and without any linesearch procedure or additional projections as well as does not need to the prior knowledge of the Lipschitz constant and the sequentially weakly continuity of the variational inequality mapping.A strong convergence is established for the proposed method to a solution of a variational inequality problem under certain mild assumptions.Finally,we give some numerical experiments illustrating the performance of the proposed method for variational inequality problems.
基金the National Natural Science Foundation of China (10572021 ,10372053)Basic Research Foundation of Beijing Institute of Tech-nology (BIT-UBF-200507A4206)
文摘A computational method of constraint stabilization and correction is introduced. The method is based on the Baumgart＇s one-step method. Constraint conditions are addressed to stabilize and correct the solution. Two examples are given to illustrate the results of the method.
