Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Muon density measurements with the KASCADE central detector

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.     

Revisiting Sustainable Development of Dry Valleys in Hengduan Mountains Region

Dry valleys are a striking geographic landscape in Hengduan Mountains Region and are characterized by low rainfall, desert type of vegetation and fragile environment. Past efforts and resources have been concentrated mainly on rehabilitation of degraded ecosystem and fragile environment, particularly reforestation, while socio-economic development has been largely overlooked. Despite successes in pocket areas, the overall trend of unsustainability and environmental deterioration are continuing. It is important to understand that uplift of the Tibetan Plateau is the root cause of development of dry valleys, and development and formation of dry valleys is a natural process. Human intervention has played a secondary role in development of dry valleys and degradation of dry valleys though human intervention in many cases has speeded up environmental degradation of the dry valleys. It is important to understand that dry valleys are climatic enclaves and an integrated approach that combines rehabilitation of degraded ecosystems and socio-economic development should be adopted if the overall goal of sustainable development of dry valleys is to be achieved. Promotion of niche-based cash crops, rural energy including hydropower, solar energy, biogas and fuelwood plantation is recommended as the priority activities.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

有关边界层问题中一类奇异积分方程新的存在性结果

对边界层理论新结果中出现的一类奇异积分方程w（t）=∫1 t（1-s）（λ＋λs＋s）/w（s）ds＋（1-t）∫t 0s/w（s）ds,t∈（0,1）进行讨论，并得出了上述方程在λ∈（-1/2,0）上正解存在性的新结果。     

关于Diophantine方程的解

对于正整数a ,设S(a)是Smarandache函数。证明了 :方程S(1·2 ) +S(2·3) +… +S(x(x +1) ) =S(x(x +1) (x +2 ) /3)仅有正整数解x =1。     

用比电导法研究活性炭的吸附特性

用比电导法研究两种活性炭自稀水溶液中吸附强电解质硫酸铬和硫酸铜的吸附平衡特性。结果表明,在本文的研究条件下,两种活性炭对硫酸铬和硫酸铜都有吸附作用。同时,活性炭吸附硫酸铜的吸附平衡特性可以用 Freundlich 方程式来描述。研究的结果对固体在溶液中的吸附基础理论以及处理工业废水具有一定的意义。     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

The efifect of wave breaking on surface wave imaging by Synthetic Aperture Radar

The effect of ocean wave breaking as a non-Bragg mechanism on backscattering cross-section and modulation transfer functions (MTF) of radar was investigated based on Bragg resonance theory and parametric method. The result showed that the additional effect of wave breaking on backscattering cross-section is not more than 20% except for the small incident angle of VV polarized electromagnetic (e.m.) wave but is significant for HH polarized e.m. wave. Breaking waves lead to increase in the modulus of tilt modulation MTF and the larger the wind speed, the faster the increase. For large incident angle, the modulus of tilt modulation MTF with wave breaking decreases quickly with incident angle for HH polarization and approach to that without wave breaking for VV polarization. The hydrodynamic MTF increases 30%-60% when considering wave breaking and the increase is larger for HH polarization than for VV polarization.