带源项的变系数非线性反应扩散方程的精确解

本文利用广义条件对称方法对带源项的变系数非线性反应扩散方程 f(x)u_t=(g(x)D(u)u_x)_x+h(x)P(u)u_x+q(x)Q(u)进行研究. 当扩散项D(u)取u^m (m≠-1,0,1)和e^u两种重要情形时, 对该方程进行对称约化,得到了具有广义泛函分离变量形式的精确解. 这些精确解包含了该方程对应常系数情况下的解. 关键词： 广义条件对称 精确解 非线性反应扩散方程     

径向积分的解析计算

本文阐明的方法能方便地求出径向积分ⁿ₁^l₁F_Lⁿ₂^l₂(q)≡∫₀^∞> (R_n1^l₁(αr)jL(Kr)R_n2^l₂(αr)r²dr(q≡(K)/(α))的解析表达式。文中给出了L=0,1,2情况下各个径向积分的解析式,同时还以多重求和的形式络出了对各种量子数L,n₁,l₁,n₂,l₂适用的径向积分表达式。     

二维完全阻挫$lt;i$gt;XY$lt;/i$gt;模型的动力学指数

采用大规模动力学蒙特卡罗模拟方法,对二维完全阻挫XY模型的Kosterlitz-Thouless（KT）型相变展开数值研究.系统从有序初始态出发演化到高于KT相变的温度,以普适的动力学标度形式为基础,通过测量磁化和Binder累积量,得出动力学关联时间和平衡态空间关联长度,确定出更精确的动力学指数z.特别是建议并证实了一种在KT相变温度以上（T>T_KT）,独立判断动力学指数z的方法.模拟结果表明,动力学指数z≈2,这与在相变温度以下（T<T_KT）测量的结果一致. 关键词： 蒙特卡罗法 动力学指数 Kosterlitz-Thouless相变 XY模型')" href="#">二维完全阻挫XY模型     

用二维含时Fokker-Planck方程分析Tokamak中α粒子的输运

按分布函数的定义不同,描述高能带电粒子在等离子体中输运的-Planck方程有不同的形式。从数值计算的观点出发对两种不同形式的Fokker-Planck方程作了比较和评价,并指出Fokker-Planck碰撞项可解释为速度空间的对流扩散项。在此基础上用有限差分方法求解二维(速度一维,几何一维)含时Fokker-Planck方程,编制了计算程序CAPT,并将其应用于α粒子的输运研究。最后计算了典型的Tokamak D-T聚变堆参数下α粒子的损失,并给出了堆内α粒子的分布及损失α粒子的速度分布。     

低Re数非线性绕流问题的一种算法

本文对低Re数非线性绕流问題提出了一种算法。做法是以Oseen线化方程为基础通过迭代修正来考虑非线性的惯性项的影响,即先将Navier-Stokes方程写作Oseen方程外加一"强迫函数"项的形式,然后用基本解‘Oseenlet’的积分形式给出其解式,它是一个含物面积分和流动空间积分的非线性积分微分方程。设沿边界上的‘Oseenlet’分布强度为待求量,由边条件加以确定。空间体积分由于有赖于流场,而它在求出解前是未知的,故采用迭代修正的做法进行处理。迭代过程从线性流场开始进行,直到算出达到规定精度要求的收敛解为止。作为算例,对圆柱绕流问題进行了计算,给出了圆柱阻力随Re数的变化规律,并同实验及有关计算资料做了比较。结果表明,本文方法是令人满意的。     

外加直流电场作用下高阶弱非线性复合介质的电势分布

利用同伦分析方法, 研究了一类由柱形杂质随机嵌入基质所形成的、电场和电流密度满足J = σ E + χ |E|²E + η|E|⁴E 形式本构关系的高阶弱非线性复合介质在外加直流电场作用下的电势分布问题. 首先利用模函数展开法, 将本构方程及边界条件化成了一系列非线性常微分方程的边值问题; 再利用同伦分析方法进行求解, 给出了电势在基质和杂质区域的渐近解析解. 关键词： 高阶弱非线性复合介质 模函数展开法 同伦分析方法 电势分布     

室内可吸入颗粒物运动特性的格子 Boltzmann方法模拟

为探究室内可吸入颗粒物的运动特性,使用格子Boltzmann方法,在颗粒物运动概率模型的基础上考虑布朗力对颗粒物的作用,利用改良的LB-CA（Lattice Boltzmann-cellular automata）模型,模拟了Re数分别为400、1 000和2 000时粒径为0.01 μm、0.1 μm和1 μm的颗粒物在上送上回和上送侧回两种回风形式中的运动特性.结果显示：颗粒物的空间分布范围随着Re数的增大而增大,小粒径的颗粒物受到气流湍动和扩散作用的影响更明显;颗粒物的均方位移（mean square displacement,MSD）与Re数、颗粒物粒径的大小成反比,而同样的Re数下,颗粒物在上送侧回的回风形式中其MSD较大.总体上看上送侧回的回风形式具有较低的悬浮颗粒数和更高的室内空气品质.     

(G′/G)展开法和(2+1)维非对称Nizhnik-Novikov-Veselov系统的新精确解

通过引入并扩展 （G′/G）展开法, 构造出（2+1）维非对称Nizhnik-Novikov-Veselov系统的三种形式的新精确通解: 双曲函数通解, 三角函数通解, 有理函数通解. 当双曲函数通解中的常数取特定值时, 通解变为相应孤立波解. 关键词： G′/G）展开法')" href="#">（G′/G）展开法 （2+1）维非对称Nizhnik-Novikov-Veselov系统 精确解 孤立波解     

计算自旋-s算子幺正演化矩阵ds(t)的新方法及其应用

提出了一种严格求解任意自旋-s算子幺正演化矩阵的方法,该方法不同于群论的方法和直接计算的方法,是一种间接的算法.方法的核心是利用两个系统表示的等价性：即自旋-s算子Hamiltonian量H_s=S_x与Heisenberg XX开链带相互作用J_n=n(N-n)的Hamiltonian量的等价性,由于存在这种等价性,自旋-s算子幺正演化矩阵的计算可通过Heisenberg XX开链中态的演化来实现.采用该方法计算了s=3/2,s=2和s=5/2时对应的幺正演化矩阵.由于初始态|sm〉在算子e^-itS_x下的演化实质上相当于对态|sm〉进行一个绕x轴转角为β＝t的转动,演化矩阵元d^s_m'm(t)=〈sm′|e^-itS_x|sm〉就是转动后的态e^-itS_x|sm〉在|sm′〉态上的投影值,所以在t＝π时刻的演化矩阵刚好对应Heisenberg XX开链上量子态的理想传输. 关键词： s算子')" href="#">自旋-s算子 幺正演化矩阵 量子态传输     

初始压入位置对Ni基单晶合金纳米压痕影响研究

针对Ni基单晶合金建立初始压入γ 相的γ /γ' 模型和初始压入γ'相的γ'/γ 模型, 采用分子动力学方法模拟金刚石压头压入两种模型的纳米压痕过程, 计算两种模型[001]晶向硬度. 采用中心对称参数分析两种模型(001)相界面错配位错对纳米压痕过程的影响. 结果显示: 弛豫后, 两种模型(001)相界面错配位错形式不同, 其中γ'/γ 模型(001)相界面错配位错以面角位错形式存在; 压入深度在0.930 nm 之前, 两种模型(001)相界面错配位错变化不大, 压入载荷-压入深度及硬度-压入深度曲线较符合; 压入深度在0.930 nm之后, γ'/γ 模型(001)相界面错配位错长大很多, 导致相同压入深度时γ'/γ 模型比γ /γ'模型压入载荷和硬度计算结果小; 压入深度在2.055 nm之后, γ /γ'模型(001)相界面错配位错对γ 相中位错进入γ'相有阻碍作用, 但仍有部分位错越过(001) 相界面进入γ' 相中, γ'/γ 模型(001)相界面处面角位错对γ' 相中位错进入γ 相有更明显的阻碍作用, 几乎无位错越过(001) 相界面进入γ 相中, 面角位错的强化作用更明显, 所以γ'/γ 模型比γ /γ'模型压入载荷上升速度快.     

Projective and affine connections on S and integrable systems

It is known that the Korteweg–de Vries (KdV) equation is a geodesic flow of an L² metric on the Bott–Virasoro group. This can also be interpreted as a flow on the space of projective connections on S¹. The space of differential operators Δ⁽ⁿ⁾=∂ⁿ+u₂∂ⁿ⁻²++u_n form the space of extended or generalized projective connections. If a projective connection is factorizable Δ⁽ⁿ⁾=(∂−((n+1)/2−1)p₁)(∂+(n−1)/2p_n) with respect to quasi primary fields p_i’s, then these fields satisfy ∑_i=1ⁿ((n+1)/2−i)p_i=0. In this paper we discuss the factorization of projective connection in terms of affine connections. It is shown that the Burgers equation and derivative non-linear Schrödinger (DNLS) equation or the Kaup–Newell equation is the Euler–Arnold flow on the space of affine connections.     

The evolution of parton distributions at small x

We investigate parton distributions at small x using moments of the Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution equation with respect to x. In this representation the kernel of the GLAP equation contains singularities in the moment variable ω at ω = 0. We show that the importance of these singularities at small x depends on the form of the starting distributions. We examine the range of x in which the GLAP equation is valid. The restrictions on the range of x depend on the form of the starting distributions. We investigate whether the GLAP equation can be used to interpolate data in the HERA region. Results are given for the structure function F₂ at small x. A possible method of determining the gluon distribution from F₂ is discussed.     

Integrable geodesic flows and super polytropic gas equations

The polytropic gas equations are shown to be the geodesic flows with respect to an L² metric on the semidirect product space Diff(S¹)C^∞(S¹), where Diff(S¹) is the group of orientation preserving diffeomorphisms of the circle. We also show that the N=1 supersymmetric polytropic gas equation constitute an integrable geodesic flow on the extended Neveu–Schwarz space. Recently other kinds of supersymmetrizations have been studied vigorously in connection with superstring theory and are called supersymmetric-B (SUSY-B) extension. In this paper we also show that the SUSY-B extension of the polytropic gas equation form a geodesic flow on the extension of the Neveu–Schwarz space.     

Riemannian curvature on the group of area-preserving diffeomorphisms (motions of fluid) of 2-sphere

Euler's equation for an incompressible fluid filled in a Riemannian manifold D is regarded as a geodesic equation on the group of volume-preserving diffeomorphisms of D provided with a one-sided invariant metric. A negative sectional curvature implies instability of the geodesic with respect to the corresponding flow and perturbation. The exponential growth of the perturbation is estimated from the values of the sectional curvatures.

This paper presents the expression of the components of Riemannian curvature tensor of the group of area-preserving diffeomorphisms of a 2-sphere in explicit formulas through 3 − j coefficients.     

信息保存方法在微尺度变温流场中的应用研究

在信息保存法(IP)方法的基础上，提出了一个简单易行的温度模型，该模型可以有效模拟计算有温度变化的流动.此外研究发现，在用IP方法模拟计算的流场局部密度变化很大时，如果采用质量守恒方程的二阶中心差分格式更新信息密度时会使计算发散，因此建议采用一阶迎风格式更新信息密度. 关键词： IP方法 低速流动 微尺度 变温流动     

Baryon flow at SIS energies and beyond

We calculate the baryon flow P_x/A(y) in the energy range from 0.25 to 4.0A GeV in a relativistic transport model for Ni + Ni and Au + Au collisions employing various models for the baryon self-energies. We find that to describe the flow data of the FOPI Collaboration above 1A GeV the strength of the vector potential has to be reduced at high relative momentum or at high density. The latter phenomenon leads to a softening of the nuclear equation of state at high density.     

对流扩散方程的涡区分离解法

对流扩散方程是流体计算中一个基本方程,常用的数值方法导至解一个高阶的代数方程组,要求较大的存贮量和较长的计算时间。本文提出一种涡区分离解法,它利用对流扩散方程的迎风性质,把涡区从对流支配区分离出来,仅在各个涡区建立代数方程组并求解。而在对流支配区,则充分利用其抛物性,只需采用显式格式进行计算。由于在各涡区建立的这些方程组阶数和带宽都较小,因此要求存贮量较小,计算速度较快。对于雷诺数较大,涡区范围较小的问题,该方法特别有效。     

Multiple lump molecules and interaction solutions of the Kadomtsev–Petviashvili I equation

In this paper, a modified version of the solution in form of a Gramian formula is employed to investigate a new type of multiple lump molecule solution of the Kadomtsev–Petviashvili I equation. The high-order multiple lump molecules consisting of M N-lump molecules are constructed by means of the Mth-order determinant and the non-homogeneous polynomial in the degree of 2N. The interaction solutions describing P line solitons radiating P of the M N-lump molecules are constructed. The dynamic behaviors of some specific solutions are analyzed through numerical simulation. All the results will enrich our understanding of the multiple lump waves of the Kadomtsev–Petviashvili I equation.     

Relation between conductivity and fluid flow permeability in porous alumina ceramics

Alumina ceramics with various porosity have been used to test the theoretical predictions of a recent model of Wong et al. relating the electrical conductivity and fluid flow permeability of sedimentary rocks. The exponent of Archie's law and modified Kozeny equation are found to be, respectively, m = 1.91 and m′ = 3.82, in agreement with the theoretical value of m′ = 2m. The values of these exponents correspond to a non uniform pore distribution which is indeed what is observed using electron microscopy.     

On the numerical solution of the biharmonic equation in the plane

The biharmonic equation arises in a variety of problems in applied mathematics, most notably in plane elasticity and in viscous incompressible flow. Integral equation methods are natural candidates for the numerical solution of such problems, since they discritize the boundary alone, are easy to apply in the case of free or moving boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains, regardless of shape. In this paper, we follow the work of Mayo and Greenbaum and make use of the Sherman-Lauricella integral equation which is a Fredholm equation with bounded kernel. We describe a fast algorithm for the evaluation of the integral operators appearing in that equation. When combined with a conjugate gradient like algorithm, we are able to solve the discretized integral equation in an amount of time proportional to N, where N is the number of nodes in the discretization of the boundary.