Micromechanics analysis of particulate-reinforced composites and their failure mechanisms

Stress fields and failure mechanisms have been investigated in composites with particles either surface treated or untreated under uniaxial tension. Previous experimental observation of failure mechanisms in a composite with untreated particles showed that tensile cracks occurred mostly at the polar region of the particle and grew into interfacial debonding. In a composite with surface-treated particles, however, shear yielding and shear cracking proceeded along the interphase-matrix interface at the polar area of the matrix and thus may improve the mechanical behaviour of the material. The finite element calculations showed that octahedral shear stress at the polar and longitudinal areas of the particle treated by coupling agents is much larger than that of materials with untreated particles, and the shear stress distribution around the interface is sensitive to the interphase property. The results suggest that a th ree-phase model can describe the composites with surface-treated fillers.

The influence of momentum addition on the acceleration of the solar and stellar winds

The influence of the momentum addition, which may be associated with the average or fluctuation transverse component of the magnetic field or others, on the acceleration the solar wind or stellar wind is studied in a local streamtube. The results show that the larger the momentum addition the stronger the acceleration of the wind. For example, if the typical transverse magnetic field is about 0.1 of the longitudinal field, the velocity of the solar wind at 1 AU may be increased by 40%. The coronal hole may be considered as a streamtube, the presence of a high stream from the coronal hole may be explained by the existence of an average or fluctuation transverse magnetic field in the streamtube. A similar conclusion may be applied to the polar region, where the velocity of the solar wind will be larger than elsewhere as if there is a transverse component of magnetic field, as well as to the stellar wind. The influence of other parameters on the acceleration of the solar wind is also discussed. From the viewpoint of the solar wind mechanism, the present paper shows that the momentum addition in the subsonic flow region can increase the velocity of the solar wind at 1 AU.

2D NMR and FT-Raman spectroscopic studies on the interaction of lanthanide ions and Ln-DTPA with phospholipid bilayers

The interactions of lanthanide ions and the Ln-DTPA (DTPA = diethylenetriaminepentaacetate) complex with di palmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylethanolamine (DPPE) bilayers are studied by 2D NOESY and FT-Raman spectroscopy. Proton NMR spectroscopic results show that lanthanide ions combine with phosphate groups in the polar region of the outer layer of DPPC liposomes, leading to the separation in chemical shift of the proton signal of N(CH3)(3) The conformational change of the O-C-C-N+ backbone from the gauche conformer to the trans one is not found; i.e., the orientation of the polar headgroup is still parallel to the surface of the bilayers. The Ln-DTPA complex at low concentration in a pH 7.4 solution localizes far away from bilayers and thereby has little effect on the structure of bilayers. The FT-Raman spectroscopic results indicate that lanthanide ions affect strongly the fluidity of acyl chains of DPPE bilayers while the Ln-DTPA complex affects it slightly.

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

AE指数的可靠性研究

Auroral electrojet index (AE) are usually used to quantitatively describe the activity of the geomagnetic field in the polar region. AE is a means to identify the level of a substorm as well. The auroral electrojet indices (AU, AL, and AE) have served well for more than four decades as measures of magnetospheric substorm activity. However, as substorm studies have progressed considerably during the past several years, the accuracy of the present electrojet indices have become an important issue. Thus it is fortunate for us to reexamine and evaluate the accuracy of the present electrojet indices and improve them if necessary. For a better use of the present indices and for future improvement we examine the limitations of the auroral electrojet indices as an accurate quantitative measure of the auroral electrojets and magnetospheric substorms. Some of the limitations arise from the data availability and also from the present simplified scheme in deriving them, but some of them originate in the definition themselves. In the present paper, we analyze the characteristics of aurora, ionospheric current system and the AE index. It is noted that for the pseudo-breakup events the variations of AE ( > 500 nT) clearly show the sudden increase and slow decrease phases. However, aurora does not expand to poleward or equatorward, and the ionospheric currents presents the features of the magnetic convection. We mainly focus on investigating what time the AE could be used to identify the subtorm.

Cyclotron resonance mass of magnetopolaron in a polar crystal slab at finite temperatures

Using the Green function method, we have studied the cyclotron resonance of an electron interacting with bulk longitudinal optical(BO) phonons as well as surface optical(SO) phonons in a polar crystal slab at finite temperatures. It is found that the temperature dependence of magnetopolaron depends strongly on the strength of the magnetic field. The numerical results show that the cyclotron resonance mass of polaron in a slab is an increasing or decreasing function of temperature when the magnetic field is lower or higher than the resonant magnetic field region, respectively. The magnetic field and slab width dependence of cyclotron resonance mass are also studied in this paper. (C) 1999 Elsevier Science B.V. All rights reserved.

Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Flow Characteristics of Non-Polar Organic Liquids with Small Molecules in a Microchannel

用去离子水及有机液体在内径约为25μm的石英圆管内进行了流量特性实验.液体分子量范围为18～160,动力黏性系数的范围为0.5～1 mPa.s.实验雷诺数范围为Re<8.所用有机液体为:四氯化碳、乙基苯及环己烷都是非极性液体,其分子结构尺度小于1 nm.实验结果表明,在定常层流条件下,圆管内的液体流量与两端压力差成正比,其压力-流量关系仍符合经典的Hagen-Poiseuille流动.这说明非极性小分子有机液体在本实验所用微米尺度管道中其流动规律仍符合连续介质假设.鉴于微尺度流动实验的特殊性,文中还介绍了微流动实验装置,分析了微尺度流动测量误差来源及提高测量精度的措施.

A generalized wave action conservative equation for the dissipative dynamical system in the nearshore region

To describe the various complex mechanisms of the dissipative dynamical system between waves, currents, and bottoms in the nearshore region that induce typically the wave motion on large-scale variation of ambient currents, a generalized wave action equation for the dissipative dynamical system in the nearshore region is developed by using the mean-flow equations based on the Navier-Stokes equations of viscous fluid, thus raising two new concepts: the vertical velocity wave action and the dissipative wave action, extending the classical concept, wave action, from the ideal averaged flow conservative system into the real averaged flow dissipative system (that is, the generalized conservative system). It will have more applications.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

On the magnetospheric structure of pulsars

The magnetospheric structure of a pulsar is discussed for a non-force-free magnetic field. The local solution to the axisymmetric equations of the pulsar is obtained by the method of expanding in the polar angle. Particular attention is given to the solutions near the polar axis and the equator. Near the pulsar surface, the magnetic field energy density is found to be larger than the other energy components; the gravitational potential and the kinetic energy are relatively larger far away from the pulsar surface. It is shown that these relations influence the mass distribution in the pulsar magnetosphere. The results also show that the plasma rotation may be nonrigid and, hence, that a corotational region with a closed magnetic field may not exist.