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.     

A penalty method for numerically handling dispersive equations with incompatible initial and boundary data

This article is the numerical counterpart of a theoretical work in progress Qin and Teman, Applicable Anal (2011), 1–19, related to the approximation of evolution hyperbolic equations with incompatible data. The Korteweg‐de Vries and Schrödinger equations with incompatible initial and boundary data are considered here. For hyperbolic equations, the lack of regularity (compatibility) is known to produce large numerical errors which propagate throughout the spatial domain, destroying convergence. In this article, we numerically test the effectiveness of the penalty‐based method proposed in Qin and Teman, Applicable Anal (2011), 1–19, which replaces the hyperbolic equations with incompatible data by a system with compatible data. We observe that convergence is increased. As explained in the text, in the case of the Schrödinger equation, the impact of incompatible (nonregular) data is most severe, and the authors are not aware of any other method that can handle such severe incompatible data. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

3D finite element simulation of explosive welding of three-layer plates

A 3D finite element model of the explosive welding process of three-layer plates with materials of steel-copper-copper is established. Based on the presented model, the bonding mechanism is simulated and analyzed, different detonation modes are also comparatively studied to indicate the driving force spread in few microseconds. The results show that the three layer plates bond together after many times of impact between the flyers and the base driven by detonation wave, which is damping rapidly at each impa...     

Zur kepplimg eines exakten verfahrens mit einem heuristiscieb verfahren für die lösung ganzzahliger linearer optimierungsprobleme

Combining an exact method with a heuristic approach possibilities for solving linear mixed integer optimization problems are investigated. For the considered exact method numerical results with problems from the practice are given. Proper heuristic methods are the interior path methods [2] for which numerical experiences are well-known or the so-called geometric approach [8], Deriving of sufficient conditions for the existence of feasible solutions is possible.     

量子阱层和垒层具有不同Al组分的270/290/330nm AlGaN基深紫外LED光电性能

为了研究AlGaN量子阱层和垒层中Al组分不同对AlGaN基深紫外发光二极管(LED)光电性能的影响,本文利用MOCVD生长、光刻和干法刻蚀工艺制备了AlGaN量子阱层和垒层具有不同Al组分的270/290/330nm深紫外LED,通过实验和数值模拟计算方法发现,量子阱层和垒层中具有低Al组分紫外LED的AlGaN材料具有较低的位错密度、较高的光输出功率和外量子效率。通过电流-电压(I-V)曲线拟合出的较大的理想因子(3.5)和能带结构图表明,AlGaN深紫外LED的电流产生是隧穿机制占据主导作用,这是因为高Al组分AlGaN量子阱中强极化场造成了有源层区域较大的能带弯曲和电势降。     

Resolution estimates for simple one-dimensional flames

The accurate resolution of flame structure is critically important to the direct numerical simulation of reacting flows. While grid-resolution estimates are readily available for cold flows, premixed flames appear to have received relatively little attention. In this paper, a premixed flame characterised by single-step chemistry at moderate-to-high Zel'dovitch numbers (β) is analysed, and its structure is used to provide estimates for sufficient grid resolution. It is found that the critical region of the flame exhibits a weak inverse Zel'dovitch number dependence, and that heuristic methods of resolution estimation based on flame thickness grow relatively less meaningful as β → ∞. Resolution estimates for second- and fourth-order finite-difference schemes are provided.     

Issue Information

Three‐dimensional direct numerical simulation results of flow past a circular cylinder are influenced by numerical aspects, for example the spanwise domain length and the lateral boundary condition adopted for the simulation. It is found that inappropriate numerical set‐up, which restricts the development of intrinsic wake structure, leads to an over‐prediction of the onset point of the secondary wake instability (Re_cr). A best practice of the numerical set‐up is presented for the accurate prediction of Re_cr by direct numerical simulation while minimizing the computational cost. The cylinder span length should be chosen on the basis of the intrinsic wavelength of the wake structure to be simulated, whereas a long span length is not necessary. For the wake transitions above Re_cr, because the wake structures no longer follow particular wavelengths but become disordered and chaotic, a span length of more than 10 cylinder diameters (approximately three times the intrinsic wavelength) is recommended for the simulations to obtain wake structures and hydrodynamic forces that are not strongly restricted by the numerical set‐up. The performances of the periodic and symmetry lateral boundary conditions are compared and discussed. The symmetry boundary condition is recommended for predicting Re_cr, while the periodic boundary condition is recommended for simulating the wake structures above Re_cr. The general conclusions drawn through a circular cylinder are expected to be applicable to other bluff body configurations. Copyright © 2017 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.