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.展开更多
This paper presents a dual-nonlinear branch linearizer for reducing the corrected amplitude overshoot of conventional single nonlinear branch linearizer. Theoretical analysis is carried out, the analysis is verified b...This paper presents a dual-nonlinear branch linearizer for reducing the corrected amplitude overshoot of conventional single nonlinear branch linearizer. Theoretical analysis is carried out, the analysis is verified by simulation, and a prototype of Ka band 25.28~26.08 GHz dual nonlinear branch linearizer is achieved. It indicates that the corrected amplitude overshoot is less than 0.5 dB, the C/I3 improvement is more than 10 dB related to a single carrier IBO 9 dB, when it is linked and tested for 50 W spacebrone Travelling Wave Tube Amplifier(TWTA).展开更多
A study is made on the overshoot phenomena in wind-generated waves. The surface displacements of time-growing. waves are measured at four fetches in a wind wave channel. The evolution of high frequency waves is displa...A study is made on the overshoot phenomena in wind-generated waves. The surface displacements of time-growing. waves are measured at four fetches in a wind wave channel. The evolution of high frequency waves is displayed with wavelet transform. The results are compared with Sutherland's. It is found that high frequency wave components experience much stronger energy overshoot in the evolution. The energy of high frequency waves decreases greatly after overshoot.展开更多
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.展开更多
Stokes' first problem has been investigated for a Maxwell fluid in a porous half-space for gaining insight into the effect of viscoelasticity on the start-up flow in a porous medium. An exact solution was obtained by...Stokes' first problem has been investigated for a Maxwell fluid in a porous half-space for gaining insight into the effect of viscoelasticity on the start-up flow in a porous medium. An exact solution was obtained by using the Fourier sine transform. It was found that at large values of the relaxation time the velocity overshoot occurs obviously and the system exhibits viscoelastic behaviours. On the other hand, for short relaxation time the velocity overshoot disappears and the system exhibits viscous behaviours. A critical value of the relaxation time was obtained for the emergence of the velocity overshoot. Furthermore, it was found that the velocity overshoot is caused by both the viscoelasticity of the Maxwell fluid and the Darcy resistance resulting from the structure of the micropore in the porous medium.展开更多
Recently, it is proved in the literature that for a given controllable pair (A, B) with A ∈R^n×n, B ∈R^n×m, and any λ ≥ 1, a gain matrix K can be designed so that ‖e^(A+BK)t‖ ≤Mλ^Le^-λt, where ...Recently, it is proved in the literature that for a given controllable pair (A, B) with A ∈R^n×n, B ∈R^n×m, and any λ ≥ 1, a gain matrix K can be designed so that ‖e^(A+BK)t‖ ≤Mλ^Le^-λt, where M and L are constants independent of λ. Here, we show that M and L can be chosen much smaller than that proposed above. As a consequence, the estimation on overshoot of a transition matrix can be bounded more precisely. This can be regarded as a complement to the existing result.展开更多
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.展开更多
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.展开更多
Pulse inductively coupled plasma has been widely used in the microelectronics industry,but the existence of overshoot phenomenon may affect the uniformity of plasma and generate high-energy ions,which could damage the...Pulse inductively coupled plasma has been widely used in the microelectronics industry,but the existence of overshoot phenomenon may affect the uniformity of plasma and generate high-energy ions,which could damage the chip.The overshoot phenomenon at various spatial locations in pulsed inductively coupled Ar and Ar/CF_(4) discharges is studied in this work.The electron density,effective electron temperature,relative light intensity,and electron energy probability function(EEPF) are measured by using a time-resolved Langmuir probe and an optical probe,as a function of axial and radial locations.At the initial stage of pulse,both electron density and relative light intensity exhibit overshoot phenomenon,i.e.,they first increase to a peak value and then decrease to a convergent value.The overshoot phenomenon gradually decays,when the probe moves away from the coils.Meanwhile,a delay appears in the variation of the electron densities,and the effective electron temperature decreases,which may be related to the reduced strength of electric field at a distance,and the consequent fewer high-energy electrons,inducing limited ionization and excitation rate.The overshoot phenomenon gradually disappears and the electron density decreases,when the probe moves away from reactor centre.In Ar/CF_(4) discharge,the overshoot phenomenon of electron density is weaker than that in the Ar discharge,and the plasma reaches a steady density within a much shorter time,which is probably due to the more ionization channels and lower ionization thresholds in the Ar/CF_(4) plasma.展开更多
The semiclassical transport equations are used to study the high frequency performance of AlGaAs / GaAs HBTs. Electron velocity overshoot effect and its influence on the cut off frequency characteristics of AlGaAs / G...The semiclassical transport equations are used to study the high frequency performance of AlGaAs / GaAs HBTs. Electron velocity overshoot effect and its influence on the cut off frequency characteristics of AlGaAs / GaAs HBTs with different collector design parameters are analyzed and discussed.展开更多
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).展开更多
Background: Frequent outbreaks of insects and diseases have been recorded in the native forests of western North America during the last few decades, but the distribution of these outbreaks has been far from uniform....Background: Frequent outbreaks of insects and diseases have been recorded in the native forests of western North America during the last few decades, but the distribution of these outbreaks has been far from uniform. In some cases, recent climatic variations may explain some of this spatial variation along with the presence of expansive forests composed of dense, older trees. Forest managers and policy makers would benefit if areas especially prone to disturbance could be recognized so that mitigating actions could be taken. Methods: We use two ponderosa pine-dominated sites in western Montana, U.S.A. to apply a modeling approach that couples information acquired via remote sensing, soil surveys, and local weather stations to assess where bark beetle outbreaks might first occur and why. Although there was a general downward trend in precipitation for both sites over the period between 1998 and 2010 (slope =-1.3, R2 = 0.08), interannual variability was high. Some years showed large increases followed by sharp decreases. Both sites had similar topography and fire histories, but bark beetle activity occurred earlier (circa 2000 to 2001) and more severely on one site than on the other. The initial canopy density of the two sites was also similar, with leaf area indices ranging between 1.7-2.0 m2. m-2. We wondered if the difference in bark beetle activity was related to soils that were higher in clay content at site I than at site II. To assess this possibility, we applied a process-based stand growth model (3-PG) to analyze the data and evaluate the hypotheses. Results: We found that when wet years were followed by drier years, the simulated annual wood production per unit of leaf area, a measure of tree vigor, dropped below a critical threshold on site I but not on site II. Conclusion: We concluded that the difference in vulnerability of the two stands to beetle outbreaks can be explained largely by differences in gross photosynthesis attributed to the fact that an equivalent amount of stored water in the rooting zone (100 mm) is extracted less efficiently from fine-textured soils than from coarse-textured ones.展开更多
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 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.
文摘This paper presents a dual-nonlinear branch linearizer for reducing the corrected amplitude overshoot of conventional single nonlinear branch linearizer. Theoretical analysis is carried out, the analysis is verified by simulation, and a prototype of Ka band 25.28~26.08 GHz dual nonlinear branch linearizer is achieved. It indicates that the corrected amplitude overshoot is less than 0.5 dB, the C/I3 improvement is more than 10 dB related to a single carrier IBO 9 dB, when it is linked and tested for 50 W spacebrone Travelling Wave Tube Amplifier(TWTA).
文摘A study is made on the overshoot phenomena in wind-generated waves. The surface displacements of time-growing. waves are measured at four fetches in a wind wave channel. The evolution of high frequency waves is displayed with wavelet transform. The results are compared with Sutherland's. It is found that high frequency wave components experience much stronger energy overshoot in the evolution. The energy of high frequency waves decreases greatly after overshoot.
基金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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10372007 and 10572006) and the New , Century Training Programme Foundation for the Talents by Chinese Ministry of Education.
文摘Stokes' first problem has been investigated for a Maxwell fluid in a porous half-space for gaining insight into the effect of viscoelasticity on the start-up flow in a porous medium. An exact solution was obtained by using the Fourier sine transform. It was found that at large values of the relaxation time the velocity overshoot occurs obviously and the system exhibits viscoelastic behaviours. On the other hand, for short relaxation time the velocity overshoot disappears and the system exhibits viscous behaviours. A critical value of the relaxation time was obtained for the emergence of the velocity overshoot. Furthermore, it was found that the velocity overshoot is caused by both the viscoelasticity of the Maxwell fluid and the Darcy resistance resulting from the structure of the micropore in the porous medium.
基金This work was supported by the National Natural Science Foundation of China (No. 60604032, 10601050).
文摘Recently, it is proved in the literature that for a given controllable pair (A, B) with A ∈R^n×n, B ∈R^n×m, and any λ ≥ 1, a gain matrix K can be designed so that ‖e^(A+BK)t‖ ≤Mλ^Le^-λt, where M and L are constants independent of λ. Here, we show that M and L can be chosen much smaller than that proposed above. As a consequence, the estimation on overshoot of a transition matrix can be bounded more precisely. This can be regarded as a complement to the existing result.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11675039,11875101,11935005,and 12075049)the Fundamental Research Funds for the Central Universities of China(Grant Nos.DUT18TD06 and DUT20LAB201)。
文摘Pulse inductively coupled plasma has been widely used in the microelectronics industry,but the existence of overshoot phenomenon may affect the uniformity of plasma and generate high-energy ions,which could damage the chip.The overshoot phenomenon at various spatial locations in pulsed inductively coupled Ar and Ar/CF_(4) discharges is studied in this work.The electron density,effective electron temperature,relative light intensity,and electron energy probability function(EEPF) are measured by using a time-resolved Langmuir probe and an optical probe,as a function of axial and radial locations.At the initial stage of pulse,both electron density and relative light intensity exhibit overshoot phenomenon,i.e.,they first increase to a peak value and then decrease to a convergent value.The overshoot phenomenon gradually decays,when the probe moves away from the coils.Meanwhile,a delay appears in the variation of the electron densities,and the effective electron temperature decreases,which may be related to the reduced strength of electric field at a distance,and the consequent fewer high-energy electrons,inducing limited ionization and excitation rate.The overshoot phenomenon gradually disappears and the electron density decreases,when the probe moves away from reactor centre.In Ar/CF_(4) discharge,the overshoot phenomenon of electron density is weaker than that in the Ar discharge,and the plasma reaches a steady density within a much shorter time,which is probably due to the more ionization channels and lower ionization thresholds in the Ar/CF_(4) plasma.
文摘The semiclassical transport equations are used to study the high frequency performance of AlGaAs / GaAs HBTs. Electron velocity overshoot effect and its influence on the cut off frequency characteristics of AlGaAs / GaAs HBTs with different collector design parameters are analyzed and discussed.
基金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).
文摘Background: Frequent outbreaks of insects and diseases have been recorded in the native forests of western North America during the last few decades, but the distribution of these outbreaks has been far from uniform. In some cases, recent climatic variations may explain some of this spatial variation along with the presence of expansive forests composed of dense, older trees. Forest managers and policy makers would benefit if areas especially prone to disturbance could be recognized so that mitigating actions could be taken. Methods: We use two ponderosa pine-dominated sites in western Montana, U.S.A. to apply a modeling approach that couples information acquired via remote sensing, soil surveys, and local weather stations to assess where bark beetle outbreaks might first occur and why. Although there was a general downward trend in precipitation for both sites over the period between 1998 and 2010 (slope =-1.3, R2 = 0.08), interannual variability was high. Some years showed large increases followed by sharp decreases. Both sites had similar topography and fire histories, but bark beetle activity occurred earlier (circa 2000 to 2001) and more severely on one site than on the other. The initial canopy density of the two sites was also similar, with leaf area indices ranging between 1.7-2.0 m2. m-2. We wondered if the difference in bark beetle activity was related to soils that were higher in clay content at site I than at site II. To assess this possibility, we applied a process-based stand growth model (3-PG) to analyze the data and evaluate the hypotheses. Results: We found that when wet years were followed by drier years, the simulated annual wood production per unit of leaf area, a measure of tree vigor, dropped below a critical threshold on site I but not on site II. Conclusion: We concluded that the difference in vulnerability of the two stands to beetle outbreaks can be explained largely by differences in gross photosynthesis attributed to the fact that an equivalent amount of stored water in the rooting zone (100 mm) is extracted less efficiently from fine-textured soils than from coarse-textured ones.
基金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.