Stress-distribution around a rigid line in dissimilar media

The elastic plane problem of a rigid line inclusion between two dissimilar media was considered. By solving the Riemann-Hilbert problem, the closed-form solution was obtained and the stress distribution around the rigid line was investigated. It was found that the modulus of the singular behavior of the stress remains proportional to the inverse square root of the distance from the rigid line end, but the stresses possess a pronounced oscillatory character as in the case of an interfacial crack tip.

The plane problem of collinear rigid lines under arbitrary loads

The elastic plane problem of collinear rigid lines under arbitrary loads is dealt with. Applying the Riemann-Schwarz symmetry principle integrated with the analysis of the singularity of complex stress functions, the general formulation is presented, and the closed-form solutions to several problems of practical importance are given, which include some published results as the special cases. Lastly the stress distribution in the immediate vicinity of the rigid line end is examined.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

A unified law of motion of lithospheric plates and the driving mechanism

From observed data on lithospheric plates, a unified empirical law for plate motion,valid for continental as well as oceanic plates, is obtained in the following form: The speedof plate motion U depends linearly on a geometric parameter T_d, ratio of the sum of effectiveridge length and trench arc length to the sum of area of continental part of plate and total areaof cold sinking slab. Based on this unified law, a simple mechanical analysis shows that, themain driving forces for lithospheric plates come from push along the mid-ocean ridge andpull by the cold sinking slab, while the main drag forces consist of the viscous traction beneaththe continental part of plate and over both faces of the sinking slab. Moreover, the specific-push along ridge and pull by slab are found to be of equal magnitude.

CMP润滑方程的Chebyshev加速超松弛法求解

化学机械抛光（chemicalmechanicalpolishing，CMP）是一项融合了化学分解和机械力学的工艺技术，其中包含了流体动力润滑的作用．文中通过分析流场效应，建立了Sum在化学机械抛光实验基础上获得的润滑（Reynolds）方程，利用Chebyshev加速超松弛技术对润滑方程进行求解，得到了新的抛光液压力分布模型，给出了一些新的载荷及转矩在不同参数下的变化情况．

IP Simulation of Micro Gas Flows under 3-D Head Sliders

Gas film lubrication of a three-dimensional flat read-write head slider is calculated using the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method, respectively. The pressure distributions on the head slider surface at different velocities and flying heights obtained by the two methods are in excellent agreement. IP method is also employed to deal with head slider with three-dimensional complex configuration. The pressure distribution on the head slider surface and the net lifting force obtained by the IP method also agree well with those of DSMC method. Much less (of the order about 10(2) less) computational time (the sum of the time used to reach a steady stage and the time used in sampling process) is needed by the IP method than the DSMC method and such an advantage is more remarkable as the gas velocity decreases.

计算流体力学

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

The effect of the incident relative phase on the four-wave mixing field and the electromagnetically induced transparency

In a A-type system employing a two-photon pump field, a four-wave mixing field can be generated simultaneously and, hence, a closed-loop system forms. We study theoretically the effect of the relative phase between the two incident fields on the generated four-wave mixing field and the electromagnetically induced transparency. It is found that the phase of the generated four-wave mixing field is the sum of the incident relative phase and a fixed phase that is irrelative to the incident relative phase. Hence, the total phase of the closed-loop system is independent of the incident relative phase. As a result, the incident relative phase has no effect on the electromagnetically induced transparency, which is different from the case of a A-type loop system closed by a third incident field. (c) 2005 Pleiades Publishing, Inc.

Sideband manipulation of population inversion in a three-level Λ atomic system

Sideband manipulation of population inversion in a three-level A atomic configuration is investigated theoretically. Compared with the case of a nearly monochromatic field, a population inversion between an excited state and a ground state has been found in a wide sideband intensity range by increasing the difference in frequency between three components. Furthermore, the population inversion can be controlled by the sum of the relative phases of the sideband components of the trichromatic pump field with respective to the phase of the central component. Changing the sum phase from 0 to pi, the population inversion between the excited state and the ground state can increase within nearly half of the sideband intensity range. At the same time, the sideband intensity range that corresponds to the system exhibiting inversion rho(00) > rho 11 also becomes wider evidently.

Trichromatic manipulation of Kerr nonlinearity in a three-level lambda atomic system

Trichromatic manipulation of Kerr nonlinearity in a three-level A atomic configuration is investigated theoretically. It is shown that for a weak monochromatic probe field, the enhanced Kerr nonlinearity can be achieved in multiple separate transparent windows due to interference effect of multiple two-photon Raman channels. Furthermore, the property of Kerr nonlinearity can be controlled by the sum of the relative phases of the sideband components of the trichromatic pump field compared to the central component.

Coherent enhancement of broadband frequency up-conversion in BBO crystal by shaping femtosecond laser pulses

An optimal feedback control of broadband frequency up-conversion in BBO crystal is experimentally demonstrated by shaping femto-second laser pulses based on genetic algorithm, and the frequency up-conversion efficiency can be enhanced by similar to 16%. SPIDER results show that the optimal laser pulses have shorter pulse-width with the little negative chirp than the original pulse with the little positive chirp. By modulating the fundamental spectral phase with periodic square distribution on SLM-256, the frequency up-conversion can be effectively controlled by the factor of about 17%. The experimental results indicate that the broadband frequency up-conversion efficiency is related to both of second harmonic generation (SHG) and sum frequency generation (SFG), where the former depends on the fundamental pulse intensity, and the latter depends on not only the fundamental pulse intensity but also the fundamental pulse spectral phase. (c) 2006 Elsevier B.V. All rights reserved.