Comparison study of several numerical integration schemes used in calculations of density functional theory

Several numerical integration schemes for the evaluation of matrix elements in density functional theory calculations have been studied and compared by computational practice. The best scheme was found to be the combination of the atomic partition function proposed by Becke with the scaled generalized Gauss-Laguerre quadrature formula for radial integration suggested by Yang, which achieve the highest convergence rate to the numerical integration. With the same number of integration points, the accuracy of the calculated results by this scheme is higher by 1 to 2 orders of magnitudes than that by other schemes. The reason for achieving higher accuracy by this scheme has been proposed preliminarily.     

强流电子束泵浦XeCl准分子激光动力学模型

强流电子束泵浦XeCl准分子激光动力学模型由三部分组成,即电子束能量沉积的计算;电子温度、电子反应速率的计算和化学/激光动力学。这个模型可以准确地预报小信号增益、吸收等激光特征量的时间变化规律。 该模型是在文献[1]、[2]、[3]报导的动力学模型的基础上提出的。采用四阶龙格_库塔法在VAX—11/780机器上进行数值求解。 计算给出XeCl准分子激光反应过程中各种粒子浓度、小信号增益、吸收、输出光强、激发速率以及平均电子能量随时间的变化规律。计算结果表明本征效率是激发速率,电流密度和工作气体(Ne/Xe/HCl)的各分压比的函数。 该模型可为高功率准分子激光器的研制提供设计参数和最佳实验条件。     

Numerical solution of minimax problems of optimal control,part 1

This paper contains general transformation techniques useful to convert minimax problems of optimal control into the Mayer-Bolza problem of the calculus of variations [Problem (P)]. We consider two types of minimax problems: minimax problems of Type (Q), in which the minimax function depends on the state and does not depend on the control; and minimax problems of Type (R), in which the minimax function depends on both the state and the control. Both Problem (Q) and Problem (R) can be reduced to Problem (P).For Problem (Q), we exploit the analogy with a bounded-state problem in combination with a transformation of the Jacobson type. This requires the proper augmentation of the state vectorx(t), the control vectoru(t), and the parameter vector , as well as the proper augmentation of the constraining relations. As a result of the transformation, the unknown minimax value of the performance index becomes a component of the parameter vector being optimized.For Problem (R), we exploit the analogy with a bounded-control problem in combination with a transformation of the Valentine type. This requires the proper augmentation of the control vectoru(t) and the parameter vector , as well as the proper augmentation of the constraining relations. As a result of the transformation, the unknown minimax value of the performance index becomes a component of the parameter vector being optimized.In a subsequent paper (Part 2), the transformation techniques presented here are employed in conjunction with the sequential gradient-restoration algorithm for solving optimal control problems on a digital computer; both the single-subarc approach and the multiple-subarc approach are discussed.This research was supported by the National Science Foundation, Grant No. ENG-79-18667, and by Wright-Patterson Air Force Base, Contract No. F33615-80-C3000. This paper is a condensation of the investigations reported in Refs. 1–7. The authors are indebted to E. M. Coker and E. M. Sims for analytical and computational assistance.     

Back-and-forth shooting method for solving two-point boundary-value problems

The paper proposes a special iterative method for a nonlinear TPBVP of the form (t)=f(t, x(t),p(t)), (t)=g(t, x(t),p(t)), subject toh(x(0),p(0))=0,e(x(T),p(T))=0. Certain stability properties of the above differential equations are taken into consideration in the method, so that the integration directions associated with these equations respectively are opposite to each other, in contrast with the conventional shooting methods. Via an embedding and a Riccati-type transformation, the TPBVP is reduced to consecutive initial-value problems of ordinary differential equations. A preliminary numerical test is given by a simple example originating in an optimal control problem.     

Convergence of the Finite Difference Scheme for a Fast Diffusion Equation in Porous Media

We are concerned with an implicit scheme for the finite difference solution to a nonlinear parabolic equation with a multivalued coefficient that describes the fast diffusion in a porous medium. The boundary conditions contain the multivalued function as well. We prove the stability and the convergence of the scheme, emphasizing the precise nature of convergence in this specific case, and compute the error level of the approximating solution. The method is aimed to simplify the numerical computations for the solutions to equations of this type, without performing an approximation of the multivalued function. The theory is illustrated by numerical results.     

Numerical study on soot removal in partial oxidation of methane to syngas reactors

The serious carbon deposition existing in catalytic partial oxidation of methane (CPOM) to syngas process is one of the key problems that impede its industrialization. In this study, 3-dimensional unsteady numerical simulations of the soot formation and oxidation in oxidation section in a heat coupling reactor were carried out by computational fluid dynamics (CFD) approach incorporating the Moss-Brookes model for soot formation. The model has been validated and proven to be in good agreement with experiment results. Effects of nozzle type, nozzle convergence angle, channel spacing, number of channels, radius/height ratio, oxygen/carbon ratio, preheat temperature and additional introduction of steam on the soot formation were simulated. Results show that the soot formation in oxidation section of the heat coupling reactor depends on both nozzle structures and operation conditions, and the soot concentration can be greatly reduced by optimization with the maximum mass fraction of soot inside the oxidation reactor from 2.28% to 0.0501%, and so that the soot mass fraction at the exit reduces from 0.74% to 0.03%.     

Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods

In this paper, we present a finite element method with a residual‐based artificial viscosity for simulation of turbulent compressible flow, with adaptive mesh refinement based on a posteriori error estimation with sensitivity information from an associated dual problem. The artificial viscosity acts as a numerical stabilization, as shock capturing, and as turbulence capturing for large eddy simulation of turbulent flow. The adaptive method resolves parts of the flow indicated by the a posteriori error estimates but leaves shocks and turbulence under‐resolved in a large eddy simulation. The method is tested for examples in 2D and 3D and is validated against experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

Accurate integration over atomic regions bounded by zero‐flux surfaces

The approach for the integration over a region covered by zero‐flux surface is described. This approach based on the surface triangulation technique is efficiently realized in a newly developed program TWOE . The elaborated method is tested on several atomic properties including the source function. TWOE results are compared with those produced by using well‐known existing programs. Absolute errors in computed atomic properties are shown to range usually from 10^?6 to 10^?5 au. The demonstrative examples prove that present realization has perfect convergence of atomic properties with increasing size of angular grid and allows to obtain highly accurate data even in the most difficult cases. It is believed that the developed program can be bridgehead that allows to implement atomic partitioning of any desired molecular property with high accuracy. © 2012 Wiley Periodicals, Inc.     

��β�ӳ��������ģ�����Ʋ����ʵĽ�����

In this paper, for a kind of risk models with heavy-tailed and delayed claims, we derive the asymptotics of the infinite-time ruin probability and the uniform asymptotics of the finite-time ruin probability. The numerical simulation results are also presented. The results of theoretical analysis and numerical simulation show that the influence of the delay for the claim payment is nearly negligible to the ruin probability when the initial capital and running-time are all large.     

Moving grid method for numerical simulation of stratified flows

We develop a second‐order accurate Navier–Stokes solver based on r‐adaptivity of the underlying numerical discretization. The motion of the mesh is based on the fluid velocity field; however, certain adjustments to the Lagrangian velocities are introduced to maintain quality of the mesh. The adjustments are based on the variational approach of energy minimization to redistribute grid points closer to the areas of rapid solution variation. To quantify the numerical diffusion inherent to each method, we monitor changes in the background potential energy, computation of which is based on the density field. We demonstrate on a standing interfacial gravity wave simulation how using our method of grid evolution decreases the rate of increase of the background potential energy compared with using the same advection scheme on the stationary grid. To further highlight the benefit of the proposed moving grid method, we apply it to the nonhydrostatic lock‐exchange flow where the evolution of the interface is more complex than in the standing wave test case. Naive grid evolution based on the fluid velocities in the lock‐exchange flow leads to grid tangling as Kelvin–Helmholtz billows develop at the interface. This is remedied by grid refinement using the variational approach. Copyright © 2012 John Wiley & Sons, Ltd.