弱形式方程中分片函数奇异性的光滑化分析

研究弱形式方程中的分片光滑函数因求导数出现的奇异性问题，意在为得到离散格式时正确处理这一问题提供依据。在不改变拟合性质的情况下，将分片函数光滑化，则分步积分公式仍然有效。通过光滑后的函数讨论其积分的性质，总结出了几条计算规则。文中的典型例子说明这些计算规则是合理而有效的。     

Some Results on Dynamics of Special Double-Loop ΣΔ-Modulators

The dynamic behaviour of a specific two-dimensional state space model with discontinuity is studied. This model arises from the study of double-loop -modulators with constant input. Using mathematical tools we explain certain simulation results, and some properties are derived. Simulations based on time-varying input are also provided.     

On Thermal Contact of Two Axially Symmetric Elastic Solids

A contact problem of two elastic convex and axially symmetric solids heated (or cooled) to temperatures of different values is considered. Pertinent formulae have been derived for relations between the contact pressure, geometrical characteristics of the solids and distributions of heat flux over the contacting region. We have analysed: 1. The problem of the loss of the contact between two solids pressed together with active heat fluxes. We discuss the cases for which the contact of the axially symmetric solids can take the form of a circle, or an annulus. 2. The problem of a paradox when the mathematically well posed contact problem of thermoelasticity leads to a physically unacceptable solution with a region of overlapping materials. Here we discuss a generalization of the cooled sphere paradox. The heat flux functions are continuously differentiable, of constant sign. The conditions have been derived for the cases when the paradox can be avoided.     

Analysis of Newton’s Method to Compute Travelling Waves in Discrete Media

We present a variant of Newton’s method for computing travelling wave solutions to scalar bistable lattice differential equations. We prove that the method converges to a solution, obtain existence and uniqueness of solutions to such equations with a small second order term and study the limiting behaviour of such solutions as this second order term tends to zero. The robustness of the algorithm will be discussed using numerical examples. These results will also be used to illustrate phenomena like propagation failure, which are encountered when studying lattice differential equations. We finish by discussing the broad application range of the method and illustrate that higher dimensional systems exhibit richer behaviour than their scalar counterparts.     

Perturbation methods for nonlinear autonomous discrete-time dynamical systems

Two perturbation methods for nonlinear autonomous discrete-time dynamical systems are presented. They generalize the classical Lindstedt-Poincaré and multiple scale perturbation methods that are valid for continuous-time systems. The Lindstedt-Poincaré method allows determination of the periodic or almost-periodic orbits of the nonlinear system (limit cycles), while the multiple scale method also permits analysis of the transient state and the stability of the limit cycles. An application to the discrete Van der Pol equation is also presented, for which the asymptotic solution is shown to be in excellent agreement with the exact (numerical) solution. It is demonstrated that, when the sampling step tends to zero the asymptotic transient and steady-state discrete-time solutions correctly tend to the asymptotic continuous-time solutions.     

A numerical exploration of secondary separation in starting flow around a flat plate by the discrete vortex model

In the present work, a further numerical simulation of the starting flow around a flat plate normal to the direction of motion in a uniform fluid has been made by means of the discrete vortex method. The secondary separation occurring at rear surface of the plate is explored, and predicted approximately using Thwait's method. The calculated results show that in the early stages of the flow secondary separation does occur. The evolution of flow field, the vortex growing process and the characteristics of secondary vortices have been described. The time dependent drag coefficients, the vorticity shed from the edges and rear surface, and the separation positions are calculated as well as the distributions of velocity and pressure on the plate. In the case of flow normal to the plate, the calculated secondary vortices are weak. Their existence will change the local velocity distributions and affect pressure distributions. However, the effect on drag coefficient is negligible.     

Further study on vortex shedding flow by a circular cylinder

In the present paper the mechanism involved in vortex shedding flows is investigated in detail. In the early stage of the unsteady separated flow the interaction between secondary vortices and primary vortices is quite strong. In the later stage of the flow, corresponding to the vortex shedding the recirculating flow region on each side of the aft body goes through such a cycle: growth-contraction-growth, the secondary separation occurs periodically rather than continuously. The reduction of circulation is taken into account in three cases with different decay factors to study its influence on the prediction of main flow characteristics. Results show that to simulate vortex shedding flow it is necessary to include the reduction of circulation to bring the calculated results into good agreement with experiments. An improved discrete vortex model is suggested in which both the secondary separation and the reduction are incorporated. The processes of vortex shedding, the forces prediction and other flow characteristics are given and some discussions are made. Porject is supported by National Natural Science Foundation of China.     

Direct Statistical Simulation of a Supersonic Flow of a Binary Mixture of Rarefied Gases around a Transversely Positioned Cylinder

A supersonic flow of a binary mixture of gases in a wide range of rarefaction (from a flow with a Knudsen number K _n = 0.1 to a free-molecular flow) around a cylinder is studied by means of direct statistical Monte Carlo simulations (DSMC method). The influence of a small fraction of heavy particles in a light gas flow on the region of significant nonequilibrium near the cylinder and on the heat flux is considered.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 53–59, September–October, 2005.     

Numerical analysis of funnel formation during unsteady fluid outflow from a rotating cylindrical vessel

The main features of the unsteady outflow of a fluid from a cylindrical vessel rotating together with it at a constant angular velocity through a central drain hole in the bottom are studied. The software package STAR-CD tested on the results of experiments with water in the absence of rotation is used. Certain important features of the phenomenon under consideration related with the formation of vortex funnels in the fluid are established. The effect of the main control parameters of the problem, namely, angular velocity, viscosity, initial depth, etc., is analyzed.