The mechanism of energy storage by shearing the solar magnetic field

In this paper, a complete set of MHD equations have been solved by numerical calculations in an attempt to study the dynamical evolutionary processes of the initial equilibrium configuration and to discuss the energy storage mechanism of the solar atmosphere by shearing the magnetic field. The initial equilibrium configuration with an arch bipolar potential field obtained from the numerical solution is similar to the configuration in the vicinity of typical solar flare before its eruption. From the magnetic induction equation in the set of MHD equations and dealing with the non-linear coupling effects between the flow field and magnetic field, the quantitative relationship has been derived for their dynamical evolution. Results show that plasma shear motion at the bottom of the solar atmosphere causes the magnetic field to shear; meanwhile the magnetic field energy is stored in local regions. With the increase of time the local magnetic energy increases and it may reach an order of 4×10^25 J during a day. Thus the local storage of magnetic energy is large enough to trigger a big solar flare and can be considered as the energy source of solar flares. The energy storage mechanism by shearing the magnetic field can well explain the slow changes in solar active regions.

Optimizations for the uniformity of the non-volatile volume grating in doubly doped LiNbO3 crystals

The formation of the non-uniformity of the non-volatile volume grating in doubly doped LiNbO3 crystals is studied in detail. We find that the non-uniformity of the grating is mainly caused by strong ultraviolet light absorption, and the average saturation space-charge field is small and the diffraction efficiency is low as a result of the non-uniformity of the grating. In order to optimize the uniformity of the grating, we propose the recording scheme by using two sensitizing beams simultaneously from the two opposite sides of the crystals. Theoretical simulations and experimental verifications are performed. Results show that the well uniformed grating with high diffraction efficiency can be obtained by using this optimization scheme. (c) 2004 Elsevier B.V. All rights reserved.

XPM-induced pulse compression in birefringent dispersion decreasing fiber

The pulse compression induced by cross-phase modulation in birefringent dispersion decreasing fiber is discussed theoretically by solving the coupled Schrodinger equations which include the contribution of the high-order non-linear effects, and third-order dispersion. In particular, it is found that a high quality compressed signal pulse can be obtained by a pump pulse of low intense through the technique. The dependence of optimum compression on the non-linear factor N, time delay tau(d) and the dispersive ratio f is also discussed in detail. (c) 2005 Elsevier B.V. All rights reserved.

Temperature-Dependent Growth and Life Cycle of Nemoura sichuanensis (Plecoptera: Nemouridae) in a Chinese Mountain Stream

Plecoptera constitute a numerically and ecologically significant component in mountain streams all over the world, but little is known of their life cycles in Asia. The life cycle of Nemoura sichuanensis and its relationship to water temperature was investigated during a 4-year study in a headwater stream (known as the Jiuchong torrent) of the Xiangxi River in Central China. Size structure histograms suggest that the life cycle was univoltine, and the relationships between the growth of Nemoura sichuanensis, physiological time, and effective accumulated water temperature were described using logistic regressions. The growth pattern was generally similar within year classes but growth rates did vary between year-classes. Our field data suggest a critical thermal threshold for emergence in Nemoura sichuanensis, that was close to 9 degrees C. The total number of physiological days required for completing larval development was 250 days. The effective accumulated water temperature was 2500 degree-days in the field. Development during the life cycle increased somewhat linearly with the physiological time and the effective accumulated water temperature, but some non-linear relationships were best developed by logistic equations.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.