Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On solving energy-dependent partitioned eigenvalue problem by genetic algorithm: The case of real symmetric Hamiltonian matrices

An energy-dependent partitioning scheme is explored for extracting a small number of eigenvalues of a real symmetric matrix with the help of genetic algorithm. The proposed method is tested with matrices of different sizes (30 × 30 to 1000 × 1000). Comparison is made with Löwdin’s strategy for solving the problem. The relative advantages and disadvantages of the GA-based method are analyzed     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

�ɲ����˶�����ά��������������һ�������ڵ����Զ������������������

??We study the linear quadratic optimal stochastic control problem which is jointly driven by Brownian motion and L\'{e}vy processes. We prove that the new affine stochastic differential adjoint equation exists an inverse process by applying the profound section theorem. Applying for the Bellman's principle of quasilinearization and a monotone iterative convergence method, we prove the existence and uniqueness of the solution of the backward Riccati differential equation. Finally, we prove that the optimal feedback control exists, and the value function is composed of the initial value of the solution of the related backward Riccati differential equation and the related adjoint equation.     

An existence result for a two‐phase Stefan problem arising in metal casting

We present a model arising from the thermal modelling of two metal casting processes. We consider an enthalpy formulation for this two‐phase Stefan problem in a time varying three‐dimensional domain and consider convective heat transfer in the liquid phase. Then, we introduce a weak formulation in a fixed domain, by means of a suitable transformation. Existence of solution is obtained by applying an abstract theorem. The proof of this theorem is done by taking an implicit discretization in time together with a regularization. By passing to the limit in the regularization parameter and in the time step, we obtain the existence of solution of the continuous problem. Copyright © 2005 John Wiley & Sons, Ltd.     

极大单调算子的一个邻近点算法

设E是自反的Banach空间,T∶E→2E是极大单调算子.T-10≠.令x0∈E,yn=(J λnT)-1xn en,xn 1=J-1(αnJxn (1-αn)Jyn),n≥0,λn>0,αn∈[0,1],本文研究了{xn}收敛性.     

Monitoring the Aggregation of Dansyl Chloride in Acetone through Fluorescence Measurements

The aggregation of dansyl chloride (DNS-Cl) in acetone has been studied in detail by steady-state fluorescence techniques.It has been demonstrated that DNS-Cl is stable in acetone during purification and aggregation study processes.The aggregates are not solvolyzed in acetone,and do not take part n any chemical reactions either.It has been found that DNS-Cl tends to aggregate even when its concentration is much lower than its solubility in acetone.The aggregation is reversible,and both the aggregation and the deaggregation are very slow processes.Introduction of SDS has a positive effect upon the formation and stabilization of the aggregates.     

一类二阶两点边值问题的单调迭代方法

通过改进经典的单调迭代方法对于一类二阶两点边值的问题的正解建立了单调迭代程序。这些迭代程序都是从常值函数开始的，因而是可行并且有效的。     

The Generalized KKM Map and Its Applications to Variational Inequalities

In this paper,applying the concept of generalized KKM map,we study problems ofvariational inequalities.We weaken convexity(concavity)conditions for a functional of two variables■(x,y)in the general variational inequalities.Last,we show a proof of non-topological degree meth-od of acute angle principle about monotone operator as an application of these results.     

分析非轴对称腔的复光线数值迭代法

将像散光束的复光线表示法用于分析非轴对称腔，并对有复杂像散的折迭腔作了数值迭代计算以说明方法的应用．