一类离散的非线性爆炸方程平衡点的稳定性

本所研究的非线性爆炸方程实质上是由可数无穷多个彼此相互关联的非线性常微分方程所组成的自治系统，它刻划了在只有基本粒子与i-粒子(i≥1)进行碰撞反应的系统里，粒子增长过程中密度随时间变化规律。本证明了如果系数满足一定的假设，那么在爆炸占优的条件下，这一系列的平衡点在Lyapunov意义下是稳定的．     

带有碰撞爆炸的离散的凝结方程密度守恒解的唯一性

在条件aij≤Amin{i，j}，i，J≥1(其中A是非负常数)成立的情况下，证明了带有碰撞爆炸的离散的凝结方程的密度守恒解是唯一的，所获得的结论改进了已有的结果．     

液滴碰撞和聚合模型研究

建立了描述液滴碰撞和聚合过程的数学模型,设计了一种高效的计算液滴碰撞对搜索算法,并在已有研究成果的基础上对液滴碰撞的结果做了进一步的发展．借助平滑粒子流体动力学方法两种实现方式各自的优点,将液滴间的相互作用局限在其周围一定数目的液滴之间,并采用积分核函数定义了液滴间碰撞的概率,通过数值模拟探讨了模型的特性．结果表明,所建模型对所采用的计算网格没有明显的依赖性,具有较高的计算精度和计算效率,不但能很好地维持系统动量的守恒性,而且对液滴初始尺寸分布没有明显的依赖性．     

相互碰撞的圆粒子在竖直通道中沉降的数值研究

采用格子Boltzmann方法,数值模拟了考虑相互碰撞的两圆粒子的沉降,分析了Re数、初始相对位移及通道宽度对粒子沉降的影响．结果表明,在0.1＜Re＜20范围内, 粒子沉降具有周期性．Re数越大,两圆粒子的相互作用越强,粒子横向位移的幅度也越大．在大Re数时,沉降的过程是两个粒子交替领先;在较小Re数时,当后面的粒子接近的时候,领先的粒子被向右侧推了一段后仍会继续领先;对中等Re数,原先在后面的粒子在第一次加速后将取得领先位置并一直保持下去．粒子的初始分布位置对沉降的形态影响不大．管道宽度变化时,粒子总的沉降特性不变,而周期改变,管道越宽,周期越长．     

双自由度非定点斜碰撞振动系统的动力学分析

对两个单摆组成的双自由度、非定点、斜碰撞振动系统的动力学行为进行了详细研究．揭示了在双自由度、非定点、斜碰撞过程中恢复系数、摩擦系数、系统参数和碰撞前后系统状态之间的关系．基于Poincaré映射方法和非定点斜碰撞关系推导出该系统单碰周期n次谐运动存在性判据．根据Floquet理论分析了该系统次谐运动周期解的稳定性问题,给出了Floquet特征乘子的计算公式．通过数值仿真证实了该方法的有效性,同时分析了非定点、斜碰撞系统碰撞点位置的概率分布情况．     

Brown凝并中两个不同直径纳米颗粒的碰撞系数

对Brown凝并中两个不同直径纳米颗粒的碰撞系数进行了研究,通过求解碰撞方程,获得了在van der Waals力和弹性变形力作用下,直径为100 nm至750 nm的邻苯二甲酸二辛酯纳米颗粒的碰撞系数．发现碰撞系数总体上随着颗粒直径和2个颗粒半径比的增加而减小;当颗粒直径为550 nm时,碰撞系数有一个突然的增加．最后给出了具有不同直径2个纳米颗粒碰撞系数的新表达式．     

基于耦合纯无网格方法时间分数阶下孤立子波碰撞过程的数值模拟研究北大核心CSCD

为数值预测时间分数阶耦合非线性Schrödinger(TF-CNLS)方程描述的孤立子波非弹性碰撞过程,首次发展了一种耦合纯无网格有限点集法(coupled finite pointset method,CFPM).其构造过程为:1)对时间分数阶Caputo导数项采用一种高精度的差分格式;2)对空间导数采用基于Taylor展开和加权最小二乘法的有限粒子法(FPM)离散格式;3)对区域进行局部加密和采用稳定性好的双曲余弦核函数以提高数值精度.数值研究中,首先,运用CFPM对有解析解的一维TF-CNLS方程进行求解,分析了节点均匀分布或局部加密情况下的误差和收敛阶,表明给出的耦合无网格法具有近似二阶精度和易局部加密求解的灵活性;其次,运用CFPM对无解析解一维TF-CNLS方程描述的孤立子波非弹性碰撞过程进行了数值预测,其出现的波塌缩现象与整数阶下出现的多波现象截然不同;最后,与有限差分结果作对比,表明CFPM数值预测时间分数阶下孤立子波非弹性碰撞过程的复杂传播现象是可靠的.     

窄带随机噪声激励下线性碰撞系统的响应

研究了单自由度线性单边碰撞系统在窄带随机噪声激励下的次共振响应问题．用Zhuravlev变换将碰撞系统转化为连续的非碰撞系统,然后用随机平均法得到了关于慢变量的随机微分方程．在约束距离为0时,用矩方法给出了系统响应幅值二阶矩的解析表达式．在约束距离不为0时,近似地得到了系统响应幅值二阶矩的解析表达式．讨论了系统阻尼项、窄带随机噪声的带宽和中心频率以及碰撞恢复系数等参数对于系统响应的影响．理论计算和数值模拟表明,系统响应幅值将在激励频率接近于次共振频率时达到最大,而当激励频率逐渐偏离次共振频率时,系统响应迅速衰减．数值模拟表明提出的方法是有效的．     

158 A Ge V p-Pb Pb-Pb碰撞中的单奇异和多奇异粒子产生

用描写相对论性核 核碰撞的LUCIAE模型及相应的MonteCarlo事例产生器分析了WA97最近发表的 1 5 8AGeVp Pb和Pb Pb碰撞中单奇异粒子 (Λ ,Λ)和多奇异粒子 (Ξ-,Ξ-,Ω-,Ω-)的多重数和横质量分布 .由于LUCIAE模型包含了弦碎裂微观过程中奇异夸克压低因子随碰撞体系的能量、中心度和质量的增大而增大 ,以及由它导致的相对论性核 核碰撞宏观过程中奇异粒子产额随上述三因素之增而增的物理机制 ;LUCIAE模型计算结果能较好描写WA97实验所揭示的相对论性核 核碰撞中奇异粒子产额随中心度之增而增和奇异粒子增强随奇异夸克数之增而增的实验事实 ,以及单奇异和多奇异粒子的横质量分布 .     

相对论性核与乳胶(Ag)碰撞中灰粒子的系统研究 *

用描写相对论性核 核碰撞的LUCIAE模型和相应的事件产生器系统地研究了 1 4.6 ,6 0和 2 0 0AGeV的O核、2 0 0AGeV的Si和S核以及 1 1 .6AGeV的Au核与乳胶 (Ag)碰撞中灰粒子产生的平均多重数、多重数分布以及角分布 3个物理量同入射能量、射弹质量及碰撞中心度间的关系 ,还研究了再散射在灰粒子产生机制中的作用 .LUCIAE的这些研究结果与相应的EMU0 1乳胶实验结果都相一致 .     

Dark solitonic excitations and collisions from a fourth-order dispersive nonlinear Schrödinger model for the alpha helical protein

Recent protein observations motivate the dark-soliton study to explain the energy transfer in the proteins. In this paper we will investigate a fourth-order dispersive nonlinear Schrödinger equation, which governs the Davydov solitons in the alpha helical protein with higher-order effects. Painlevé analysis is performed to prove the equation is integrable. Through the introduction of an auxiliary function, bilinear forms and dark N-soliton solutions are constructed with the Hirota method and symbolic computation. Asymptotic analysis on the two-soliton solutions indicates that the soliton collisions are elastic. Decrease of the coefficient of higher-order effects can increase the soliton velocities. Graphical analysis on the two-soliton solutions indicates that the head-on collision between the two solitons, overtaking collision between the two solitons and collision between a moving soliton and a stationary one are all elastic. Collisions among the three solitons are all pairwise elastic.     

New Approaches and New Applications for Computer Simulation of N-Body Problems

In this paper we study new computer models for turbulent flow in the small and in the large, for soliton collision, and for the fundamental problem of electrostatics. These are typical for models of microdrop collision, stress of a slotted copper plate, contact angle of adhesion, biological self reorganization, the bounce of an elastic ball, the motion of a top on a smooth surface, and elastic snap-through, which are only mentioned and referenced appropriately.     

球形纳米颗粒Brown凝并碰撞效率的新表达式

研究了邻苯二甲酸二辛酯纳米颗粒在Brown凝并过程中的碰撞效率．在考虑Stokes阻力、润滑力、van der Waals力和颗粒变形恢复力的情况下推导了一组碰撞方程,通过数值求解这组方程,得到了当颗粒半径从50 nm到500 nm变化时,颗粒碰撞效率和半径之间的关系,计算得到的结果和实验结果符合较好．计算结果表明,在颗粒半径为50 nm到500 nm的情况下,颗粒的碰撞效率随颗粒半径的增加而减小．基于计算结果,提出了颗粒碰撞效率的新表达式．     

Global analysis of boundary and interior crises in an elastic impact oscillator

The crisis phenomena of a Duffing–Van der Pol oscillator with a one-side elastic constraint are studied by the composite cell coordinate system method in this paper. By computing the global properties such as attractors, basins of attraction and saddles, the vivid evolutionary process of two kinds of crises: boundary crisis and interior crisis are shown. The boundary crisis is resulted by the collision of a chaotic attractor and a periodic saddle on the basin boundary. It is observed that there are two types of interior crises. One is caused by the collision of a chaotic attractor and a chaotic saddle within the interior of basin of attraction. The other one occurs because a period attractor collides with a chaotic saddle within the interior of basin of attraction. The saddles of system play an important role in the crisis process. The results show that this method is an efficient tool to perform the global analysis of elastic impact oscillators.     

Energy conservation property of MW-CRF deterministic particle method

A deterministic particle method for kinetic equation in several dimensions was presented by Motta and Wick some years ago. In the present paper, we formulate the method in general coordinates and apply this formulation to an elastic collision kernel. We show that the method preserves exactly the energy conservation property of the elastic kernel.     

A general mathematical model for the collision between a free-fall hammer of a pile-driver and an elastic pile: Continuous dependence and low-frequency asymptotic expansion

The main goal of this paper is to study with continuous dependence and a low-frequency asymptotic expansion of a unique weak solution to an initial-boundary value problem. This problem is arisen from a general mathematical model for the collision between the free-fall hammer of a pile-driver and an elastic pile whose ends are furnished with a bearing.     

Bright-soliton collisions with shape change by intensity redistribution for the coupled Sasa–Satsuma system in the optical fiber communications

By means of symbolic computation and Darboux transformation, analytically and numerically investigated in this paper is a two-coupled Sasa–Satsuma system, which can describe the pulse propagation in birefringent fibers, so as to increase the bit rate in optical fibers, or achieve wavelength-division multiplexing. Analytical bright N-soliton solution of the system is firstly derived. Based on the bright one- and two-soliton solutions, numerical simulation and figure illustration are carried out on through the multi-parametric management, i.e., different choices among eight parameters in the two-soliton solutions. The interaction mechanisms for the bright two-solitons are revealed in three aspects: Separating evolution behaviors, elastic collision behaviors and inelastic collision behaviors. There exist three different cases for the inelastic collision for the two-soliton, which reflect correspondingly different energy transfer mechanisms (by intensity redistribution) between the two components: Manakov-typed collision; a near-elastic collision and another completely inelastic collision between the two components; and four single-solitons in two components undergo shape changes (inelastic and elastic) due to intensity redistribution, where one single-soliton keeps invariant and the other three single-solitons change during the collision. The collision mechanisms may be viewed as the two-solitons interact in a waveguide supporting propagation of two nonlinear waves simultaneously. In general, partial suppression (enhancement) of intensity between the components is dependent on the values of the soliton parameters.     

Soliton solutions via auxiliary function method for a coherently-coupled model in the optical fiber communications

A coherently-coupled nonlinear Schrödinger system in the optical fiber communications, with the mixed self-phase modulation (SPM), cross-phase modulation (XPM) and positive coherent coupling terms, is studied through the bilinear method with an auxiliary function. Solutions for that system are found to be of two types: singular and non-singular ones, and the latter appear as the soliton-typed. Vector bright one- and two-solitons are derived with the corresponding phase-shift parameter constraints. In virtue of computerized symbolic computation and asymptotic behavior analysis, elastic collision mechanisms of such vector solitons are investigated. With the aid of graphical simulation, vector solitons are displayed to be of the single- or double-hump profiles. The formation and collision mechanisms of the vector bright solitons for that system are generated based on the combined effects of SPM, XPM and coherent coupling. Only elastic collisions of the vector solitons occur for that system, which is a distinctive feature amid those of other coherently-coupled nonlinear Schrödinger systems.     

Mixed impedance transmission problems for vibrating layered elastic bodies

Direct scattering problems for partially coated piecewise homogenous and inhomogeneous layered obstacles in linear elasticity lead to mixed impedance transmission problems for the steady‐state elastic oscillation equations. For a piecewise homogenous isotropic composite body, we employ the potential method and reduce the mixed impedance transmission problem to an equivalent system of boundary pseudodifferential equations. We give a detailed analysis of the corresponding pseudodifferential operators, which live on the interface between the layers and on a proper submanifold of the boundary of the composite elastic body, and establish uniqueness and existence results for the original mixed impedance transmission problem for arbitrary values of the oscillation frequency parameter; this is crucial in the study of inverse elastic scattering problems for partially coated layered obstacles. We also investigate regularity properties of solutions near the collision curves, where the different boundary conditions collide, and establish almost best Hölder smoothness results. Further, we analyze the asymptotic behavior of the stress vector near the collision curve and derive explicit formulas for the stress singularity exponents. The case of Lipschitz surfaces is briefly treated separately. In the case of a composite body containing homogeneous or inhomogeneous finite anisotropic inclusions, we develop an alternative hybrid method based on the so‐called nonlocal approach and reduce the mixed transmission problem to an equivalent functional‐variational equation with a sesquilinear form that ‘lives’ on a bounded part of the layered composite body and its boundary. We show that this sesquilinear form is coercive and that the corresponding variational equation is uniquely solvable. Copyright © 2015 John Wiley & Sons, Ltd.