Intensity Oscillations Research Articles

Numerical investigation on unstable oscillation of oxygen jet condensation in cryogenic liquid rocket

Oxygen jet condensation always occurs with intense unstable oscillation in the cryogenic delivery pipe of a space launch vehicle. This phenomenon threatens the safety of the propulsion pipeline. It is essential to explore the transient physical features of oxygen jet condensation precisely, including the spectrum of pressure oscillation. Aiming at a deeper understanding of the thermal hydraulic characteristic of oxygen jet condensation, abundant simulations are carried out using a modified phase change method. In this model, the height function method is introduced to capture the interfacial curvature to reveal the mechanism of periodic fluctuation. The simulation proves that three typical flow patterns are found at variant operating conditions in oxygen jet condensation: chugging, oscillation, and swinging. The condensation regime diagram is summarized according to the flow pattern distribution. The chugging flow, along with the phenomenon of liquid oxygen sucked-back flow, shows the most intense pressure oscillation with a maximum amplitude of 133 kPa. In the oscillation flow, the oxygen vapor plume shrinks periodically with a low frequency of about 10 Hz. The swinging flow is a stable flow pattern with a slight pulsation of oxygen vapor plume wake. Its amplitude of pressure oscillation is approximately 1–3 kPa. The influence of liquid oxygen mass flow rate and subcooling on the unstable characteristic is also obtained. These findings offer ample theoretical direction for the development of the space launch vehicle.

Acoustic and optical plasmons excitation in double-channel AlGaN/GaN HEMT under asymmetric boundaries

We have simulatively studied the electrical excitation of acoustic and optical plasmons in double-channel AlGaN/GaN high electron mobility transistors (HEMT) under asymmetric boundaries. By solving the self-consistent hydrodynamic model with Maxwell’s equations, it is found that the drift plasmons are excited in the bilayer 2DEGs under asymmetric boundaries. The oscillation intensity in a bilayer system is much stronger than that in a monolayer 2DEG system because of the simultaneous excitation of the fundamental and the second harmonic components. The fundamental component always corresponds to the acoustic plasmon, while the optical plasmon can be excited when it is resonant with the second harmonic of the acoustic plasmon. Because of the out-of-phase properties of the acoustic plasmons, their radiated power is limited by the current cancellation. On the other hand, as the optical plasmons are excited, the radiated powers are much higher than those of the single channel HEMT due to the in-phase currents generated in the double layers. The effects of different parameters such as gate lengths, 2DEG spacings, and electron concentrations on plasmon excitations and THz emissions of double-channel HEMTs are analyzed in detail. The numerical results provide a more powerful terahertz radiation mechanism in double-channel HEMT compared with the Dyakonov–Shur instability in traditional single-channel HEMT.

Tunable High-Frequency Acoustoelectric Current Oscillations in Fluorine-Doped Single-Walled Carbon Nanotubes

Herein, we report on a fluorine-doped single-walled carbon nanotube (FSWCNT) phenomenon, that yields tunable high-frequency self-sustained acoustoelectric direct current (ADC) oscillations. A tractable analytical method was used in the hypersound domain, to base the calculations on carriers in the lowest miniband. Hypothetically, the energy of interaction between the carriers and the acoustic phonons is less than the energy of the typical carriers. High-order harmonics of the acoustic phonons’ effective field could be disregarded under this supposition. The ADC was observed to exhibit a nonlinearity, that resulted from the carrier distribution function’s distortion as a result of interaction with the acoustic phonons, which had strong nonlinear effects. Theoretically, we demonstrated that the dynamics of space charge instabilities, due to Bragg reflection of Bloch oscillating carriers in the FSWCNT’s miniband, were the only factors which contributed to the creation of radiation in the terahertz (THz) frequency range. The study also investigated the influence of various FSWCNT parameters such as the overlapping integrals (Δs and Δz), ac-field E1, and carrier concentration noon the behaviour of the ADC. The results showed that the intensity of the ADC oscillation Jzzae/Joae could be tuned by adjusting Δs, Δz, E1, and no.This tunability suggests that FSWCNTs could be used as an active device operating at very high frequencies, potentially reaching the submillimeter wavelength range. The study also suggests the possibility of domain suppression and acoustic Bloch gain through dynamic ADC stabilisation.

Dynamical Behavior of a Stacked Blue Laser Diode Consisting of Two Active Regions with Wide InGaN Quantum Wells and Two Tunnel Junctions

In nitride‐based laser diodes, the poor conductivity of the p‐type region, higher light absorption, and the difficult processing of an ohmic p‐contact with low resistance are major challenges. Recently interband tunnel junctions (TJs) have been used to replace much of the hole transport layers by electron transport layers, and also to allow for two n‐type contacts. In this work, GaN‐based laser diode stacks with TJs are analyzed. The devices are grown with plasma‐assisted molecular beam epitaxy and include two waveguides with active regions consisting of 25 nm wide InGaN quantum well (QW), each, which are connected through a TJ. A second TJ is placed on top to eliminate the p‐type contact. In this study, the time‐dependent spectral behavior of such structures is analyzed by streak camera measurements. The two QWs have different lasing thresholds and emission wavelength. Using the temporal resolution of the streak camera, a lasing onset of the second QW a few nanoseconds after the first one is observed. Single shot measurements show unusual intensity oscillations from both QWs, however, one is stable and the other laser diode (LD) has irregular oscillations.

A Pressure-Oscillation-Based RON Estimation Method for Spark Ignition Fuels beyond RON 100

Knock in spark ignition (SI) engines occurs when the air–fuel mixture in the combustion chamber ignites spontaneously ahead of the flame front, reducing combustion efficiency and possibly leading to engine damage if left unattended. The use of knock sensors to prevent it is common practice in modern engines. Another measure to mitigate knock is the use of higher-octane fuels. The American Society for Testing and Materials’ (ASTM) determination of the Research Octane Number (RON) and Motor Octane Number (MON) of spark ignition fuels has been based on measuring cylinder pressure rise at the onset of knock since its inception in the 1930s. This is achieved through a low-pass filtered pressure signal. Knock detection in contemporary engines, however, relies on measuring engine vibrations caused by high-frequency pressure oscillations during knock. The difference between conditions in which fuels are evaluated for their octane rating and the conditions that generate a knock intensity signal from the knock sensor suggests a potential difference between octane rating and the knock limit typically identified by a contemporary knock sensor. To address this disparity, a modified RON measurement method has been developed, incorporating pressure oscillation measurements. This test method addresses the historical lack of correlation between RON and high-frequency pressure oscillation intensity during knock. Using toluene standardization fuels (TSFs) as a reference, the obtained results demonstrate excellent high-frequency knock intensity-based RON estimations for gasoline. The method is able to differentiate between two fuels that share the same ASTM RON, associating them with a RON-like metric that is more aligned with their performance in a modern SI engine. This alternative method could potentially serve as a template for an upgrade to the existing ASTM RON method without significantly disrupting the current approach. Additionally, its capability to evaluate fuels beyond RON 100 opens the door to assessing a wider range of fuels for antiknock properties and the intensity of fuel oscillations during knocking combustion.

Iterative suppression of Kerr-induced instabilities in Bessel beams using on-axis intensity shaping

Kerr-induced instabilities in zeroth-order Bessel beams with low focusing angle prevent the formation of longitudinally uniform plasma rods in the filamentation regime. These instabilities lead to the oscillation of the beam on-axis intensity via the generation of new spatial frequencies by a first stage of spectral broadening followed by a second stage of four-wave mixing. Here, we numerically demonstrate an efficient approach to drastically reduce the instabilities due to the second stage. It is based on shaping the longitudinal intensity profile with spatio-spectral components in opposition of phase to the Kerr-generated ones via an iterative approach. Zeroth-order Bessel beams with a longitudinal flat intensity plateau can be generated in a few iterations in the nonlinear regime. This is performed in both monochromatic and pulsed femtosecond regimes.

Mechanisms of enhanced heat transfer of piston oscillating cooling in a heavy-duty diesel engine

The thermal load and reliability of pistons become major concerns with the increase of power density of internal combustion engines. The oscillating cooling technology is often adopted to enhance the heat transfer of pistons in heavy-duty diesel engines. To understand the heat transfer mechanisms, this study developed an optical oscillating cooling piston based on a real heavy-duty engine. The characteristics of two-phase oscillating flows were optically observed under different conditions. As engine speed increased from 100 to 700 r/min or oil supply pressure decreased from 0.2 to 0.8 bar, the oscillation intensity of the two-phase flow gradually increased and the oscillating flow pattern developed from slug to intermittent mist flow. Furthermore, computational fluid dynamics simulations were performed to analyze the heat transfer mechanisms of oscillating cooling and optimize the structure of the piston. The results showed that the heat transfer coefficient of oscillating cooling can be decomposed into three components, namely conventional flow, oscillation enhanced and jet flushing parts. The conventional flow part accounts for the largest proportion of the total heat transfer coefficient, whose value is positively correlated with oil flow rate. The oscillation enhanced part is optimal when the oil filling rate is between 40 % and 60 %. The jet flushing part can lead to a higher temperature gradient in the outer cooling gallery and thus increased thermal stress. The cooling capacities of inner and outer galleries are mainly influenced by the diameters of the oil outlet and distributor pipes, respectively. The jet flushing part can be controlled by changing the inclination angle of the distributor pipe. Finally, optimization achieved an improvement of 8.3 % in cooling capacity of inner cooling gallery and a reduction of 27.5 % of heat dissipation inhomogeneity in the piston.

Phenomenon of reverberation intensity oscillation observed in shallow-water reverberation experiment

Reverberation has always been a disturbance for active sonar. In the Yellow Sea reverberation experiment conducted in July 2014, the motion of the soliton internal wave was recorded by the temperature sensor attached to the vertical array. The effect of the soliton internal wave on reverberation was investigated from both the temperature chain and the reverberation static experiment. Through the processing of the data in this section, it is found that the soliton internal wave can cause the generation of clutter and the oscillation phenomenon of the average reverberation intensity, and the subsequent research can realize the acoustic telemetry to forecast the internal wave and further improve the performance of the active sonar.

Vortex dynamics and boundary layer transition in flow around a rectangular cylinder with different aspect ratios at medium Reynolds number

The numerical investigation focuses on the flow patterns around a rectangular cylinder with three aspect ratios ( $L/D=5$ , $10$ , $15$ ) at a Reynolds number of $1000$ . The study delves into the dynamics of vortices, their associated frequencies, the evolution of the boundary layer and the decay of the wake. Kelvin–Helmholtz (KH) vortices originate from the leading edge (LE) shear layer and transform into hairpin vortices. Specifically, at $L/D=5$ , three KH vortices merge into a single LE vortex. However, at $L/D=10$ and $15$ , two KH vortices combine to form a LE vortex, with the rapid formation of hairpin vortex packets. A fractional harmonic arises due to feedback from the split LE shear layer moving upstream, triggering interaction with the reverse flow. Trailing edge (TE) vortices shed, creating a Kármán-like street in the wake. The intensity of wake oscillation at $L/D=5$ surpasses that in the other two cases. Boundary layer transition occurs after the saturation of disturbance energy for $L/D=10$ and $15$ , but not for $L/D=5$ . The low-frequency disturbances are selected to generate streaks inside the boundary layer. The TE vortex shedding induces the formation of a favourable pressure gradient, accelerating the flow and fostering boundary layer relaminarization. The self-similarity of the velocity defect is observed in all three wakes, accompanied by the decay of disturbance energy. Importantly, the decrease in the shedding frequency of LE (TE) vortices significantly contributes to the overall decay of disturbance energy. This comprehensive exploration provides insights into complex flow phenomena and their underlying dynamics.

Pump-stopping-induced hydraulic oscillations in long-distance district heating system: Modelling and a comprehensive analysis of critical factors

To achieve low-carbon goals, long-distance district heating (LDH) systems are widely applied, which utilize the waste heats from renewable energy power plants and waste incineration plants. The hydraulic oscillations induced by pump-stopping in LDH systems may give rise to critical incidents, such as pipeline ruptures and heating disruptions. The systematic analysis of hydraulic oscillations induced by pump-stopping and the subsequent revelation of critical factors that mitigate their intensity are fundamental prerequisites for ensuring the safety and reliability of LDH system. In this paper, firstly, a pump-stopping hydraulic transient model of the LDH system is established based on the distributed parameter method. Then, combining the method of characteristics (MOC) and the Newton iteration method, the numerical solution scheme for this model is proposed. Finally, taking a 20 km LDH system as an example, the characteristics of hydraulic oscillations induced by pump-stopping is analyzed, and the key factors influencing the intensity of these hydraulic oscillations are investigated. In addition, the differences between hydraulic oscillations induced by pump-stopping and those induced by sudden faults on pipelines are also analyzed. The results indicate that an optimal pump inertia moment value and a properly sized bypass pipe at the pump significantly mitigate hydraulic oscillations, while the influence of constant pressure at the heat source inlet on hydraulic oscillations can be negligible. Our study has important reference value for improving the safe and reliable operation of LDH systems.

Research on the influence factors of hydraulic oscillator on drag reduction efficiency in horizontal well drilling

The application of hydraulic oscillator can solve the problems of large friction and serious pressure support in the drilling process of extended reach well and long horizontal well, but the conventional hydraulic oscillator parameter setting is not reasonable, high pressure consumption, easy to make the ground machine pump overload operation and failure. Therefore, based on the LuGre dynamic friction model theory, a calculation model for the friction between drill string and rock wall under the condition of longitudinal oscillation is established, and the influence law of hydraulic oscillator oscillation parameters on drag reduction efficiency is studied. The results show that the change of oscillation intensity, oscillation frequency and oscillation amplitude can improve the drag reduction efficiency, and the increase trend of drag reduction efficiency is first increased and then decreased. The selection of each parameter has the optimal value. According to the analysis of orthogonal test results, it is found that the reduction degree of friction between drill string and borehole wall rock by each parameter is in the order of oscillation intensity, oscillation frequency and oscillation amplitude from the largest to the smallest.

Interannual variability of mid-high-latitude intraseasonal oscillation intensity at the southern hemisphere during austral summer

Interannual variability of mid-high-latitude intraseasonal oscillation intensity at the southern hemisphere during austral summer

A computational approach to predict and enhance the sensitivity of X-ray resonant magnetic reflectometry to the magnetic behavior of deeply buried interfaces

In the present work a computational approach is applied to model and predict the results of X-ray resonant magnetic reflectometry – a non-destructive synchrotron-based technique to probe chemical composition, crystallographic environment and magnetization in multilayer epitaxial heterostructures with nanoscale depth resolution. The discussed 2D mapping approach is a step forward with respect to conventional resonant X-ray reflectometry and consists of collecting a fine step array of reflected intensity as a function of grazing angle and photon energy across the absorption edge of a particular chemical element. With the use of circularly polarized photons the method can be extended to magnetic systems to produce a map of dichroic reflectance directly related to the magnetization profile of the heterostructure. Studying the magnetic field dependence of dichroic reflectance maps can provide valuable information on the magnetization reversal of individual sublayers of a multilayer heterostructure. In the present paper modeling is performed for a bilayer system mimicking the behavior of a 30 nm ɛ-Fe2O3 thin film that is known to exhibit a pronounced two-component magnetic hysteresis. A technique to find optimal energy/angle combinations in order to sense magnetization of individual sublayers is proposed. Also discussed is the advantage of heavy-element capping, which leads to a substantial increase of the dichroic intensity oscillation contrast in the pre-edge region where the sensitivity to the magnetic behavior of the deeply buried interfaces is most pronounced.

Hybrid simulation of instabilities in capacitively coupled RF CF4/Ar plasmas driven by a dual frequency source

Instabilities in capacitively coupled Ar/CF4 plasma discharges driven by dual frequency sources are investigated using a one-dimensional fluid/electron Monte Carlo hybrid model. Periodic oscillations of the electron density and temperature on the timescale of multiple low frequency (LF) periods are observed. As the electron density increases, an intense oscillation of the electron temperature within each high frequency (HF) period is initiated. This causes a fluctuation of the electron density and results in a discharge instability. This phenomenon is consistent with the discharge behavior observed in scenarios with single-frequency (SF) sources, as reported by Dong et al (2022 Plasma Sources Sci. Technol. 31 025006). However, unlike the SF case, plasma parameters such as the electron density, electric field, electron power absorption and ionization rate exhibit not only periodic fluctuations but also a spatial asymmetry under the influence of the dual-frequency source. This spatial asymmetry leads to a non-uniform distribution of the electron density between the electrodes, which is related to a spatially asymmetric electric field, electron heating, and ionization around a region of minimum electron density (inside the bulk). This region of minimum electron density is shifted back and forth through the entire plasma bulk from one electrode to the other within multiple LF period. The above phenomena are related to superposition effect between the instabilities and the dual-frequency source. Moreover, the time averaged electric field influences the spatio-temporal evolution of ion fluxes. The ion fluxes at the electrodes, which play an important role in etching processes, are affected by both the high and LF components of the driving voltage waveform as well as the observed instabilities. As the HF increases, the electronegativity and electron temperature are reduced and the electron density increases, resulting in a gradual disappearance of the instabilities.

Identification of Intensive Oscillation Parameters of Power System Based on CS-CNN-GRU Network

Broadband oscillations caused by power electronic devices contain a large number of high harmonics and interharmonics, and the frequencies of these harmonic/interharmonic components are likely to be very close. In order to improve the accuracy of dense oscillation parameter identification, this paper proposes a cyclic gate unit-convolutional neural network (CNN-GRU) dense oscillation parameter identification method optimized by the cuckoo search algorithm (CS). The Cuckoo search algorithm combined with the cross-entropy function realizes the automatic acquisition of the parameters of the CNN-GRU network and improves the adaptive ability of the model. The analysis of the measured data shows that the method in this paper has high accuracy in identifying the dense oscillation parameters of dense frequency harmonics, interharmonics, and neighboring fundamental interharmonic pairs, etc., and the computational speed is also improved compared with the traditional deep learning network.

The Effects of Parameter Settings on Triggering Time and Climb Rate during Lean-Premixed Combustion Thermoacoustic Oscillations

This study theoretically explored the effects of parameter settings on thermoacoustic oscillations with a low-order model. Three factors were explored—combustor length, inlet gas temperature and thermal power. The research findings indicate that optimizing the parameter settings can yield better thermoacoustic oscillation suppression results. The sound pressure amplitude decreased from 3.2 × 105 Pa to 2.1 × 105 Pa as the combustor length increased from 1.2 m to 6.0 m. The triggering time increased from 0.32 s to 0.91 s when the combustion chamber length increased. The climb rate declined from 23.38 × 105 Pa/s to 3.75 × 105 Pa/s when the combustor length was elongated. The sound pressure amplitude decreased from 3.44 × 105 Pa to 2.4 × 105 Pa as the gas temperature rose from 0 to 100 °C. The triggering time and climb rate variation tendency were similar when the gas temperature increased—both declined as the gas temperature rose. The sound pressure amplitude experienced a slight fluctuation when the thermal power rose. However, the triggering time decreased from 0.26 s to 0.043 s when the thermal power improved. The climb rate increased from 18.72 × 105 Pa/s to 27.65 × 105 Pa/s when the thermal power rose. The oscillation frequency presented was completely different in three cases that had different wavelengths and oscillation intensities. The triggering time and climb rate fluctuated extensively in varying conditions, and the above two factors were interrelated and contradictory to each other when thermoacoustic oscillation was excited. This study explored parameters’ effects on triggering time and climb rate, thereby providing references for constructing a model-based control system for thermoacoustic oscillation feedback control.

Modeling and Analysis of a Thermo-Electro-Magneto-Viscoelastic Actuator

This study examines the dynamics of a smart actuator known as a thermo-electro-magneto-viscoelastic plate. This type of actuator is made of an electro-magneto-active polymer filled with suitably added fillers, which enhance its performance. An energy-based dynamic model is derived to analyze the behavior of the actuator, taking into account thermal, electrical, and mechanical factors. The study uses a well-known Euler–Lagrange equation to develop the governing dynamic equation of the actuator. In addition, the stability, periodicity, and resonant behavior of the actuator are evaluated through time-evolution diagrams, Poincaré plots, and phase portraits. The results indicate that the nonlinear behavior of the actuator changes when subjected to different modes of actuation such as DC and AC. Increasing particle reinforcement and reducing temperature change strengthens the polymer and depletes the level of deformation. The intensity of oscillation also decreases, and the excitation frequency increases as the temperature change reduces and the volume fraction increases. Lastly, this study provides important insights into the design of modern actuators, which have numerous potential applications in both engineering and medical industries.

The impact of an ultrasonic standing wave on the sorption behavior of proteins: Investigation of the role of acoustically induced non-spherical bubble oscillations

HypothesisProtein molecules adsorb on the air/liquid interface due to possessing a hydrophobic side. A full surface coverage is important in many processes such as in protein harvesting by foam fractionation. The adsorption of proteins in low concentration solutions is preceded by a relatively long time lag known as the induction period. This has been attributed to the formation of an adsorbed monolayer, which relies on the reorientation of the protein molecules. The reduction of the induction period can significantly facilitate the sorption process to reach full protein coverage. For this purpose acoustically induced non-spherical bubble oscillations can aid in the formation of the monolayer and enhance the sorption process. ExperimentIn this study, low frequency ultrasound was used to induce non-spherical oscillations on an air bubble attached to a capillary. Profile analysis tensiometry was deployed to examine the effect of these non-spherical oscillations on the sorption dynamics of different proteins. FindingsWe observed that during the initial stages of adsorption, when the bubble surface is almost empty, non-spherical oscillations occur, which were found to significantly expedite the adsorption process. However, during later stages of the adsorption process, despite the continued presence of several sonication phenomena such as the primary radiation force and acoustic streaming, no change in adsorption behavior of the proteins could be noted. The occurrence, duration, and intensity of the non-spherical bubble oscillations appeared to be the sole contributing factors for the change of the sorption dynamics of proteins.

Theoretical study on the photophysical properties of thiophene-fused-type BODIPY series molecules in fluorescence imaging and photodynamic therapy.

As a class of photosensitizers (PSs) with dual functions of photodynamic therapy (PDT) and fluorescence imaging, the relationship between the structure and dual-function of thiophene-fused-type BODIPY dyes has not been studied in depth before. We found that the thiophene-fused-type BODIPY triplet photosensitizer is produced according to the energy level matching rule and the introduction of the thiophene ring significantly reduces the energy gap ΔEST between singlet and triplet states, as revealed by our investigation of the excited state structures and energies of thieno-fused BODIPY dyes. At the same time, a tiny ΔEST also results in a greatly enhanced intersystem crossing (ISC) rate, kISC. The kISC value of MeO-BODIPY, having the highest singlet oxygen quantum yield (ΦΔ), is the largest. Substitution with a strong electron donor N,N-dimethylaminophenyl (DMA) leads to the vertical configuration in the T1 state. The small ΔE (0.0029 eV) between the HOMO and HOMO-1 triggers the photo induced electron transfer (PET) of inhibiting ISC and fluorescence. When thieno-fused BODIPYs react with pyrrole, the increase of π-conjugation and smaller ΔEHOMO-LUMO explain the redshift in emission wavelength of thieno-pyrrole-fused BODIPY. The more planar configuration of the S1 state and the stronger oscillator intensity reflect a higher fluorescence quantum yield (ΦF). The extension of π-conjugation can cause molecules to transition to higher-level singlet excited states (Sn states, n ≥ 1) after absorbing energy and reduce the energy level of the excited state, resulting in multiple channels and favoring 1O2 production for thieno-pyrrole-fused BODIPYs with electron-withdrawing groups at the para-position of the phenyl groups. Due to ΔES0-T1 < 0.980 eV, the substitution of electron-donating groups cannot produce 1O2. In this work, we have revealed the mechanism of ISC and the fluorescence emission process in the thiophene-fused-type BODIPY dye, which has provided a theoretical foundation and guidance for the future design of BODIPY-based heavy-atom-free PSs for molecular applications in PDT.

Multi-Scale Periodic Variations in Soil Moisture in the Desert Steppe Environment of Inner Mongolia, China

Uncovering the complex periodic variations in soil moisture can provide a significant reference for climate prediction and hydrological process simulation. We used wavelet analysis to quantify and identify the multi-scale periodic variations of soil moisture in the desert steppe of Inner Mongolia from 2009 to 2019 and analyzed the differences between the multi-scale periodic changes in soil moisture at the bottom (BS) and upper slope (US). The results show that the interannual variability in soil moisture at the BS and US has a significant upward trend. Moreover, the amount and volatility decrease with the increase in soil depth in the range of 0–30 cm, and the soil moisture at the BS is 26.4% higher than the US. The soil moisture has periodic changes of 0.5a, 1a, 2a, 3a and 3.5a in the desert steppe environment of Inner Mongolia. The periodic structure and intensity of different slope positions are greatly different. Soil moisture at the US has more complex multi-scale periodic changes, and the periodic oscillations of 3.5a, 3a, and 1a are dominant. The periodic oscillations of 0.5a and 1a are dominant at the BS. At the BS, the intensity of periodic oscillation of 1a after January 2015 has weakened. The weakening of soil moisture by temperature, rainfall and soil temperature caused the change in the multiple time-scale periodic oscillation of soil moisture.