Permeability of Fractal Porous Media by Monte Carlo Simulations

The permeability of the fractal porous media is simulated by Monte Carlo technique in this work. Based oil the fractal character of pore size distribution in porous media, the probability models for pore diameter and for permeability are derived. Taking the bi-dispersed fractal porous media as examples, the permeability calculations are performed by the present Monte Carlo method. The results show that the present simulations present a good agreement compared with the existing fractal analytical solution in the general interested porosity range. The proposed simulation method may have the potential in prediction of other transport properties (such as thermal conductivity, dispersion conductivity and electrical conductivity) in fractal porous media, both saturated and unsaturated.

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

Diffusion dominated process for the crystal growth of a binary alloy

The pure diffusion process has been often used to study the crystal growth of a binary alloy in the microgravity environment. In the present paper, a geometric parameter, the ratio of the maximum deviation distance of curved solidification and melting interfaces from the plane to the radius of the crystal rod, was adopted as a small parameter, and the analytical solution was obtained based on the perturbation theory. The radial segregation of a diffusion dominated process was obtained for cases of arbitrary Peclet number in a region of finite extension with both a curved solidification interface and a curved melting interface. Two types of boundary conditions at the melting interface were analyzed. Some special cases such as infinite extension in the longitudinal direction and special range of Peclet number were reduced from the general solution and discussed in detail.

Asymptotic solutions to a class of nonlinear oscillation systems

In this paper, we present an asymptotic method for the analysis of a class of strongly nonlinear oscillators, derive second-order approximate solutions to them expressed in terms of their amplitudes and phases, and obtain the equations governing the amplitudes and phases, by which the amplitudes of the corresponding limit cycles and their behaviour can be determined. As an example, we investigate the modified van der Pol oscillator and give the second-order approximate analytical solution of its limit cycle. The comparison with the numerical solutions shows that the two results agree well with each other.

大洋中脊的形成过程

本文从热传导方程出发,得到了大洋中脊下岩石层温度分布的分析表达式及数值计算结果。结果表明,软流层上涌流动所提供的热源可以使大洋中脊下岩石层逐步融化;岩石层的相对移动速度对大洋中脊岩石层温度场及融化深度影响较大。 

圆弧形沉积场地对平面瑞利波散射的解析分析

分别采用Biot饱和多孔介质动力学理论和单相介质弹性动力学理论模拟层状沉积谷场地和周围场地,利用波函数展开法,得出了层状饱和沉积场地对Rayleigh波散射的解析解.进行参量分析,研究了无量纲入射波频率、层状沉积层排列顺序及相对刚度和厚度等因素对Rayleigh波散射效应的影响

Diffusive waves in a channel with concentrated lateral inflow

In the current paper an analytical solution for diffusive wave equation with the concentrate-distributed lateral inflow is yielded. Finite-difference numerical method is also employed to validate this model. The backwater effects drawn from lateral inflow on the mainstream are examined finally.

A Breakdown Criterion of Free Molecular Flows and an Optimum Analysis of EBPVD

Two important issues in electron beam physical vapor deposition (EBPVD) are addressed. The first issue is a validity condition of the classical cosine law widely used in the engineering context. This requires a breakdown criterion of the free molecular assumption on which the cosine law is established. Using the analytical solution of free molecular effusion flow, the number of collisions (N-c) for a particle moving from an evaporative source to a substrate is estimated that is proven inversely proportional to the local Knudsen number at the evaporation surface. N-c = 1 is adopted as a breakdown criterion of the free molecular assumption, and it is verified by experimental data and DSMC results. The second issue is how to realize the uniform distributions of thickness and component over a large-area thin film. Our analysis shows that at relatively low evaporation rates the goal is easy achieved through arranging the evaporative source positions properly and rotating the substrate.

A Mechanical Model For The Quantification Of The Effect Of Laser Quenching On Ctod In Steels

The technology of laser quenching is widely used to improve the surface properties of steels in surface engineering. Generally, laser quenching of steels can lead to two important results. One is the generation of residual stress in the surface layer. In general, the residual stress varies from the surface to the interior along the quenched track depth direction, and the residual stress variation is termed as residual stress gradient effect in this work. The other is the change of mechanical properties of the surface layer, such as the increases of the micro-hardness, resulting from the changes of the microstructure of the surface layer. In this work, a mechanical model of a laser-quenched specimen with a crack in the middle of the quenched layer is developed to quantify the effect of residual stress gradient and the average micro-hardness over the crack length on crack tip opening displacement (CTOD). It is assumed that the crack in the middle of the quenched layer is created after laser quenching, and the crack can be a pre-crack or a defect due to some reasons, such as a void, cavity or a micro-crack. Based on the elastic-plastic fracture mechanics theory and using the relationship between the micro-hardness and yield strength, a concise analytical solution, which can be used to quantify the effect of residual stress gradient and the average micro-hardness over the crack length resulting from laser quenching on CTOD, is obtained. The concise analytical solution obtained in this work, cannot only be used as a means to predict the crack driving force in terms of the CTOD, but also serve as a baseline for further experimental investigation of the effect after laser-quenching treatment on fracture toughness in terms of the critical CTOD of a specimen, accounting for the laser-quenching effect. A numerical example presented in this work shows that the CTOD of the quenched can be significantly decreased in comparison with that of the unquenched. (C) 2008 Elsevier B.V. All rights reserved.

Lasing without or with inversion in an open four-level system with a phase-fluctuation driving field

The effect of exit rate and the ratio of atomic injection rate on gain behaviour has been investigated, and the effects of phase fluctuation on absorption, dispersion and population difference in an open four-level system have been analysed by using numerical simulation from the steady linear, analytical solution. The variation of the linewidth, Rabi frequency of the driving field, the exit rate or the ratio of atomic injection rate can change the lasing properties in the open system. The presence of finite linewidth due to driving-field phase fluctuation prevents the open four-level atomic system from obtaining a high refractive index along with zero absorption.

Spatiotemporal behaviors and singularity of ultrashort pulsed Elegant Hermite-Gaussian beams

Analytical propagation expressions of ultrashort pulsed Elegant Hermite-Gaussian beams are derived and spatiotemporal properties of the pulses with different transverse modes are studied. Singularity of the complex amplitude envelope solution of the pulses obtained under slowly varying envelope approximation is analyzed in detail. The rigorous analytical solution of the pulse is deduced and no singularity emerges in the solution. The obtained results indicate that the transverse mode affects not only the spatiotemporal properties but also the singularity of the pulses. Time delay of the off-axis maximum intensity is more obvious and the singularity is located nearer to the z-axis for the pulse with higher transverse modes. (C) 2007 Elsevier GmbH. All rights reserved.

Mode analysis and Q-factor enhancement due to mode coupling in rectangular resonators

Modes in rectangular resonators are analyzed and classified according to symmetry properties, and quality factor (Q-factor) enhancement due to mode coupling is observed. In the analysis, mode numbers p and q are used to denote the number of wave nodes in the direction of two orthogonal sides. The even and odd mode numbers correspond to symmetric and antisymmetric field distribution relative to the midlines of sides, respectively. Thus, the modes in a rectangle resonator can be divided into four classes according to the parity of p and q. Mode coupling between modes of different classes is forbidden; however, anti-crossing mode coupling between the modes in the same class exists and results in new modes due to the combination of the coupled modes. One of the combined modes has very low power loss and high Q-factor based on far-field emission of the analytical field distribution, which agrees well with the numerical results of the finite-difference time-domain (FDTD) simulation. Both the analytical and FDTD results show that the Q-factors of the high Q-factor combined modes are over one order larger than those of the original modes. Furthermore, the general condition required to achieve high-Q modes in the rectangular resonator is given based on the analytical solution.

Measurement of chirp parameter and modulation index of a semiconductor laser based on optical spectrum analysis

A novel and simple method for measuring the chirp parameter, frequency, and intensity modulation indexes of directly modulated lasers is proposed in a small-signal modulation scheme. A graphical approach is presented. An analytical solution to the measurement of low chirp parameters is also given. The measured results agree well with those obtained using the conventional methods.

Investigation of mode characteristics for microdisk resonators by S-matrix and three-dimensional finite-difference time-domain technique

The mode characteristics of a three-dimensional (3D) microdisk with a vertical refractive index distribution of n(2)/3.4/n(2) are investigated by the S-matrix method and 3D finite-difference time-domain (FDTD) technique. For the microdisk with a thickness of 0.2 mu m. and a radius of 1 mu m, the mode wavelengths and quality factors for the HE7,1 mode obtained by 3D FDTD simulation and the S-matrix method are in good agreement as n(2) increases from 1.0 to 2.6. But the Q factor obtained by the 3D FDTD rapidly decreases from 1.12 X 10(4) to 379 as n2 increases from 2.65 to 2.8 owing to the vertical radiation losses, which cannot be predicted by the proposed S-matrix method. The comparisons also show that quality factors obtained from the analytical solution of two-dimensional microdisks under the effective index approximation are five to seven times smaller than those of the 3D FDTD as n(2) = 1 and R = 1 mu m. (c) 2006 Optical Society of America.