Reducing Huge Gyroscopic Eigenproblems by Automated Multi-level Substructuring

Simulating numerically the sound radiation of a rolling tire requires the solution of a very large and sparse gyroscopic eigenvalue problem. Taking advantage of the automated multi-level substructuring (AMLS) method it can be projected to a much smaller gyroscopic problem, the solution of which however is still quite costly since the eigenmodes are non-real and complex arithmetic is necessary. This paper discusses the application of AMLS to huge gyroscopic problems and the numerical solution of the AMLS reduction. A numerical example demonstrates the efficiency of AMLS.     

Stability of chirped bright and dark soliton-like solutions of the cubic complex Ginzburg-Landau equation with variable coefficients

We consider an inhomogeneous optical fiber system described by the generalized cubic complex Ginzburg-Landau (CGL) equation with varying dispersion, nonlinearity, gain (loss), nonlinear gain (absorption) and the effect of spectral limitation. Exact chirped bright and dark soliton-like solutions of the CGL equation were found by using a suitable ansatz. Furthermore, we analyze the features of the solitons and consider the problem of stability of these soliton-like solutions under finite initial perturbations. It is shown by extensive numerical simulations that both bright and dark soliton-like solutions are stable in an inhomogeneous fiber system. Finally, the interaction between two chirped bright and dark soliton-like pulses is investigated numerically.     

On an Application of Newton's Method to Nonlinear Operators with w-Conditioned Second Derivative

We present a new approach to study the convergence of Newton's method in Banach spaces, which relax the conditions appearing in the usual studies. The approach is based on the bound required for the second derivative of the operator involved. An application to a nonlinear integral equation is presented.     

On a class of new and practical performance indexes for approximation of fold bifurcations of nonlinear power flow equations

Efficient measurement of the performance index (the distance of a loading parameter from the voltage collapse point) is one of the key problems in power system operations and planning and such an index indicates the severity of a power system with regard to voltage collapse. There exist many interesting methods and ideas to compute this index. However, some successful methods are not yet mathematically justified while other mathematically sound methods are often proposed directly based on the bifurcation theory and they require the initial stationary state to be too close to the unknown turning point to make the underlying methods practical.This paper first gives a survey of several popular methods for estimating the fold bifurcation point including the continuation methods, bifurcation methods and the test function methods (Seydel's direct solution methods, the tangent vector methods and the reduced Jacobian method) and discuss their relative advantages and problems. Test functions are usually based on scaling of the determinant of the Jacobian matrix and it is generally not clear how to determine the behaviour of such functions. As the underlying nonlinear equations are of a particular type, this allows us to do a new analysis of the determinants of the Jacobian and its submatrices in this paper. Following the analysis, we demonstrate how to construct a class of test functions with a predictable analytical behaviour so that a suitable index can be produced. Finally, examples of two test functions from this class are proposed. For several standard IEEE test systems, promising numerical results have been achieved.     

Stable multi-sets

In this paper we introduce a generalization of stable sets: stable multi-sets. A stable multi-set is an assignment of integers to the vertices of a graph, such that specified bounds on vertices and edges are not exceeded. In case all vertex and edge bounds equal one, stable multi-sets are equivalent to stable sets.  For the stable multi-set problem, we derive reduction rules and study the associated polytope. We state necessary and sufficient conditions for the extreme points of the linear relaxation to be integer. These conditions generalize the conditions for the stable set polytope. Moreover, the classes of odd cycle and clique inequalities for stable sets are generalized to stable multi-sets and conditions for them to be facet defining are determined.  The study of stable multi-sets is initiated by optimization problems in the field of telecommunication networks. Stable multi-sets emerge as an important substructure in the design of optical networks. Received: February 14, 2001/Revised version: September 7, 2001     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

Generalized quasilinearization method for reaction-diffusion equations under nonlinear and nonlocal flux conditions

In this paper we consider an initial boundary value problem for a reaction-diffusion equation under nonlinear and nonlocal Robin type boundary condition. Assuming the existence of an ordered pair of upper and lower solutions we establish a generalized quasilinearization method for the problem under consideration whose characteristic feature consists in the construction of monotone sequences converging to the unique solution within the interval of upper and lower solutions, and whose convergence rate is quadratic. Thus this method provides an efficient iteration technique that produces not only improved approximations due to the monotonicity of its iterates, but yields also a measure of the convergence rate.     

烷氧基取代聚对苯乙炔三阶非线性光学性能

通过双醚化反应、氯甲基化反应以及在强碱性条件下进行的脱氯化氢反应制备聚（2-甲氧基-5-丁氧基）对苯乙炔（PMOBOPV）、聚（2-甲氧基-5-己氧基）对苯乙炔（PMOHOPV）、聚（2,5-二丁氧基）对苯乙炔（PDBOPV）和聚（2,5-二己氧基）对苯乙炔（PDHOPV）等四种可溶性聚对苯乙炔（PPV）衍生物，通过紫外-可见吸收光谱对产物分子结构进行表征.结果显示，PMOBOPV、PMOHOPV、PDBOPV和PDHOPV的共轭π电子发生π → π* 跃迁的吸收峰分别位于491 nm、495nm、504nm和510nm处，相应的光学禁带宽度分别为2.23eV、2.18eV、2.12eV和2.07eV.利用简并四波混频技术测量PPV衍生物的三阶非线性光学性能，探讨了分子结构对PPV衍生物三阶非线性极化率（χ(3)）的影响.研究发现，激发波长为532 nm时，PMOBOPV、PMOHOPV、PDBOPV和PDHOPV的共振χ(3)值分别为3.45×10－10、5.13×10－10、7.15×10－10和9.61×10－10 esu；激发波长为1064 nm时，它们的非共振χ(3)值分别为1.09×10－11、1.42×10－11、1.62×10－11和2.14×10－11 esu.     

Simulation of the nonlinear vibration of a string using the Cellular Automata based on the reflection rule

In this study, the nonlinear dynamic responses of a string are simulated using the Cellular Automata method based on the reflection rule. In the case of nonlinear systems, the velocity of wave propagation is not constant and depends on the amplitude. A new treatment of the dynamic time step is proposed for the Cellular Automata method considering the effect of the propagation velocity. As numerical examples, first, the dynamic responses of a string with linear characteristic are simulated using the Cellular Automata method. A typical resonance curve can be obtained. Second, the dynamic responses of a string with nonlinear characteristic are simulated using the proposed method. Some characteristic types of vibration can be obtained. It is concluded that the linear and nonlinear dynamic responses of a string may be obtained by simulation using the Cellular Automata method.     

Time delay as a key to apoptosis induction in the p53 network

A feedback mechanism that involves the proteins p53 and mdm2, induces cell death as a controlled response to severe DNA damage. A minimal model for this mechanism demonstrates that the response may be dynamic and connected with the time needed to translate the mdm2 protein. The response takes place if the dissociation constant k between p53 and mdm2 varies from its normal value. Although it is widely believed that it is an increase in k that triggers the response, we show that the experimental behaviour is better described by a decrease in the dissociation constant. The response is quite robust upon changes in the parameters of the system, as required by any control mechanism, except for few weak points, which could be connected with the onset of cancer. Received 8 May 2002 / Received in final form 9 July 2002 Published online 17 September 2002