Local and Parallel Finite Element Algorithms for Eigenvalue Problems

Abstract Some new local and parallel finite element algorithms are proposed and analyzed in this paper foreigenvalue problems.With these algorithms, the solution of an eigenvalue problem on a fine grid is reduced tothe solution of an eigenvalue problem on a relatively coarse grid together with solutions of some linear algebraicsystems on fine grid by using some local and parallel procedure.A theoretical tool for analyzing these algorithmsis some local error estimate that is also obtained in this paper for finite element approximations of eigenvectorson general shape-regular grids.     

Second-order random wave solutions for interfacial internal waves in N-layer density-stratified fluid

This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins （1963） and Dean （1979） in the study of random surface waves and by Song （2004） in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts： the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song （2004） for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.     

On a stochastic nonlinear equation in one-dimensional viscoelasticity

In this paper we discuss an initial-boundary value problem for a stochastic nonlinear equation arising in one-dimensional viscoelasticity. We propose to use a new direct method to obtain a solution. This method is expected to be applicable to a broad class of nonlinear stochastic partial differential equations.

     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Generic replica symmetric field-theory for short range Ising spin glasses

Symmetry considerations and a direct, Hubbard-Stratonovich type, derivation are used to construct a replica field-theory relevant to the study of the spin glass transition of short range models in a magnetic field. A mean-field treatment reveals that two different types of transitions exist, whenever the replica number n is kept larger than zero. The Sherrington-Kirkpatrick critical point in zero magnetic field between the paramagnet and replica magnet (a replica symmetric phase with a nonzero spin glass order parameter) separates from the de Almeida-Thouless line, along which replica symmetry breaking occurs. We argue that for studying the de Almeida-Thouless transition around the upper critical dimension d = 6, it is necessary to use the generic cubic model with all the three bare masses and eight cubic couplings. The critical role n may play is also emphasized. To make perturbative calculations feasible, a new representation of the cubic interaction is introduced. To illustrate the method, we compute the masses in one-loop order. Some technical details and a list of vertex rules are presented to help future renormalisation-group calculations. Received 9 October 2001     

Diffusion from sessile droplets through membranes

We analyze diffusion from a periodic array of hemispherical droplets through a membrane. We find that the multiple sources do not interact strongly, even when the droplets are closely spaced, so that the flux through the membrane appears nearly additive.     

仓库容量有限条件下的随机存贮管理问题的评注

给出2005年全国部分高校研究生数学建模竞赛D题的研究背景,对参赛论文作了简要评注,并提供了一种参考解答.     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Inequalities for the Moments and Distribution of the Ladder Height of a Random Walk

We obtain upper bounds for the tail distribution of the first nonnegative sum of a random walk and for the moments of the overshoot over an arbitrary nonnegative level if the expectation of jumps is positive and close to zero. In addition, we find an estimate for the expectation of the first ladder epoch.     

电感耦合等离子体原子发射光谱法测定中药方剂煎煮液中的微量元素

本文采用HNO3－H2O2对中药方剂煎煮液进行处理，拟定了电感耦合等离子体原子发射光谱法同时测定试样中12个微量元素的方法，选择了仪器的最佳工作条件，考察了基本效应的影响。在选定条件下，12个微量元素的检出限为0．10－5．0ng/g，相对标准偏差为0.10%-9.5%(n=10),加标回收率为85％－103％。该方法简单、快速，已用于大量样品的检测，分析结果满意。