Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

共轭A-调和张量的双权Hardy-Littlewood不等式

本文证明了共轭A-调和张量的局部双权积分不等式,此结果类似于共轭调和函数的古典Hardy-Littlewood不等式.作为局部结果的应用,还证明了John域上的共轭A-调和张量的全局双权积分不等式.     

A priori error analysis of the mortar finite element method for variational inequalities

The mortar finite element method is a special domain decomposition method, which can handle the situation where meshes on different subdomains need not align across the interface. In this article, we will apply the mortar element method to general variational inequalities of free boundary type, such as free seepage flow, which may show different behaviors in different regions. We prove that if the solution of the original variational inequality belongs to H²(D), then the mortar element solution can achieve the same order error estimate as the conforming P₁ finite element solution. Application of the mortar element method to a free surface seepage problem and an obstacle problem verifies not only its convergence property but also its great computational efficiency. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Reorientational transition and stripe domains in Co films

We have grown, by sputtering, Co films on (1 1 1) Si substrate, with the aim to find the critical thickness of the reorientational transition of the magnetization from the plane to out of the plane as the hcp axis is oriented perpendicular to the film. Stripe Domains (SD) by Magnetic Force Microscopy, characteristic signature of Perpendicular Magnetic Anisotropy, have been found only in samples grown in some series while samples grown with the same growth parameters do not show stripe domains, indicating an in-plane orientation of the magnetization. These apparently controversial results will be explained in terms of the system fundamental parameters, magnetization, exchange stiffness constant and perpendicular magnetic anisotropy of the sample, which play a crucial role especially in Co films.     

Modeling and simulation of transient responses of a flexible beam floating in finite depth water under moving loads

The dynamic response of a free–free flexible beam floating in an unbounded water domain under the effect of moving loads is numerically analyzed. The water is assumed compressible and inviscid. The surface disturbance satisfies a linear free surface wave condition and an undisturbed condition at infinity. In the present work, a finite element procedure was developed directly in time domain and implemented to solve the two-dimensional problem of the transient behavior of an elastic beam floating on the surface of finite deep water under the passage of a moving force with uniform speed. The presented data demonstrates the applicability of the proposed mathematical model and numerical approach. The influences on the dynamic responses of floating beam of some factors were studied.     

一类指数项级数的若干性质

证明了一类指数项级数有结构简单的收敛域,其和函数的性质与幂级数的和函数性质相似.     

An overlapping domain decomposition preconditioner for a class of discontinuous Galerkin approximations of advection-diffusion problems

We consider a scalar advection-diffusion problem and a recently proposed discontinuous Galerkin approximation, which employs discontinuous finite element spaces and suitable bilinear forms containing interface terms that ensure consistency. For the corresponding sparse, nonsymmetric linear system, we propose and study an additive, two-level overlapping Schwarz preconditioner, consisting of a coarse problem on a coarse triangulation and local solvers associated to a family of subdomains. This is a generalization of the corresponding overlapping method for approximations on continuous finite element spaces. Related to the lack of continuity of our approximation spaces, some interesting new features arise in our generalization, which have no analog in the conforming case. We prove an upper bound for the number of iterations obtained by using this preconditioner with GMRES, which is independent of the number of degrees of freedom of the original problem and the number of subdomains. The performance of the method is illustrated by several numerical experiments for different test problems using linear finite elements in two dimensions.

     

光脉冲传输模拟的快速数值差分递推算法及其应用研究

利用零邻域的马克劳林展开对归一化非线性薛定谔方程的频域差分形式进行分析，得出一个能够同时考虑光学媒质的色散作用以及非线性克尔作用的时域快速数值差分递推公式.选取若干算例，将运用该公式的数值计算结果与已知的解析结果和传统分步傅里叶方法所得的计算结果做了相应的比较，结果表明这种快速数值差分递推公式不但拥有相当快的计算速度，也有很高的计算精度，而且在物理上符合光脉冲在光学媒质中传输时色散和非线性作用同时施加影响的客观实际，这说明它是研究非线性光学媒质中光脉冲传输的一种科学、合理而快速有效的数值计算方法. 关键词： 快速数值差分递推公式 零邻域马克劳林展开 非线性薛定谔方程 光脉冲传输     

Radiation from a Line Source Placed in Two-dimensional Photonic Crystals

A novel formulation of radiation from a localized line source placed in two-dimensional photonic crystals consisting of layered periodic arrays of parallel circular cylinders is presented. The method employs the spectral domain approach. The spectral response of the photonic crystals to the line source excitation is calculated using the lattice sums, the T-matrix of a circular cylinder, and the generalized reflection and transmission matrices of the layered system. The far-zone radiated field is obtained using the conventional asymptotic method to the spectral response. The radiation patterns of the localized line source sandwiched by two photonic crystals are numerically studied.     

Subextensive Singularity in the 2D ± <Emphasis Type="Italic">J</Emphasis> Ising Spin Glass

The statistics of low energy states of the 2D Ising spin glass with +1 and −1 bonds are studied for L × L square lattices with L≤ 48, and p = 0.5, where p is the fraction of negative bonds, using periodic and/or antiperiodic boundary conditions. The behavior of the density of states near the ground state energy is analyzed as a function of L, in order to obtain the low temperature behavior of the model. For large finite L there is a range of T in which the heat capacity is proportional to T ^{5.33 ± 0.12}. The range of T in which this behavior occurs scales slowly to T = 0 as L increases. Similar results are found for p = 0.25. Our results indicate that this model probably obeys the ordinary hyperscaling relation d ν = 2 − α, even though T _c = 0. The existence of the subextensive behavior is attributed to long-range correlations between zero-energy domain walls, and evidence of such correlations is presented. PACS numbers: 75.10.Nr, 75.40.Mg, 75.60.Ch, 05.50.+q