热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

在同位旋SU(2)不变耦合模型中研究 p pbar -> n nbar 核子反应的重整化无穷阶链图效应

本文采用描述中性及荷电 pion 介子（pion0,pion+,pion-）与核子-反核子强相互作用的同位旋SU(2)不变耦合模型,计算出在pion0 及 pion+, pion- 重整化链图传播下 p pbar -> n nbar 核子反应微分截面的“精确”解析结果;并且将此结果与在 pion0 及 pion+, pion- 树图传播下的微分截面作了对比分析,得到相应的辐射修正重要信息。本文完成的工作对进一步深入研究三种荷电状态 pion 介子与核子-反核子强相互作用的同位旋SU(2)不变耦合模型以及深入理论探讨质子-反质子对撞实验,都将提供某些理论研究的参考价值。     

求解Euler方程的区域分解方法与并行算法

将复杂形状区域划分成多块子区域,研究发展了一种多块区域之间迎风守恒型的内边界耦合方法,实现相邻子区域解的光滑过渡,使多区耦合得到总体流场的数值解。对二维翼型跨音速流动和圆弧形隆起物超音速流动等进行了分区数值计算,并将计算结果与单区计算结果和实验结果作了比较。并行分区计算引入"先进先出"的同步控制等待机制,实现了高效率并行计算,还分析了影响并行效率的主要因素。     

A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

应用多光路主动差分光学吸收光谱仪观测大气污染物的空间分布

研制了一套具有4条测量光路的主动差分光学吸收光谱仪(DOAS)系统,以一盏短孤氙灯为4台集成一体的发射/接收望远镜提供光源,可同时发射4条相互独立的测量光路,4台望远镜的角度可分别独立调整以对准不同位置的角反射镜,光路之间的最大发散角可达28°.为测试仪器性能,使4条光路测量同一区域内的NO₂浓度. 对4条光路的测量结果进行相互对比,任意两条光路的相关系数均在0.97以上,测得NO₂的体积分数平均差异小于10^-9, 表明结果具有良好的一致性 关键词： 差分光学吸收光谱 2')" href="#">NO₂ 2')" href="#">SO₂ 大气污染     

Velocity–vorticity–helicity formulation and a solver for the Navier–Stokes equations

For the three-dimensional incompressible Navier–Stokes equations, we present a formulation featuring velocity, vorticity and helical density as independent variables. We find the helical density can be observed as a Lagrange multiplier corresponding to the divergence-free constraint on the vorticity variable, similar to the pressure in the case of the incompressibility condition for velocity. As one possible practical application of this new formulation, we consider a time-splitting numerical scheme based on an alternating procedure between vorticity–helical density and velocity–Bernoulli pressure systems of equations. Results of numerical experiments include a comparison with some well-known schemes based on pressure–velocity formulation and illustrate the competitiveness on the new scheme as well as the soundness of the new formulation.     

Stability and bifurcation analysis in the delay-coupled van der Pol oscillators

In this paper, the dynamics of a system of two van der Pol equations with a finite delay are investigated. We show that there exist the stability switches and a sequence of Hopf bifurcations occur at the zero equilibrium when the delay varies. Using the theory of normal form and the center manifold theorem, the explicit expression for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived.     

求解二维Poisson方程的重心插值配点法

采用重心Lagrange插值配点法计算了二维Poisson方程.采用重心Lagrange插值法构造近似函数,由配点法离散Poisson方程及其边界条件.数值算例表明方法具有理论简单、计算精度高的特点.     

One-dimensional photonic crystal design

In this article we present a method to determine the band spectrum, band gaps, and discrete energy levels, of a one-dimensional photonic crystal with localized impurities. For one-dimensional crystals with piecewise constant refractive indices we develop an algorithm to recover the refractive index distribution from the period map. Finally, we derive the relationship between the period map and the scattering matrix containing the information on the localized modes.     

Traveling waves for a reaction–diffusion–advection predator–prey model

In this paper we study a reaction–diffusion–advection predator–prey model in a river. The existence of predator-invasion traveling wave solutions and prey-spread traveling wave solutions in the upstream and downstream directions is established and the corresponding minimal wave speeds are obtained. While some crucial improvements in theoretical methods have been established, the proofs of the existence and nonexistence of such traveling waves are based on Schauder’s fixed-point theorem, LaSalle’s invariance principle and Laplace transform. Based on theoretical results, we investigate the effect of the hydrological and biological factors on minimal wave speeds and hence on the spread of the prey and the invasion of the predator in the river. The linear determinacy of the predator–prey Lotka–Volterra system is compared with nonlinear determinacy of the competitive Lotka–Volterra system to investigate the mechanics of linear and nonlinear determinacy.