原子尺度毛细凝聚的开尔文方程

实验表明宏观尺度建立的开尔文方程经适当改造后依然可描述原子尺度水的毛细凝聚现象.      

量子化学在凝聚炸药研究中的应用

详细地介绍了量子化学计算分子结构的方法，并例举了在研究凝聚炸药中的各种应用。     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

“凝聚物质冲击压缩”研究的若干论题

讨论了凝聚物质冲击压缩研究领域包括的范围和各类问题所占的比例，对1999年和2001年举行的两次APS凝聚物质冲击压缩专题会议发表的论文进行了分类统计，并对会议反映的最新动态和我所关心的热点问题的意义作了简要的说明。     

对流-扩散相互作用结构的不变性

本文提出并证明了不可压缩剪切层流中对流-扩散相互作用结构不变性诸定理:即二难剪切层流与其线性化及非线性扰动存在同一的对流-扩散相互作用结构,且物理尺度(指时间、空间和速度尺度)相同。给出十个推论,例如:对流-扩散相互作用可在剪切层流及其扰动场内“激发“快时间尺度和小空间尺度结构,线性化稳定性原理的约定对剪切流体系统成立等。应用题例导出计及时间-空间尺度效应和非平行流效应的广义Orr-Sommerfeld(GOS)方程,证实它有两个粘性解:阻尼层解和干扰层解;经典OS方程及其两个粘性解:边界层解和Heisenberg临界层解,Triple-deck稳定性理论基本方程及其两个粘性解,均是本文GOS方程及其两个粘性解的特例。     

凝聚炸药爆轰研究中的光电测试技术进展

着重介绍了８０年代末以来凝聚炸药爆轰研究中，光电测试技术的发展状况，主要包括引爆机理研究的谱方法和传播爆轰测试所用的ＶＩＳＡＲ、ＦＰＩ、ＥＭＶ、阵铜汁、冲击极化汁、散离探针技术等，对两种典型的数学化调整ＣＣＤ相机系统地作了简单的介绍。     

连续体结构拓扑优化应力约束凝聚化的ICM方法

为克服应力约束下拓扑优化问题约束数目多、应力敏度计算量大的困难，提出 了应力约束化凝聚化的ICM方法. 在利用Mises强度理论将应力约束转换成应变能约束后， 提出了应力约束凝聚化的两条途径：其一为应力全局化的方法，其二为应力约束集成化的方 法. 由此建立了多工况下以重量为目标、以凝聚化应变能为约束的连续体结构优化模型，并 利用对偶理论对优化模型进行了求解. 4个数值算例表明：该方法具有较高的计算效率，得 到的拓扑结构比较合理，不仅适用于二维连续体结构，也适用于三维连续体结构.     

充液或部分充液转子的动力稳定性

本文讨论了具有内外阻尼的高速充液转子的动力稳定性。首先通过对旋转流体的平面流场的求解，导出充液转子作简谐运动时流体对转子的动压力，由此导出转子的运动方程；讨论了充液转子的动力稳定性，给出了稳定性解析判据和稳定性边界。结果表明，存在转速门槛值，低于该转速时，充液转子可存在稳定区；当高于该转速时，系统永远失稳，这一结论复盖了已有文献的结果。     

两自由度可调非线性减振器

研究了可调非线性减振器的优化设计.基于哈密尔顿最小势能原理建立非线性动力学模型,系统局域参数内,实现非线性系统幅值优化.利用平均法求解可调非线性减振器频响方程.分析系统解的稳定性,优化系统参数,降低系统幅值响应.     

基于四叉树法的有限元网格自动生成及凝聚方法

本文介绍了用于有限元网格自动划分的四叉树方法及相应的非均匀网格的生成方法──网格凝聚法。此方法可以用于处理任意形状的单连通或多连通的平面结构，其边界以折线及二次曲线描述。相应的网格生成器用户界面友好，极少需要人工干预。所生成的单元大多为四边形元，在边界处理时用了少量三角元，所有单元性态良好。非均匀网格可以实现多处、多重加密。     

On the local stability analysis of the approximate harmonic balance solutions

Periodic response of nonlinear oscillators is usually determined by approximate methods. In the "steady state" type methods, first an approximate solution for the steady state periodic response is determined, and then the local stability of this solution is determined by analyzing the equation of motion linearized about this predicted "solution". An exact stability analysis of this linear variational equation can provide erroneous stability type information about the approximate solutions. It is shown that a consistent stability type information about these solutions can be obtained only when the linearized variational equation is analyzed by approximate methods, and the level of accuracy of this analysis is consistent with that of the approximate solutions. It is demonstrated that these consistent stability results do not imply that the approximate solution is qualitatively correct. It is also shown that the difference between an approximate and the next higher order stability analysis can be used to "guess" the role of higher harmonics in the periodic response. This trial and error procedure can be used to ensure the qualitatively correct and numerically accurate nature of the approximate solutions and the corresponding stability analysis.     

Dynamic analysis of mathematical model of ethanol fermentation with gas stripping

In this paper, a mathematical model for ethanol fermentation with gas stripping is investigated. Firstly, the model with continuous substrate input is taken. We study the existence and local stability of two equilibrium points. According to Poincare–Bendixson Theorem, the sufficient condition for the globally asymptotical stability of positive equilibrium point is obtained, which implies that we can get stable ethanol product. Secondly, we study the model with impulsive substrate input and obtain the sufficient condition for the local stability of cell-free periodic solution by using the Floquet’s theory of impulsive differential equation and small-amplitude perturbation skills. In a certain limiting case, it is shown that a nontrivial periodic solution emerges via a supercritical (subcritical) bifurcation. Finally, our results are confirmed by means of numerical simulation.     

局部彼得洛夫-伽辽金法分析各向异性板屈曲

基于Kirchhoff板理论和对挠度函数采用移动最小二乘近似函数进行插值，进一步研究无网格局部Petrov-Galerkin(MLPG)方法在各向异性板稳定问题中的应用．分析中，本质边界条件采用罚因子法施加，离散的特征值方程由板稳定控制方程的局部积分对称弱形式中得到．通过数值算例并与其他方法的结果进行比较，表明MLPG法求解各向异性薄板稳定问题具有收敛性好、精度高等一系列优点．     

Wave scattering from a rough surface: Solution by an iterative method

Two iterative solutions of the Helmholtz equation for a scalar field in above a rough surface that admits the Dirichlet boundary condition are derived. The bases for the two iterative methods are two different boundary integral equations that represent the solution. The first integral equation is classified as a Fredholm integral equation of the first kind. The second is classified as a Fredholm integral equation of the second kind. This classification suggests that it is easier to find stable solution methods to the second equation. In both methods, the boundary integral was separated into a major part which is easy to calculate and a local residual part. The major part is a convolution and thus can be calculated using FFT in complexity O(N log N), where N is the number of surface points in which the surface height and its first derivatives together with the incoming wave and its normal derivative are all known. The residual element of the equations can be approximated efficiently only for surfaces where their amplitude is less than the wavelength of the incoming wave. The iterative schemes were tested numerically against a reference solution in order to examine the applicability range, the error estimation and the stability of the schemes. All tests supported the superiority of the second method. In particular the error estimation and stability tests indicated good performance for surfaces with slope up to 1. Yet, being an equation in the scattered field alone, makes the first method useful as a benchmark solution in its domain of applicability.     

Stability and Hopf bifurcation in a three-species system with feedback delays

A kind of three-species system with Holling type II functional response and feedback delays is introduced. By analyzing the associated characteristic equation, its local stability and the existence of Hopf bifurcation are obtained. We derive explicit formulas to determine the direction of the Hopf bifurcation and the stability of periodic solution bifurcated out by using the normal-form method and center manifold theorem. Numerical simulations confirm our theoretical findings.     

Asymptotic Stability of N-Solitary Waves of the FPU Lattices

We study stability of N-solitary wave solutions of the Fermi-Pasta-Ulam (FPU) lattice equation. Solitary wave solutions of the FPU lattice equation cannot be characterized as critical points of conservation laws due to the lack of infinitesimal invariance in the spatial variable. In place of standard variational arguments for Hamiltonian systems, we use an exponential stability property of the linearized FPU equation in a weighted space which is biased in the direction of motion. The dispersion of the linearized FPU equation balances the potential term for low frequencies, whereas the dispersion is superior for high frequencies.We approximate the low frequency part of a solution of the linearized FPU equation by a solution to the linearized Korteweg-de Vries (KdV) equation around an N-soliton solution. We prove an exponential stability property of the linearized KdV equation around N-solitons by using the linearized Bäcklund transformation and use the result to analyze the linearized FPU equation.     

Stability analysis for two kinds of equations in two-species population dynamics

In this paper we study the stability for equilibrium points of equations in two-population dynamics.We discuss two predator-prey-patch models.Model1 is described by a differential equation.Model2 is described by an integral differential equation.We obtain the conditions for the stability of their equilibrium points.The results show that the overall population stability despite local extinction is realizable.     

Propagation of Perturbations in Plane Poiseuille Flow between Walls of Nonuniform Compliance

The evolution of small perturbations in longitudinally nonuniform flows is studied with reference to the problem of the propagation of flow perturbations in a plane channel with walls of variable elasticity. Using the solution of the problem of the receptivity of the flow to local vibrations of the walls, the problem considered can be reduced to the solution of an integral equation for a single function, namely, the complex vibration amplitude of the walls. A numerical method for solving this equation on the basis of a piecewise-linear approximation of the unknown function is proposed. It is shown that the instability wave amplitude changes discontinuously at the junction of the rigid and elastic channel sections. A second method of investigating the process of propagation of perturbations in the flow considered is proposed. This method is based on laws of evolution of perturbations in nonuniform flows and an analytic solution of the problem of perturbation scattering on the junction of walls with different compliance. On the basis of this method the classical stability theory is generalized to include the case of nonuniform flows.     

Propagation of a Tollmien-Schlichting wave over the junction between rigid and compliant surfaces

The propagation of an instability wave over the junction region between rigid and compliant panels is studied theoretically. The problem is investigated using three different methods with reference to flow in a plane channel containing sections with elastic walls. Within the framework of the first approach, using the solution of the problem of flow receptivity to local wall vibration, the problem considered is reduced to the solution of an integro-differential equation for the complex wall oscillation amplitude. It is shown that at the junction of rigid and elastic channel walls the instability-wave amplitude changes stepwise. For calculating the step value, another, analytical, method of investigating the perturbation propagation process, based on representing the solution as a superposition of modes of the locally homogeneous problem, is proposed. This method is also applied to calculating the flow stability characteristics in channels containing one or more elastic sections or consisting of periodically alternating rigid and compliant sections. The third method represents the unknown solution as the sum of a local forced solution and a superposition of orthogonal modes of flow in a channel with rigid walls. The latter method can be used for calculating the eigenvalues and eigenfunctions of the stability problem for flow in a channel with uniformly compliant walls.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 31–48. Original Russian Text Copyright © 2004 by Manuilovich.