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.

Adiabatic compression of plasma with purely toroidal rotation

A variational method is developed to find approximate solutions to the generalized Grad-Shafranov equations for an adiabatic compression of the plasma with toroidal rotation, via the expansion in Fourier series in poloidal angle of the flux surface coordinates. The numerical results, which are carried out by the present method and by the usual two-dimensional method for a static equilibrium state, agree well.

静动结合的化学腐蚀法制备近场光学探针的理论分析

理论上分析了静动结合的化学腐蚀法制备探针的具体机理及过程。在静态腐蚀的过程中, 利用流体力学Young-Laplace方程的一级近似解获得了光纤插入到HF酸中形成的新月形高度。在动态腐蚀过程中, 详细分析了当静态腐蚀时间和动态腐蚀时间分别取不同值时, 光纤移动速度对光纤探针结构的影响。利用此法可制备出尖端锐利、大锥角或多锥体等各种结构的光纤探针。这为实验上制备出性能优良的探针, 为拓宽扫描近场光学显微镜的应用范围奠定基础。将上述理论分析的结果与本文实验中所得初步结果进行了比较, 所得结果一致。

Kinetics of ytterbium(III) extraction with Cyanex 272 using a constant interfacial cell with laminar flow

Studies have been made on the kinetics of ytterbium(III) with bis-(2,4,4-trimethylpentyl) phosphinic acid (Cyanex 272, HA) in n-heptane using a constant interfacial cell with laminar flow. The stiochiometry and the equilibrium constant of the extracted complex formation reaction between Yb3+ and Cyanex 272 are determined. The extraction rate is dependent of the stirring rate. This fact together with the Ea value suggests that the mass transfer process is a mixed chemical reaction-diffusion controlled at lower temperature, whereas it is entirely diffusion controlled at higher temperature. The rate equations for the ytterbium extraction with Cyanex 272 have been obtained. The rate-determining step is also made by predictions derived from interfacial reaction models, and through the approximate solutions of the flux equation, diffusion parameters and thickness of the diffusion film have been calculated.

Mass transfer kinetics of yttriium(III) using a constant interfacial cell with laminar flow. Part I. Extraction with Cyanex 302

The interfacial tension is measured for Cyanex 302 in heptane and adsorption parameters are calculated according to Gibbs equation and Szyskowski isotherm. The results indicate that Cyanex 302 has a high interfacial activity, allowing easy extraction reaction to take place at the liquid-liquid interface. The extraction kinetics of yttrium(III) with Cyanex 302 in heptane are investigated by a constant interfacial cell with laminar flow. The effects of stirring rate, temperature and specific interfacial area on the extraction rate are discussed. The results suggest that the extraction kinetics is a mixed regime with film diffusion and an aqueous one-step chemical reaction proposed to be the rate-controlling step. Assuming the mass transfer process can be formally treated as a pseudo-first-order reversible reaction with respect to the metal cation, the rate equation for the extraction reaction of yttrium(III) with Cyanex 302 at pH <5 is obtained as follows:R-f = 10(-7.85)[Y(OH)(2)(+)]((a))[H(2)A(2)]((o))(1.00)[H+]((a))(-1.00)Diffusion parameters and rate constants are calculated through approximate solutions of the flux equation.

Kinetics and mechanism of Yb(III) extraction and separation from Y(III) with mixtures of bis(2,4,4-trimethylpentyl)phosphinic acid and 2-ethylhexyl phosphonic acid mono-2-ethylhexyl ester

The ytterbium(III) extraction kinetics and mechanism with mixtures of bis(2,4,4-trimethylpentyl)phosphinic acid (Cyanex272) and 2-ethylhexyl phosphonic acid mono-2-ethylhexyl ester (P507) dissolved in heptane have been investigated by constant interfacial cell with laminar flow. The effects of the stirring rate, temperature, extractant concentration, and pH on the extraction with mixtures of Cyanex272 and P507 have been studied. The results are compared with those of the system with Cyanex272 or P507 alone. It is concluded that the Yb(III) extraction rate is enhanced with mixtures extractant of Cyanex272 and P507 according to their values of the extraction rate constant, which is due to decreasing the activation energy of the mixtures. At the same time, the mixtures exhibits no synergistic effects for Y(III), which provides better possibilities for Yb(III) and Y(III) separations at a proper conditions than anyone alone. Moreover, thermodynamic extraction separation Yb(III) and Y(III) by the mixtures has been discussed, which agrees with kinetics results. Extraction rate equations have also been obtained, and through the approximate solutions of the flux equation, diffusion parameters and thickness of the diffusion film have been calculated.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

ω循环型预条件矩阵及其应用

In this paper, we apply the preconditioned conjugate gradient method to the solution of positive-definite Toeplitz systems, especially we introduce a new kind of co-circulant preconditioners Pn[ca] by the use of embedding method. We have also discussed the properties of these new preconditioners and proved that many of former preconditioners can be considered as some special cases of Pn[co\. Because of the introduction of co-circulant preconditioners pn[a>], we can greatly overcome the singularity caused by circulant preconditioners. We have discussed the oo-circulant series and functions. We compare the ordinary circularity with the co-circularity, showing that the latter one can be considered as the extended form of the former one; correspondingly, many methods and theorems of the ordinary circularity can be extended. Furthermore, we present the co-circulant decompositional method. By the use of this method, we can divide any co-circulant signal into a summation of many sub-signals; especially among those sub-signals, there are many subseries of which their period is just equal to 1, which are actually the frequency elements of the original co-circulant signal. In this way, we can establish the relationship between the signal and its frequency elements, that is, the frequency elements hi the frequency domain are actually signals with the period of 1 in the spatial domain. We have also proved that the co-circulant has already existed in the traditional Fourier theory. By the use of different criteria for constructing preconditioners, we can get many different preconditioned systems. From the preconditioned systems PN[approximate solutions. Theoretical analysis and actually numerical experiments have already proved the effectiveness of this approximation. By combining the co-circulant with the boundary condition, we can overcome end effects efficiently. We have also presented the co-circulant boundary condition and compared it with the helix boundary condition. Furthermore, we can also get another new boundary condition: the mixed boundary condition. Numerical experiments have proved that the co-circulant boundary condition is better than the mixed one while the helix boundary condition is the worst among them.

Approximate theory of natural-convection with small grashof number

To gain some insight into the behaviour of low-gravity flows in the material processing in space, an approximate theory has been developed for the convective motion of fluids with a small Grashof number Gr. The expansion of the variables into a series of Gr reduces the Boussinesq equation to a system of weakly coupled linearly inhomogeneous equations. Moreover, the analogy concept is proposed and utilized in the study of the plate bending problems in solid mechanics. Two examples are investigated in detail, i. e. the 2-dimensional steady flows in either circular or square infinite closed cylinder, which is horizontally imposed at a specified temperature of linear distribution on the boundaries. The results for stream function ψ, velocity u and temperature T are provided. The analysis of the influences of some parameters such as the Grashof number Gr and the Prandtl number Pr, on motions will lead to several interesting conclusions. The theory seems to be useful for seeking for an analytical solutions. At least, it will greatly simplify the complicated problems originally governed by the Navier-Stokes equation including buoyancy. It is our hope that the theory might be applicable to unsteady or 3-dimensional cases in future.

Approximate analytical description for fundamental-mode fields of graded-index fibers: Beyond the Gaussian approximation

An approximate analytical description for fundamental-mode fields of graded-index fibers is explicitly presented by use of the power-series expansion method, the maximum-value condition at the fiber axis, the decay properties of fundamental-mode fields at large distance from the fiber axis, and the approximate modal parameters U obtained from the Gaussian approximation. This analytical description is much more accurate than the Gaussian approximation and at the same time keep the simplicity of the latter. As two special examples, we present the approximate analytical formulas for the fundamental-mode fields of a step profile fiber and a Gaussian profile fiber, and we find that they are both highly accurate in the single-mode range by comparing them with the corresponding exact solutions.

General Solution to Two-Dimensional Nonslipping JKR Model with a Pulling Force in an Arbitrary Direction

In this paper, a generalized JKR model is investigated, in which an elastic cylinder adhesively contacts with an elastic half space and the contact region is assumed to be perfect bonding. An external pulling force is acted on the cylinder in an arbitrary direction. The contact area changes during the pull-off process, which can be predicted using the dynamic Griffith energy balance criterion as the contact edge shifts. Full coupled solution with an oscillatory singularity is obtained and analyzed by numerical calculations. The effect of Dundurs' parameter on the pull-off process is analyzed, which shows that a nonoscillatory solution can approximate the general one under some conditions, i.e., larger pulling angle (pi/2 is the maximum value), smaller a/R or larger nondimensional parameter value of Delta gamma/E*R. Relations among the contact half width, the external pulling force and the pulling angle are used to determine the pull-off force and pull-off contact half width explicitly. All the results in the present paper as basic solutions are helpful and applicable for experimenters and engineers.

Numerical Solutions for the Transient Flow in the Homogenous Closed Circle Reservoirs

There are many fault block fields in China. A fault block field consists of fault pools. The small fault pools can be viewed as the closed circle reservoirs in some case. In order to know the pressure change of the developed formation and provide the formation data for developing the fault block fields reasonably, the transient flow should be researched. In this paper, we use the automatic mesh generation technology and the finite element method to solve the transient flow problem for the well located in the closed circle reservoir, especially for the well located in an arbitrary position in the closed circle reservoir. The pressure diffusion process is visualized and the well-location factor concept is first proposed in this paper. The typical curves of pressure vs time for the well with different well-location factors are presented. By comparing numerical results with the analytical solutions of the well located in the center of the closed circle reservoir, the numerical method is verified.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.