借助Matlab遗传算法反演火山区压力源参数

首先分析了Mogi模型的经典最小二乘反演由于出现矩阵奇异或接近奇异而发散的局限性,决定采用具有全局高效搜索性和优化计算时不需要求导的遗传算法来反演火山区压力源参数。反演程序调用Matlab7.0 Release14工具箱中的选择、交叉和变异等相关程序。最后应用该程序和长白山天池火山区水准和GPS观测资料联合反演火山区压力源参数,经过比较遗传算法和阻尼最小二乘方法的反演结果,可看出2种方法结果近似,可见此方法方便快捷地实现了目标函数的最小化。     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

对求解非线性方程方法的探索

BaiCL在2001年提出了双曲函数法，并把解设成具体的双曲函数，即sinhw和coshw的线性组合．本文把解设成是由2个函数f(ξ)和g(ξ)组成的线性组合，它们可以有多种取法，从而使方程更直接、更有效．借助数学软件Maple，用改进后的方法和吴文俊消元法，求解Bai C L文中的一个例子，获得了包含Bai文结果在内的更为丰富、精确的行波解．     

真空冷冻干燥机凝结器的改进

分析了真空冷冻干燥机水汽凝结器存在容易造成"冰塞"的原因,并提出了安装气流扩散器的改进方法,对改进前后的凝结器捕水性能进行了比较.由于气流扩散器的安装,捕水周期得以延长,捕水能力及工作效率显著提高,实用性较好,而且扩散器的制作费用低,性价比高,有推广价值.     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

潮流计算中的数学方法——牛顿—拉夫逊法

本文介绍了一种解非线性代数方程组的数学方法——牛顿-拉夫逊法,并将其应用于电力系统潮流计算中。采用本文方法的计算结果表明该教学方法在电力系统潮流问题中是收敛的;同时计算结果与电力专业软件的计算结果保持了一致,从而验证了结果的正确性。     

对面向对象分析法进行比较评估的一种新方法

提出一种新的比较框架，其中包含了一系列面向对象分析（OOA）方法共同的基本特性．以5种著名的面向对象分析方法为例，阐述了如何通过它们的特性细节，对面向对象分析方法进行分析、比较和评估．     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

共轭斜量法在某些特征向量上收敛速度的估计

设 Ａ是对称正定矩阵，λ_1是 Ａ 的最大或最小特征值，χ_1是对应的特征向量．｛zk｝是用共轭斜量法求解方程组 Αχ＝ｂ时的近似解序列，ｅi＝Ａ～(_1)b-zi,本文给出了|x_1～Tei|较合理的上界估计式。从而为分析预处理共轭斜量法提供了进一步的理论基础。     

Detecting instabilities in flows of viscoelastic fluids

There is a growing interest in developing numerical tools to investigate the onset of physical instabilities observed in experiments involving viscoelastic flows, which is a difficult and challenging task as the simulations are very sensitive to numerical instabilities. Following a recent linear stability analysis carried out in order to better understand qualitatively the origin of numerical instabilities occurring in the simulation of flows viscoelastic fluids, the present paper considers a possible extension for more complex flows. This promising method could be applied to track instabilities in complex (i.e. essentially non‐parallel) flows. In addition, results related to transient growth mechanism indicate that it might be responsible for the development of numerical instabilities in the simulation of viscoelastic fluids. Copyright © 2003 John Wiley & Sons, Ltd.