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

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

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.     

Augmented Lagrangian Algorithms Based on the Spectral Projected Gradient Method for Solving Nonlinear Programming Problems

The spectral projected gradient method SPG is an algorithm for large-scale bound-constrained optimization introduced recently by Birgin, Martínez, and Raydan. It is based on the Raydan unconstrained generalization of the Barzilai-Borwein method for quadratics. The SPG algorithm turned out to be surprisingly effective for solving many large-scale minimization problems with box constraints. Therefore, it is natural to test its perfomance for solving the sub-problems that appear in nonlinear programming methods based on augmented Lagrangians. In this work, augmented Lagrangian methods which use SPG as the underlying convex-constraint solver are introduced (ALSPG) and the methods are tested in two sets of problems. First, a meaningful subset of large-scale nonlinearly constrained problems of the CUTE collection is solved and compared with the perfomance of LANCELOT. Second, a family of location problems in the minimax formulation is solved against the package FFSQP.     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

Stability Conditions of the MMAP [ K ]/ G [ K ]/1/ LCFS Preemptive Repeat Queue

In this paper, we study the stability conditions of the MMAP[K]/G[K]/1/LCFS preemptive repeat queue. We introduce an embedded Markov chain of matrix M/G/1 type with a tree structure and identify conditions for the Markov chain to be ergodic. First, we present three conventional methods for the stability problem of the queueing system of interest. These methods are either computationally demanding or do not provide accurate information for system stability. Then we introduce a novel approach that develops two linear programs whose solutions provide sufficient conditions for stability or instability of the queueing system. The new approach is numerically efficient. The advantages and disadvantages of the methods introduced in this paper are analyzed both theoretically and numerically.     

晶体微观结构的数值计算

晶体微观结构是晶体材料在特定物理条件下其多个能量极小平衔态在空间形成的某种微尺度的规则分布．几何非线性的连续介质力学理论可以用能量极小化原理来解释晶体微观结构的形成，并用Young测度来刻画平衡态各变体在空间的概率分布．定性的理解与定量地分析和计算晶体材料的微观结构对于发展和改进高级晶体功能材料，如形状记忆合金、铁电体、磁至伸缩材料等，有重要的意义．本文回顾了近年来晶体微观结构数值计算方面的最新进展．介绍了计算晶体微观结构的几种数值方法及有关的数值分析结果。     

美式期权定价中非局部问题的有限元方法

在本文中 ,我们关心的是美式期权的有限元方法 .首先 ,根据 [4 ]我们对所讨论的问题引进一个新奇的实用的方法 ,它涉及到对原问题重新形成准确的数学公式 ,使得数值解的计算可以在非常小的区域上进行 ,从而该算法计算速度快精度高 .进而 ,我们利用超逼近分析技术得到了有限元解关于 L2 -模的最优估计 .     

Laguerre approximation of stable manifolds with application to connecting orbits

We present an algorithm, based on approximation by Laguerre polynomials, for computing a point on the stable manifold of a stationary solution of an autonomous system. A superconvergence phenomenon means that the accuracy of our results is much higher than the usual spectral accuracy. Both the theory and the implementation of the method are considered. Finally, as an application of the algorithm, we describe a fully spectral approximation of homo- and heteroclinic orbits.

     

Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.