To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition （EMD）, an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the ...To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition （EMD）, an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with ＂not-a-knot＂ conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie＇ s derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.展开更多
High-power O-mode radio waves can excite artificial instabilities in the F region,according to experiments conducted at the European Incoherent Scatter Science Association(EISCAT)heating facility.The main instabilitie...High-power O-mode radio waves can excite artificial instabilities in the F region,according to experiments conducted at the European Incoherent Scatter Science Association(EISCAT)heating facility.The main instabilities include the parametric decay instability(PDI),oscillating two-stream instability(OTSI),and thermal parametric instability(TPI).The PDI and OTSI not only compete with each other,but also compete with the TPI,leading to a two-stage overshoot phenomenon:a miniovershoot occurs on a millisecond time scale after pump-on,followed by the main overshoot.We gain insight into the miniovershoot via a generalized Zakharov model,whereas the main overshoot can be observed as an enhanced plasma line overshoot phenomenon in incoherent scatter radar spectra.We can also observe that the zero-frequency ion line exists only in the initial heating period after a cold start and that the upshifted and downshifted ion lines behave irregularly in the spectra.The simulation results show that competition between the PDI and OTSI leads to an initial peak,which we named the pre-miniovershoot.The following processes,namely ion density caviton generation,and collapse and cascade in the development of the PDI,contribute to the miniovershoot phenomenon.展开更多
There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. T...There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. The cubic spline fitting has been analysed, and the causes of producing the overshoot and undershoot phenomenon and the end swing phenomenon have been pointed out in this paper. Two new methods of cubic spline fitting and sine spline fitting and the new technique of handling the end points of the original data curve can completely re- move the overshoot and undershoot phenomenon and the end swing phenomenon on the condition of unchanging original data, and have the advantages of the continuous fitting functions and its continuous one order derivative, the simple and convenient calculations, the small calculation amount and the easy work on it.展开更多
Based on the turbulent convection model （TCM） of Li ＆ Yang, we have studied the characteristics of turbulent convection in the envelopes of 2 and 5M⊙ stars at the red giant branch and asymptotic giant branch phase...Based on the turbulent convection model （TCM） of Li ＆ Yang, we have studied the characteristics of turbulent convection in the envelopes of 2 and 5M⊙ stars at the red giant branch and asymptotic giant branch phases. The TCM has been successfully applied over the entire convective envelopes, including the convective unstable zone and the overshooting regions. We find that the convective motions become progressively stronger when the stellar models are located farther up along the Hayashi line. In the convective unstable zone, we find that the turbulent correlations are proportional to functions of a common factor （V - V^d）T, which explains similar distributions in those correlations. For the TCM we find that if the obtained stellar temperature structure is close to that of the mixing length theory （MLT）, the convective motion will have a much larger velocity and thus be more violent. However, if the turbulent velocity is adjusted to be close to that of the MLT, the superadiabatic convection zone would be much more extended inward, which would lead to a lower effective temperature of the stellar model. For the overshooting distance, we find that the e-folding lengths of the turbulent kinetic energy k in both the top and bottom overshooting regions decrease as the stellar model is progressively located farther up along the Hayashi line, but both the extents of the decrease are not obvious. The overshooting distances of the turbulent correlation /u＇rT＂ are almost the same for the different stellar models with the same set of TCM parameters. For the decay modes of the kinetic energy k, we find that they are very similar for different stellar models based on the same set of TCM parameters, and there is a nearly linear relationship between lg k and In P for different stellar models. When Cs or α increases while the other parameters are fixed, the obtained linearly decaying distance will become longer.展开更多
Based on Multifunctional Transport Satellite data and the infrared window-texture detection algorithm, the level of overshooting top（OT） activity within a tropical cyclone（TC）, which is defined as the hourly mean ...Based on Multifunctional Transport Satellite data and the infrared window-texture detection algorithm, the level of overshooting top（OT） activity within a tropical cyclone（TC）, which is defined as the hourly mean number of OT occurrence,was statistically investigated in the western North Pacific basin for the period 2005–12. Based on the level of OT activity,the samples were divided into OT and non-OT cases or high-activity-OT（HA-OT） and low-activity-OT（LA-OT） cases. The differences in large-scale environmental variables between OT（HA-OT） and non-OT（LA-OT） cases were examined 12 hours prior to the OT occurrence. Statistical analysis showed that environmental differences did exist between the OT and non-OT cases. The OTs were more skewed towards the early stage of the TC life cycle, and mostly concentrated in low latitudes.Meanwhile, a sufficiently deep warm-water layer, large temperature difference between the upper- and lower-level troposphere, large humidity at the middle and upper levels, and large atmospheric instability, were favorable for OT occurrence.The differences in large-scale environmental characteristics between HA-OTs and LA-OTs were not as significant as those between OTs and non-OTs, but the HA-OT samples tended to occur when the vertical shear was weak and the TC intensity was low. Finally, statistical models were designed to predict the OT and HA-OT. When at least three OT（HA-OT） predictor thresholds were satisfied, the Peirce skill score reached a maximum value of 0.49（0.30）.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
基金the Ministerial Level Advanced Research Foundation (445030705QB0301)
文摘To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition （EMD）, an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with ＂not-a-knot＂ conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie＇ s derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.
基金the National Natural Science Foundation of China(NSFC Grant Nos.42104150,42074187,41774162,and 41704155)the Foundation of the National Key Laboratory of the Electromagnetic Environment(Grant No.6142403200303)+2 种基金the Chinese Academy of Sciences,Key Laboratory of Geospace Environment,University of Science&Technology of China(Grant No.GE2020-01)the Fundamental Research Funds for the Central Universities(Grant No.2042021kf0020)the Excellent Youth Foundation of Hubei Provincial Natural Science Foundation(Grant No.2019CFA054).
文摘High-power O-mode radio waves can excite artificial instabilities in the F region,according to experiments conducted at the European Incoherent Scatter Science Association(EISCAT)heating facility.The main instabilities include the parametric decay instability(PDI),oscillating two-stream instability(OTSI),and thermal parametric instability(TPI).The PDI and OTSI not only compete with each other,but also compete with the TPI,leading to a two-stage overshoot phenomenon:a miniovershoot occurs on a millisecond time scale after pump-on,followed by the main overshoot.We gain insight into the miniovershoot via a generalized Zakharov model,whereas the main overshoot can be observed as an enhanced plasma line overshoot phenomenon in incoherent scatter radar spectra.We can also observe that the zero-frequency ion line exists only in the initial heating period after a cold start and that the upshifted and downshifted ion lines behave irregularly in the spectra.The simulation results show that competition between the PDI and OTSI leads to an initial peak,which we named the pre-miniovershoot.The following processes,namely ion density caviton generation,and collapse and cascade in the development of the PDI,contribute to the miniovershoot phenomenon.
基金The Foundation Research and Development Programs of China (2004CB418404).
文摘There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. The cubic spline fitting has been analysed, and the causes of producing the overshoot and undershoot phenomenon and the end swing phenomenon have been pointed out in this paper. Two new methods of cubic spline fitting and sine spline fitting and the new technique of handling the end points of the original data curve can completely re- move the overshoot and undershoot phenomenon and the end swing phenomenon on the condition of unchanging original data, and have the advantages of the continuous fitting functions and its continuous one order derivative, the simple and convenient calculations, the small calculation amount and the easy work on it.
基金Supported by the National Natural Science Foundation of China(Grant Nos. 10973035 and 10673030)the Knowledge Innovation Key Program of the Chinese Academy of Sciences under Grant No. KJCX2-YW-T24
文摘Based on the turbulent convection model （TCM） of Li ＆ Yang, we have studied the characteristics of turbulent convection in the envelopes of 2 and 5M⊙ stars at the red giant branch and asymptotic giant branch phases. The TCM has been successfully applied over the entire convective envelopes, including the convective unstable zone and the overshooting regions. We find that the convective motions become progressively stronger when the stellar models are located farther up along the Hayashi line. In the convective unstable zone, we find that the turbulent correlations are proportional to functions of a common factor （V - V^d）T, which explains similar distributions in those correlations. For the TCM we find that if the obtained stellar temperature structure is close to that of the mixing length theory （MLT）, the convective motion will have a much larger velocity and thus be more violent. However, if the turbulent velocity is adjusted to be close to that of the MLT, the superadiabatic convection zone would be much more extended inward, which would lead to a lower effective temperature of the stellar model. For the overshooting distance, we find that the e-folding lengths of the turbulent kinetic energy k in both the top and bottom overshooting regions decrease as the stellar model is progressively located farther up along the Hayashi line, but both the extents of the decrease are not obvious. The overshooting distances of the turbulent correlation /u＇rT＂ are almost the same for the different stellar models with the same set of TCM parameters. For the decay modes of the kinetic energy k, we find that they are very similar for different stellar models based on the same set of TCM parameters, and there is a nearly linear relationship between lg k and In P for different stellar models. When Cs or α increases while the other parameters are fixed, the obtained linearly decaying distance will become longer.
基金partially supported by the National Fundamental Research(973)Program of China(Grant Nos.2013CB430100 and 2015CB452800)the National Natural Science Foundation of China(Grant Nos.41275057 and 41505086)
文摘Based on Multifunctional Transport Satellite data and the infrared window-texture detection algorithm, the level of overshooting top（OT） activity within a tropical cyclone（TC）, which is defined as the hourly mean number of OT occurrence,was statistically investigated in the western North Pacific basin for the period 2005–12. Based on the level of OT activity,the samples were divided into OT and non-OT cases or high-activity-OT（HA-OT） and low-activity-OT（LA-OT） cases. The differences in large-scale environmental variables between OT（HA-OT） and non-OT（LA-OT） cases were examined 12 hours prior to the OT occurrence. Statistical analysis showed that environmental differences did exist between the OT and non-OT cases. The OTs were more skewed towards the early stage of the TC life cycle, and mostly concentrated in low latitudes.Meanwhile, a sufficiently deep warm-water layer, large temperature difference between the upper- and lower-level troposphere, large humidity at the middle and upper levels, and large atmospheric instability, were favorable for OT occurrence.The differences in large-scale environmental characteristics between HA-OTs and LA-OTs were not as significant as those between OTs and non-OTs, but the HA-OT samples tended to occur when the vertical shear was weak and the TC intensity was low. Finally, statistical models were designed to predict the OT and HA-OT. When at least three OT（HA-OT） predictor thresholds were satisfied, the Peirce skill score reached a maximum value of 0.49（0.30）.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.
