Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.

     

Tomographic reconstruction using heuristic Monte Carlo methods

A tomographic reconstruction method based on Monte Carlo random searching guided by the information contained in the projections of radiographed objects is presented. In order to solve the optimization problem, a multiscale algorithm is proposed to reduce computation. The reconstruction is performed in a coarse-to-fine multigrid scale that initializes each resolution level with the reconstruction of the previous coarser level, which substantially improves the performance. The method was applied to a real case reconstructing the internal structure of a small metallic object with internal components, showing excellent results.     

ITERATIVE APPROXIMATION OF FIXED POINTS OF NON-LIPSCHITZIAN ASYMPTOTICALLY PSEUDOCONTRACTIVE MAPPINGS

The purpose of this paper is to investigate the problem of approximating fixed points of non-Lipschitizian asymptotically pseudocontractive mappings in an arbitrary real Banach space by the modified Ishikawa iterative sequences with errors.     

激光诱导间质肿瘤热疗的数值模拟和实验研究

本文在考虑生物组织物性动态变化的情况下建立了激光诱导间质肿瘤热疗(LITT)的物理数学模型,采用MonteCarlo方法数值模拟了LITT中激光能量在生物组织内的传输过程,基于Pennes生物传热方程和Arrhenius方程数值求解了组织内的温度分布和热损伤体积的变化,分析了热物性及血液灌注率的动态变化对LITT过程的影响,并与相应的离体实验结果进行了对比。数值模拟结果表明,组织的热物性及血液灌注率的动态变化对于热损伤体积的变化具有重要的影响。因此在激光诱导间质肿瘤热疗的数值模拟中应该考虑热物性及血液灌注率的动态变化以期为临床治疗方案的制定提供更为准确的依据。     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

晶格失配对异质外延超薄膜生长中成核特性的影响

利用动力学蒙特卡罗方法模拟了异质外延超薄膜生长中的成核过程.研究了薄膜与衬底的晶格失配对超薄膜生长中成核密度、平均核尺寸、标度关系及生长模式的影响.结果发现产生压（张）应变的晶格负（正）失配使生长过程更早（迟）从成核区进入过渡区，失配越大，这一效应越明显.在相同的沉积条件下，负失配导致超薄膜形成较低的成核密度与较大的平均核尺寸，而正失配则相反.成核密度满足标度关系Ns≈(F/D)χ，随着失配度从-0.04增加到0.02，标度系数χ从0.37逐渐减小到0.33，对应超薄膜生长过程从包含二聚体扩散模式转变到无 关键词： 薄膜生长 成核 晶格失配 蒙特卡罗模拟     

闪烁光纤在γ射线辐照下部分特性的蒙特卡罗模拟

 闪烁光纤在射线成像方面的应用越来越广泛。为了进一步了解闪烁光纤在射线辐照下的基本特性，基于蒙特卡罗方法，利用计算机模拟分析了γ射线在闪烁光纤中的好事例率与光纤长度及射线能量的关系，能量沉积效率与光纤长度及射线能量的关系。此项工作对于闪烁光纤阵列在射线成像，剂量场测量等方面的研究很有价值。     

Studies on the Adsorption of Asymmetrical Triblock Copolymers by Scheutjens-Fleer Theory and Monte Carlo Simulation

The adsorption of asymmetrical triblock copolymers from a non-selective solvent on solid surface has been studied by using Scheutjens-Fleer mean-field theory and Monte Carlo simulation method on lattice model. The main aim of this paper is to provide detailed computer simulation data, taking A8-kB20Ak as a key example, to study the influence of the structure of copolymer on adsorption behavior and make a comparison between MC and SF results. The simulated results show that the size distribution of various configurations and density-profile are dependent on molecular structure and adsorption energy. The molecular structure will lead to diversity of adsorption behavior. This discrepancy between different structures would be enlarged for the surface coverage and adsorption amount with increasing of the adsorption energy. The surface coverage and the adsorption amount as well as the bound fraction will become larger as symmetry of the molecular structure becomes gradually worse. The adsorption layer becomes thicker with increasing of symmetry of the molecule when adsorption energy is smaller but it becomes thinner when adsorption energy is higher. It is shown that SF theory can reproduce the adsorption behavior of asymmetrical triblock copolymers. However, systematic discrepancy between the theory and simulation still exists.The approximations inherited in the mean-filed theory such as random mixing and the allowance of direct back folding may be responsible for those deviations.     

On General Mixed Quasivariational Inequalities

In this paper, we suggest and analyze several iterative methods for solving general mixed quasivariational inequalities by using the technique of updating the solution and the auxiliary principle. It is shown that the convergence of these methods requires either the pseudomonotonicity or the partially relaxed strong monotonicity of the operator. Proofs of convergence is very simple. Our new methods differ from the existing methods for solving various classes of variational inequalities and related optimization problems. Various special cases are also discussed.     

A monotone iterative scheme for a nonlinear second order equation based on a generalized anti–maximum principle

In this paper, a generalized anti–maximum principle for the second order differential operator with potentials is proved. As an application, we will give a monotone iterative scheme for periodic solutions of nonlinear second order equations. Such a scheme involves the L^p norms of the growth, 1 ≤ p ≤ ∞, while the usual one is just the case p = ∞.