带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

特征值问题的预变换方法(I): 杨辉三角阵变换与二阶PDE 特征多项式

本文提出一类求解特征值问题的下三角预变换方法, 目标是通过相似变换后矩阵下三角元素平方和明显减少、且变换后的特征值及其特征向量较易求解, 使变换后的对角线可作为全体特征值很好的一组初值, 其作用如同对于解方程组找到好的预条件子, 加速迭代收敛. 以二阶PDE 数值计算为例,对于以Laplace 方程为代表的特征波向量组及正交多项式组有广泛的应用前景.
杨辉三角是我国古代数学家的一项重要成就. 本文引入杨辉三角矩阵作为预变换子, 给出一般矩阵用杨辉三角矩阵作为左、右预变换子时变为上三角矩阵的充要条件, 给出了元素为行指标二次多项式的两个矩阵类(三对角线阵与五对角线阵) 中特征值何时保持二次多项式的充要条件, 并应用于构造新的二元PDE 正交多项式.     

一类推广后的Feigenbaum函数方程的光滑解

考虑一类推广后的Feigenbaum函数方程｛g（0）=1,-1≤g（x）≤1,x∈[-1,1],h（g（x））=g（g（h（x）））其中h（x）是[-1,1]上的递减光滑奇函数且满足h（0）=0,-1〈h’（x）〈0,z∈[-1,1]．利用构造性方法讨论上述方程的光滑解的存在性及唯一性．     

Next‐to‐next‐to‐leading order post‐Newtonian spin(1)‐spin(2) Hamiltonian for self‐gravitating binaries

We present the next‐to‐next‐to‐leading order post‐Newtonian (PN) spin(1)‐spin(2) Hamiltonian for two self‐gravitating spinning compact objects. If both objects are rapidly rotating, then the corresponding interaction is comparable in strength to a 4PN effect. The Hamiltonian is checked via the global Poincaré algebra with the center‐of‐mass vector uniquely determined by an ansatz.     

Improvements on open and traction boundary conditions for Navier–Stokes time-splitting methods

We present in this paper a numerical scheme for incompressible Navier–Stokes equations with open and traction boundary conditions, in the framework of pressure-correction methods. A new way to enforce this type of boundary condition is proposed and provides higher pressure and velocity convergence rates in space and time than found in the present state of the art. We illustrate this result by computing some numerical and physical tests. In particular, we establish reference solutions of a laminar flow in a geometry where a bifurcation takes place and of the unsteady flow around a square cylinder.     

A fast directional algorithm for high-frequency electromagnetic scattering

This paper is concerned with the fast solution of high-frequency electromagnetic scattering problems using the boundary integral formulation. We extend the O(N log N) directional multilevel algorithm previously proposed for the acoustic scattering case to the vector electromagnetic case. We also detail how to incorporate the curl operator of the magnetic field integral equation into the algorithm. When combined with a standard iterative method, this results in an almost linear complexity solver for the combined field integral equations. In addition, the butterfly algorithm is utilized to compute the far field pattern and radar cross section with O(N log N) complexity.     

Artificial boundary conditions for the numerical solution of the Euler equations by the discontinuous galerkin method

We present artificial boundary conditions for the numerical simulation of compressible flows using high-order accurate discretizations with the discontinuous Galerkin (DG) finite element method. The construction of the proposed boundary conditions is based on characteristic analysis and applied for boundaries with arbitrary shape and orientation. Numerical experiments demonstrate that the proposed boundary treatment enables to convect out of the computational domain complex flow features with little distortion. In addition, it is shown that small-amplitude acoustic disturbances could be convected out of the computational domain, with no significant deterioration of the overall accuracy of the method. Furthermore, it was found that application of the proposed boundary treatment for viscous flow over a cylinder yields superior performance compared to simple extrapolation methods.     

Numerical modeling of transient two-dimensional viscoelastic waves

This paper deals with the numerical modeling of transient mechanical waves in linear viscoelastic solids. Dissipation mechanisms are described using the generalized Zener model. No time convolutions are required thanks to the introduction of memory variables that satisfy local-in-time differential equations. By appropriately choosing the relaxation parameters, it is possible to accurately describe a large range of materials, such as solids with constant quality factors. The evolution equations satisfied by the velocity, the stress, and the memory variables are written in the form of a first-order system of PDEs with a source term. This system is solved by splitting it into two parts: the propagative part is discretized explicitly, using a fourth-order ADER scheme on a Cartesian grid, and the diffusive part is then solved exactly. Jump conditions along the interfaces are discretized by applying an immersed interface method. Numerical experiments of wave propagation in viscoelastic and fluid media show the efficiency of this numerical modeling for dealing with challenging problems, such as multiple scattering configurations.     

Random laser action in bovine semen

Experiments using bovine semen reveal that the addition of a high-gain water soluble dye results in random laser action when excited by a Q-switched, frequency doubled, Nd:Yag laser. The data shows that the linewidth collapse of the emission is correlated to the sperm count of the individual samples, potentially making this a rapid, low sample volume approach to count determination.     

Sixth moment optimization of multi-channel Bragg filters

The design of multichannel optical filter is studied theoretically. New efficient two-step method based on combination of 6-th moment minimization and straightforward numerical minimization is developed. It is applied to find a fiber Bragg grating profile with minimal refractive index variation for 5 ÷ 32 channels. Obtained results are closer to the theoretical limit than previously reported.